From d03af8d8c8d01d5e417c23a9b2ec8638d5eee5f1 Mon Sep 17 00:00:00 2001 From: abn Date: Wed, 7 Feb 2024 11:09:10 +0100 Subject: [PATCH] [TetraIntersect] Settting TransformedTriangle::THRESHOLD to 100 + see previous commit : we can not keep 20 because of P1P1 intersector which uses barycentric computation. A small threshold produces sometimes more intersection points, which disturbs the (improper) barycentric computation. --- src/INTERP_KERNEL/TransformedTriangleMath.cxx | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/INTERP_KERNEL/TransformedTriangleMath.cxx b/src/INTERP_KERNEL/TransformedTriangleMath.cxx index bac06d5cd..30ecbc621 100644 --- a/src/INTERP_KERNEL/TransformedTriangleMath.cxx +++ b/src/INTERP_KERNEL/TransformedTriangleMath.cxx @@ -68,7 +68,7 @@ namespace INTERP_KERNEL const double TransformedTriangle::MULT_PREC_F = 4.0 * TransformedTriangle::MACH_EPS; /// Threshold for resetting double and triple products to zero; ( F / f in Grandy ) - const double TransformedTriangle::THRESHOLD_F = 20.0; + const double TransformedTriangle::THRESHOLD_F = 100.0; /// Threshold for what is considered a small enough angle to warrant correction of triple products by Grandy, [57] const double TransformedTriangle::TRIPLE_PRODUCT_ANGLE_THRESHOLD = 0.1; -- 2.39.2