Boost C++ Libraries: Ticket #12961: boost::math::cyl_bessel_i crashes program.

attachment set

theivorytower@… — Wed, 12 Apr 2017 05:07:26 GMT

attachment → bessel.cpp

Sample code that crashes.

John Maddock — Thu, 13 Apr 2017 09:50:30 GMT

I'm unable to reproduce this, either on Windows (MSVC) or on Ubuntu Linux with GCC.

What's the OS?

You can help me a little by stepping through your sample program in the debugger - you should get to:

bessel.hpp:221 return bessel_i_small_z_series(v, x, pol);

And then stepping into that should lead to the rather trivial evaluation in:

bessel_ik.hpp:78 return prefix * result;

With

prefix = 0.40373565439688797 and result = 1.0000408337880704

theivorytower@… — Thu, 13 Apr 2017 17:03:40 GMT

Thanks for the reply. I'm using Arch Linux.

When I get to bessel_ik.hpp:78 return prefix * result, the variables are prefix = inf and result = 1.000040833788070502945698192487356.

When I step in a little further, I find the reason prefix was evaluated to be infinite was because pow(x / 2, v) was evaluated to be infinite! So this is must be a bug with libstdc++. Sorry this is not a bug with boost. FYI with my libstdc++ I find std::pow(x, y) ALWAYS returns infinite whenever y is long double. I will dig further into libstdc++...

John Maddock — Sat, 15 Apr 2017 11:20:06 GMT

As a workaround you could define BOOST_MATH_NO_LONG_DOUBLE_MATH_FUNCTIONS which would disable all use of the long double std lib functions internally.

status changed; resolution set

John Maddock — Tue, 31 Jul 2018 17:39:20 GMT

status new → closed
resolution → invalid