From fc174d29da7690e820c6efcfd86ac28cfe3cf1b5 Mon Sep 17 00:00:00 2001 From: Andreas Kloeckner Date: Thu, 10 Aug 2017 15:12:13 -0500 Subject: [PATCH 1/2] Use infinity norm to measure FMM error --- test/test_fmm.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/test/test_fmm.py b/test/test_fmm.py index c85759fe..62210e6a 100644 --- a/test/test_fmm.py +++ b/test/test_fmm.py @@ -225,7 +225,7 @@ def test_sumpy_fmm(ctx_getter, knl, local_expn_class, mpole_expn_class): pot = pot.get() ref_pot = ref_pot.get() - rel_err = la.norm(pot - ref_pot) / la.norm(ref_pot) + rel_err = la.norm(pot - ref_pot, np.inf) / la.norm(ref_pot, np.inf) logger.info("order %d -> relative l2 error: %g" % (order, rel_err)) pconv_verifier.add_data_point(order, rel_err) -- GitLab From 1c3a76f4e22c840776398e3e47d72883a9b8acc2 Mon Sep 17 00:00:00 2001 From: Andreas Kloeckner Date: Thu, 10 Aug 2017 20:38:34 -0500 Subject: [PATCH 2/2] Mildly increase FMM test scale to address artificially low errors --- test/test_fmm.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/test/test_fmm.py b/test/test_fmm.py index 62210e6a..0d8f36d8 100644 --- a/test/test_fmm.py +++ b/test/test_fmm.py @@ -111,8 +111,8 @@ def test_sumpy_fmm(ctx_getter, knl, local_expn_class, mpole_expn_class): ctx = ctx_getter() queue = cl.CommandQueue(ctx) - nsources = 500 - ntargets = 50 + nsources = 1000 + ntargets = 300 dtype = np.float64 from boxtree.tools import ( -- GitLab