\log \left(1 + x\right)
\begin{array}{l}
\mathbf{if}\;1 + x \le 1.0000031303406751:\\
\;\;\;\;\log 1 + \left(1 \cdot x + x \cdot \left(x \cdot \frac{\frac{-1}{2}}{1 \cdot 1}\right)\right)\\
\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.0000031303406751)) {
VAR = ((double) (((double) log(1.0)) + ((double) (((double) (1.0 * x)) + ((double) (x * ((double) (x * ((double) (-0.5 / ((double) (1.0 * 1.0))))))))))));
} else {
VAR = ((double) log(((double) (1.0 + x))));
}
return VAR;
}




Bits error versus x
Results
| Original | 39.0 |
|---|---|
| Target | 0.2 |
| Herbie | 0.2 |
if (+ 1.0 x) < 1.0000031303406751Initial program 59.1
Taylor expanded around 0 0.3
Simplified0.3
rmApplied distribute-lft-in0.3
Simplified0.3
if 1.0000031303406751 < (+ 1.0 x) Initial program 0.1
Final simplification0.2
herbie shell --seed 2020184
(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)))