bos #26461 [EDF] (2021) MedCoupling: debug tetra/tetra 3d intersector algorithm

diff --git a/src/INTERP_KERNEL/TransformedTriangleInline.hxx b/src/INTERP_KERNEL/TransformedTriangleInline.hxx
index 5f35327818af093d7793bb634d548046ab52d5f1..cc00d407463f74c101a61069ec40f00bcddfd8fe 100644 (file)
--- a/src/INTERP_KERNEL/TransformedTriangleInline.hxx
+++ b/src/INTERP_KERNEL/TransformedTriangleInline.hxx
@@ -188,7 +188,8 @@ inline bool TransformedTriangle::testEdgeIntersectsTriangle(const TetraEdge edge
 
   //? should equality with zero use epsilon?
   LOG(5, "testEdgeIntersectsTriangle : t1 = " << t1 << " t2 = " << t2 );
-  return (t1*t2 <= 0.0) && (t1 - t2 != 0.0);
+  return (t1*t2 <= 0.0) && !epsilonEqual(t1 - t2, 0.0);
+  //return (t1*t2 <= 0.0) && (t1 - t2 != 0.0);
 }
 
 inline bool TransformedTriangle::testFacetSurroundsSegment(const TriSegment seg, const TetraFacet facet) const