diff options
Diffstat (limited to 'libm/e_exp.c')
| -rw-r--r-- | libm/e_exp.c | 25 | 
1 files changed, 25 insertions, 0 deletions
diff --git a/libm/e_exp.c b/libm/e_exp.c index ebf2e5af0..abfb17675 100644 --- a/libm/e_exp.c +++ b/libm/e_exp.c @@ -155,3 +155,28 @@ double attribute_hidden __ieee754_exp(double x)	/* default IEEE double exp */  	    return y*twom1000;  	}  } + +/* + * wrapper exp(x) + */ +#ifndef _IEEE_LIBM +double exp(double x) +{ +	static const double o_threshold =  7.09782712893383973096e+02; /* 0x40862E42, 0xFEFA39EF */ +	static const double u_threshold = -7.45133219101941108420e+02; /* 0xc0874910, 0xD52D3051 */ + +	double z = __ieee754_exp(x); +	if (_LIB_VERSION == _IEEE_) +		return z; +	if (isfinite(x)) { +		if (x > o_threshold) +			return __kernel_standard(x, x, 6); /* exp overflow */ +		if (x < u_threshold) +			return __kernel_standard(x, x, 7); /* exp underflow */ +	} +	return z; +} +#else +strong_alias(__ieee754_exp, exp) +#endif +libm_hidden_def(exp)  | 
