0028458: Some faces are missing when loading a specific CAD. This is a regression...

Function ACos and ASin are modified to avoid exception when argument is +-(1. + Epsilon(1.)) because of "numerical noise"

diff --git a/src/Standard/Standard_Real.cxx b/src/Standard/Standard_Real.cxx
index fdc4220c2c1ba19ac5472dd939336a8ad5f4c29a..106b7ac4cf06b128f5cd5d47e0ccb3309f357b36 100644 (file)
--- a/src/Standard/Standard_Real.cxx
+++ b/src/Standard/Standard_Real.cxx
@@ -20,6 +20,8 @@
 #include <Standard_Stream.hxx>
 #include <Standard_OStream.hxx>
 
+static const Standard_Real ACosLimit = 1. + Epsilon(1.);
+
 // ------------------------------------------------------------------
 // Hascode : Computes a hascoding value for a given real
 // ------------------------------------------------------------------
@@ -43,10 +45,18 @@ Standard_Integer HashCode(const Standard_Real me, const Standard_Integer Upper)
 //-------------------------------------------------------------------
 Standard_Real ACos (const Standard_Real Value) 
 { 
-  if ( (Value < -1.) || (Value > 1.) ){
+  if ((Value < -ACosLimit) || (Value > ACosLimit)){
     throw Standard_RangeError();
-  } 
-  return acos(Value); 
+  }
+  else if (Value > 1.)
+  {
+    return 0.; //acos(1.)
+  }
+  else if (Value < -1.)
+  {
+    return M_PI; //acos(-1.)
+  }
+  return acos(Value);
 }
 
 //-------------------------------------------------------------------
@@ -93,10 +103,18 @@ Standard_Real ACosApprox (const Standard_Real Value)
 //-------------------------------------------------------------------
 Standard_Real ASin (const Standard_Real Value) 
 { 
-  if ( Value < -1 || Value > 1 ){
+  if ((Value < -ACosLimit) || (Value > ACosLimit)){
     throw Standard_RangeError();
   }
-  return asin(Value); 
+  else if (Value > 1.)
+  {
+    return M_PI_2; //asin(1.)
+  }
+  else if (Value < -1.)
+  {
+    return -M_PI_2; //asin(-1.)
+  }
+  return asin(Value);
 }
 
 //-------------------------------------------------------------------