Protection of acos against values out of [-1.,1.]

diff --git a/src/INTERP_KERNEL/TransformedTriangleMath.cxx b/src/INTERP_KERNEL/TransformedTriangleMath.cxx
index 4dd23859cb58022737228ee627aa1c0bd325e28f..047c59bf95608364a1d72110593957e8dd3c36f4 100644 (file)
--- a/src/INTERP_KERNEL/TransformedTriangleMath.cxx
+++ b/src/INTERP_KERNEL/TransformedTriangleMath.cxx
@@ -354,7 +354,10 @@ namespace INTERP_KERNEL
     
     //? is this more stable? -> no subtraction
     //    return asin( dotProd / ( lenNormal * lenEdgeVec ) ) + 3.141592625358979 / 2.0;
-    return atan(1.0)*4.0 - acos( dotProd / ( lenNormal * lenEdgeVec ) );
+    double tmp=dotProd / ( lenNormal * lenEdgeVec );
+    tmp=std::max(tmp,-1.);
+    tmp=std::min(tmp,1.);
+    return atan(1.0)*4.0 - acos(tmp);
 
   }