From 4bf68b6f9f2b92b11150e6148b6c03b3a0e2bf07 Mon Sep 17 00:00:00 2001
From: Andreas Kloeckner <inform@tiker.net>
Date: Wed, 27 Sep 2017 11:10:41 -0500
Subject: [PATCH] Fix curve-pot example [ci skip]
---
examples/curve-pot.py | 4 +++-
1 file changed, 3 insertions(+), 1 deletion(-)
diff --git a/examples/curve-pot.py b/examples/curve-pot.py
index f4212208..4b0a5d44 100644
--- a/examples/curve-pot.py
+++ b/examples/curve-pot.py
@@ -161,6 +161,7 @@ def draw_pot_figure(aspect_ratio,
evt, (y,) = lpot(queue, native_curve.pos, ovsmp_curve.pos,
centers,
[xovsmp * ovsmp_curve.speed * ovsmp_weights],
+ expansion_radii=np.ones(centers.shape[1]),
**lpot_kwargs)
return y
@@ -185,6 +186,7 @@ def draw_pot_figure(aspect_ratio,
evt, (curve_pot,) = lpot(queue, native_curve.pos, ovsmp_curve.pos,
centers,
[ovsmp_density * ovsmp_curve.speed * ovsmp_weights],
+ expansion_radii=np.ones(centers.shape[1]),
**lpot_kwargs)
# }}}
@@ -263,7 +265,7 @@ def draw_pot_figure(aspect_ratio,
if __name__ == "__main__":
- draw_pot_figure(aspect_ratio=1, nsrc=100, novsmp=100, helmholtz_k=(15+4j)*0.3,
+ draw_pot_figure(aspect_ratio=1, nsrc=100, novsmp=100, helmholtz_k=(35+4j)*0.3,
what_operator="D", what_operator_lpot="D", force_center_side=1)
# pt.savefig("eigvals-ext-nsrc100-novsmp100.pdf")
--
GitLab