From cd66b5b54d64752b5878cf335ef36c48d27eda70 Mon Sep 17 00:00:00 2001
From: Andreas Kloeckner <inform@tiker.net>
Date: Thu, 30 Jul 2020 16:00:50 -0500
Subject: [PATCH] Eager inf-norm: actually use inf norm in recursion
---
grudge/eager.py | 4 +++-
1 file changed, 3 insertions(+), 1 deletion(-)
diff --git a/grudge/eager.py b/grudge/eager.py
index a1385685..566f3cd1 100644
--- a/grudge/eager.py
+++ b/grudge/eager.py
@@ -174,7 +174,9 @@ class EagerDGDiscretization(DGDiscretizationWithBoundaries):
self.norm(vec[idx])**2
for idx in np.ndindex(vec.shape))**0.5
elif p == np.inf:
- return max(self.norm(vec[idx]) for idx in np.ndindex(vec.shape))
+ return max(
+ self.norm(vec[idx], np.inf)
+ for idx in np.ndindex(vec.shape))
else:
raise ValueError("unsupported norm order")
--
GitLab