From f0dfd4fc9fa4023e938c4ca91356cc40e9c13be4 Mon Sep 17 00:00:00 2001
From: Isuru Fernando <isuruf@gmail.com>
Date: Mon, 12 Mar 2018 10:09:53 -0500
Subject: [PATCH] Remove use of ID
---
sumpy/expansion/__init__.py | 11 +++--------
1 file changed, 3 insertions(+), 8 deletions(-)
diff --git a/sumpy/expansion/__init__.py b/sumpy/expansion/__init__.py
index 5d8c25d6..afdd9248 100644
--- a/sumpy/expansion/__init__.py
+++ b/sumpy/expansion/__init__.py
@@ -350,13 +350,8 @@ class LinearRecurrenceBasedDerivativeWrangler(DerivativeWrangler):
pde_mat = np.array(pde_mat, dtype=np.float64)
n = nullspace(pde_mat, atol=tol)
idx = self.get_reduced_coeffs()
- if idx is None or len(idx) < n.shape[1]:
- from scipy.linalg.interpolative import interp_decomp
- k, idx, proj = interp_decomp(n.T, tol)
- s = np.hstack([np.eye(k), proj])[:, np.argsort(idx)]
- idx = idx[:k]
- else:
- s = np.linalg.solve(n.T[:, idx], n.T)
+ assert len(idx) >= n.shape[1]
+ s = np.linalg.solve(n.T[:, idx], n.T)
stored_identifiers = [mis[i] for i in idx]
else:
s = np.eye(len(mis))
@@ -393,7 +388,7 @@ class LinearRecurrenceBasedDerivativeWrangler(DerivativeWrangler):
indicates the reduced coefficients are unknown and will be calculated
using an Interpolative Decomposition of the PDE matrix
"""
- return None
+ raise NotImplementedError
def get_derivative_taker(self, expr, var_list):
from sumpy.tools import NewLinearRecurrenceBasedMiDerivativeTaker
--
GitLab