\log \left(1 + x\right)
\begin{array}{l}
\mathbf{if}\;1 + x \le 1.0000014320063078:\\
\;\;\;\;\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 r62172 = 1.0;
double r62173 = x;
double r62174 = r62172 + r62173;
double r62175 = log(r62174);
return r62175;
}
double f(double x) {
double r62176 = 1.0;
double r62177 = x;
double r62178 = r62176 + r62177;
double r62179 = 1.0000014320063078;
bool r62180 = r62178 <= r62179;
double r62181 = r62176 * r62177;
double r62182 = log(r62176);
double r62183 = r62181 + r62182;
double r62184 = 0.5;
double r62185 = 2.0;
double r62186 = pow(r62177, r62185);
double r62187 = pow(r62176, r62185);
double r62188 = r62186 / r62187;
double r62189 = r62184 * r62188;
double r62190 = r62183 - r62189;
double r62191 = log(r62178);
double r62192 = r62180 ? r62190 : r62191;
return r62192;
}




Bits error versus x
Results
| Original | 38.6 |
|---|---|
| Target | 0.3 |
| Herbie | 0.2 |
if (+ 1.0 x) < 1.0000014320063078Initial program 59.2
Taylor expanded around 0 0.3
if 1.0000014320063078 < (+ 1.0 x) Initial program 0.1
Final simplification0.2
herbie shell --seed 2020018
(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)))