From ddce690b77898fd36494997025c7877ce8aee362 Mon Sep 17 00:00:00 2001 From: eap Date: Mon, 17 Jan 2022 14:39:32 +0300 Subject: [PATCH] bos #26461 [EDF] (2021) MedCoupling: debug tetra/tetra 3d intersector algorithm --- src/INTERP_KERNEL/TransformedTriangleInline.hxx | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/src/INTERP_KERNEL/TransformedTriangleInline.hxx b/src/INTERP_KERNEL/TransformedTriangleInline.hxx index 5f3532781..cc00d4074 100644 --- 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 -- 2.39.2