doc/examples/grouping_elements_ex05.py

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