Salome HOME
23189: EDF 11603 - Dyssymetry in meshing
[modules/smesh.git] / doc / salome / examples / filters_ex01.py
index 86be381826b9b04c262cf4bbe392bbb0b939052a..88305c97f9523b52ad8809b37456bf7bf1d031a1 100644 (file)
@@ -2,7 +2,33 @@
 
 # create mesh
 from SMESH_mechanic import *
-# get faces with aspect ratio > 6.5
-filter = smesh.GetFilter(SMESH.FACE, SMESH.FT_AspectRatio, SMESH.FT_MoreThan, 6.5)
+
+# get faces with aspect ratio > 1.5
+filter = smesh.GetFilter(SMESH.FACE, SMESH.FT_AspectRatio, SMESH.FT_MoreThan, 1.5)
 ids = mesh.GetIdsFromFilter(filter)
-print "Number of faces with aspect ratio > 6.5:", len(ids)
+print "Number of faces with aspect ratio > 1.5:", len(ids)
+
+# copy the faces with aspect ratio > 1.5 to another mesh;
+# this demostrates that a filter can be used where usually a group or sub-mesh is acceptable
+filter.SetMesh( mesh.GetMesh() )
+mesh2 = smesh.CopyMesh( filter, "AR > 1.5" )
+print "Number of copied faces with aspect ratio > 1.5:", mesh2.NbFaces()
+
+# create a group (Group on Filter) of faces with Aspect Ratio < 1.5
+group = mesh.MakeGroup("AR < 1.5", SMESH.FACE, SMESH.FT_AspectRatio, '<', 1.5)
+print "Number of faces with aspect ratio < 1.5:", group.Size()
+
+# combine several criteria to Create a Group of only Triangular faces with Aspect Ratio < 1.5;
+# note that contents of a GroupOnFilter is dynamically updated as the mesh changes
+crit = [ smesh.GetCriterion( SMESH.FACE, SMESH.FT_AspectRatio, '<', 1.5, BinaryOp=SMESH.FT_LogicalAND ),
+         smesh.GetCriterion( SMESH.FACE, SMESH.FT_ElemGeomType,'=', SMESH.Geom_TRIANGLE ) ]
+triaGroup = mesh.MakeGroupByCriteria( "Tria AR < 1.5", crit )
+print "Number of triangles with aspect ratio < 1.5:", triaGroup.Size()
+
+# get range of values of Aspect Ratio of all faces in the mesh
+aspects = mesh.GetMinMax( SMESH.FT_AspectRatio )
+print "MESH: Min aspect = %s, Max aspect = %s" % ( aspects[0], aspects[1] )
+
+# get max value of Aspect Ratio of faces in triaGroup
+grAspects = mesh.GetMinMax( SMESH.FT_AspectRatio, triaGroup )
+print "GROUP: Max aspect = %s" % grAspects[1]