From 49adfd9a1deddf73e9896c8c1564cf296d621da3 Mon Sep 17 00:00:00 2001 From: Matthew Scroggs Date: Sat, 1 Feb 2025 10:41:56 +0000 Subject: [PATCH] Correct points for edge adjacent triangle and quad (#14) * eta2, not eta3 * correct another tpyo in rule --- src/duffy/triangle_quadrilateral.rs | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/duffy/triangle_quadrilateral.rs b/src/duffy/triangle_quadrilateral.rs index be9ea69..e9aea74 100644 --- a/src/duffy/triangle_quadrilateral.rs +++ b/src/duffy/triangle_quadrilateral.rs @@ -66,7 +66,7 @@ fn tri_quad_edge_points( // Second part quadrilateral_output_points.push((1.0 - xi) * eta3); - quadrilateral_output_points.push(xi * eta3); + quadrilateral_output_points.push(xi * eta2); triangle_output_points.push(xi * (1.0 - eta3) + eta3); triangle_output_points.push(xi * eta1); output_weights.push(weight0); @@ -97,7 +97,7 @@ fn tri_quad_edge_points( // Sixth part - quadrilateral_output_points.push(xi * eta1 * (1.0 - eta2 - eta3) + eta1); + quadrilateral_output_points.push(xi * eta1 * (1.0 - eta2 - eta3) + eta3); quadrilateral_output_points.push(xi); triangle_output_points.push(xi * eta1 * (1.0 - eta3) + eta3); triangle_output_points.push(xi * eta1);