\log \left(1 + x\right)
\begin{array}{l}
\mathbf{if}\;1 + x \le 1.00002052412234321:\\
\;\;\;\;\left(1 \cdot x + 0.333333333333333315 \cdot {x}^{3}\right) - 0.5 \cdot {x}^{2}\\
\mathbf{else}:\\
\;\;\;\;\log \left(1 + x\right)\\
\end{array}double code(double x) {
return ((double) log(((double) (1.0 + x))));
}
double code(double x) {
double VAR;
if ((((double) (1.0 + x)) <= 1.0000205241223432)) {
VAR = ((double) (((double) (((double) (1.0 * x)) + ((double) (0.3333333333333333 * ((double) pow(x, 3.0)))))) - ((double) (0.5 * ((double) pow(x, 2.0))))));
} else {
VAR = ((double) log(((double) (1.0 + x))));
}
return VAR;
}




Bits error versus x
Results
| Original | 39.0 |
|---|---|
| Target | 0.3 |
| Herbie | 0.2 |
if (+ 1.0 x) < 1.00002052412234321Initial program 59.1
rmApplied flip-+59.1
Applied log-div59.1
Taylor expanded around 0 0.2
Simplified0.2
Taylor expanded around 0 0.2
Simplified0.2
if 1.00002052412234321 < (+ 1.0 x) Initial program 0.1
Final simplification0.2
herbie shell --seed 2020150
(FPCore (x)
:name "ln(1 + x)"
:precision binary64
:herbie-target
(if (== (+ 1.0 x) 1.0) x (/ (* x (log (+ 1.0 x))) (- (+ 1.0 x) 1.0)))
(log (+ 1.0 x)))