\frac{e^{x} - 1}{x}\begin{array}{l}
\mathbf{if}\;x \le -1.03925223585416463 \cdot 10^{-4}:\\
\;\;\;\;\frac{\frac{\sqrt[3]{\log \left(e^{\mathsf{fma}\left(-1, 1, e^{x + x}\right)}\right)} \cdot \sqrt[3]{\log \left(e^{\mathsf{fma}\left(-1, 1, e^{x + x}\right)}\right)}}{\sqrt[3]{e^{x} + 1} \cdot \sqrt[3]{e^{x} + 1}}}{\frac{x \cdot \sqrt[3]{e^{x} + 1}}{\sqrt[3]{\mathsf{fma}\left(-1, 1, e^{x + x}\right)}}}\\
\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 temp;
if ((x <= -0.00010392522358541646)) {
temp = (((cbrt(log(exp(fma(-1.0, 1.0, exp((x + x)))))) * cbrt(log(exp(fma(-1.0, 1.0, exp((x + x))))))) / (cbrt((exp(x) + 1.0)) * cbrt((exp(x) + 1.0)))) / ((x * cbrt((exp(x) + 1.0))) / cbrt(fma(-1.0, 1.0, exp((x + x))))));
} else {
temp = fma(0.16666666666666666, pow(x, 2.0), fma(0.5, x, 1.0));
}
return temp;
}




Bits error versus x
Results
| Original | 39.6 |
|---|---|
| Target | 40.1 |
| Herbie | 0.3 |
if x < -0.00010392522358541646Initial program 0.1
rmApplied flip--0.1
Simplified0.0
rmApplied add-log-exp0.1
rmApplied add-cube-cbrt0.1
Applied add-cube-cbrt0.1
Applied times-frac0.1
Applied associate-/l*0.1
Simplified0.1
if -0.00010392522358541646 < x Initial program 60.0
Taylor expanded around 0 0.4
Simplified0.4
Final simplification0.3
herbie shell --seed 2020056 +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))