From b4539e3e6fe3b283117b9edcf753857a42fd5f94 Mon Sep 17 00:00:00 2001
From: Matthew Smith <mjsmith6@illinois.edu>
Date: Tue, 15 Dec 2020 12:47:27 -0600
Subject: [PATCH] change default dd_out value in MassOperatorBase
---
grudge/symbolic/operators.py | 8 ++++----
1 file changed, 4 insertions(+), 4 deletions(-)
diff --git a/grudge/symbolic/operators.py b/grudge/symbolic/operators.py
index 5155cbb6..b6db10eb 100644
--- a/grudge/symbolic/operators.py
+++ b/grudge/symbolic/operators.py
@@ -486,7 +486,7 @@ class MassOperatorBase(Operator):
if dd_in is None:
dd_in = prim.DD_VOLUME
if dd_out is None:
- dd_out = dd_in
+ dd_out = prim.DD_VOLUME
super().__init__(dd_in, dd_out)
@@ -705,7 +705,7 @@ def integral(arg, dd=None):
return NodalSum(dd)(
arg * prim.cse(
- MassOperator(dd_in=dd)(prim.Ones(dd)),
+ MassOperator(dd_in=dd, dd_out=dd)(prim.Ones(dd)),
"mass_quad_weights",
prim.cse_scope.DISCRETIZATION))
@@ -764,7 +764,7 @@ def h_max_from_volume(ambient_dim, dim=None, dd=None):
return NodalMax(dd_in=dd)(
ElementwiseSumOperator(dd)(
- MassOperator(dd_in=dd)(prim.Ones(dd))
+ MassOperator(dd_in=dd, dd_out=dd)(prim.Ones(dd))
)
)**(1.0/dim)
@@ -785,7 +785,7 @@ def h_min_from_volume(ambient_dim, dim=None, dd=None):
return NodalMin(dd_in=dd)(
ElementwiseSumOperator(dd)(
- MassOperator(dd_in=dd)(prim.Ones(dd))
+ MassOperator(dd_in=dd, dd_out=dd)(prim.Ones(dd))
)
)**(1.0/dim)
--
GitLab