Lines Matching refs:expm1
29 /* expm1(x)
41 * 2. Approximating expm1(r) by a special rational function on
67 * expm1(r) = exp(r)-1 is then computed by the following
71 * expm1(r) = r + --- + --- * [--------------------]
75 * expm1(r+c) = expm1(r) + c + expm1(r)*c
76 * ~ expm1(r) + c + r*c
78 * expm1(r+c). Now rearrange the term to avoid optimization
82 * expm1(r+c)~r - ({r*(--- * [--------------------]-c)-c} - --- )
87 * 3. Scale back to obtain expm1(x):
89 * expm1(x) = either 2^k*[expm1(r)+1] - 1
90 * = or 2^k*[expm1(r) + (1-2^-k)]
94 * (B). To achieve maximum accuracy, we compute expm1(x) by
105 * expm1(INF) is INF, expm1(NaN) is NaN;
106 * expm1(-INF) is -1, and
107 * for finite argument, only expm1(0)=0 is exact.
115 * if x > 7.09782712893383973096e+02 then expm1(x) overflow
140 expm1 (double x)
292 } /* expm1 */