doc/salome/examples/grouping_elements_ex05.py

   1 # Union of groups
   2
   3 import SMESH_mechanic
   4
   5 smesh  = SMESH_mechanic.smesh
   6 mesh   = SMESH_mechanic.mesh
   7 salome = SMESH_mechanic.salome
   8
   9 # Criterion : AREA > 20
  10 aFilter = smesh.GetFilter(smesh.FACE, smesh.FT_Area, smesh.FT_MoreThan, 20.)
  11
  12 anIds = mesh.GetIdsFromFilter(aFilter)
  13
  14 print "Criterion: Area > 20, Nb = ", len( anIds )
  15
  16 # create a group by adding elements with area > 20
  17 aGroup1 = mesh.CreateEmptyGroup(smesh.FACE, "Area > 20")
  18 aGroup1.Add(anIds)
  19
  20 # Criterion : AREA = 20
  21 aFilter = smesh.GetFilter(smesh.FACE, smesh.FT_Area, smesh.FT_EqualTo, 20.)
  22
  23 anIds = mesh.GetIdsFromFilter(aFilter)
  24
  25 print "Criterion: Area = 20, Nb = ", len( anIds )
  26
  27 # create a group by adding elements with area = 20
  28 aGroup2 = mesh.CreateEmptyGroup( smesh.FACE, "Area = 20" )
  29
  30 aGroup2.Add(anIds)
  31
  32 # create union group : area >= 20
  33 aGroup3 = mesh.UnionListOfGroups([aGroup1, aGroup2], "Area >= 20")
  34 print "Criterion: Area >= 20, Nb = ", len(aGroup3.GetListOfID())
  35 # Please note that also there is UnionGroups() method which works with two groups only
  36
  37 # Criterion : AREA < 20
  38 aFilter = smesh.GetFilter(smesh.FACE, smesh.FT_Area, smesh.FT_LessThan, 20.)
  39
  40 anIds = mesh.GetIdsFromFilter(aFilter)
  41
  42 print "Criterion: Area < 20, Nb = ", len(anIds)
  43
  44 # create a group by adding elements with area < 20
  45 aGroup4 = mesh.CreateEmptyGroup(smesh.FACE, "Area < 20")
  46 aGroup4.Add(anIds)
  47
  48 # create union group : area >= 20 and area < 20
  49 aGroup5 = mesh.UnionListOfGroups([aGroup3, aGroup4], "Any Area")
  50 print "Criterion: Any Area, Nb = ", len(aGroup5.GetListOfID())
  51
  52 salome.sg.updateObjBrowser(1)