[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Bug-gsl] Possible bug in error handling of gsl_sf_bessel_Knu_e
From: |
Johannes Marbach |
Subject: |
[Bug-gsl] Possible bug in error handling of gsl_sf_bessel_Knu_e |
Date: |
Fri, 29 Jul 2011 13:47:21 +0200 |
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
address@hidden gsl]$ gcc -lgsl -lgsl -lgslcblas -lm bessel.c
address@hidden 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] Possible bug in error handling of gsl_sf_bessel_Knu_e,
Johannes Marbach <=