X-Git-Url: http://git.salome-platform.org/gitweb/?a=blobdiff_plain;f=doc%2Fexamples%2Fcreating_meshes_ex03.py;h=88b3d0b1afee1886ca1cfaef928cc5f116c8579a;hb=796122663b68ea1cee8c90033398a34d93518d39;hp=5ad4e10a014575d178ca45ca17339e928af4658b;hpb=4cf07a14111e98e8889620ee7e6371574c31a50c;p=modules%2Fsmesh.git diff --git a/doc/examples/creating_meshes_ex03.py b/doc/examples/creating_meshes_ex03.py index 5ad4e10a0..88b3d0b1a 100644 --- a/doc/examples/creating_meshes_ex03.py +++ b/doc/examples/creating_meshes_ex03.py @@ -9,50 +9,43 @@ 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"]) +box = geom_builder.MakeBoxDXDYDZ(200, 200, 200) +[Face_1, Face_2, Face_3, Face_4, Face_5, Face_6] = geom_builder.SubShapeAllSorted(box, geom_builder.ShapeType["FACE"]) # create Mesh object on Box shape -Mesh_1 = smesh_builder.Mesh(Box_1) +mesh = smesh_builder.Mesh(box) # 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) +mesh.Segment().NumberOfSegments(20) +mesh.Triangle().MaxElementArea(1200) +mesh.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() +mesh.Segment(geom=Face_1).NumberOfSegments(4) +mesh.Triangle(geom=Face_1) # 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() +mesh.Segment(geom=Face_2).NumberOfSegments(8) +mesh.Triangle(geom=Face_2) # 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()) +mesh.Segment(geom=Face_3).NumberOfSegments(12) +mesh.Triangle(geom=Face_3) + +# get existing sub-mesh priority order: F1 -> F2 -> F3 +[[SubMesh_F1, SubMesh_F3, SubMesh_F2]] = mesh.GetMeshOrder() +isDone = mesh.Compute() +if not isDone: raise Exception("Error when computing Mesh") +print("Nb elements at initial order of sub-meshes:", mesh.NbElements()) + +# set new sub-mesh order: F2 -> F1 -> F3 +isDone = mesh.SetMeshOrder([[SubMesh_F2, SubMesh_F1, SubMesh_F3]]) +isDone = mesh.Compute() +if not isDone: raise Exception("Error when computing Mesh") +print("Nb elements at new order of sub-meshes:", mesh.NbElements()) + +# compute with other sub-mesh order: F3 -> F2 -> F1 +isDone = mesh.SetMeshOrder([[SubMesh_F3, SubMesh_F2, SubMesh_F1]]) +isDone = mesh.Compute() +if not isDone: raise Exception("Error when computing Mesh") +print("Nb elements at another order of sub-meshes:", mesh.NbElements())