-anElemType = SMESH.ALL;
-print "anElemType =", anElemType
-#anIds = CheckBelongToGeomFilter(smesh,mesh,box,subShapeList[1],SMESH.FACE)
-anIds = CheckBelongToGeomFilter(smesh,mesh,box,box,SMESH.FACE)
-print "anIds = ", anIds
+## Current style
+def CheckBelongToGeomFilter(theMesh, theShape, theSubShape, theElemType):
+ if theShape != theSubShape:
+ aName = str(theSubShape)
+ geompy.addToStudyInFather(theShape,theSubShape,aName)
+
+ theMesh.Compute()
+ aFilter = smesh.GetFilter(theElemType, SMESH.FT_BelongToGeom, theSubShape)
+ return aFilter.GetElementsId(theMesh.GetMesh())
+
+
+anElemType = SMESH.FACE;
+print("anElemType =", anElemType)
+#anIds = CheckBelongToGeomFilter(mesh,box,subShapeList[1],anElemType)
+anIds = CheckBelongToGeomFilter(mesh,box,box,anElemType)
+print("Number of ids = ", len(anIds))
+print("anIds = ", anIds)
+## Check old version
+#anIds = CheckBelongToGeomFilterOld(smesh,mesh.GetMesh(),box,box,anElemType)
+#print "anIds = ", anIds