fem: explicitly set use_canonical_quadrature_point_ordering=False for…

… special node types
firedrakeproject · Nov 22, 2024 · 9c88515 · 9c88515
1 parent 76d8daa
commit 9c88515
Showing 1 changed file with 13 additions and 0 deletions.
diff --git a/tests/regression/test_integral_hex.py b/tests/regression/test_integral_hex.py
@@ -53,3 +53,16 @@ def test_integral_hex_interior_facet_solve(mesh_from_file):
     solve(a == L, sol, bcs=[bc])
     err = assemble((sol - f)**2 * dx)**0.5
     assert err < 1.e-14
+
+
+def test_integral_hex_interior_facet_solve():
+    mesh = UnitCubeMesh(2, 2, 2, hexahedral=True)
+    a = CellSize(mesh)('+') * dS
+    A = assemble(a)
+    print(A)
+    a = CellVolume(mesh)('+') * dS
+    A = assemble(a)
+    print(A)
+    a = FacetArea(mesh)('+') * dS
+    A = assemble(a)
+    print(A)