src/GEOM_SWIG/GEOM_shellSolid.py

   1 #  Copyright (C) 2005  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
   2 #  CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
   3 #
   4 #  This library is free software; you can redistribute it and/or
   5 #  modify it under the terms of the GNU Lesser General Public
   6 #  License as published by the Free Software Foundation; either
   7 #  version 2.1 of the License.
   8 #
   9 #  This library is distributed in the hope that it will be useful,
  10 #  but WITHOUT ANY WARRANTY; without even the implied warranty of
  11 #  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  12 #  Lesser General Public License for more details.
  13 #
  14 #  You should have received a copy of the GNU Lesser General Public
  15 #  License along with this library; if not, write to the Free Software
  16 #  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  17 #
  18 # See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  19 #
  20 #==============================================================================
  21 #  Info.
  22 #  Bug (from script, bug)   : shellSolid_modified.py, PAL6611
  23 #  Modified                 : 25/11/2004
  24 #  Author                   : Kovaltchuk Alexey
  25 #  Project                  : PAL/SALOME
  26 #==============================================================================
  27
  28 # Imports
  29 # -------
  30
  31 import salome
  32 import SALOMEDS
  33
  34 import geompy
  35
  36 # Interface with geometry
  37 # -----------------------
  38
  39 #geom = salome.lcc.FindOrLoadComponent("FactoryServer", "GEOM")
  40 geomgui = salome.ImportComponentGUI("GEOM")
  41
  42 geom = geompy
  43
  44 # Interface with study
  45 # --------------------
  46
  47 def addToStudy(shape, name):
  48     i = geompy.addToStudy(shape, name)
  49     salome.sg.updateObjBrowser(0)
  50 #    geomgui.createAndDisplayGO(i)
  51     return i
  52
  53 def setName(obj):
  54     obj.SetName(salome.orb.object_to_string(obj))
  55
  56
  57 def MakeVertex(x, y, z):
  58     v = geom.MakeVertex(x, y, z)
  59     setName(v)
  60     return v
  61
  62 def MakeEdge(v1, v2):
  63     e = geom.MakeEdge(v1, v2)
  64     setName(e)
  65     return e
  66
  67 def MakeArc(f, p, t):
  68     e = geom.MakeArc(f, p, t)
  69     setName(e)
  70     return e
  71
  72 def MakeArcCenter(c, f, t):
  73     p = fkl(c, f, t)
  74     e = geom.MakeArc(f, p, t)
  75     setName(e)
  76     return e
  77
  78 # Face with 4 edges
  79 # -----------------
  80
  81 # cf. makeTiledSurfByContour
  82
  83 def MakeQuadFace(e1, e2, e3, e4):
  84     l = []
  85     l.append(e1)
  86     l.append(e2)
  87     l.append(e3)
  88     l.append(e4)
  89     w = geom.MakeWire(l)
  90     f = geom.MakeFace(w, 1)
  91     setName(f)
  92     return f
  93
  94 # cf. makeTiledSurfBetweenTwolines
  95
  96 def MakeQuadFaceEdge(e1, e2):
  97     return fkl
  98
  99 # cf. makeTiledSurfByPointsOfContour
 100 # cf. makeMesh
 101
 102 def MakeQuadFaceVertex(v1, v2, v3, v4):
 103     l = []
 104     e = MakeEdge(v1, v2)
 105     l.append(e)
 106     e = MakeEdge(v2, v3)
 107     l.append(e)
 108     e = MakeEdge(v3, v4)
 109     l.append(e)
 110     e = MakeEdge(v4, v1)
 111     l.append(e)
 112     w = geom.MakeWire(l)
 113     f = geom.MakeFace(w, 1)
 114     setName(f)
 115     return f
 116
 117 # Solid with 6 faces
 118 # ------------------
 119
 120 # cf. makePavedVolWithContour
 121 # cf. makeVolWithEnvolope
 122
 123 def MakeHexaSolid(f1, f2, f3, f4, f5, f6):
 124     l = []
 125     l.append(f1)
 126     l.append(f2)
 127     l.append(f3)
 128     l.append(f4)
 129     l.append(f5)
 130     l.append(f6)
 131     s = geom.MakeShell(l)
 132 # FKL: bad shell
 133     l = []
 134     l.append(s)
 135     s = geom.MakeSolid(l)
 136 # FKL: MakeSolid crash
 137     setName(s)
 138     return s
 139
 140 # cf. makeVolBetweenTwoSurfaces
 141
 142 def MakeHexaSolidFace(f1, f2):
 143     return fkl
 144
 145 # Face (with 4 edges) and solid (with 6 faces)
 146 # --------------------------------------------
 147
 148 # makeTiledSurfByRotationOfLine
 149
 150 def MakeRevolution(g, pt, ve, angle):
 151     """
 152     Creates a face  by rotation of an edge with an angle around an axis defined by a point and a vector or
 153     creates a solid by rotation of a  face with an angle around an axis defined by a point and a vector
 154     """
 155     axis = geom.MakeAxisStruct(pt.x, pt.y, pt.z, ve.x, ve.y, ve.z)
 156     s = geom.MakeRevolution(g, axis, angle)
 157     setName(s)
 158     return s
 159
 160 # cf. makeSurfSewingTwolines
 161
 162 def MakeSewing(ls):
 163     """
 164     Creates a face  by sewing common edges between a list of faces or
 165     Creates a solid by sewing common faces between a list of solids
 166     """
 167     s = geom.MakeSewing(ls)
 168     setName(s)
 169     return s
 170
 171 # cf. makeIntersectionBetweenTwoMeshes
 172
 173 def MakeCommon(s1, s2):
 174     """
 175     Creates the common face  between 2 faces or
 176     Creates the common solid between 2 solids
 177     """
 178     s = geom.MakeBoolean(s1, s2, 1)
 179     setName(s)
 180     return s
 181
 182 # cf. makeMerge
 183
 184 def MakeFuse(s1, s2):
 185     """
 186     Fuses 2 faces or 2 solids
 187     """
 188     s = geom.MakeBoolean(s1, s2, 3)
 189     setName(s)
 190     return s
 191
 192 # Additional
 193
 194 def MakeCut(s1, s2):
 195     """
 196     Cuts 2 faces or 2 solids
 197     """
 198     s = geom.MakeBoolean(s1, s2, 2)
 199     setName(s)
 200     return s
 201
 202 # Edge, face (with 4 edges) and solid (with 6 faces)
 203 # --------------------------------------------------
 204
 205 # cf. makeSurfByTranslation
 206 # cf. makeVolByTranslation
 207
 208 def MakePrism(s, v):
 209     """
 210     Creates a edge  defined by a vertex and along a vector
 211     Creates a face  defined by a edge   and along a vector
 212     Creates a solid defined by a face   and along a vector
 213     """
 214     r = geom.MakePrism(s, point(0, 0, 0), v)
 215     setName(r)
 216     return r
 217
 218 # cf. homothety
 219
 220 def MakeScaleTransform(s, center, factor):
 221     """
 222     Creates a homothety of a geometric object
 223     """
 224     r = geom.MakeScaleTransform(s, center, factor)
 225     setName(r)
 226     return r
 227
 228 # Vertex, edge, face (with 4 edges) and solid (with 6 faces)
 229 # ----------------------------------------------------------
 230
 231 # cf. translate
 232
 233 def MakeTranslation(s, vx, vy, vz):
 234     """
 235     Translates a vertex, an edge, a face or a solid
 236     """
 237     r = geom.MakeTranslation(s, vx, vy, vz)
 238     setName(r)
 239     return r
 240
 241 # cf. rotation
 242
 243 def MakeRotation(g, pt, ve, angle):
 244     """
 245     Creates a rotation of the geometric object with an angle around an axis defined by a point and a vector
 246     """
 247     axis = geom.MakeAxisStruct(pt.x, pt.y, pt.z, ve.x, ve.y, ve.z)
 248     s = geom.MakeRotation(g, axis, angle)
 249     setName(s)
 250     return s
 251
 252 # cf. symmetryByPoint
 253
 254 def MakeMirrorByPoint(s, point):
 255     """
 256     Creates a symmetric object by point symetry
 257     """
 258     return fkl
 259
 260 # cf. symmetryByLine
 261
 262 def MakeMirrorByLine(s, line):
 263     """
 264     Creates a symmetric object by line symetry
 265     """
 266     return fkl
 267
 268 # cf. symmetryByPlane
 269
 270 def MakeMirrorByPlane(s, pt, ve):
 271     """
 272     Creates a symmetric object by plane symetry defined by a point and a normal vector
 273     """
 274     p = plane(pt, ve)
 275     r = geom.MakeMirrorByplane(s, p)
 276     setName(r)
 277     return r
 278
 279 # Compound
 280 # --------
 281
 282 # Additional
 283
 284 def BlockMakeCompound(ls):
 285     """
 286     Creates a compound defined by a list and glue faces
 287     """
 288     c = MakeCompound(ls)
 289     r = MakeGlueFaces(c, 1.e-5)
 290     setName(r)
 291     return r
 292
 293 # Boolean operators
 294 # -----------------
 295
 296 def mergePointsNearBy(l, v):
 297     """
 298     Modifies a list of shape by replacing all vertex near a new vextex
 299     """
 300     return fkl
 301
 302 # Acces operators
 303 # ---------------
 304
 305 def getCoords(v):
 306     """
 307     Gets the coordinates of the vertex
 308     """
 309     return fkl
 310
 311 def extractBeginPoint(e):
 312     """
 313     Gets the first point of an edge
 314     """
 315     return fkl
 316
 317 def extractEndPoint(e):
 318     """
 319     Gets the second point of an edge
 320     """
 321     return fkl
 322
 323 # fkl more functionalities
 324
 325 # Addionnal
 326 # =========
 327
 328 # Check geometry
 329 # --------------
 330
 331 def CheckQuadFace(f):
 332     """
 333     Checks if a face has 4 edges
 334     """
 335     return fkl
 336
 337 def CheckHexaSolid(s):
 338     """
 339     Checks if a solid has 6 faces and each faces have 4 edges
 340     """
 341     return fkl
 342
 343 # fkl more if other kind of solid like: tetra, penta, etc
 344
 345 # point for building
 346 # ------------------
 347
 348 def point(x, y, z):
 349     p = geom.MakeVertex(x, y, z)
 350     return p
 351
 352 def plane(pt, dir):
 353     p = geom.MakePlane(pt, d, 100)
 354     return p
 355
 356 # Solid
 357 # -----
 358
 359 def MakeCylinder(center, dir, radius, height):
 360     s = geom.MakeCylinder(center, dir, radius, height)
 361     setName(s)
 362     return s
 363
 364 def MakeBox(p1, p2):
 365     s = geom.MakeBoxTwoPnt(p1, p2)
 366     setName(s)
 367     return s
 368
 369 # Compound
 370 # --------
 371
 372 ShapeType = {"COMPOUND":0, "COMPSOLID":1, "SOLID":2, "SHELL":3, "FACE":4, "WIRE":5, "EDGE":6, "VERTEX":7, "SHAPE":8}
 373
 374 def MakePartitionList(solids, tools):
 375     """
 376     Creates a list of shape by a partition of a list of solids by a list of tools
 377     """
 378     p = geom.Partition(solids, tools, [], [], ShapeType["SHAPE"]);
 379     l = geom.SubShapeAll(p, ShapeType["SHAPE"])
 380     return l
 381
 382 def MakePartition(solids, tools):
 383     """
 384     Creates a partition of a list of solids by a list of tools
 385     """
 386     l = MakePartitionList(solids, tools)
 387     c = MakeCompound(l)
 388     setName(c)
 389     return c
 390
 391 def BlockMakeMultiTranslation1D(shape, dir, step, times):
 392     m = geom.MakeMultiTranslation1D(shape, dir, step, times)
 393     c  = MakeGlueFaces(m, 1.e-5)
 394     setName(c)
 395     return c
 396
 397 def BlockMakeMultiTranslation2D(shape, dir1, step1, times1, dir2, step2, times2):
 398     m  = geom.MakeMultiTranslation2D(shape, dir1, step1, times1, dir2, step2, times2)
 399     c  = MakeGlueFaces(m, 1.e-5)
 400     setName(c)
 401     return c
 402
 403 def BlockMakeMultiRotation1D(shape, dir, point, times):
 404     m = geom.MakeMultiRotation1D(shape, dir, point, times)
 405     c  = MakeGlueFaces(m, 1.e-5)
 406     setName(c)
 407     return c
 408
 409 def BlockMakeMultiRotation2D(shape, dir, point, angle, times1, step, times2):
 410     m = geom.MakeMultiRotation2D(shape, dir, point, angle, times1, step, times2)
 411     c  = MakeGlueFaces(m, 1.e-5)
 412     setName(c)
 413     return c
 414
 415 # internal
 416 # --------
 417
 418 def MakeCompound(ls):
 419     """
 420     Creates a compound defined by a list
 421     """
 422     c = geom.MakeCompound(ls)
 423     return c
 424
 425 def MakeSewingShape(s, eps):
 426     """
 427     Creates a shape fully sewed
 428     """
 429     r = geom.MakeSewingShape(s, eps)
 430     return r
 431
 432 def MakeGlueFaces(s, eps):
 433     """
 434     Touched faces are replaced by one
 435     """
 436     r = geom.MakeGlueFaces(s, eps)
 437     return r
 438
 439
 440 # Geometry
 441 # ========
 442
 443 # A small cube centered and put on a great cube
 444
 445 # Points
 446 # ------
 447
 448 greatPoint111 = point( 0,  0,  0)
 449 greatPoint211 = point(10,  0,  0)
 450 greatPoint311 = point(20,  0,  0)
 451 greatPoint411 = point(30,  0,  0)
 452
 453 greatPoint121 = point( 0, 10,  0)
 454 greatPoint221 = point(10, 10,  0)
 455 greatPoint321 = point(20, 10,  0)
 456 greatPoint421 = point(30, 10,  0)
 457
 458 greatPoint112 = point( 0,  0, 10)
 459 greatPoint212 = point(10,  0, 10)
 460 greatPoint312 = point(20,  0, 10)
 461 greatPoint412 = point(30,  0, 10)
 462
 463 greatPoint122 = point( 0, 10, 10)
 464 greatPoint222 = point(10, 10, 10)
 465 greatPoint322 = point(20, 10, 10)
 466 greatPoint422 = point(30, 10, 10)
 467
 468 greatPoint113 = point( 0,  0, 20)
 469 greatPoint213 = point(10,  0, 20)
 470 greatPoint313 = point(20,  0, 20)
 471 greatPoint413 = point(30,  0, 20)
 472
 473 greatPoint123 = point( 0, 10, 20)
 474 greatPoint223 = point(10, 10, 20)
 475 greatPoint323 = point(20, 10, 20)
 476 greatPoint423 = point(30, 10, 20)
 477
 478 greatPoint114 = point( 0,  0, 30)
 479 greatPoint214 = point(10,  0, 30)
 480 greatPoint314 = point(20,  0, 30)
 481 greatPoint414 = point(30,  0, 30)
 482
 483 greatPoint124 = point( 0, 10, 30)
 484 greatPoint224 = point(10, 10, 30)
 485 greatPoint324 = point(20, 10, 30)
 486 greatPoint424 = point(30, 10, 30)
 487
 488
 489 smallPoint111 = greatPoint222
 490 smallPoint211 = greatPoint322
 491 smallPoint121 = point(10, 20, 10)
 492 smallPoint221 = point(20, 20, 10)
 493
 494 smallPoint112 = greatPoint223
 495 smallPoint212 = greatPoint323
 496 smallPoint122 = point(10, 20, 20)
 497 smallPoint222 = point(20, 20, 20)
 498
 499 # Edges
 500 # -----
 501
 502 smallEdgeX11 = MakeEdge(smallPoint111, smallPoint211)
 503 smallEdgeX21 = MakeEdge(smallPoint121, smallPoint221)
 504 smallEdgeX12 = MakeEdge(smallPoint112, smallPoint212)
 505 smallEdgeX22 = MakeEdge(smallPoint122, smallPoint222)
 506
 507 smallEdgeY11 = MakeEdge(smallPoint111, smallPoint121)
 508 smallEdgeY21 = MakeEdge(smallPoint211, smallPoint221)
 509 smallEdgeY12 = MakeEdge(smallPoint112, smallPoint122)
 510 smallEdgeY22 = MakeEdge(smallPoint212, smallPoint222)
 511
 512 smallEdgeZ11 = MakeEdge(smallPoint111, smallPoint112)
 513 smallEdgeZ21 = MakeEdge(smallPoint211, smallPoint212)
 514 smallEdgeZ12 = MakeEdge(smallPoint121, smallPoint122)
 515 smallEdgeZ22 = MakeEdge(smallPoint221, smallPoint222)
 516
 517
 518 greatEdgeX111 = MakeEdge(greatPoint111, greatPoint211)
 519 greatEdgeX211 = MakeEdge(greatPoint211, greatPoint311)
 520 greatEdgeX311 = MakeEdge(greatPoint311, greatPoint411)
 521 greatEdgeX121 = MakeEdge(greatPoint121, greatPoint221)
 522 greatEdgeX221 = MakeEdge(greatPoint221, greatPoint321)
 523 greatEdgeX321 = MakeEdge(greatPoint321, greatPoint421)
 524
 525 greatEdgeX112 = MakeEdge(greatPoint112, greatPoint212)
 526 greatEdgeX212 = MakeEdge(greatPoint212, greatPoint312)
 527 greatEdgeX312 = MakeEdge(greatPoint312, greatPoint412)
 528 greatEdgeX122 = MakeEdge(greatPoint122, greatPoint222)
 529 greatEdgeX222 = smallEdgeX11
 530 greatEdgeX322 = MakeEdge(greatPoint322, greatPoint422)
 531
 532 greatEdgeX113 = MakeEdge(greatPoint113, greatPoint213)
 533 greatEdgeX213 = MakeEdge(greatPoint213, greatPoint313)
 534 greatEdgeX313 = MakeEdge(greatPoint313, greatPoint413)
 535 greatEdgeX123 = MakeEdge(greatPoint123, greatPoint223)
 536 greatEdgeX223 = smallEdgeX12
 537 greatEdgeX323 = MakeEdge(greatPoint323, greatPoint423)
 538
 539 greatEdgeX114 = MakeEdge(greatPoint114, greatPoint214)
 540 greatEdgeX214 = MakeEdge(greatPoint214, greatPoint314)
 541 greatEdgeX314 = MakeEdge(greatPoint314, greatPoint414)
 542 greatEdgeX124 = MakeEdge(greatPoint124, greatPoint224)
 543 greatEdgeX224 = MakeEdge(greatPoint224, greatPoint324)
 544 greatEdgeX324 = MakeEdge(greatPoint324, greatPoint424)
 545
 546 greatEdgeY11 = MakeEdge(greatPoint111, greatPoint121)
 547 greatEdgeY21 = MakeEdge(greatPoint211, greatPoint221)
 548 greatEdgeY31 = MakeEdge(greatPoint311, greatPoint321)
 549 greatEdgeY41 = MakeEdge(greatPoint411, greatPoint421)
 550
 551 greatEdgeY12 = MakeEdge(greatPoint112, greatPoint122)
 552 greatEdgeY22 = MakeEdge(greatPoint212, greatPoint222)
 553 greatEdgeY32 = MakeEdge(greatPoint312, greatPoint322)
 554 greatEdgeY42 = MakeEdge(greatPoint412, greatPoint422)
 555
 556 greatEdgeY13 = MakeEdge(greatPoint113, greatPoint123)
 557 greatEdgeY23 = MakeEdge(greatPoint213, greatPoint223)
 558 greatEdgeY33 = MakeEdge(greatPoint313, greatPoint323)
 559 greatEdgeY43 = MakeEdge(greatPoint413, greatPoint423)
 560
 561 greatEdgeY14 = MakeEdge(greatPoint114, greatPoint124)
 562 greatEdgeY24 = MakeEdge(greatPoint214, greatPoint224)
 563 greatEdgeY34 = MakeEdge(greatPoint314, greatPoint324)
 564 greatEdgeY44 = MakeEdge(greatPoint414, greatPoint424)
 565
 566 greatEdgeZ111 = MakeEdge(greatPoint111, greatPoint112)
 567 greatEdgeZ211 = MakeEdge(greatPoint211, greatPoint212)
 568 greatEdgeZ311 = MakeEdge(greatPoint311, greatPoint312)
 569 greatEdgeZ411 = MakeEdge(greatPoint411, greatPoint412)
 570
 571 greatEdgeZ121 = MakeEdge(greatPoint121, greatPoint122)
 572 greatEdgeZ221 = MakeEdge(greatPoint221, greatPoint222)
 573 greatEdgeZ321 = MakeEdge(greatPoint321, greatPoint322)
 574 greatEdgeZ421 = MakeEdge(greatPoint421, greatPoint422)
 575
 576 greatEdgeZ112 = MakeEdge(greatPoint112, greatPoint113)
 577 greatEdgeZ212 = MakeEdge(greatPoint212, greatPoint213)
 578 greatEdgeZ312 = MakeEdge(greatPoint312, greatPoint313)
 579 greatEdgeZ412 = MakeEdge(greatPoint412, greatPoint413)
 580
 581 greatEdgeZ122 = MakeEdge(greatPoint122, greatPoint123)
 582 greatEdgeZ222 = smallEdgeZ11
 583 greatEdgeZ322 = smallEdgeZ21
 584 greatEdgeZ422 = MakeEdge(greatPoint422, greatPoint423)
 585
 586 greatEdgeZ113 = MakeEdge(greatPoint113, greatPoint114)
 587 greatEdgeZ213 = MakeEdge(greatPoint213, greatPoint214)
 588 greatEdgeZ313 = MakeEdge(greatPoint313, greatPoint314)
 589 greatEdgeZ413 = MakeEdge(greatPoint413, greatPoint414)
 590
 591 greatEdgeZ123 = MakeEdge(greatPoint123, greatPoint124)
 592 greatEdgeZ223 = MakeEdge(greatPoint223, greatPoint224)
 593 greatEdgeZ323 = MakeEdge(greatPoint323, greatPoint324)
 594 greatEdgeZ423 = MakeEdge(greatPoint423, greatPoint424)
 595
 596 # Faces
 597 # -----
 598
 599 smallFaceX1 = MakeQuadFace(smallEdgeY11, smallEdgeZ11, smallEdgeY12, smallEdgeZ12)
 600 smallFaceX2 = MakeQuadFace(smallEdgeY21, smallEdgeZ21, smallEdgeY22, smallEdgeZ22)
 601 smallFaceY1 = MakeQuadFace(smallEdgeX11, smallEdgeZ11, smallEdgeX12, smallEdgeZ21)
 602 smallFaceY2 = MakeQuadFace(smallEdgeX21, smallEdgeZ12, smallEdgeX22, smallEdgeZ22)
 603 smallFaceZ1 = MakeQuadFace(smallEdgeX11, smallEdgeY11, smallEdgeX21, smallEdgeY21)
 604 smallFaceZ2 = MakeQuadFace(smallEdgeX12, smallEdgeY12, smallEdgeX22, smallEdgeY22)
 605
 606
 607 greatFaceX11 = MakeQuadFace(greatEdgeY11, greatEdgeZ111, greatEdgeY12, greatEdgeZ121)
 608 greatFaceX21 = MakeQuadFace(greatEdgeY21, greatEdgeZ211, greatEdgeY22, greatEdgeZ221)
 609 greatFaceX31 = MakeQuadFace(greatEdgeY31, greatEdgeZ311, greatEdgeY32, greatEdgeZ321)
 610 greatFaceX41 = MakeQuadFace(greatEdgeY41, greatEdgeZ411, greatEdgeY42, greatEdgeZ421)
 611
 612 greatFaceX12 = MakeQuadFace(greatEdgeY12, greatEdgeZ112, greatEdgeY13, greatEdgeZ122)
 613 greatFaceX22 = MakeQuadFace(greatEdgeY22, greatEdgeZ212, greatEdgeY23, greatEdgeZ222)
 614 greatFaceX32 = MakeQuadFace(greatEdgeY32, greatEdgeZ312, greatEdgeY33, greatEdgeZ322)
 615 greatFaceX42 = MakeQuadFace(greatEdgeY42, greatEdgeZ412, greatEdgeY43, greatEdgeZ422)
 616
 617 greatFaceX13 = MakeQuadFace(greatEdgeY13, greatEdgeZ113, greatEdgeY14, greatEdgeZ123)
 618 greatFaceX23 = MakeQuadFace(greatEdgeY23, greatEdgeZ213, greatEdgeY24, greatEdgeZ223)
 619 greatFaceX33 = MakeQuadFace(greatEdgeY33, greatEdgeZ313, greatEdgeY34, greatEdgeZ323)
 620 greatFaceX43 = MakeQuadFace(greatEdgeY43, greatEdgeZ413, greatEdgeY44, greatEdgeZ423)
 621
 622 greatFaceY111 = MakeQuadFace(greatEdgeX111, greatEdgeZ111, greatEdgeX112, greatEdgeZ211)
 623 greatFaceY211 = MakeQuadFace(greatEdgeX211, greatEdgeZ211, greatEdgeX212, greatEdgeZ311)
 624 greatFaceY311 = MakeQuadFace(greatEdgeX311, greatEdgeZ311, greatEdgeX312, greatEdgeZ411)
 625 greatFaceY121 = MakeQuadFace(greatEdgeX121, greatEdgeZ121, greatEdgeX122, greatEdgeZ221)
 626 greatFaceY221 = MakeQuadFace(greatEdgeX221, greatEdgeZ221, greatEdgeX222, greatEdgeZ321)
 627 greatFaceY321 = MakeQuadFace(greatEdgeX321, greatEdgeZ321, greatEdgeX322, greatEdgeZ421)
 628
 629 greatFaceY112 = MakeQuadFace(greatEdgeX112, greatEdgeZ112, greatEdgeX113, greatEdgeZ212)
 630 greatFaceY212 = MakeQuadFace(greatEdgeX212, greatEdgeZ212, greatEdgeX213, greatEdgeZ312)
 631 greatFaceY312 = MakeQuadFace(greatEdgeX312, greatEdgeZ312, greatEdgeX313, greatEdgeZ412)
 632 greatFaceY122 = MakeQuadFace(greatEdgeX122, greatEdgeZ122, greatEdgeX123, greatEdgeZ222)
 633 greatFaceY222 = smallFaceY1
 634 greatFaceY322 = MakeQuadFace(greatEdgeX322, greatEdgeZ322, greatEdgeX323, greatEdgeZ422)
 635
 636 greatFaceY113 = MakeQuadFace(greatEdgeX113, greatEdgeZ113, greatEdgeX114, greatEdgeZ213)
 637 greatFaceY213 = MakeQuadFace(greatEdgeX213, greatEdgeZ213, greatEdgeX214, greatEdgeZ313)
 638 greatFaceY313 = MakeQuadFace(greatEdgeX313, greatEdgeZ313, greatEdgeX314, greatEdgeZ413)
 639 greatFaceY123 = MakeQuadFace(greatEdgeX123, greatEdgeZ123, greatEdgeX124, greatEdgeZ223)
 640 greatFaceY223 = MakeQuadFace(greatEdgeX223, greatEdgeZ223, greatEdgeX224, greatEdgeZ323)
 641 greatFaceY323 = MakeQuadFace(greatEdgeX323, greatEdgeZ323, greatEdgeX324, greatEdgeZ423)
 642
 643 greatFaceZ11 = MakeQuadFace(greatEdgeX111, greatEdgeY11, greatEdgeX121, greatEdgeY21)
 644 greatFaceZ21 = MakeQuadFace(greatEdgeX211, greatEdgeY21, greatEdgeX221, greatEdgeY31)
 645 greatFaceZ31 = MakeQuadFace(greatEdgeX311, greatEdgeY31, greatEdgeX321, greatEdgeY41)
 646
 647 greatFaceZ12 = MakeQuadFace(greatEdgeX112, greatEdgeY12, greatEdgeX122, greatEdgeY22)
 648 greatFaceZ22 = MakeQuadFace(greatEdgeX212, greatEdgeY22, greatEdgeX222, greatEdgeY32)
 649 greatFaceZ32 = MakeQuadFace(greatEdgeX312, greatEdgeY32, greatEdgeX322, greatEdgeY42)
 650
 651 greatFaceZ13 = MakeQuadFace(greatEdgeX113, greatEdgeY13, greatEdgeX123, greatEdgeY23)
 652 greatFaceZ23 = MakeQuadFace(greatEdgeX213, greatEdgeY23, greatEdgeX223, greatEdgeY33)
 653 greatFaceZ33 = MakeQuadFace(greatEdgeX313, greatEdgeY33, greatEdgeX323, greatEdgeY43)
 654
 655 greatFaceZ14 = MakeQuadFace(greatEdgeX114, greatEdgeY14, greatEdgeX124, greatEdgeY24)
 656 greatFaceZ24 = MakeQuadFace(greatEdgeX214, greatEdgeY24, greatEdgeX224, greatEdgeY34)
 657 greatFaceZ34 = MakeQuadFace(greatEdgeX314, greatEdgeY34, greatEdgeX324, greatEdgeY44)
 658
 659 # Solids
 660 # ------
 661
 662 smallBlock   = MakeHexaSolid(smallFaceX1, smallFaceX2, smallFaceY1, smallFaceY2, smallFaceZ1, smallFaceZ2)
 663
 664 greatBlock11 = MakeHexaSolid(greatFaceX11, greatFaceX21, greatFaceY111, greatFaceY121, greatFaceZ11, greatFaceZ12)
 665 greatBlock21 = MakeHexaSolid(greatFaceX21, greatFaceX31, greatFaceY211, greatFaceY221, greatFaceZ21, greatFaceZ22)
 666 greatBlock31 = MakeHexaSolid(greatFaceX31, greatFaceX41, greatFaceY311, greatFaceY321, greatFaceZ31, greatFaceZ32)
 667
 668 greatBlock12 = MakeHexaSolid(greatFaceX12, greatFaceX22, greatFaceY112, greatFaceY122, greatFaceZ12, greatFaceZ13)
 669 greatBlock22 = MakeHexaSolid(greatFaceX22, greatFaceX32, greatFaceY212, greatFaceY222, greatFaceZ22, greatFaceZ23)
 670 greatBlock32 = MakeHexaSolid(greatFaceX32, greatFaceX42, greatFaceY312, greatFaceY322, greatFaceZ32, greatFaceZ33)
 671
 672 greatBlock13 = MakeHexaSolid(greatFaceX13, greatFaceX23, greatFaceY113, greatFaceY123, greatFaceZ13, greatFaceZ14)
 673 greatBlock23 = MakeHexaSolid(greatFaceX23, greatFaceX33, greatFaceY213, greatFaceY223, greatFaceZ23, greatFaceZ24)
 674 greatBlock33 = MakeHexaSolid(greatFaceX33, greatFaceX43, greatFaceY313, greatFaceY323, greatFaceZ33, greatFaceZ34)
 675
 676 # Compound
 677 # --------
 678
 679 l = []
 680 l.append(smallBlock)
 681 l.append(greatBlock11)
 682 l.append(greatBlock21)
 683 l.append(greatBlock31)
 684 l.append(greatBlock12)
 685 l.append(greatBlock22)
 686 l.append(greatBlock32)
 687 l.append(greatBlock13)
 688 l.append(greatBlock23)
 689 l.append(greatBlock33)
 690
 691 piece = MakeCompound(l)
 692
 693 # Add in study
 694 # ------------
 695
 696 idpiece = addToStudy(piece, "Cubes2pyGibi")
 697 #geom.InsertOp.Export(piece,"piece.brep", "BREP")