[Bug-gsl] problem with gsl_sf_exp()

[Thomas Weber]

Hi, 

it seems gsl_sf_exp can't handle large negative arguments:

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

int main(void)
{
  double a = 1125000.0;
  double b = 0.0, c = 0.0;

  b = exp(-a);
  printf("b: %f\n", b);
  
  c = gsl_sf_exp(-a);
  printf("c: %f\n", c);

  return 0;
}


Output:
b: 0.000000
gsl: exp.c:114: ERROR: underflow
Default GSL error handler invoked.
Aborted

system: Debian unstable (with gsl 1.6.92)

Please CC me, I'm not subscribed to the list

Regards
        Thomas