Average Error: 39.5 → 0.2
Time: 21.9s
Precision: 64
Internal Precision: 1408
\[\log \left(1 + x\right)\]
↓
\[\begin{array}{l}
\mathbf{if}\;\log \left(1 + x\right) \le 2.686244946081019 \cdot 10^{-06}:\\
\;\;\;\;x + \left(x \cdot x\right) \cdot \left(x \cdot \frac{1}{3} - \frac{1}{2}\right)\\
\mathbf{else}:\\
\;\;\;\;\log \left(1 + x\right)\\
\end{array}\]
Target
| Original | 39.5 |
|---|
| Target | 0.2 |
|---|
| Herbie | 0.2 |
|---|
\[\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)) < 2.686244946081019e-06
Initial program 59.1
\[\log \left(1 + x\right)\]
Taylor expanded around 0 0.2
\[\leadsto \color{blue}{\left(\frac{1}{3} \cdot {x}^{3} + x\right) - \frac{1}{2} \cdot {x}^{2}}\]
Applied simplify0.2
\[\leadsto \color{blue}{x + \left(x \cdot x\right) \cdot \left(x \cdot \frac{1}{3} - \frac{1}{2}\right)}\]
if 2.686244946081019e-06 < (log (+ 1 x))
Initial program 0.1
\[\log \left(1 + x\right)\]
- Recombined 2 regimes into one program.
Runtime
herbie shell --seed '#(1070258749 1877548225 2229079127 1588002776 3179087814 1886870650)'
(FPCore (x)
:name "ln(1 + x)"
:herbie-target
(if (== (+ 1 x) 1) x (/ (* x (log (+ 1 x))) (- (+ 1 x) 1)))
(log (+ 1 x)))