doc/examples/grouping_elements_ex06.py

   1 # Intersection of groups
   2
   3 from mechanic import *
   4
   5 # Criterion : AREA > 20
   6 aFilter = smesh_builder.GetFilter(SMESH.FACE, SMESH.FT_Area, SMESH.FT_MoreThan, 20.)
   7
   8 anIds = mesh.GetIdsFromFilter(aFilter)
   9
  10 print("Criterion: Area > 20, Nb = ", len(anIds))
  11
  12 # create a group by adding elements with area > 20
  13 aGroup1 = mesh.CreateEmptyGroup(SMESH.FACE, "Area > 20")
  14 aGroup1.Add(anIds)
  15
  16 # Criterion : AREA < 60
  17 aFilter = smesh_builder.GetFilter(SMESH.FACE, SMESH.FT_Area, SMESH.FT_LessThan, 60.)
  18
  19 anIds = mesh.GetIdsFromFilter(aFilter)
  20
  21 print("Criterion: Area < 60, Nb = ", len(anIds))
  22
  23 # create a group by adding elements with area < 60
  24 aGroup2 = mesh.CreateEmptyGroup(SMESH.FACE, "Area < 60")
  25 aGroup2.Add(anIds)
  26
  27 # create an intersection of groups : 20 < area < 60
  28 aGroup3 = mesh.IntersectListOfGroups([aGroup1, aGroup2], "20 < Area < 60")
  29 print("Criterion: 20 < Area < 60, Nb = ", len(aGroup3.GetListOfID()))
  30 # Please note that also there is IntersectGroups() method which works with two groups only
  31