Skip to content
Snippets Groups Projects
Commit 75dc33f3 authored by Isuru Fernando's avatar Isuru Fernando
Browse files

Fix tests

parent 82103baf
No related branches found
No related tags found
1 merge request!68Specify the PDE symbolically
...@@ -400,7 +400,7 @@ class LinearRecurrenceBasedExpansionTermsWrangler(ExpansionTermsWrangler): ...@@ -400,7 +400,7 @@ class LinearRecurrenceBasedExpansionTermsWrangler(ExpansionTermsWrangler):
n = nullspace(pde_mat) n = nullspace(pde_mat)
n = n[offset*iexpr:offset*(iexpr+1), :] n = n[offset*iexpr:offset*(iexpr+1), :]
idx = self.get_reduced_coeffs(n) idx = self.get_reduced_coeffs(n)
n = n[:,:len(idx)] n = n[:, :len(idx)]
s = solve_symbolic(n.T[:, idx], n.T) s = solve_symbolic(n.T[:, idx], n.T)
stored_identifiers = [mis[i] for i in idx] stored_identifiers = [mis[i] for i in idx]
else: else:
...@@ -417,7 +417,7 @@ class LinearRecurrenceBasedExpansionTermsWrangler(ExpansionTermsWrangler): ...@@ -417,7 +417,7 @@ class LinearRecurrenceBasedExpansionTermsWrangler(ExpansionTermsWrangler):
plog.done() plog.done()
logger.debug("number of Taylor coefficients was reduced from {orig} to {red}" print("number of Taylor coefficients was reduced from {orig} to {red}"
.format(orig=len(self.get_full_coefficient_identifiers()), .format(orig=len(self.get_full_coefficient_identifiers()),
red=len(stored_identifiers))) red=len(stored_identifiers)))
......
...@@ -412,8 +412,10 @@ def test_translations(ctx_getter, knl, local_expn_class, mpole_expn_class): ...@@ -412,8 +412,10 @@ def test_translations(ctx_getter, knl, local_expn_class, mpole_expn_class):
issubclass(local_expn_class, VolumeTaylorExpansionBase): issubclass(local_expn_class, VolumeTaylorExpansionBase):
# FIXME: Embarrassing--but we run out of memory for higher orders. # FIXME: Embarrassing--but we run out of memory for higher orders.
orders = [2, 3] orders = [2, 3]
else: elif isinstance(knl, StokesletKernel):
orders = [3, 4, 5] orders = [3, 4, 5]
else:
orders = [2, 3, 4]
nboxes = centers.shape[-1] nboxes = centers.shape[-1]
def eval_at(e2p, source_box_nr, rscale): def eval_at(e2p, source_box_nr, rscale):
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment