Average Error: 39.0 → 0.3
Time: 1.9m
Precision: 64
Internal Precision: 1408
\[\log \left(1 + x\right)\]
↓
\[\begin{array}{l}
\mathbf{if}\;\log \left(1 + x\right) \le 5.0712818213415073 \cdot 10^{-08}:\\
\;\;\;\;x - \left(\frac{1}{2} - x \cdot \frac{1}{3}\right) \cdot \left(x \cdot x\right)\\
\mathbf{else}:\\
\;\;\;\;\log \left(1 + x\right)\\
\end{array}\]
Target
| Original | 39.0 |
|---|
| Target | 0.2 |
|---|
| Herbie | 0.3 |
|---|
\[\begin{array}{l}
\mathbf{if}\;1 + x = 1:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;\frac{x \cdot \log \left(1 + x\right)}{\left(1 + x\right) - 1}\\
\end{array}\]
Derivation
- Split input into 2 regimes
if (log (+ 1 x)) < 5.0712818213415073e-08
Initial program 59.1
\[\log \left(1 + x\right)\]
Taylor expanded around 0 0.3
\[\leadsto \color{blue}{\left(\frac{1}{3} \cdot {x}^{3} + x\right) - \frac{1}{2} \cdot {x}^{2}}\]
Applied simplify0.3
\[\leadsto \color{blue}{x - \left(\frac{1}{2} - x \cdot \frac{1}{3}\right) \cdot \left(x \cdot x\right)}\]
if 5.0712818213415073e-08 < (log (+ 1 x))
Initial program 0.3
\[\log \left(1 + x\right)\]
- Recombined 2 regimes into one program.
Runtime
herbie shell --seed '#(1070578969 3140398606 632207097 462683394 1189254563 964980650)'
(FPCore (x)
:name "ln(1 + x)"
:herbie-target
(if (== (+ 1 x) 1) x (/ (* x (log (+ 1 x))) (- (+ 1 x) 1)))
(log (+ 1 x)))