From bd6867980f70f518feea48e27c640781d2e5dd5d Mon Sep 17 00:00:00 2001 From: unalmis Date: Tue, 17 Sep 2024 04:03:04 -0400 Subject: [PATCH] Use interior nodes only for more fair test --- tests/test_integrals.py | 11 +++++------ 1 file changed, 5 insertions(+), 6 deletions(-) diff --git a/tests/test_integrals.py b/tests/test_integrals.py index 9d157a2466..63dcc8c431 100644 --- a/tests/test_integrals.py +++ b/tests/test_integrals.py @@ -910,14 +910,13 @@ class TestBounce1DQuadrature: [ (True, tanh_sinh(40), None), (False, tanh_sinh(20), None), - # Node density near boundary is (1−x²)⁻¹. + # Node density near boundary is 1/(1−x²). (True, leggauss(25), auto_sin), (True, _chebgauss1(30), auto_sin), # Lobatto nodes - (False, leggauss_lob(10), auto_sin), - (False, chebgauss2(10), auto_sin), - # Node density near boundary is (1−x²)⁻⁰ᐧ⁵. - (False, leggauss_lob(14), None), + (False, leggauss_lob(8, interior_only=True), auto_sin), + # Node density near boundary is 1/√(1−x²). + (False, leggauss_lob(13, interior_only=True), None), (False, chebgauss2(8), None), ], ) @@ -968,7 +967,7 @@ def test_bounce_quadrature(self, is_strong, quad, automorphism): automorphism, check=True, ) - result = bounce.integrate(integrand, pitch_inv, check=True, plot=True) + result = bounce.integrate(integrand, pitch_inv, check=True, plot=False) assert np.count_nonzero(result) == 1 np.testing.assert_allclose(result.sum(), truth, rtol=1e-4)