From 83221641540c00af04d135d3395100ef6aa20c76 Mon Sep 17 00:00:00 2001 From: Addison Alvey-Blanco Date: Fri, 11 Oct 2024 00:09:31 -0500 Subject: [PATCH] actually add the high-order quadrature rule referenced in docs --- meshmode/discretization/poly_element.py | 7 +++++++ 1 file changed, 7 insertions(+) diff --git a/meshmode/discretization/poly_element.py b/meshmode/discretization/poly_element.py index 9b830e65..00b32cb3 100644 --- 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)