From d2adf02e1f2e2357049e8c7666d291eac44cf6a9 Mon Sep 17 00:00:00 2001 From: Thomas Gibson Date: Sun, 25 Apr 2021 19:53:41 -0500 Subject: [PATCH] Add debugging remark --- grudge/op.py | 7 +++++++ 1 file changed, 7 insertions(+) diff --git a/grudge/op.py b/grudge/op.py index 50b96ec9..d6513eda 100644 --- a/grudge/op.py +++ b/grudge/op.py @@ -367,6 +367,13 @@ def _apply_mass_operator(dcoll, dd_out, dd_in, vec): area_elements = area_element(actx, dcoll, dd=dd_in) data = [] + # FIXME: I think there is a bug in loopy's implementation + # of einsum; this code works when dispatching to numpy. But + # fails when using actx.einsum due to strides mismatch + # on the mass matrix... (error raised when attempting to + # invoke the loopy kernel) + # NOTE: This ONLY fails when in_grp != out_grp + # TODO: Need to debug the loopy program...? for in_grp, out_grp, ae_i, vec_i in zip( in_discr.groups, out_discr.groups, area_elements, vec): -- GitLab