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"

Conflicts:
	src/Standard/Standard_Real.cxx

diff --git a/src/Standard/Standard_Real.cxx b/src/Standard/Standard_Real.cxx
index b587af4f826fc021f746708deac92d54507f3372..7c03f07d0a5f75bfeac5da82d648c051c0c7c392 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
 // ------------------------------------------------------------------
@@ -44,9 +46,17 @@ 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)){
     Standard_RangeError::Raise();
   } 
+  else if (Value > 1.)
+  {
+    return 0.; //acos(1.)
+  }
+  else if (Value < -1.)
+  {
+    return M_PI; //acos(-1.)
+  }
   return acos(Value); 
 }
 
@@ -94,9 +104,17 @@ Standard_Real ACosApprox (const Standard_Real Value)
 //-------------------------------------------------------------------
 Standard_Real ASin (const Standard_Real Value) 
 { 
-  if ( Value < -1 || Value > 1 ){
+  if ((Value < -ACosLimit) || (Value > ACosLimit)){
     Standard_RangeError::Raise();
   }
+  else if (Value > 1.)
+  {
+    return M_PI_2; //asin(1.)
+  }
+  else if (Value < -1.)
+  {
+    return -M_PI_2; //asin(-1.)
+  }
   return asin(Value); 
 }