[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.
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.