This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
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?
> > > >