Re: [Help-gsl] GSL 1.16: gsl_sf_bessel_jl failure

[Patrick Alken]

Hi, looks like a bug, possibly a duplicate of bug #37209. Could you send 
along a simple test program which demonstrates the problem? Thanks

Patrick

On 09/16/2013 07:27 AM, JE Campagne wrote:
> Dear developpers
>
> I am running GSL 1.16 with such architecture
>
> Linux opera3.lal.in2p3.fr 2.6.32-220.13.1.el6.x86_64 #1 SMP Tue Apr 17
> 15:16:22 CDT 2012 x86_64 x86_64 x86_64 GNU/Linux
>
> and  gcc compilor version gcc (GCC) 4.4.5 20110214
>
> I have a problem both with
> gsl_sf_bessel_jl_e(int l, double x, gsl_sf_result* r)
> or
> gsl_sf_bessel_jl (int l, double x)
>
> routine that fails to compute with l>~256 and produce such message
>
> gsl: gamma.c:1467: ERROR: underflow
> Default GSL error handler invoked.
> Abort (core dumped)
>
> Is there a way to avoid such features or is it a bug?
>
> Best
> J.E Campagne