This is the mail archive of the gsl-discuss@sources.redhat.com mailing list for the GSL project.

Index Nav: [Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav: [Date Prev] [Date Next] [Thread Prev] [Thread Next]
Other format: [Raw text]

Re: Bessel function bug or limitation?

Rodolphe Conan writes:
 >   cout << "besselj(2,15000)="  << gsl_sf_bessel_Jn(2,15000) << endl;
 > and obtained:
 > gsl: bessel.c:531: ERROR: error
 > Abort (core dumped)


Thanks, I can reproduce that error. 
It looks like one of the continued fraction
computations is breaking down there.

regards

Brian Gough

 > >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?
 > > > >
Index Nav: [Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav: [Date Prev] [Date Next] [Thread Prev] [Thread Next]