--- /dev/null
+# Check patch face functionality
+
+def CheckFacesOnHoles( faces ):
+ """
+ check that none of the faces in the list have holes
+ param faces - list of faces
+ return bool - *True* if all faces not have hole; *False* otherwise
+ """
+ for face in faces:
+ newFaces = geompy.PatchFace(face)
+ if len(newFaces) != 1:
+ return False
+ return True
+
+def CompareAreaAfterPatchFace( originalFace, patchFaceResult, eps=1e-06 ):
+ """
+ Specific function for check result of PatchFace operation.
+ Check that area of original face equal area first face from list (face by outer wire of original face)
+ minus the areas of the rest elements from list (holes)
+ param originalFace - original face
+ param patchFaceResult - list of faces result PatchFace operation.
+ First element - face from outer wire of original face,
+ other elements - faces from holes on original face
+ param eps - defines tolerance for comparison
+ return bool - *True* if area of original face is equal of area face[0] minus area of holes; *False* otherwise
+ """
+ areaOfHoles = 0.
+ for index in range(1, len(patchFaceResult)):
+ areaOfHoles += geompy.BasicProperties(patchFaceResult[index])[1]
+ return geompy.BasicProperties(originalFace)[1] - (geompy.BasicProperties(patchFaceResult[0])[1] - areaOfHoles) <= eps
+
+import math
+import salome
+salome.salome_init_without_session()
+import GEOM
+from salome.geom import geomBuilder
+geompy = geomBuilder.New()
+
+# Create shape
+vz = geompy.MakeVectorDXDYDZ(0, 0, 1)
+pyz = geompy.MakeVertex(0, -150, 100)
+
+geomObj = geompy.MakeMarker(0, 0, 0, 1, 0, 0, 0, 1, 0)
+sk = geompy.Sketcher2D()
+sk.addPoint(30.0000000, 50.0000000)
+sk.addArcAbsolute(70.0000000, 20.0000000)
+sk.addSegmentAbsolute(0.0000000, 0.0000000)
+sk.close()
+sketch = sk.wire(geomObj)
+face_1 = geompy.MakeFaceWires([sketch], 1)
+multiTranslation = geompy.MakeMultiTranslation1D(face_1, None, 105, 3)
+[face_2,face_3,face_4] = geompy.ExtractShapes(multiTranslation, geompy.ShapeType["FACE"], True)
+scale_1 = geompy.MakeScaleTransform(face_3, None, 0.3)
+scale_2 = geompy.MakeScaleTransform(face_4, None, 0.3)
+
+translation_1 = geompy.MakeTranslation(scale_1, -10, 25, 0)
+translation_2 = geompy.MakeTranslation(scale_2, -25, 20, 0)
+
+rotation_1 = geompy.MakeRotation(translation_1, vz, -19*math.pi/180.0)
+rotation_2 = geompy.MakeRotation(translation_2, vz, 15*math.pi/180.0)
+
+cut = geompy.MakeCutList(face_1, [rotation_2, rotation_1], True)
+
+#Perform oepration
+faces = geompy.PatchFace(cut)
+
+# Check, that result faces haven't holes
+assert(CheckFacesOnHoles(faces))
+
+#Check area
+assert(CompareAreaAfterPatchFace(cut, faces))
