Average Error: 40.0 → 0.2
Time: 48.7s
Precision: 64
Internal Precision: 1344
\[\frac{e^{x} - 1}{x}\]
↓
\[\begin{array}{l}
\mathbf{if}\;\frac{1}{6} \cdot {x}^{2} + \left(1 + \frac{1}{2} \cdot x\right) \le 1.000036663958746:\\
\;\;\;\;\frac{1}{6} \cdot {x}^{2} + \left(1 + \frac{1}{2} \cdot x\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{e^{x}}{x} - \frac{1}{x}\\
\end{array}\]
Target
| Original | 40.0 |
|---|
| Target | 39.1 |
|---|
| Herbie | 0.2 |
|---|
\[\begin{array}{l}
\mathbf{if}\;x \lt 1 \land x \gt -1:\\
\;\;\;\;\frac{e^{x} - 1}{\log \left(e^{x}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{e^{x} - 1}{x}\\
\end{array}\]
Derivation
- Split input into 2 regimes
if (+ (* 1/6 (pow x 2)) (+ 1 (* 1/2 x))) < 1.000036663958746
Initial program 60.4
\[\frac{e^{x} - 1}{x}\]
Taylor expanded around 0 0.2
\[\leadsto \color{blue}{\frac{1}{6} \cdot {x}^{2} + \left(1 + \frac{1}{2} \cdot x\right)}\]
if 1.000036663958746 < (+ (* 1/6 (pow x 2)) (+ 1 (* 1/2 x)))
Initial program 0.1
\[\frac{e^{x} - 1}{x}\]
- Using strategy
rm Applied div-sub0.1
\[\leadsto \color{blue}{\frac{e^{x}}{x} - \frac{1}{x}}\]
- Recombined 2 regimes into one program.
Runtime
herbie shell --seed '#(1072330854 3074818769 591214268 3603999196 3863745332 3332387116)'
(FPCore (x)
:name "Kahan's exp quotient"
:herbie-target
(if (and (< x 1) (> x -1)) (/ (- (exp x) 1) (log (exp x))) (/ (- (exp x) 1) x))
(/ (- (exp x) 1) x))