Lines Matching defs:expm1

30 #pragma weak __expm1 = expm1
34  * expm1(x)
46  *   2. Approximating expm1(r) by a special rational function on
72  *	expm1(r) = exp(r)-1 is then computed by the following
76  *	      expm1(r) = r + --- + --- * [--------------------]
80  *		expm1(r+c) = expm1(r) + c + expm1(r)*c
81  *			   ~ expm1(r) + c + r*c
83  *	expm1(r+c). Now rearrange the term to avoid optimization
87  *	 expm1(r+c)~r - ({r*(--- * [--------------------]-c)-c} - --- )
92  *   3. Scale back to obtain expm1(x):
94  *	   expm1(x) = either 2^k*[expm1(r)+1] - 1
95  *		    = or     2^k*[expm1(r) + (1-2^-k)]
99  *	(B). To achieve maximum accuracy, we compute expm1(x) by
110  *	expm1(INF) is INF, expm1(NaN) is NaN;
111  *	expm1(-INF) is -1, and
112  *	for finite argument, only expm1(0)=0 is exact.
120  *	    if x >  7.09782712893383973096e+02 then expm1(x) overflow
141 /* scaled coefficients related to expm1 */
162 expm1(double x) {