From a9b45e35ec5498b6807992b59d24a0fa7edfb41f Mon Sep 17 00:00:00 2001
From: Andreas Kloeckner <inform@tiker.net>
Date: Fri, 8 Jan 2016 15:38:58 -0600
Subject: [PATCH] Sketchy face jacobian sign flip
---
grudge/symbolic/mappers/__init__.py | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/grudge/symbolic/mappers/__init__.py b/grudge/symbolic/mappers/__init__.py
index d7f3e71c..745d56be 100644
--- a/grudge/symbolic/mappers/__init__.py
+++ b/grudge/symbolic/mappers/__init__.py
@@ -524,7 +524,7 @@ class GlobalToReferenceMapper(CSECachingMapperMixin, IdentityMapper):
1/jac_in * self.rec(expr.field))
elif isinstance(expr.op, op.FaceMassOperator):
- jac_in_surf = sym.area_element(self.ambient_dim, self.dim - 1,
+ jac_in_surf = - sym.area_element(self.ambient_dim, self.dim - 1,
dd=expr.op.dd_in)
return op.RefFaceMassOperator(expr.op.dd_in, expr.op.dd_out)(
jac_in_surf * self.rec(expr.field))
--
GitLab