1 # SPDX-License-Identifier: GPL-2.0-or-later
3 from .vefm_271
import mesh
, vertex
, edge
, face
4 from math
import pi
, acos
, sin
, cos
, atan
, tan
, fabs
, sqrt
7 def check_contains(cl
, name
, print_value
=False):
10 if el
.startswith("_"):
15 print(name
, " contains ==>", el
, " value = ", tmp
)
17 print(name
, " contains ==>", el
)
18 print("\ncheck_contains finished\n\n")
25 self
.PKHG_parameters
= None
29 self
.skeleton
= [] # List of verts in the full skeleton edges.
30 self
.vertskeleton
= [] # config needs this member
31 self
.edgeskeleton
= [] # config needs this member
32 self
.sphericalverts
= []
41 self
.eccentricity
= None
58 self
.panelpoints
= None
59 self
.paneledges
= None
60 self
.reversepanel
= None
61 self
.edgelength
= None
65 def setparameters(self
, params
):
66 parameters
= self
.PKHG_parameters
= params
67 self
.frequency
= parameters
[0] # How many subdivisions - up to 20.
68 self
.eccentricity
= parameters
[1] # Elliptical if >1.0.
69 self
.squish
= parameters
[2] # Flattened if < 1.0.
70 self
.radius
= parameters
[3] # Exactly what it says.
71 self
.square
= parameters
[4] # Controls amount of superellipse in X/Y plane.
72 self
.squarez
= parameters
[5] # Controls amount of superellipse in Z dimension.
73 self
.cart
= parameters
[6] # Cuts out sphericalisation step.
74 self
.shape
= parameters
[7] # Full sphere, dome, flatbase.
75 self
.baselevel
= parameters
[8] # Where the base is cut on a flatbase dome.
76 self
.faceshape
= parameters
[9] # Triangular, hexagonal, tri-hex.
77 self
.dualflag
= parameters
[10]
78 self
.rotxy
= parameters
[11]
79 self
.rotz
= parameters
[12]
80 self
.klass
= parameters
[13]
81 self
.sform
= parameters
[14]
82 self
.super = 0 # Toggles superellipse.
83 if self
.square
!= 2.0 or self
.squarez
!= 2.0:
85 self
.odd
= 0 # Is the frequency odd. It matters for dome building.
86 if self
.frequency
% 2 != 0:
89 def makegeodesic(self
):
90 self
.vertedgefacedata() # PKHG only a pass 13okt11
91 self
.config() # Generate all the configuration information.
94 if self
.faceshape
== 1:
95 self
.hexify() # Hexagonal faces
96 elif self
.faceshape
== 2:
97 self
.starify() # Hex and Triangle faces
101 self
.sphericalize() # Convert x,y,z positions into spherical u,v.
102 self
.sphere2cartesian() # Convert spherical uv back into cartesian x,y,z for final shape.
103 for i
in range(len(self
.verts
)):
104 self
.verts
[i
].index
= i
105 for edg
in self
.edges
:
108 def vertedgefacedata(self
):
112 for i
in range(len(self
.vertskeleton
)):
113 self
.vertskeleton
[i
].index
= i
114 for edges
in self
.edgeskeleton
:
115 s
= skeletonrow(self
.frequency
, edges
, 0, self
) # self a geodesic
116 self
.skeleton
.append(s
)
117 for i
in range(len(self
.verts
)):
118 self
.verts
[i
].index
= i
119 for i
in range(len(self
.panelpoints
)):
120 a
= self
.vertsdone
[self
.panelpoints
[i
][0]][1]
121 b
= self
.vertsdone
[self
.panelpoints
[i
][1]][1]
122 c
= self
.vertsdone
[self
.panelpoints
[i
][2]][1]
123 panpoints
= [self
.verts
[a
],
126 panedges
= [self
.skeleton
[self
.paneledges
[i
][0]],
127 self
.skeleton
[self
.paneledges
[i
][1]],
128 self
.skeleton
[self
.paneledges
[i
][2]]]
130 for flag
in self
.reversepanel
:
133 p
= panel(panpoints
, panedges
, reverseflag
, self
)
135 def sphericalize(self
):
138 for vert
in(self
.verts
):
144 u
= self
.usphericalise(x
, y
, z
)
145 v
= self
.vsphericalise(x
, y
, z
)
146 self
.sphericalverts
.append([u
, v
])
148 def sphere2cartesian(self
):
149 for i
in range(len(self
.verts
)):
152 x
= self
.verts
[i
].vector
.x
* self
.radius
* self
.eccentricity
153 y
= self
.verts
[i
].vector
.y
* self
.radius
154 z
= self
.verts
[i
].vector
.z
* self
.radius
* self
.squish
156 u
= self
.sphericalverts
[i
][0]
157 v
= self
.sphericalverts
[i
][1]
158 if self
.squish
!= 1.0 or self
.eccentricity
> 1.0:
159 scalez
= 1 / self
.squish
160 v
= self
.ellipsecomp(scalez
, v
)
161 u
= self
.ellipsecomp(self
.eccentricity
, u
)
163 r1
= self
.superell(self
.square
, u
, self
.rotxy
)
164 r2
= self
.superell(self
.squarez
, v
, self
.rotz
)
170 r1
= r1
* self
.superform(self
.sform
[0], self
.sform
[1],
171 self
.sform
[2], self
.sform
[3],
172 self
.sform
[14] + u
, self
.sform
[4],
173 self
.sform
[5], self
.sform
[16] * v
)
175 r2
= r2
* self
.superform(self
.sform
[6], self
.sform
[7],
176 self
.sform
[8], self
.sform
[9],
177 self
.sform
[15] + v
, self
.sform
[10],
178 self
.sform
[11], self
.sform
[17] * v
)
179 x
, y
, z
= self
.cartesian(u
, v
, r1
, r2
)
181 self
.verts
[i
] = vertex((x
, y
, z
))
183 def usphericalise(self
, x
, y
, z
):
197 if x
< 0.0 and y
< 0.0:
198 theta
= theta
+ self
.a180
199 elif x
< 0.0 and y
> 0.0:
200 theta
= theta
+ self
.a180
204 def vsphericalise(self
, x
, y
, z
):
208 rho
= sqrt(x
** 2 + y
** 2 + z
** 2)
213 def ellipsecomp(self
, efactor
, theta
):
214 if theta
== self
.a90
:
216 elif theta
== self
.a270
:
219 result
= atan(tan(theta
) / efactor
**0.5)
222 y
= self
.a180
+ result
223 if fabs(x
- theta
) <= fabs(y
- theta
):
228 x
= self
.a180
+ result
231 if fabs(x
- theta
) <= fabs(y
- theta
):
237 def cutbasecomp(self
):
240 def cartesian(self
, u
, v
, r1
, r2
):
241 x
= r1
* cos(u
) * r2
* sin(v
) * self
.radius
* self
.eccentricity
242 y
= r1
* sin(u
) * r2
* sin(v
) * self
.radius
243 z
= r2
* cos(v
) * self
.radius
* self
.squish
248 def __init__(self
, count
, anchor
, leftindex
, rightindex
, stepvector
, endflag
, parentgeo
):
251 # Make a row of evenly spaced points.
252 for i
in range(count
+ 1):
254 self
.points
.append(leftindex
)
255 elif i
== count
and not endflag
:
256 self
.points
.append(rightindex
)
257 else: # PKHG Vectors added!
258 newpoint
= anchor
+ (stepvector
* i
)
259 vertcount
= len(parentgeo
.verts
)
260 self
.points
.append(vertcount
)
261 newpoint
.index
= vertcount
262 parentgeo
.verts
.append(newpoint
)
263 for i
in range(count
):
264 a
= parentgeo
.verts
[self
.points
[i
]]
265 b
= parentgeo
.verts
[self
.points
[i
+ 1]]
267 self
.edges
.append(len(parentgeo
.edges
))
268 parentgeo
.edges
.append(line
)
272 def __init__(self
, count
, skeletonedge
, shortflag
, parentgeo
):
275 self
.vect
= skeletonedge
.vect
276 self
.step
= skeletonedge
.vect
/ float(count
)
277 # Make a row of evenly spaced points.
278 for i
in range(count
+ 1):
279 vert1
= skeletonedge
.a
280 vert2
= skeletonedge
.b
282 if parentgeo
.vertsdone
[vert1
.index
][0]:
283 self
.points
.append(parentgeo
.vertsdone
[vert1
.index
][1])
285 newpoint
= vertex(vert1
.vector
)
286 vertcount
= len(parentgeo
.verts
)
287 self
.points
.append(vertcount
)
288 newpoint
.index
= vertcount
289 parentgeo
.vertsdone
[vert1
.index
] = [1, vertcount
]
290 parentgeo
.verts
.append(newpoint
)
293 if parentgeo
.vertsdone
[vert2
.index
][0]:
294 self
.points
.append(parentgeo
.vertsdone
[vert2
.index
][1])
296 newpoint
= vertex(vert2
.vector
)
297 vertcount
= len(parentgeo
.verts
)
298 self
.points
.append(vertcount
)
299 newpoint
.index
= vertcount
300 parentgeo
.vertsdone
[vert2
.index
] = [1, vertcount
]
301 parentgeo
.verts
.append(newpoint
)
303 newpoint
= vertex(vert1
.vector
+ (self
.step
* i
)) # must be a vertex!
304 vertcount
= len(parentgeo
.verts
)
305 self
.points
.append(vertcount
)
306 newpoint
.index
= vertcount
307 parentgeo
.verts
.append(newpoint
)
308 for i
in range(count
):
309 a
= parentgeo
.verts
[self
.points
[i
]]
310 b
= parentgeo
.verts
[self
.points
[i
+ 1]]
312 self
.edges
.append(len(parentgeo
.edges
))
313 parentgeo
.edges
.append(line
)
317 def __init__(self
, upper
, lower
, reverseflag
, parentgeo
, finish
):
318 for i
in range(finish
):
319 a
, b
, c
= upper
.points
[i
], lower
.points
[i
+ 1], lower
.points
[i
]
321 upface
= face([parentgeo
.verts
[a
], parentgeo
.verts
[c
], parentgeo
.verts
[b
]])
323 upface
= face([parentgeo
.verts
[a
], parentgeo
.verts
[b
], parentgeo
.verts
[c
]])
324 parentgeo
.faces
.append(upface
)
328 d
= upper
.points
[i
+ 1]
330 downface
= face([parentgeo
.verts
[b
], parentgeo
.verts
[d
], parentgeo
.verts
[a
]])
332 downface
= face([parentgeo
.verts
[b
], parentgeo
.verts
[a
], parentgeo
.verts
[d
]])
333 line
= edge(parentgeo
.verts
[a
], parentgeo
.verts
[b
])
334 line2
= edge(parentgeo
.verts
[d
], parentgeo
.verts
[b
])
335 parentgeo
.faces
.append(downface
)
336 parentgeo
.edges
.append(line
)
337 parentgeo
.edges
.append(line2
)
341 def __init__(self
, points
, edges
, reverseflag
, parentgeo
):
342 self
.cardinal
= points
[0]
343 self
.leftv
= points
[1]
344 self
.rightv
= points
[2]
345 self
.leftedge
= edges
[0]
346 self
.rightedge
= edges
[1]
347 self
.baseedge
= edges
[2]
349 self
.orient(parentgeo
, edges
)
350 self
.createrows(parentgeo
)
351 self
.createfaces(parentgeo
, reverseflag
)
353 def orient(self
, parentgeo
, edges
):
354 if self
.leftedge
.points
[0] != self
.cardinal
.index
:
355 self
.leftedge
.points
.reverse()
356 self
.leftedge
.vect
.negative()
358 if self
.rightedge
.points
[0] != self
.cardinal
.index
:
359 self
.rightedge
.points
.reverse()
360 self
.rightedge
.vect
.negative()
362 if self
.baseedge
.points
[0] != self
.leftv
.index
:
364 self
.baseedge
.points
.reverse()
365 self
.baseedge
.vect
.negative()
367 def createrows(self
, parentgeo
):
368 for i
in range(len(self
.leftedge
.points
)):
369 if i
== parentgeo
.frequency
:
370 newrow
= self
.baseedge
372 newrow
= edgerow(i
, parentgeo
.verts
[self
.leftedge
.points
[i
]], self
.leftedge
.points
[i
],
373 self
.rightedge
.points
[i
], self
.baseedge
.step
, 0, parentgeo
)
374 self
.rows
.append(newrow
)
376 def createfaces(self
, parentgeo
, reverseflag
):
377 for i
in range(len(self
.leftedge
.points
) - 1):
378 facefill(self
.rows
[i
], self
.rows
[i
+ 1], reverseflag
, parentgeo
, len(self
.rows
[i
].points
))
381 # for point on top? YES!
382 class tetrahedron(geodesic
, mesh
):
383 def __init__(self
, parameter
):
384 geodesic
.__init
__(mesh
)
385 geodesic
.setparameters(self
, parameter
)
386 self
.set_vert_edge_skeleons()
388 def set_vert_edge_skeleons(self
):
389 self
.vertskeleton
= [vertex((0.0, 0.0, 1.73205080757)),
390 vertex((0.0, -1.63299316185, -0.577350269185)),
391 vertex((1.41421356237, 0.816496580927, -0.57735026919)),
392 vertex((-1.41421356237, 0.816496580927, -0.57735026919))]
394 self
.edgeskeleton
= [edge(self
.vertskeleton
[0], self
.vertskeleton
[1]),
395 edge(self
.vertskeleton
[0], self
.vertskeleton
[2]),
396 edge(self
.vertskeleton
[0], self
.vertskeleton
[3]),
397 edge(self
.vertskeleton
[1], self
.vertskeleton
[2]),
398 edge(self
.vertskeleton
[2], self
.vertskeleton
[3]),
399 edge(self
.vertskeleton
[1], self
.vertskeleton
[3])]
401 self
.panelpoints
= [[0, 1, 2], [0, 2, 3], [0, 1, 3], [1, 2, 3]]
402 self
.paneledges
= [[0, 1, 3], [1, 2, 4], [0, 2, 5], [3, 5, 4]]
403 self
.reversepanel
= [2, 3]
405 self
.vertsdone
= [[0, 0]] * len(self
.vertskeleton
)
408 # for edge on top? YES
409 class tetraedge(geodesic
):
410 def __init__(self
, parameter
):
411 geodesic
.__init
__(mesh
)
412 geodesic
.setparameters(self
, parameter
)
413 self
.set_vert_edge_skeleons()
415 def set_vert_edge_skeleons(self
):
416 self
.vertskeleton
= [vertex((0.0, -1.41421356237, 1.0)),
417 vertex((0.0, 1.41421356237, 1.0)),
418 vertex((1.41421356237, 0.0, -1.0)),
419 vertex((-1.41421356237, 0.0, -1.0))]
421 self
.edgeskeleton
= [edge(self
.vertskeleton
[0], self
.vertskeleton
[1]),
422 edge(self
.vertskeleton
[0], self
.vertskeleton
[2]),
423 edge(self
.vertskeleton
[0], self
.vertskeleton
[3]),
424 edge(self
.vertskeleton
[1], self
.vertskeleton
[3]),
425 edge(self
.vertskeleton
[1], self
.vertskeleton
[2]),
426 edge(self
.vertskeleton
[2], self
.vertskeleton
[3])]
428 for i
in range(len(self
.vertskeleton
)):
429 self
.vertskeleton
[i
].index
= i
431 self
.panelpoints
= [[0, 1, 2], [1, 2, 3], [0, 1, 3], [0, 2, 3]]
432 self
.paneledges
= [[0, 1, 4], [4, 3, 5], [0, 2, 3], [1, 2, 5]]
433 self
.reversepanel
= [0, 3]
435 self
.vertsdone
= [[0, 0]] * len(self
.vertskeleton
)
438 # for face on top? YES
439 class tetraface(geodesic
):
440 def __init__(self
, parameter
):
441 geodesic
.__init
__(mesh
)
442 geodesic
.setparameters(self
, parameter
)
443 self
.set_vert_edge_skeleons()
445 def set_vert_edge_skeleons(self
):
446 self
.vertskeleton
= [vertex((-1.41421356237, -0.816496580927, 0.57735026919)),
447 vertex((1.41421356237, -0.816496580927, 0.57735026919)),
448 vertex((0.0, 1.63299316185, 0.577350269185)),
449 vertex((0.0, 0.0, -1.73205080757))]
451 self
.edgeskeleton
= [edge(self
.vertskeleton
[0], self
.vertskeleton
[1]),
452 edge(self
.vertskeleton
[2], self
.vertskeleton
[1]),
453 edge(self
.vertskeleton
[2], self
.vertskeleton
[0]),
454 edge(self
.vertskeleton
[0], self
.vertskeleton
[3]),
455 edge(self
.vertskeleton
[1], self
.vertskeleton
[3]),
456 edge(self
.vertskeleton
[2], self
.vertskeleton
[3])
458 self
.panelpoints
= [[2, 0, 1], [0, 1, 3], [2, 1, 3], [2, 0, 3]]
460 self
.paneledges
= [[2, 1, 0], [0, 3, 4], [1, 5, 4], [2, 5, 3]]
461 self
.reversepanel
= [1, 3]
463 self
.vertsdone
= [[0, 0]] * len(self
.vertskeleton
)
466 class octahedron(geodesic
):
467 def __init__(self
, parameter
):
468 geodesic
.__init
__(mesh
)
469 geodesic
.setparameters(self
, parameter
)
470 self
.set_vert_edge_skeleons()
472 def set_vert_edge_skeleons(self
):
473 self
.vertskeleton
= [vertex((0.0, 0.0, 1.0)),
474 vertex((0.0, 1.0, 0.0)),
475 vertex((-1.0, 0.0, 0.0)),
476 vertex((0.0, -1.0, 0.0)),
477 vertex((1.0, 0.0, 0.0)),
478 vertex((0.0, 0.0, -1.0))]
480 for i
in range(len(self
.vertskeleton
)):
481 self
.vertskeleton
[i
].index
= i
482 self
.edgeskeleton
= [edge(self
.vertskeleton
[0], self
.vertskeleton
[1]),
483 edge(self
.vertskeleton
[0], self
.vertskeleton
[2]),
484 edge(self
.vertskeleton
[0], self
.vertskeleton
[3]),
485 edge(self
.vertskeleton
[0], self
.vertskeleton
[4]),
486 edge(self
.vertskeleton
[1], self
.vertskeleton
[2]),
487 edge(self
.vertskeleton
[2], self
.vertskeleton
[3]),
488 edge(self
.vertskeleton
[3], self
.vertskeleton
[4]),
489 edge(self
.vertskeleton
[4], self
.vertskeleton
[1]),
490 edge(self
.vertskeleton
[1], self
.vertskeleton
[5]),
491 edge(self
.vertskeleton
[2], self
.vertskeleton
[5]),
492 edge(self
.vertskeleton
[3], self
.vertskeleton
[5]),
493 edge(self
.vertskeleton
[4], self
.vertskeleton
[5])]
495 self
.panelpoints
= [[0, 1, 2], [0, 2, 3], [0, 3, 4], [0, 4, 1], [1, 2, 5],
496 [2, 3, 5], [3, 4, 5], [4, 1, 5]]
497 self
.paneledges
= [[0, 1, 4], [1, 2, 5], [2, 3, 6], [3, 0, 7], [4, 8, 9],
498 [5, 9, 10], [6, 10, 11], [7, 11, 8]]
499 self
.reversepanel
= [4, 5, 6, 7]
501 self
.vertsdone
= [[0, 0]] * len(self
.vertskeleton
)
504 class octaedge(geodesic
):
505 def __init__(self
, parameter
):
506 geodesic
.__init
__(mesh
)
507 geodesic
.setparameters(self
, parameter
)
508 self
.set_vert_edge_skeleons()
510 def set_vert_edge_skeleons(self
):
511 self
.vertskeleton
= [vertex((0.0, -0.707106781187, 0.707106781187)),
512 vertex((0.0, 0.707106781187, 0.707106781187)),
513 vertex((1.0, 0.0, 0.0)),
514 vertex((-1.0, 0.0, 0.0)),
515 vertex((0.0, -0.707106781187, -0.707106781187)),
516 vertex((0.0, 0.707106781187, -0.707106781187))]
518 self
.edgeskeleton
= [edge(self
.vertskeleton
[0], self
.vertskeleton
[1]),
519 edge(self
.vertskeleton
[0], self
.vertskeleton
[4]),
520 edge(self
.vertskeleton
[0], self
.vertskeleton
[2]),
521 edge(self
.vertskeleton
[1], self
.vertskeleton
[2]),
522 edge(self
.vertskeleton
[1], self
.vertskeleton
[5]),
523 edge(self
.vertskeleton
[1], self
.vertskeleton
[3]),
524 edge(self
.vertskeleton
[0], self
.vertskeleton
[3]),
525 edge(self
.vertskeleton
[2], self
.vertskeleton
[4]),
526 edge(self
.vertskeleton
[2], self
.vertskeleton
[5]),
527 edge(self
.vertskeleton
[3], self
.vertskeleton
[5]),
528 edge(self
.vertskeleton
[3], self
.vertskeleton
[4]),
529 edge(self
.vertskeleton
[4], self
.vertskeleton
[5])]
531 self
.panelpoints
= [[0, 1, 2], [0, 1, 3], [0, 2, 4], [1, 2, 5], [1, 3, 5],
532 [0, 3, 4], [2, 4, 5], [3, 4, 5]]
533 self
.paneledges
= [[0, 2, 3], [0, 6, 5], [2, 1, 7], [3, 4, 8], [5, 4, 9],
534 [6, 1, 10], [7, 8, 11], [10, 9, 11]]
536 self
.reversepanel
= [0, 2, 4, 7]
538 self
.vertsdone
= [[0, 0]] * len(self
.vertskeleton
)
541 class octaface(geodesic
):
542 def __init__(self
, parameter
):
543 geodesic
.__init
__(mesh
)
544 geodesic
.setparameters(self
, parameter
)
545 self
.set_vert_edge_skeleons()
547 def set_vert_edge_skeleons(self
):
548 self
.vertskeleton
= [vertex((0.408248458663, -0.707106781187, 0.577350150255)),
549 vertex((0.408248458663, 0.707106781187, 0.577350150255)),
550 vertex((-0.816496412728, 0.0, 0.577350507059)),
551 vertex((-0.408248458663, -0.707106781187, -0.577350150255)),
552 vertex((0.816496412728, 0.0, -0.577350507059)),
553 vertex((-0.408248458663, 0.707106781187, -0.577350150255))]
555 self
.edgeskeleton
= [edge(self
.vertskeleton
[0], self
.vertskeleton
[1]),
556 edge(self
.vertskeleton
[2], self
.vertskeleton
[1]),
557 edge(self
.vertskeleton
[2], self
.vertskeleton
[0]),
558 edge(self
.vertskeleton
[0], self
.vertskeleton
[3]),
559 edge(self
.vertskeleton
[0], self
.vertskeleton
[4]),
560 edge(self
.vertskeleton
[1], self
.vertskeleton
[4]),
561 edge(self
.vertskeleton
[1], self
.vertskeleton
[5]),
562 edge(self
.vertskeleton
[2], self
.vertskeleton
[5]),
563 edge(self
.vertskeleton
[2], self
.vertskeleton
[3]),
564 edge(self
.vertskeleton
[3], self
.vertskeleton
[4]),
565 edge(self
.vertskeleton
[4], self
.vertskeleton
[5]),
566 edge(self
.vertskeleton
[3], self
.vertskeleton
[5])]
568 self
.panelpoints
= [[2, 0, 1], [0, 3, 4], [0, 1, 4], [1, 4, 5],
569 [2, 1, 5], [2, 3, 5], [2, 0, 3], [3, 4, 5]]
571 self
.paneledges
= [[2, 1, 0], [3, 4, 9], [0, 4, 5], [5, 6, 10],
572 [1, 7, 6], [8, 7, 11], [2, 8, 3], [9, 11, 10]]
574 self
.reversepanel
= [2, 5, 6, 7]
576 self
.vertsdone
= [[0, 0]] * len(self
.vertskeleton
)
579 class icosahedron(geodesic
):
580 def __init__(self
, parameter
):
581 geodesic
.__init
__(mesh
)
582 geodesic
.setparameters(self
, parameter
)
583 self
.set_vert_edge_skeleons()
585 def set_vert_edge_skeleons(self
):
586 self
.vertskeleton
= [vertex((0.0, 0.0, 0.587785252292)),
587 vertex((0.0, -0.525731096637, 0.262865587024)),
588 vertex((0.5, -0.162459832634, 0.262865565628)),
589 vertex((0.309016994375, 0.425325419658, 0.262865531009)),
590 vertex((-0.309016994375, 0.425325419658, 0.262865531009)),
591 vertex((-0.5, -0.162459832634, 0.262865565628)),
592 vertex((0.309016994375, -0.425325419658, -0.262865531009)),
593 vertex((0.5, 0.162459832634, -0.262865565628)),
594 vertex((0.0, 0.525731096637, -0.262865587024)),
595 vertex((-0.5, 0.162459832634, -0.262865565628)),
596 vertex((-0.309016994375, -0.425325419658, -0.262865531009)),
597 vertex((0.0, 0.0, -0.587785252292))]
599 self
.edgeskeleton
= [edge(self
.vertskeleton
[0], self
.vertskeleton
[1]),
600 edge(self
.vertskeleton
[0], self
.vertskeleton
[2]),
601 edge(self
.vertskeleton
[0], self
.vertskeleton
[3]),
602 edge(self
.vertskeleton
[0], self
.vertskeleton
[4]),
603 edge(self
.vertskeleton
[0], self
.vertskeleton
[5]),
604 edge(self
.vertskeleton
[1], self
.vertskeleton
[2]),
605 edge(self
.vertskeleton
[2], self
.vertskeleton
[3]),
606 edge(self
.vertskeleton
[3], self
.vertskeleton
[4]),
607 edge(self
.vertskeleton
[4], self
.vertskeleton
[5]),
608 edge(self
.vertskeleton
[5], self
.vertskeleton
[1]),
609 edge(self
.vertskeleton
[1], self
.vertskeleton
[6]),
610 edge(self
.vertskeleton
[2], self
.vertskeleton
[6]),
611 edge(self
.vertskeleton
[2], self
.vertskeleton
[7]),
612 edge(self
.vertskeleton
[3], self
.vertskeleton
[7]),
613 edge(self
.vertskeleton
[3], self
.vertskeleton
[8]),
614 edge(self
.vertskeleton
[4], self
.vertskeleton
[8]),
615 edge(self
.vertskeleton
[4], self
.vertskeleton
[9]),
616 edge(self
.vertskeleton
[5], self
.vertskeleton
[9]),
617 edge(self
.vertskeleton
[5], self
.vertskeleton
[10]),
618 edge(self
.vertskeleton
[1], self
.vertskeleton
[10]),
619 edge(self
.vertskeleton
[6], self
.vertskeleton
[7]),
620 edge(self
.vertskeleton
[7], self
.vertskeleton
[8]),
621 edge(self
.vertskeleton
[8], self
.vertskeleton
[9]),
622 edge(self
.vertskeleton
[9], self
.vertskeleton
[10]),
623 edge(self
.vertskeleton
[10], self
.vertskeleton
[6]),
624 edge(self
.vertskeleton
[6], self
.vertskeleton
[11]),
625 edge(self
.vertskeleton
[7], self
.vertskeleton
[11]),
626 edge(self
.vertskeleton
[8], self
.vertskeleton
[11]),
627 edge(self
.vertskeleton
[9], self
.vertskeleton
[11]),
628 edge(self
.vertskeleton
[10], self
.vertskeleton
[11])]
630 self
.panelpoints
= [[0, 1, 2], [0, 2, 3], [0, 3, 4], [0, 4, 5], [0, 5, 1], [1, 2, 6],
631 [2, 6, 7], [2, 3, 7], [3, 7, 8], [3, 4, 8], [4, 8, 9], [4, 5, 9],
632 [5, 9, 10], [5, 1, 10], [1, 10, 6], [6, 7, 11], [7, 8, 11],
633 [8, 9, 11], [9, 10, 11], [10, 6, 11]]
635 self
.paneledges
= [[0, 1, 5], [1, 2, 6], [2, 3, 7], [3, 4, 8], [4, 0, 9], [5, 10, 11],
636 [11, 12, 20], [6, 12, 13], [13, 14, 21], [7, 14, 15], [15, 16, 22],
637 [8, 16, 17], [17, 18, 23], [9, 18, 19], [19, 10, 24], [20, 25, 26],
638 [21, 26, 27], [22, 27, 28], [23, 28, 29], [24, 29, 25]]
640 self
.reversepanel
= [5, 7, 9, 11, 13, 15, 16, 17, 18, 19]
642 self
.vertsdone
= [[0, 0]] * len(self
.vertskeleton
)
645 class icoedge(geodesic
):
646 def __init__(self
, parameter
):
647 geodesic
.__init
__(mesh
)
648 geodesic
.setparameters(self
, parameter
)
649 self
.set_vert_edge_skeleons()
651 def set_vert_edge_skeleons(self
):
652 self
.vertskeleton
= [vertex((0, 0.309016994375, 0.5)),
653 vertex((0, -0.309016994375, 0.5)),
654 vertex((-0.5, 0, 0.309016994375)),
655 vertex((0.5, 0, 0.309016994375)),
656 vertex((-0.309016994375, -0.5, 0)),
657 vertex((0.309016994375, -0.5, 0)),
658 vertex((0.309016994375, 0.5, 0)),
659 vertex((-0.309016994375, 0.5, 0)),
660 vertex((-0.5, 0, -0.309016994375)),
661 vertex((0.5, 0, -0.309016994375)),
662 vertex((0, 0.309016994375, -0.5)),
663 vertex((0, -0.309016994375, -0.5))]
665 self
.edgeskeleton
= [edge(self
.vertskeleton
[0], self
.vertskeleton
[1]),
666 edge(self
.vertskeleton
[0], self
.vertskeleton
[7]),
667 edge(self
.vertskeleton
[0], self
.vertskeleton
[2]),
668 edge(self
.vertskeleton
[1], self
.vertskeleton
[2]),
669 edge(self
.vertskeleton
[1], self
.vertskeleton
[4]),
670 edge(self
.vertskeleton
[1], self
.vertskeleton
[5]),
671 edge(self
.vertskeleton
[1], self
.vertskeleton
[3]),
672 edge(self
.vertskeleton
[0], self
.vertskeleton
[3]),
673 edge(self
.vertskeleton
[0], self
.vertskeleton
[6]),
674 edge(self
.vertskeleton
[2], self
.vertskeleton
[7]),
675 edge(self
.vertskeleton
[2], self
.vertskeleton
[8]),
676 edge(self
.vertskeleton
[2], self
.vertskeleton
[4]),
677 edge(self
.vertskeleton
[4], self
.vertskeleton
[5]),
678 edge(self
.vertskeleton
[3], self
.vertskeleton
[5]),
679 edge(self
.vertskeleton
[3], self
.vertskeleton
[9]),
680 edge(self
.vertskeleton
[3], self
.vertskeleton
[6]),
681 edge(self
.vertskeleton
[6], self
.vertskeleton
[7]),
682 edge(self
.vertskeleton
[7], self
.vertskeleton
[10]),
683 edge(self
.vertskeleton
[7], self
.vertskeleton
[8]),
684 edge(self
.vertskeleton
[4], self
.vertskeleton
[8]),
685 edge(self
.vertskeleton
[4], self
.vertskeleton
[11]),
686 edge(self
.vertskeleton
[5], self
.vertskeleton
[11]),
687 edge(self
.vertskeleton
[5], self
.vertskeleton
[9]),
688 edge(self
.vertskeleton
[6], self
.vertskeleton
[9]),
689 edge(self
.vertskeleton
[6], self
.vertskeleton
[10]),
690 edge(self
.vertskeleton
[8], self
.vertskeleton
[10]),
691 edge(self
.vertskeleton
[8], self
.vertskeleton
[11]),
692 edge(self
.vertskeleton
[9], self
.vertskeleton
[11]),
693 edge(self
.vertskeleton
[9], self
.vertskeleton
[10]),
694 edge(self
.vertskeleton
[10], self
.vertskeleton
[11])]
696 self
.panelpoints
= [[0, 1, 2], [0, 1, 3], [0, 2, 7], [1, 2, 4], [1, 4, 5],
697 [1, 3, 5], [0, 3, 6], [0, 6, 7], [2, 7, 8], [2, 4, 8],
698 [3, 5, 9], [3, 6, 9], [7, 8, 10], [4, 8, 11], [4, 5, 11],
699 [5, 9, 11], [6, 9, 10], [6, 7, 10], [8, 10, 11], [9, 10, 11]]
701 self
.paneledges
= [[0, 2, 3], [0, 7, 6], [2, 1, 9], [3, 4, 11], [4, 5, 12], [6, 5, 13],
702 [7, 8, 15], [8, 1, 16], [9, 10, 18], [11, 10, 19], [13, 14, 22],
703 [15, 14, 23], [18, 17, 25], [19, 20, 26], [12, 20, 21], [22, 21, 27],
704 [23, 24, 28], [16, 24, 17], [25, 26, 29], [28, 27, 29]]
706 self
.reversepanel
= [0, 2, 5, 9, 11, 12, 14, 15, 17, 19]
708 self
.vertsdone
= [[0, 0]] * len(self
.vertskeleton
)
711 class icoface(geodesic
):
712 def __init__(self
, parameter
):
713 geodesic
.__init
__(mesh
)
714 geodesic
.setparameters(self
, parameter
)
715 self
.set_vert_edge_skeleons()
717 def set_vert_edge_skeleons(self
):
718 self
.vertskeleton
= [vertex((-0.17841104489, 0.309016994375, 0.46708617948)),
719 vertex((-0.17841104489, -0.309016994375, 0.46708617948)),
720 vertex((0.35682208977, 0.0, 0.467086179484)),
721 vertex((-0.57735026919, 0.0, 0.110264089705)),
722 vertex((-0.288675134594, -0.5, -0.11026408971)),
723 vertex((0.288675134594, -0.5, 0.11026408971)),
724 vertex((0.57735026919, 0.0, -0.110264089705)),
725 vertex((0.288675134594, 0.5, 0.11026408971)),
726 vertex((-0.288675134594, 0.5, -0.11026408971)),
727 vertex((-0.35682208977, 0.0, -0.467086179484)),
728 vertex((0.17841104489, -0.309016994375, -0.46708617948)),
729 vertex((0.17841104489, 0.309016994375, -0.46708617948))]
731 self
.edgeskeleton
= [edge(self
.vertskeleton
[0], self
.vertskeleton
[1]),
732 edge(self
.vertskeleton
[2], self
.vertskeleton
[1]),
733 edge(self
.vertskeleton
[2], self
.vertskeleton
[0]),
734 edge(self
.vertskeleton
[0], self
.vertskeleton
[3]),
735 edge(self
.vertskeleton
[1], self
.vertskeleton
[3]),
736 edge(self
.vertskeleton
[1], self
.vertskeleton
[4]),
737 edge(self
.vertskeleton
[1], self
.vertskeleton
[5]),
738 edge(self
.vertskeleton
[2], self
.vertskeleton
[5]),
739 edge(self
.vertskeleton
[2], self
.vertskeleton
[6]),
740 edge(self
.vertskeleton
[2], self
.vertskeleton
[7]),
741 edge(self
.vertskeleton
[0], self
.vertskeleton
[7]),
742 edge(self
.vertskeleton
[0], self
.vertskeleton
[8]),
743 edge(self
.vertskeleton
[3], self
.vertskeleton
[9]),
744 edge(self
.vertskeleton
[3], self
.vertskeleton
[4]),
745 edge(self
.vertskeleton
[5], self
.vertskeleton
[4]),
746 edge(self
.vertskeleton
[5], self
.vertskeleton
[10]),
747 edge(self
.vertskeleton
[5], self
.vertskeleton
[6]),
748 edge(self
.vertskeleton
[7], self
.vertskeleton
[6]),
749 edge(self
.vertskeleton
[7], self
.vertskeleton
[11]),
750 edge(self
.vertskeleton
[7], self
.vertskeleton
[8]),
751 edge(self
.vertskeleton
[3], self
.vertskeleton
[8]),
752 edge(self
.vertskeleton
[4], self
.vertskeleton
[9]),
753 edge(self
.vertskeleton
[4], self
.vertskeleton
[10]),
754 edge(self
.vertskeleton
[6], self
.vertskeleton
[10]),
755 edge(self
.vertskeleton
[6], self
.vertskeleton
[11]),
756 edge(self
.vertskeleton
[8], self
.vertskeleton
[11]),
757 edge(self
.vertskeleton
[8], self
.vertskeleton
[9]),
758 edge(self
.vertskeleton
[9], self
.vertskeleton
[10]),
759 edge(self
.vertskeleton
[11], self
.vertskeleton
[10]),
760 edge(self
.vertskeleton
[11], self
.vertskeleton
[9])]
762 self
.panelpoints
= [[2, 0, 1], [0, 1, 3], [2, 1, 5], [2, 0, 7], [1, 3, 4], [1, 5, 4],
763 [2, 5, 6], [2, 7, 6], [0, 7, 8], [0, 3, 8], [3, 4, 9], [5, 4, 10],
764 [5, 6, 10], [7, 6, 11], [7, 8, 11], [3, 8, 9], [4, 9, 10],
765 [6, 11, 10], [8, 11, 9], [11, 9, 10]]
767 self
.paneledges
= [[2, 1, 0], [0, 3, 4], [1, 7, 6], [2, 9, 10], [4, 5, 13], [6, 5, 14],
768 [7, 8, 16], [9, 8, 17], [10, 11, 19], [3, 11, 20], [13, 12, 21],
769 [14, 15, 22], [16, 15, 23], [17, 18, 24], [19, 18, 25], [20, 12, 26],
770 [21, 22, 27], [24, 23, 28], [25, 26, 29], [29, 28, 27]]
772 self
.reversepanel
= [1, 3, 5, 7, 9, 10, 12, 14, 17, 19]
774 self
.vertsdone
= [[0, 0]] * len(self
.vertskeleton
)
777 # PKHG TODO this does not work yet ...
778 def creategeo(geo
, polytype
, orientation
, parameters
):
780 if polytype
== 'Tetrahedron':
781 if orientation
== 'PointUp':
782 my_tetrahedron
= tetrahedron(geodesic
)
783 my_tetrahedron
.set_vert_edge_skeleons()
784 my_tetrahedron
.config()
785 check_contains(my_tetrahedron
, "my_tetra", True)
787 elif orientation
== 'EdgeUp':
788 geo
= tetraedge(parameters
)
789 else: # orientation == 2:
790 geo
= tetraface(parameters
)
791 elif polytype
== 'Octahedron':
792 if orientation
== 'PointUp':
793 geo
= octahedron(parameters
)
794 elif orientation
== 'EdgeUp':
795 geo
= octaedge(parameters
)
796 else: # if orientation == 2:
797 geo
= octaface(parameters
)
798 elif polytype
== 'Icosahedron':
799 if orientation
== 'PointUp':
800 geo
= icosahedron(parameters
)
801 elif orientation
== 'EdgeUp':
802 geo
= icoedge(parameters
)
803 else: # if orientation == 2:
804 geo
= icoface(parameters
)