Salome HOME
bos #29171 Refactor testing procedure
[modules/smesh.git] / doc / examples / creating_meshes_ex03.py
diff --git a/doc/examples/creating_meshes_ex03.py b/doc/examples/creating_meshes_ex03.py
new file mode 100644 (file)
index 0000000..5ad4e10
--- /dev/null
@@ -0,0 +1,58 @@
+# Change priority of sub-meshes in Mesh
+
+import salome
+salome.salome_init_without_session()
+
+from salome.geom import geomBuilder
+from salome.smesh import smeshBuilder
+
+geom_builder = geomBuilder.New()
+smesh_builder = smeshBuilder.New()
+
+Box_1 = geom_builder.MakeBoxDXDYDZ(200, 200, 200)
+[Face_1,Face_2,Face_3,Face_4,Face_5,Face_6] = geom_builder.SubShapeAllSorted(Box_1, geom_builder.ShapeType["FACE"])
+
+# create Mesh object on Box shape
+Mesh_1 = smesh_builder.Mesh(Box_1)
+
+# assign mesh algorithms and hypotheses
+Regular_1D = Mesh_1.Segment()
+Nb_Segments_1 = Regular_1D.NumberOfSegments(20)
+MEFISTO_2D = Mesh_1.Triangle()
+Max_Element_Area_1 = MEFISTO_2D.MaxElementArea(1200)
+Tetrahedron = Mesh_1.Tetrahedron()
+Max_Element_Volume_1 = Tetrahedron.MaxElementVolume(40000)
+
+# create sub-mesh and assign algorithms on Face_1
+Regular_1D_1 = Mesh_1.Segment(geom=Face_1)
+Nb_Segments_2 = Regular_1D_1.NumberOfSegments(4)
+MEFISTO_2D_1 = Mesh_1.Triangle(algo=smeshBuilder.MEFISTO,geom=Face_1)
+SubMesh_1 = MEFISTO_2D_1.GetSubMesh()
+
+# create sub-mesh and assign algorithms on Face_2
+Regular_1D_2 = Mesh_1.Segment(geom=Face_2)
+Nb_Segments_3 = Regular_1D_2.NumberOfSegments(8)
+MEFISTO_2D_2 = Mesh_1.Triangle(algo=smeshBuilder.MEFISTO,geom=Face_2)
+SubMesh_2 = MEFISTO_2D_2.GetSubMesh()
+
+# create sub-mesh and assign algorithms on Face_3
+Regular_1D_3 = Mesh_1.Segment(geom=Face_3)
+Nb_Segments_4 = Regular_1D_3.NumberOfSegments(12)
+MEFISTO_2D_3 = Mesh_1.Triangle(algo=smeshBuilder.MEFISTO,geom=Face_3)
+SubMesh_3 = MEFISTO_2D_3.GetSubMesh()
+
+# check existing sub-mesh priority order
+[ [ SubMesh_1, SubMesh_3, SubMesh_2 ] ] = Mesh_1.GetMeshOrder()
+isDone = Mesh_1.Compute()
+print("Nb elements at initial order of sub-meshes:", Mesh_1.NbElements())
+
+# set new sub-mesh order
+isDone = Mesh_1.SetMeshOrder( [ [ SubMesh_1, SubMesh_2, SubMesh_3 ] ])
+# compute mesh
+isDone = Mesh_1.Compute()
+print("Nb elements at new order of sub-meshes:", Mesh_1.NbElements())
+
+# compute with other sub-mesh order
+isDone = Mesh_1.SetMeshOrder( [ [ SubMesh_2, SubMesh_1, SubMesh_3 ] ])
+isDone = Mesh_1.Compute()
+print("Nb elements at another order of sub-meshes:", Mesh_1.NbElements())