Hi.
When I'm calling gsl_sf_bessel_Knu_e with rather large arguments the program
immediately invokes the default error handler without allowing me to check
the return code. As far as I understood this isn't the intended behaviour,
or is it?
Here's my sample program's source code
#include <stdio.h>
#include <gsl/gsl_sf_bessel.h>
int main () {
gsl_sf_result res;
int rc = gsl_sf_bessel_Knu_e (2.5, 2000, &res);
printf ("%i\n", rc);
printf ("%e\n", res.val);
printf ("%e\n", res.err);
return 0;
}
After compiling and running the program I receive
[jm@elroy gsl]$ gcc -lgsl -lgsl -lgslcblas -lm bessel.c
[jm@elroy gsl]$ ./a.out
gsl: exp.c:257: ERROR: underflow
Default GSL error handler invoked.
Aborted
I'm using GSL 1.15 on Arch Linux.
Regards
Johannes
_______________________________________________
Bug-gsl mailing list
[email protected]
https://lists.gnu.org/mailman/listinfo/bug-gsl
