doc/salome/examples/grouping_elements_ex06.py

   1 # Intersection 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 < 60
  21 aFilter = smesh.GetFilter(smesh.FACE, smesh.FT_Area, smesh.FT_LessThan, 60.)
  22
  23 anIds = mesh.GetIdsFromFilter(aFilter)
  24
  25 print "Criterion: Area < 60, Nb = ", len(anIds)
  26
  27 # create a group by adding elements with area < 60
  28 aGroup2 = mesh.CreateEmptyGroup(smesh.FACE, "Area < 60")
  29 aGroup2.Add(anIds)
  30
  31 # create an intersection of groups : 20 < area < 60
  32 aGroup3 = mesh.IntersectListOfGroups([aGroup1, aGroup2], "20 < Area < 60")
  33 print "Criterion: 20 < Area < 60, Nb = ", len(aGroup3.GetListOfID())
  34 # Please note that also there is IntersectGroups() method which works with two groups only
  35
  36 salome.sg.updateObjBrowser(1)