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




Bits error versus x
Results
| Original | 39.3 |
|---|---|
| Target | 0.2 |
| Herbie | 0.3 |
if (+ 1.0 x) < 1.0000000000568305Initial program 59.3
rmApplied flip3-+59.4
Applied log-div59.4
Taylor expanded around 0 0.2
Simplified0.2
if 1.0000000000568305 < (+ 1.0 x) Initial program 0.5
Final simplification0.3
herbie shell --seed 2020057
(FPCore (x)
:name "ln(1 + x)"
:precision binary64
:herbie-target
(if (== (+ 1 x) 1) x (/ (* x (log (+ 1 x))) (- (+ 1 x) 1)))
(log (+ 1 x)))