[Help-gsl] function gsl_sf_exprel_e

[Daniel Rohe]

The gsl documentation says the following:

Function:

int gsl_sf_exprel_e (double x, gsl_sf_result * result)
These routines compute the quantity (\exp(x)-1)/x using an algorithm that is 
accurate for small x.
For small x the algorithm is based on the expansion (\exp(x)-1)/x = 1 + x/2 + 
x^2/(2*3) + x^3/(2*3*4) + \dots.


My question is: when I use this function, shall I dustinguish between small and non-small values of x myself, or does the function handle this 
internally, i.e. does it call the usual exp() function if x is larger than a certain value?

thanks for any help,

daniel