From cc2e8660281f0ac4801e74803f8fd78201f55437 Mon Sep 17 00:00:00 2001
From: Thomas Gibson <gibsonthomas1120@hotmail.com>
Date: Mon, 10 May 2021 14:42:55 -0500
Subject: [PATCH] Fix face mass application
---
grudge/op.py | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/grudge/op.py b/grudge/op.py
index 415b1435..4e80ed46 100644
--- a/grudge/op.py
+++ b/grudge/op.py
@@ -718,8 +718,8 @@ def _apply_face_mass_operator(dcoll, dd, vec):
0<=idof<nvol_nodes and
0<=j<nface_nodes}""",
"""
- result[iel,idof] = sum(f, sum(j, mat[idof, f, j] \
- * jac_surf[f, iel, j] \
+ result[iel,idof] = sum(f, sum(j, mat[idof, f, j] \
+ * jac_surf[f, iel, j] \
* vec[f, iel, j]))
""",
name="face_mass"
@@ -737,7 +737,7 @@ def _apply_face_mass_operator(dcoll, dd, vec):
vol_element_group=volgrp,
dtype=dtype
)
- input_view = result[afgrp.index].reshape(
+ input_view = vec[afgrp.index].reshape(
nfaces, volgrp.nelements, afgrp.nunit_dofs
)
jac_surf = surf_area_elements[afgrp.index].reshape(
--
GitLab