1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
|
#############################################################
## ##
## Copyright (c) 2003-2011 by The University of Queensland ##
## Earth Systems Science Computational Centre (ESSCC) ##
## http://www.uq.edu.au/esscc ##
## ##
## Primary Business: Brisbane, Queensland, Australia ##
## Licensed under the Open Software License version 3.0 ##
## http://www.opensource.org/licenses/osl-3.0.php ##
## ##
#############################################################
"""Defines the L{Cylinder} base class"""
from esys.lsm.vis import core
from Modifier import Modifiable
class Cylinder(core.Cylinder, Modifiable):
def __init__(self, endPt1, endPt2, radius):
core.Cylinder.__init__(self, endPt1, endPt2, radius)
Modifiable.__init__(self)
def writeSdl(self, f):
f.write("cylinder {")
f.write(
"\n<"
+
",".join(
map(
str,
core.Vec3(self.getEndPt1())-core.Vec3(self.getCenter())
)
)
+
"> "
)
f.write(
"\n<"
+
",".join(
map(
str,
core.Vec3(self.getEndPt2())-core.Vec3(self.getCenter())
)
)
+
"> "
)
f.write(str(self.getRadius()))
f.write(" ")
Modifiable.writeSdl(self, f)
f.write(" translate <")
f.write(",".join(map(str,self.getCenter())))
f.write(">")
f.write("}")
class Disk(Cylinder):
"""
Objects of this class represent cylinders.
"""
def __init__(self, center, radius, height=None, direction=None):
if (height == None):
height = radius*0.01
if (direction == None):
direction = core.Vec3(0,0,1)
direction = core.Vec3(direction)
nrm = direction.norm()
if (nrm == 0.0):
direction = core.Vec3(0,0,1)
nrm = 1.0
center = core.Vec3(center)
direction = direction/nrm
Cylinder.__init__(
self,
endPt1 = center - (direction*(height*0.5)),
endPt2 = center + (direction*(height*0.5)),
radius = radius
)
|
