
Boost Users : 
Subject: Re: [Boostusers] Boost_1_54_0: boost::math::sph_bessel(n, x) crash at x=0 n>0
From: JE Campagne (campagne_at_[hidden])
Date: 20131018 06:22:10
Your fix pass all my personnal tests. Thanks.
Sorry for the trouble, according to the limit close to zero.
Bye the way do you know when 1_55 version will be publicly available?
JeanEric
Le 18/10/2013 11:01, John Maddock a écrit :
>> 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.
>
> There are no changes to operator()(), the n variable appears in the
> next part of the diff in
>
> inline T sph_bessel_j_imp(unsigned n, T x, const Policy& pol)
>
>> 2) is there a way to get your current version of the bessel.hpp file ?
>
> Attached.
>
> Note I had to change the previous fix somewhat as it failed for types
> with extended exponent ranges (I realised that just as I was dozing
> off to sleep last night!)
>
>> That's true that for x close to zero but positive there are also
>> problems to the use of power series close to 0 is certainly the best
>> thing to do. I am astonished that this not yet used.
>
> Sorry I don't understand: the problem was occurring inside the power
> series for small x.
>
> John.
>
>
> _______________________________________________
> Boostusers mailing list
> Boostusers_at_[hidden]
> http://lists.boost.org/mailman/listinfo.cgi/boostusers
Boostusers list run by williamkempf at hotmail.com, kalb at libertysoft.com, bjorn.karlsson at readsoft.com, gregod at cs.rpi.edu, wekempf at cox.net