From 8e3731b3ed3872ab796f28c51a372897ec939cfc Mon Sep 17 00:00:00 2001 From: haakoek Date: Mon, 4 May 2026 14:03:28 +0200 Subject: [PATCH] Precision of int exp(-a*x^2)x^k dx seems to depend on SciPy version. Lower atol value. --- tests/test_quadrature.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/tests/test_quadrature.py b/tests/test_quadrature.py index 007939d..a78696b 100644 --- a/tests/test_quadrature.py +++ b/tests/test_quadrature.py @@ -87,7 +87,7 @@ def test_integrate_xk_exp_minus_a_x2(a): integral_quadrature = np.dot(w, f) integral_exact = 0.5 * a ** (-(k + 1) / 2) * gamma((k + 1) / 2) np.testing.assert_allclose( - integral_quadrature, integral_exact, rtol=0.0, atol=1e-12 + integral_quadrature, integral_exact, rtol=0.0, atol=1e-10 ) test_integrate_xk_exp_minus_a_x2(1.0)