From 6f70f21b104cba5a38c6642d314d313aaa1450d8 Mon Sep 17 00:00:00 2001
From: Isuru Fernando <isuruf@gmail.com>
Date: Tue, 7 Jul 2020 09:21:02 -0500
Subject: [PATCH] Comment on expansion order less than PDE order case
---
sumpy/expansion/__init__.py | 2 ++
1 file changed, 2 insertions(+)
diff --git a/sumpy/expansion/__init__.py b/sumpy/expansion/__init__.py
index ea70b0eb..eeaaa650 100644
--- a/sumpy/expansion/__init__.py
+++ b/sumpy/expansion/__init__.py
@@ -416,6 +416,8 @@ class LinearPDEBasedExpansionTermsWrangler(ExpansionTermsWrangler):
pde_dict = self.get_pde().eq
for ident in pde_dict.keys():
if ident not in coeff_ident_enumerate_dict:
+ # Order of the expansion is less than the order of the PDE.
+ # In that case, the compression matrix is the identity matrix
coeff_matrix = defaultdict(list)
reconstruct_matrix = []
for i in range(len(mis)):
--
GitLab