\log \left(1 + x\right)
\begin{array}{l}
\mathbf{if}\;1 + x \le 1:\\
\;\;\;\;\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 r88199 = 1.0;
double r88200 = x;
double r88201 = r88199 + r88200;
double r88202 = log(r88201);
return r88202;
}
double f(double x) {
double r88203 = 1.0;
double r88204 = x;
double r88205 = r88203 + r88204;
bool r88206 = r88205 <= r88203;
double r88207 = r88203 * r88204;
double r88208 = log(r88203);
double r88209 = r88207 + r88208;
double r88210 = 0.5;
double r88211 = 2.0;
double r88212 = pow(r88204, r88211);
double r88213 = pow(r88203, r88211);
double r88214 = r88212 / r88213;
double r88215 = r88210 * r88214;
double r88216 = r88209 - r88215;
double r88217 = log(r88205);
double r88218 = r88206 ? r88216 : r88217;
return r88218;
}




Bits error versus x
Results
| Original | 39.3 |
|---|---|
| Target | 0.3 |
| Herbie | 0.6 |
if (+ 1.0 x) < 1.0Initial program 59.6
Taylor expanded around 0 0.3
if 1.0 < (+ 1.0 x) Initial program 1.2
Final simplification0.6
herbie shell --seed 2020046
(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)))