diff --git a/grudge/symbolic/operators.py b/grudge/symbolic/operators.py index 217fcf02ed262c2091f8188dfb952f4a70799e3d..7009cdd835f442fa155784c92e992cba3575d618 100644 --- a/grudge/symbolic/operators.py +++ b/grudge/symbolic/operators.py @@ -551,10 +551,10 @@ def norm(p, arg, dd=None): return sym.CFunction("sqrt")(norm_squared) elif p == np.Inf: - result = sym.NodalMax()(sym.CFunction("fabs")(arg)) + result = sym.NodalMax(dd_in=dd)(sym.CFunction("fabs")(arg)) from pymbolic.primitives import Max - if isinstance(norm_squared, np.ndarray): + if isinstance(result, np.ndarray): from functools import reduce result = reduce(Max, result)