boost_mirror/math
1
0
Fork 0
mirror of https://github.com/boostorg/math.git synced 2025-05-11 21:33:52 +00:00
Code Issues Packages Projects Releases Wiki Activity

Coverage: Improve gamma.hpp.

Browse Source
This commit is contained in:
jzmaddock 2025-04-22 10:13:08 +01:00
parent ca6e976716
commit 4a3d3bc66e
2 changed files with 3 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
include/boost/math/special_functions/gamma.hpp
View File

@ -131,11 +131,6 @@ BOOST_MATH_GPU_ENABLED T gamma_imp_final(T z, const Policy& pol, const Lanczos&
if(z <= 0)
{
if(floor(z) == z)
{
return policies::raise_pole_error<T>(function, "Evaluation of tgamma at a negative integer %1%.", z, pol);
}
// shift z to > 1:
while(z < 0)
{
@ -225,7 +220,8 @@ BOOST_MATH_GPU_ENABLED BOOST_MATH_FORCEINLINE T gamma_imp(T z, const Policy& pol
{
if ((fabs(result) < 1) && (tools::max_value<T>() * fabs(result) < boost::math::constants::pi<T>()))
{
return policies::raise_overflow_error<T>(function, nullptr, pol);
return policies::raise_overflow_error<T>(function, nullptr, pol); // LCOV_EXCL_LINE MP only.
}
}
else
@ -1302,12 +1298,6 @@ template <class T, class Policy>
BOOST_MATH_GPU_ENABLED T gamma_incomplete_imp_final(T a, T x, bool normalised, bool invert,
const Policy& pol, T* p_derivative)
{
constexpr auto function = "boost::math::gamma_p<%1%>(%1%, %1%)";
if(a <= 0)
return policies::raise_domain_error<T>(function, "Argument a to the incomplete gamma function must be greater than zero (got a=%1%).", a, pol);
if(x < 0)
return policies::raise_domain_error<T>(function, "Argument x to the incomplete gamma function must be >= 0 (got x=%1%).", x, pol);
BOOST_MATH_STD_USING
typedef typename lanczos::lanczos<T, Policy>::type lanczos_type;

1
test/Jamfile.v2
View File

@ -575,6 +575,7 @@ test-suite special_fun :
[ run test_tgamma_for_issue396_part1.cpp /boost/test//boost_unit_test_framework : : : $(float128_type) <toolset>gcc-mingw:<cxxflags>-Wa,-mbig-obj <debug-symbols>off <toolset>msvc:<cxxflags>/bigobj [ check-target-builds ../config//is_cygwin_run "Cygwin CI run" : <build>no ] ]
[ run test_tgamma_for_issue396_part2.cpp /boost/test//boost_unit_test_framework : : : <toolset>gcc-mingw:<cxxflags>-Wa,-mbig-obj <debug-symbols>off <toolset>msvc:<cxxflags>/bigobj [ check-target-builds ../config//is_cygwin_run "Cygwin CI run" : <build>no ] ]
[ run test_tgamma_ratio.cpp test_instances//test_instances pch_light /boost/test//boost_unit_test_framework ]
[ run test_tgamma_eh.cpp /boost/test//boost_unit_test_framework ]
[ run test_trig.cpp test_instances//test_instances pch_light /boost/test//boost_unit_test_framework ]
[ run test_zeta.cpp /boost/test//boost_unit_test_framework test_instances//test_instances pch_light ]
[ run test_sinc.cpp /boost/test//boost_unit_test_framework pch_light ]