Correct points for edge adjacent triangle and quad (#14)

* eta2, not eta3 * correct another tpyo in rule
bempp · Feb 1, 2025 · 49adfd9 · 49adfd9
1 parent 48abec5
commit 49adfd9
Showing 1 changed file with 2 additions and 2 deletions.
diff --git 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);