Fix gen-quad reference simplex

firedrakeproject · Jul 19, 2024 · 14a5948 · 14a5948
1 parent afc4bf9
commit 14a5948
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/FIAT/quadrature_schemes.py b/FIAT/quadrature_schemes.py
@@ -36,7 +36,7 @@
 # FIAT
 from FIAT.reference_element import (HEXAHEDRON, QUADRILATERAL, TENSORPRODUCT,
                                     TETRAHEDRON, TRIANGLE, UFCTetrahedron,
-                                    UFCTriangle, ufc_simplex, symmetric_simplex)
+                                    UFCTriangle, symmetric_simplex)
 
 
 def create_quadrature(ref_el, degree, scheme="default"):
@@ -334,8 +334,8 @@ def _kmv_lump_scheme(ref_el, degree):
 
 
 def gen_quad_simplex(dim):
-    ref_el = ufc_simplex(dim)
-    ref_el.vertices = ((0., 0., 0.), (1., 0., 0.), (0., 1., 1.), (1., 1., 1.))
+    ref_el = symmetric_simplex(dim)
+    ref_el.vertices = tuple(tuple(float(i > j) for j in range(dim)) for i in range(dim+1))
     return ref_el