From b2cbd22fa38c930bec721cce55b4b3d9023924aa Mon Sep 17 00:00:00 2001
From: Thomas Gibson <gibsonthomas1120@hotmail.com>
Date: Sun, 25 Apr 2021 15:54:49 -0500
Subject: [PATCH] Push ref mass matrix to the array context
---
grudge/op.py | 4 +++-
1 file changed, 3 insertions(+), 1 deletion(-)
diff --git a/grudge/op.py b/grudge/op.py
index 92de47b9..f401be68 100644
--- a/grudge/op.py
+++ b/grudge/op.py
@@ -348,7 +348,9 @@ def reference_mass_matrix(actx, out_element_group, in_element_group):
def get_ref_mass_mat(out_grp, in_grp):
if out_grp == in_grp:
from meshmode.discretization.poly_element import mass_matrix
- return mass_matrix(in_grp)
+ return actx.freeze(
+ actx.from_numpy(mass_matrix(in_grp))
+ )
from modepy import vandermonde
basis = out_grp.basis_obj()
--
GitLab