From f14fa3fe4f208419a3123ec80b9ddd20a9ddf4b7 Mon Sep 17 00:00:00 2001 From: Thomas Gibson Date: Wed, 14 Apr 2021 19:45:43 -0500 Subject: [PATCH] Always check quadrature accuracy (not just for unisolvent case) --- grudge/symbolic/operators.py | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/grudge/symbolic/operators.py b/grudge/symbolic/operators.py index 07c41383..35a2d33a 100644 --- a/grudge/symbolic/operators.py +++ b/grudge/symbolic/operators.py @@ -682,12 +682,13 @@ class RefFaceMassOperator(ElementwiseLinearOperator): face ) + assert face_quadrature.exact_to >= n + m + # If the group has a nodal basis and is unisolvent, # we use the basis on the face to compute the face mass matrix if (isinstance(afgrp, ElementGroupWithBasis) and afgrp.space.space_dim == afgrp.nunit_dofs): - assert face_quadrature.exact_to >= n + m face_basis = afgrp.basis_obj() matrix[:, iface, :] = mp.nodal_mass_matrix_for_face( -- GitLab