boost/math
2
0
Fork 0
mirror of https://github.com/boostorg/math.git synced 2026-01-19 04:22:09 +00:00
Code Issues Projects Releases Wiki Activity

Merge pull request #805 from mckib2/bessel_ik-overflow

Browse Source
This commit is contained in:
Matt Borland
2022-07-20 09:38:52 -07:00
committed by GitHub
parent 8f5520cea5 b0096e649f
commit d4aacdd30f
1 changed files with 9 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

10
include/boost/math/special_functions/detail/bessel_ik.hpp
View File

@@ -377,7 +377,15 @@ int bessel_ik(T v, T x, T* I, T* K, int kind, const Policy& pol)
for (k = 1; k <= n; k++) // forward recurrence for K
{
T fact = 2 * (u + k) / x;
if((tools::max_value<T>() - fabs(prev)) / fact < fabs(current))
// Check for overflow: if (max - |prev|) / fact > max, then overflow
// (max - |prev|) / fact > max
// max * (1 - fact) > |prev|
// if fact < 1: safe to compute overflow check
// if fact >= 1: won't overflow
const bool will_overflow = (fact < 1)
? tools::max_value<T>() * (1 - fact) > fabs(prev)
: false;
if(!will_overflow && ((tools::max_value<T>() - fabs(prev)) / fact < fabs(current)))
{
prev /= current;
scale /= current;