\log \left(1 + x\right)
\begin{array}{l}
\mathbf{if}\;1 + x \le 1.00000008116859207:\\
\;\;\;\;\left(1 \cdot x + \log 1\right) - \frac{1}{2} \cdot \frac{{x}^{2}}{{1}^{2}}\\
\mathbf{else}:\\
\;\;\;\;\log \left(1 + x\right)\\
\end{array}double f(double x) {
double r62535 = 1.0;
double r62536 = x;
double r62537 = r62535 + r62536;
double r62538 = log(r62537);
return r62538;
}
double f(double x) {
double r62539 = 1.0;
double r62540 = x;
double r62541 = r62539 + r62540;
double r62542 = 1.000000081168592;
bool r62543 = r62541 <= r62542;
double r62544 = r62539 * r62540;
double r62545 = log(r62539);
double r62546 = r62544 + r62545;
double r62547 = 0.5;
double r62548 = 2.0;
double r62549 = pow(r62540, r62548);
double r62550 = pow(r62539, r62548);
double r62551 = r62549 / r62550;
double r62552 = r62547 * r62551;
double r62553 = r62546 - r62552;
double r62554 = log(r62541);
double r62555 = r62543 ? r62553 : r62554;
return r62555;
}




Bits error versus x
Results
| Original | 38.9 |
|---|---|
| Target | 0.2 |
| Herbie | 0.3 |
if (+ 1.0 x) < 1.000000081168592Initial program 59.0
Taylor expanded around 0 0.4
if 1.000000081168592 < (+ 1.0 x) Initial program 0.3
Final simplification0.3
herbie shell --seed 2020045
(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)))