Re: [Boost-users] Boost_1_54_0: boost::math::sph_bessel(n, x) crash at x=0 n>0

Subject: Re: [Boost-users] Boost_1_54_0: boost::math::sph_bessel(n, x) crash at x=0 n>0
From: Paul A. Bristow (pbristow_at_[hidden])
Date: 2013-10-18 05:03:08

> -----Original Message-----
> From: Boost-users [mailto:boost-users-bounces_at_[hidden]] On Behalf Of JE Campagne
> Sent: Friday, October 18, 2013 9:09 AM
> To: boost-users_at_[hidden]
> Subject: Re: [Boost-users] Boost_1_54_0: boost::math::sph_bessel(n, x) crash at x=0 n>0
> 
> Hi
> 
> 1) in your patch for T operator()() part, what's the "n" variable, looking at 1_54_0 version there
is no
> mention to this variable.

Ask the organ grinder ;-)

> 2) is there a way to get your current version of the bessel.hpp file ?

http://svn.boost.org/svn/boost/trunk/boost/math/special_functions/bessel.hpp

in case you want to try this now.

(There might be other dependencies but I guess not).

Revision: 86346
Author: johnmaddock
Date: 18 October 2013 09:56:42
Message:
Previous commit failed for types with an extended exponent range - use logarithms rather than
assuming the result is zero.
----
Modified : /trunk/boost/math/special_functions/bessel.hpp
Revision: 86343
Author: johnmaddock
Date: 17 October 2013 19:21:56
Message:
Fix for sph_bessel when v is large and x is small.
----
Modified : /trunk/boost/math/special_functions/bessel.hpp
Modified : /trunk/libs/math/test/test_bessel_j.hpp
Paul
---
Paul A. Bristow,
Prizet Farmhouse, Kendal LA8 8AB  UK
[phone]  07714330204
pbristow_at_[hidden]