doc/salome/examples/filters_ex16.py

   1 # Double edges, Double faces, Double volumes
   2
   3 from smesh import *
   4 # make a mesh on a box
   5 box = geompy.MakeBoxDXDYDZ(100,100,100)
   6 mesh = Mesh( box, "Box" )
   7 mesh.Segment().NumberOfSegments(10)
   8 mesh.Quadrangle()
   9 mesh.Hexahedron()
  10 mesh.Compute()
  11 # copy all elements with translation and Merge nodes
  12 mesh.TranslateObject( mesh, MakeDirStruct( 10,0,0), Copy=True )
  13 mesh.MergeNodes( mesh.FindCoincidentNodes(1e-7) )
  14 # create filters to find equal elements
  15 equalEdgesFilter   = GetFilter(SMESH.EDGE, FT_EqualEdges)
  16 equalFacesFilter   = GetFilter(SMESH.FACE, FT_EqualFaces)
  17 equalVolumesFilter = GetFilter(SMESH.VOLUME, FT_EqualVolumes)
  18 # get equal elements
  19 print "Number of equal edges:",   len( mesh.GetIdsFromFilter( equalEdgesFilter ))
  20 print "Number of equal faces:",   len( mesh.GetIdsFromFilter( equalFacesFilter ))
  21 print "Number of equal volumes:", len( mesh.GetIdsFromFilter( equalVolumesFilter ))