\frac{e^{x} - 1}{x}\begin{array}{l}
\mathbf{if}\;x \le -1.31978793712926563 \cdot 10^{-4}:\\
\;\;\;\;\frac{\frac{\mathsf{fma}\left(-1, 1, e^{x + x}\right)}{e^{x} \cdot e^{x} - 1 \cdot 1}}{\frac{x}{e^{x} - 1}}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\frac{1}{6}, {x}^{2}, \mathsf{fma}\left(\frac{1}{2}, x, 1\right)\right)\\
\end{array}double code(double x) {
return ((exp(x) - 1.0) / x);
}
double code(double x) {
double VAR;
if ((x <= -0.00013197879371292656)) {
VAR = ((fma(-1.0, 1.0, exp((x + x))) / ((exp(x) * exp(x)) - (1.0 * 1.0))) / (x / (exp(x) - 1.0)));
} else {
VAR = fma(0.16666666666666666, pow(x, 2.0), fma(0.5, x, 1.0));
}
return VAR;
}




Bits error versus x
Results
| Original | 38.8 |
|---|---|
| Target | 39.1 |
| Herbie | 0.3 |
if x < -0.00013197879371292656Initial program 0.0
rmApplied flip--0.0
Simplified0.0
rmApplied flip-+0.0
Applied associate-/r/0.0
Applied associate-/l*0.0
if -0.00013197879371292656 < x Initial program 59.8
Taylor expanded around 0 0.5
Simplified0.5
Final simplification0.3
herbie shell --seed 2020102 +o rules:numerics
(FPCore (x)
:name "Kahan's exp quotient"
:precision binary64
:herbie-target
(if (and (< x 1) (> x -1)) (/ (- (exp x) 1) (log (exp x))) (/ (- (exp x) 1) x))
(/ (- (exp x) 1) x))