[Bug-gsl] Spherical Bessel function failure

[Jakub Benda]

The spherical Bessel function returns NaN for some (large) inputs without
setting the error code.
Example program:

--------------------------------------------------------------------------------
#include <stdio.h>
#include <gsl/gsl_sf.h>

int main(void)
{
gsl_sf_result res;
int err = gsl_sf_bessel_jl_e(364, 36.62, &res);

printf("err = %d\n", err);
printf("res.val = %g\n", res.val);
printf("res.err = %g\n", res.err);

return 0;
}
--------------------------------------------------------------------------------

Output:

--------------------------------------------------------------------------------
err = 0
res.val = -nan
res.err = nan
--------------------------------------------------------------------------------

Jakub