[Bug-gsl] Problem with function gsl_sf_exp in gsl-1.12

[Pierrick Bruneau]

Hello,

I had little trouble with the function (double) gsl_sf_exp(double) :  
when a pretty low value is passed as parameter (empirically : lower  
than -750.0), this function causes a crash in the calling program.

This is rather surprising as exp is theorically defined for arbitrary  
low values.
I overcame the problem with a little patch block,

if(val < -700.0) {
        val = 0.0;
} else {
        val = gsl_sf_exp(val);
}

But I thought it was relevant to notice it.

Pierrick Bruneau

pbruneau.vcf
Description: Text Data