diff --git a/meshmode/discretization/poly_element.py b/meshmode/discretization/poly_element.py index c15a1b18a32bd2ec0698eccdc538ae49e3d94243..bbaff473e56362e18f7a9a82e933062a208cff36 100644 --- a/meshmode/discretization/poly_element.py +++ b/meshmode/discretization/poly_element.py @@ -25,13 +25,9 @@ THE SOFTWARE. import numpy as np #import numpy.linalg as la -from pytools import memoize_method, memoize_method_nested -from meshmode.discretization import Discretization +from pytools import memoize_method from meshmode.mesh import SimplexElementGroup as _MeshSimplexElementGroup -import pyopencl as cl - -import loopy as lp import modepy as mp __doc__ = """ @@ -124,7 +120,9 @@ class PolynomialWarpAndBlendElementGroup(PolynomialSimplexElementGroupBase): @property @memoize_method def weights(self): - raise NotImplementedError() + return np.dot( + mp.mass_matrix(self.basis(), self.unit_nodes), + np.ones(len(self.basis()))) # }}}