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)

diff --git a/src/BRepClass3d/BRepClass3d_SClassifier.cxx b/src/BRepClass3d/BRepClass3d_SClassifier.cxx
index 4d01b2021b8bf148003886ac2cbb57ead4dfc0c7..14a34e9e773c6415f3c6f25a797071a0b619e823 100644 (file)
--- a/src/BRepClass3d/BRepClass3d_SClassifier.cxx
+++ b/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();
-  if (aSE.MoreShell())
+  // Collect faces in sequence to iterate
+  std::vector<TopoDS_Face> aFaces;
+  for (aSE.InitShell(); aSE.MoreShell(); aSE.NextShell())
+  {
+    for (aSE.InitFace(); aSE.MoreFace(); aSE.NextFace())
+    {
+      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++)
   {
-    aSE.InitFace();
-    if (aSE.MoreFace())
+    for (std::vector<TopoDS_Face>::iterator iFace = aFaces.begin(); iFace != aFaces.end(); ++iFace)
     {
-      TopoDS_Face aF = aSE.CurrentFace();
+      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)
-          return;
+      if (!bFound || !FaceNormal(aF, aU, aV, aDN))
+        continue;
 
-        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,11 +173,9 @@ void BRepClass3d_SClassifier::PerformInfinitePoint(BRepClass3d_SolidExplorer& aS
             return;
           }
         }
-        aF = MinFace;
+    } // iteration by faces
+  } // iteration by points
       }
-    } //if (aSE.MoreFace())
-  } //if (aSE.MoreShell())
-}
 
 //=======================================================================
 //function : Perform

diff --git a/tests/boolean/volumemaker/A3 b/tests/boolean/volumemaker/A3
index 5b5b2c37b30ea751034909444174a41e85eb847c..9aeac529d7f63bb930257243c8a723c05d7fd40a 100644 (file)
--- a/tests/boolean/volumemaker/A3
+++ b/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.
 

diff --git a/tests/boolean/volumemaker/B5 b/tests/boolean/volumemaker/B5
index 64f680918c6f4b4f533b28044b7af8436b9ee6be..06903f3eb5f12d1d1a21d59fc96a7bd2dd2ed3fc 100644 (file)
--- a/tests/boolean/volumemaker/B5
+++ b/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.

diff --git a/tests/boolean/volumemaker/B7 b/tests/boolean/volumemaker/B7
index 795efc44d0c6d9bf7078500190a4d4d2894cbbe2..a6e96335410d3dd7a81ae46c29bb63f36f364eba 100644 (file)
--- a/tests/boolean/volumemaker/B7
+++ b/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

diff --git a/tests/bugs/modalg_6/bug26627 b/tests/bugs/modalg_6/bug26627
new file mode 100644 (file)
index 0000000..b33ce91
--- /dev/null
+++ b/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!"
+}