bug-gsl
[Top][All Lists]
Advanced

[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


reply via email to

[Prev in Thread] Current Thread [Next in Thread]