test/SMESH_create_dual_mesh.py

   1 #!/usr/bin/env python
   2
   3 ###
   4 ### This file is generated automatically by SALOME v9.9.0 with dump python functionality
   5 ###
   6
   7 import sys
   8 import salome
   9
  10 salome.salome_init()
  11 import salome_notebook
  12 notebook = salome_notebook.NoteBook()
  13 sys.path.insert(0, r'/home/B61570/work_in_progress/dual_mesh')
  14
  15 ###
  16 ### GEOM component
  17 ###
  18
  19 import GEOM
  20 from salome.geom import geomBuilder
  21 import math
  22 import SALOMEDS
  23
  24
  25 geompy = geomBuilder.New()
  26
  27 O = geompy.MakeVertex(0, 0, 0)
  28 OX = geompy.MakeVectorDXDYDZ(1, 0, 0)
  29 OY = geompy.MakeVectorDXDYDZ(0, 1, 0)
  30 OZ = geompy.MakeVectorDXDYDZ(0, 0, 1)
  31 Sphere_1 = geompy.MakeSphereR(100)
  32 [geomObj_1] = geompy.ExtractShapes(Sphere_1, geompy.ShapeType["FACE"], True)
  33 geompy.addToStudy( O, 'O' )
  34 geompy.addToStudy( OX, 'OX' )
  35 geompy.addToStudy( OY, 'OY' )
  36 geompy.addToStudy( OZ, 'OZ' )
  37 geompy.addToStudy( Sphere_1, 'Sphere_1' )
  38
  39 ###
  40 ### SMESH component
  41 ###
  42
  43 import  SMESH, SALOMEDS
  44 from salome.smesh import smeshBuilder
  45
  46 smesh = smeshBuilder.New()
  47
  48 NETGEN_3D_Parameters_1 = smesh.CreateHypothesisByAverageLength( 'NETGEN_Parameters', 'NETGENEngine', 34.641, 0 )
  49 Mesh_1 = smesh.Mesh(Sphere_1,'Mesh_1')
  50 status = Mesh_1.AddHypothesis( Sphere_1, NETGEN_3D_Parameters_1 )
  51 NETGEN_1D_2D_3D = Mesh_1.Tetrahedron(algo=smeshBuilder.NETGEN_1D2D3D)
  52 isDone = Mesh_1.Compute()
  53 dual_Mesh_1 = smesh.CreateDualMesh( Mesh_1, 'dual_Mesh_1', True)
  54
  55
  56 assert(dual_Mesh_1.NbPolyhedrons() > 0)
  57 assert(dual_Mesh_1.NbTetras() == 0)
  58
  59
  60 if salome.sg.hasDesktop():
  61   salome.sg.updateObjBrowser()