Hi folks!
I'm having some trouble with gsl_sf_bessel_il_scaled(l,x) for large l / small x.
The problem manifests itself with the error message:
gsl: gamma.c:1454: ERROR: underflow
Default GSL error handler invoked.
Aborted
This results e.g. from trying to compute gsl_sf_bessel_il_scaled(100, 0.01)
Thus my question:
Can I find documentation that describes for which values of (l,x) for which the
function will compute a value successfully?
On wikipedia I found an asymptotic expression that works fine as long as
0 < x << sqrt(l+1), so I can use that instead if GSL only fails in cases when
that approximation is OK to my accuracy.
But it would be nice if I don't have to experiment to determine for which (l,x)
I can run the GSL function.
Or does anyone know of a more "robust" bessel function, e.g. in a different C /
C++ library?
Thanks!
Torquil Sørensen
_______________________________________________
Help-gsl mailing list
[email protected]
http://lists.gnu.org/mailman/listinfo/help-gsl
