Bessel function bug or limitation?

[Rodolphe Conan]

Hi.

I installed GSL 1.1.1 on Redhat 7.2
with
configure
make
make install

I compile the following program

#include <iostream>
#include <gsl/gsl_sf_bessel.h>

/*
Compiled with: g++ test_bjn.cc -o test_bjn -I/usr/local/include/ 
-L/usr/local/lib -lgsl -lgslcblas -lm

*/

int main() {

  cout << "besselj(2,15000)="  << gsl_sf_bessel_Jn(2,15000) << endl;

}

with g++ version 2.96

and obtained:
gsl: bessel.c:531: ERROR: error
Abort (core dumped)

Rodolphe

Brian Gough wrote:

>Rodolphe Conan writes:
> > > I need to compute bessel function for large arguments (for example: 
> > > gsl_sf_bessel_Jn(0,20000))
> > > and GSL crashes with the following error message:
> > > gsl: bessel.c:531: ERROR: error
> > > Abort (core dumped)
> > >
> > > I check the gsl_sf_bessel_Jn with 2 version of GSL 0.9.4 and 1.1.1:
> > > on GSL 0.9.4, GSL gave me the good value but with the wrong sign and
> > > now with the GSL 1.1.1 it crashes.
> > >
> > > Is it a bug or a known limitation of GSL?
> > >
>
>Hi,
>Please could you send more information:
>- an example program which reproduces the problem
>- the details of how you compiled/installed gsl (platform, compiler, etc).
>
>Thanks
>Brian Gough
>
>

-- 
Rodolphe Conan

European Southern Observatory
Instrumentation Division
Adaptive Optics Group
Karl-Schwarzschild-Str.  2
D-85748 Garching bei Muenchen
Germany

Tel: +49-89-32006541    Fax: +49-89-3202362
e-mail : rconan@eso.org