doc/examples/grouping_elements_ex04.py

   1 # Edit a Group
   2
   3 from mechanic import *
   4
   5 # Get ids of all faces with area > 35
   6 aFilter = smesh_builder.GetFilter(SMESH.FACE, SMESH.FT_Area, SMESH.FT_MoreThan, 35.)
   7
   8 anIds = mesh.GetIdsFromFilter(aFilter)
   9
  10 print("Criterion: Area > 35, Nb = ", len(anIds))
  11
  12 # create a group by adding elements with area > 35
  13 aGroup = mesh.CreateEmptyGroup(SMESH.FACE, "Area > 35")
  14 aGroup.Add(anIds)
  15
  16 # Get ids of all faces with area > 40
  17 aFilter = smesh_builder.GetFilter(SMESH.FACE, SMESH.FT_Area, SMESH.FT_MoreThan, 40.)
  18
  19 anIds = mesh.GetIdsFromFilter(aFilter)
  20
  21 print("Criterion: Area > 40, Nb = ", len(anIds))
  22
  23 # create a group of elements with area [35; 40] by removing elements with area > 40 from group aGroup
  24 aGroup.Remove(anIds)
  25 aGroup.SetName("35 < Area < 40")
  26
  27 # print the result
  28 aGroupElemIDs = aGroup.GetListOfID()
  29
  30 print("Criterion: 35 < Area < 40, Nb = ", len(aGroupElemIDs))
  31
  32 j = 1
  33 for i in range(len(aGroupElemIDs)):
  34   if j > 20: j = 1; print("")
  35   print(aGroupElemIDs[i], end=' ')
  36   j = j + 1
  37   pass
  38 print("")