-----Original Message----- From: Boost-users [mailto:boost-users-bounces@lists.boost.org] On Behalf Of JE Campagne Sent: Friday, October 18, 2013 9:09 AM To: boost-users@lists.boost.org 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 +44 1539 561830 07714330204 pbristow@hetp.u-net.com