\log \left(1 + x\right)
\begin{array}{l}
\mathbf{if}\;x + 1 \le 1.000052777771681:\\
\;\;\;\;x + \left(\frac{1}{3} \cdot x - \frac{1}{2}\right) \cdot \left(x \cdot x\right)\\
\mathbf{else}:\\
\;\;\;\;\log \left(\sqrt{x + 1}\right) + \left(\log \left(\sqrt{\sqrt{x + 1}}\right) + \log \left(x + 1\right) \cdot \frac{1}{4}\right)\\
\end{array}double f(double x) {
double r3745833 = 1.0;
double r3745834 = x;
double r3745835 = r3745833 + r3745834;
double r3745836 = log(r3745835);
return r3745836;
}
double f(double x) {
double r3745837 = x;
double r3745838 = 1.0;
double r3745839 = r3745837 + r3745838;
double r3745840 = 1.000052777771681;
bool r3745841 = r3745839 <= r3745840;
double r3745842 = 0.3333333333333333;
double r3745843 = r3745842 * r3745837;
double r3745844 = 0.5;
double r3745845 = r3745843 - r3745844;
double r3745846 = r3745837 * r3745837;
double r3745847 = r3745845 * r3745846;
double r3745848 = r3745837 + r3745847;
double r3745849 = sqrt(r3745839);
double r3745850 = log(r3745849);
double r3745851 = sqrt(r3745849);
double r3745852 = log(r3745851);
double r3745853 = log(r3745839);
double r3745854 = 0.25;
double r3745855 = r3745853 * r3745854;
double r3745856 = r3745852 + r3745855;
double r3745857 = r3745850 + r3745856;
double r3745858 = r3745841 ? r3745848 : r3745857;
return r3745858;
}




Bits error versus x
Results
| Original | 39.2 |
|---|---|
| Target | 0.2 |
| Herbie | 0.2 |
if (+ 1 x) < 1.000052777771681Initial program 59.1
Taylor expanded around 0 0.2
Simplified0.2
if 1.000052777771681 < (+ 1 x) Initial program 0.1
rmApplied add-sqr-sqrt0.1
Applied log-prod0.1
rmApplied add-sqr-sqrt0.1
Applied log-prod0.1
rmApplied pow1/20.1
Applied sqrt-pow10.1
Applied log-pow0.1
Final simplification0.2
herbie shell --seed 2019149
(FPCore (x)
:name "ln(1 + x)"
:herbie-target
(if (== (+ 1 x) 1) x (/ (* x (log (+ 1 x))) (- (+ 1 x) 1)))
(log (+ 1 x)))