Re: [Bug-gsl] incorrect values by bessel_steed_array

[Hu Zhan]

hypot works. Thanks.

Hu

On Fri, 2005-02-11 at 07:49, Brian Gough wrote:
> Hu Zhan writes:
>  > There is one line in gsl_sf_bessel_jl_steed_array:
>  > 
>  >     W = x_inv / sqrt(FP*FP + F*F);
>  > 
>  > When LMAX is large, FP and F can be larger than 1e160, and cause W = 0.
>  > The option -mfpmath=387 (90-bit arithmetics internally) can make it work
>  > on the opteron, but not pentium.
> 
> Thanks for the explanation, I see the problem.  Could you try the
> following patch (to replace sqrt(a^2+b^2) by hypot) and let me know if
> it works.