actually add the high-order quadrature rule referenced in docs

inducer · Oct 11, 2024 · 8322164 · 8322164
1 parent a449133
commit 8322164
Showing 1 changed file with 7 additions and 0 deletions.
diff --git a/meshmode/discretization/poly_element.py b/meshmode/discretization/poly_element.py
@@ -592,8 +592,15 @@ def __init__(self, mesh_el_group, order):
         unit_nodes_1d = legendre_gauss_lobatto_nodes(order)
         unit_nodes = mp.tensor_product_nodes([unit_nodes_1d]*mesh_el_group.dim)
 
+        self._quadrature_rule = mp.LegendreGaussLobattoTensorProductQuadrature(
+            order, mesh_el_group.dim)
+
         super().__init__(mesh_el_group, order, unit_nodes=unit_nodes)
 
+    @memoize_method
+    def quadrature_rule(self):
+        return self._quadrature_rule
+
     def discretization_key(self):
         return (type(self), self.dim, self.order)