From 8ba98b79f67c233931a4e4b28c7ff0d34ab4b1da Mon Sep 17 00:00:00 2001 From: jzmaddock Date: Thu, 1 Jan 2026 10:48:31 +0000 Subject: [PATCH] Use test values from wolframapha --- test/test_igamma.hpp | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/test/test_igamma.hpp b/test/test_igamma.hpp index 1d3be9fd4..7de179c13 100644 --- a/test/test_igamma.hpp +++ b/test/test_igamma.hpp @@ -257,11 +257,11 @@ void test_spots(T) BOOST_CHECK_EQUAL(::boost::math::gamma_p(static_cast(1770), static_cast(1e-12)), 0); // - // Check that lgamma_q returns correct values + // Check that lgamma_q returns correct values with spot values calculated via wolframalpha log(Q[a, x]) // - BOOST_CHECK_CLOSE(::boost::math::lgamma_q(static_cast(5), static_cast(100)), static_cast(log(1.6139305336977304790405739225035685228527400976549e-37L)), tolerance); - BOOST_CHECK_CLOSE(::boost::math::lgamma_q(static_cast(22.5), static_cast(2000)), static_cast(-1883.4897732037716195918619632721L), tolerance * 10); // calculated via mpmath - BOOST_CHECK_CLOSE(::boost::math::lgamma_q(static_cast(501.2), static_cast(2000)), static_cast(-810.31461624182202285737730562687L), tolerance * 10); // calculated via mpmath + BOOST_CHECK_CLOSE(::boost::math::lgamma_q(static_cast(5), static_cast(100)), static_cast(-84.71697591169848944613823640968965801339401810393519310714864307L), tolerance); + BOOST_CHECK_CLOSE(::boost::math::lgamma_q(static_cast(22.5), static_cast(2000)), static_cast(-1883.489773203771543025750308264545743305089849873060383828767138L), tolerance * 10); + BOOST_CHECK_CLOSE(::boost::math::lgamma_q(static_cast(501.2), static_cast(2000)), static_cast(-810.2453406781655559126505101822969531699112391075198076300675402L), tolerance * 10); // // Coverage: //