0024829: sprops with precision throws an exception

Reason of exception is reading nonexistent element of array.

Test case for issue CR24829

diff --git a/src/BRepGProp/BRepGProp_Sinert.cxx b/src/BRepGProp/BRepGProp_Sinert.cxx
index fcd2674..24b8877 100644 (file)
--- a/src/BRepGProp/BRepGProp_Sinert.cxx
+++ b/src/BRepGProp/BRepGProp_Sinert.cxx
@@ -216,9 +216,9 @@ static Standard_Integer LFillIntervalBounds(Standard_Real               A,
                                             const TColStd_Array1OfReal& Knots, 
                                             const Standard_Integer      NumSubs)
 {
-  Standard_Integer iEnd = Knots.Upper(), jEnd = L1->Upper();
+  Standard_Integer iEnd = MaxSubs(Knots.Upper()-1, NumSubs), jEnd = L1->Upper();
+  iEnd = Max(iEnd, Knots.Upper());
   if(iEnd - 1 > jEnd){
-    iEnd = MaxSubs(iEnd-1,NumSubs); 
     L1    = new math_Vector(1,iEnd);
     L2    = new math_Vector(1,iEnd);
     DimL  = new math_Vector(1,iEnd);
@@ -242,9 +242,9 @@ static Standard_Integer UFillIntervalBounds(Standard_Real               A,
                                             const TColStd_Array1OfReal& Knots, 
                                             const Standard_Integer      NumSubs)
 {
-  Standard_Integer iEnd = Knots.Upper(), jEnd = U1->Upper();
+  Standard_Integer iEnd = MaxSubs(Knots.Upper()-1, NumSubs), jEnd = U1->Upper();
+  iEnd = Max(iEnd, Knots.Upper());
   if(iEnd - 1 > jEnd){
-    iEnd = MaxSubs(iEnd-1,NumSubs); 
     U1   = new math_Vector(1,iEnd);
     U2   = new math_Vector(1,iEnd);
     DimU = new math_Vector(1,iEnd);

diff --git a/tests/bugs/modalg_5/bug24829 b/tests/bugs/modalg_5/bug24829
new file mode 100644 (file)
index 0000000..0400b1d
--- /dev/null
+++ b/tests/bugs/modalg_5/bug24829
@@ -0,0 +1,15 @@
+puts "========="
+puts "OCC24829"
+puts "========="
+puts ""
+###########################################################
+# sprops with precision throws an exception
+###########################################################
+
+restore [locate_data_file bug24829_faceOutput-wrong.brep] result
+
+sprops result 0.0001
+
+set square 15.1187
+
+set 3dviewer 1