0026627: [Regression] Shape Healing hangs as of OCC 6.8.0

Check for orientation of the solid corrected to ensure that cycle always finishes.

Test case added: bugs modalg_6 bug26627
Tests boolean volumemaker A3, B5, B7 corrected (improvements)

src/BRepClass3d/BRepClass3d_SClassifier.cxx

@@ -33,6 +33,9 @@
 #include <math_BullardGenerator.hxx>
 #include <BRepTopAdaptor_FClass2d.hxx>
 
+#include <vector>
+
+// modified by NIZHNY-MKK  Mon Jun 21 15:13:40 2004
 static
   Standard_Boolean FaceNormal (const TopoDS_Face& aF,
                                const Standard_Real U,
@@ -102,25 +105,32 @@ void BRepClass3d_SClassifier::PerformInfinitePoint(BRepClass3d_SolidExplorer& aS
   myFace.Nullify();
   myState=2;
 
-  aSE.InitShell();
+  // Collect faces in sequence to iterate
-  if (aSE.MoreShell())
+  std::vector<TopoDS_Face> aFaces;
+  for (aSE.InitShell(); aSE.MoreShell(); aSE.NextShell())
   {
-    aSE.InitFace();
+    for (aSE.InitFace(); aSE.MoreFace(); aSE.NextFace())
-    if (aSE.MoreFace())
     {
-      TopoDS_Face aF = aSE.CurrentFace();
+      aFaces.push_back (aSE.CurrentFace());
+    }
+  }
+
+  // iteratively try up to 10 probing points from each face
+  const int NB_MAX_POINTS_PER_FACE = 10;
+  for (int itry = 0; itry < NB_MAX_POINTS_PER_FACE; itry++)
+  {
+    for (std::vector<TopoDS_Face>::iterator iFace = aFaces.begin(); iFace != aFaces.end(); ++iFace)
+    {
+      TopoDS_Face aF = *iFace;
+
       TopAbs_State aState = TopAbs_OUT;
       IntCurveSurface_TransitionOnCurve aTransition = IntCurveSurface_Tangent;
-      TopoDS_Face MinFace = aF;
+
-      for (;;)
-      {
         aParam = 0.1 + 0.8 * aRandomGenerator.NextReal(); // random number in range [0.1, 0.9]
         bFound = aSE.FindAPointInTheFace(aF, aPoint, aU, aV, aParam);
-        if (!bFound)
+      if (!bFound || !FaceNormal(aF, aU, aV, aDN))
-          return;
-
-        if (!FaceNormal(aF, aU, aV, aDN))
         continue;
+
         gp_Lin aLin(aPoint, -aDN);
         Standard_Real parmin = RealLast();
         for (aSE.InitShell();aSE.MoreShell();aSE.NextShell()) { 
@@ -141,7 +151,6 @@ void BRepClass3d_SClassifier::PerformInfinitePoint(BRepClass3d_SolidExplorer& aS
                     parmin = Intersector3d.WParameter(imin);
                     aState = Intersector3d.State(imin);
                     aTransition = Intersector3d.Transition(imin);
-                    MinFace = CurFace;
                   }
                 }
               }
@@ -164,10 +173,8 @@ void BRepClass3d_SClassifier::PerformInfinitePoint(BRepClass3d_SolidExplorer& aS
             return;
           }
         }
-        aF = MinFace;
+    } // iteration by faces
-      }
+  } // iteration by points
-    } //if (aSE.MoreFace())
-  } //if (aSE.MoreShell())
       }
 
 //=======================================================================

tests/boolean/volumemaker/A3

@@ -2,7 +2,6 @@
 # plane unstable
 
 puts "TODO OCC26020 ALL: Error: bopcheck failed"
-puts "TODO ?OCC26019 ALL: Faulty shapes in variables faulty_1 to faulty_"
 
 # planar face 
 plane pln_f1 0 -870 -1.3877787807814457e-014 0 1 1.1102230246251565e-016
@@ -37,5 +36,5 @@ mkface f6 pln_f6 -1000000 1000000 -1000000 1000000
 # make volume operation 
 mkvolume result f1 f2 f3 f4 f5 f6
 
-set square 2.4e+013
+set square 1183220.
 

tests/boolean/volumemaker/B5

@@ -1,8 +1,6 @@
 # test script on make volume operation
 # plane
 
-puts "TODO OCC26020 ALL: Error: bopcheck failed"
-
 # planar face 
 plane pln_f1 2949.26242211 -890 552.5 1.110223024625157e-016 -1 1.1102230246251563e-016
 erase pln_f1
@@ -36,5 +34,4 @@ mkface f6 pln_f6 -1000000 1000000 -1000000 1000000
 # make volume operation 
 mkvolume result f1 f2 f3 f4 f5 f6
 
-set square 2.4e+013
+set square 12400.
-

tests/boolean/volumemaker/B7

@@ -1,8 +1,7 @@
 # test script on make volume operation
 # plane
 
-puts "TODO ?OCC26020 ALL: Faulty shapes in variables faulty_1 to faulty_"
+puts "TODO OCC26020 ALL: Error: bopcheck failed"
-puts "TODO ?OCC26020 ALL: Error: bopcheck failed"
 
 # planar face 
 plane pln_f1 18.855982726712998 17.500000000800412 0 -0.96152394764524818 -0.27472112788189063 0
@@ -52,4 +51,4 @@ mkface f9 pln_f9 -1000000 1000000 -1000000 1000000
 # make volume operation 
 mkvolume result f1 f2 f3 f4 f5 f6 f7 f8 f9
 
-set square 2.75501e+013
+set square 1.3412e+013

tests/bugs/modalg_6/bug26627

@@ -0,0 +1,19 @@
+puts "##################################################"
+puts "0026627: Shape Healing hangs as of OCC 6.8.0"
+puts "##################################################"
+
+# load and check shape
+restore [locate_data_file bug26627_fixed.brep] a
+tolerance a
+checkshape a
+
+# call fixshape -- it should finish in fraction of second
+cpulimit 10
+fixshape result a 1e-6
+
+# result should have positive volume
+if { [regexp {Mass\s*:\s*([0-9.e+-]*)} [vprops result] dummy volume] } {
+  checkreal "Volume of the solid" $volume 4332.63 0.1 0.
+} else {
+  puts "Error: cannot get volume of result!"
+}