[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [Bug-gsl] incorrect values by bessel_steed_array

From: Hu Zhan
Subject: Re: [Bug-gsl] incorrect values by bessel_steed_array
Date: Fri, 11 Feb 2005 09:53:19 -0800

hypot works. Thanks.


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.

reply via email to

[Prev in Thread] Current Thread [Next in Thread]