From e0b8df396828bcea90d58c7118821f5c62ca4b69 Mon Sep 17 00:00:00 2001 From: Isuru Fernando Date: Tue, 7 Jul 2020 09:47:07 -0500 Subject: [PATCH] Fix typo and comment on max_mi --- sumpy/expansion/__init__.py | 3 ++- sumpy/expansion/local.py | 4 ++++ 2 files changed, 6 insertions(+), 1 deletion(-) diff --git a/sumpy/expansion/__init__.py b/sumpy/expansion/__init__.py index d3dfbccb..17264e8e 100644 --- a/sumpy/expansion/__init__.py +++ b/sumpy/expansion/__init__.py @@ -196,7 +196,8 @@ class ExpansionTermsWrangler(object): if self.max_mi is None: return res - return [mi for mi in res if all(mi[i] <= self.max_mi[i])] + return [mi for mi in res if + all(mi[i] <= self.max_mi[i] for i in range(self.dim))] def copy(self, **kwargs): new_kwargs = dict( diff --git a/sumpy/expansion/local.py b/sumpy/expansion/local.py index 6080f646..3acdfea0 100644 --- a/sumpy/expansion/local.py +++ b/sumpy/expansion/local.py @@ -240,6 +240,10 @@ class VolumeTaylorLocalExpansionBase(LocalExpansionBase): from sumpy.tools import add_mi + # Calculate a elementwise maximum multi-index because the number + # of multi-indices needed is much less than + # gnitstam(src_order + tgt order) when PDE conforming expansions + # are used. For full Taylor, there's no difference. max_mi = [0]*self.dim for i in range(self.dim): max_mi[i] = max(mi[i] for mi in -- GitLab