\log \left(1 + x\right)
\begin{array}{l}
\mathbf{if}\;1 + x \le 1.000000152443675460744998417794704437256:\\
\;\;\;\;\left(x \cdot 1 + x \cdot \frac{\frac{-1}{2}}{\frac{1 \cdot 1}{x}}\right) + \log 1\\
\mathbf{else}:\\
\;\;\;\;\log \left(\sqrt{1 + x}\right) + \log \left(\sqrt{1 + x}\right)\\
\end{array}double f(double x) {
double r47680 = 1.0;
double r47681 = x;
double r47682 = r47680 + r47681;
double r47683 = log(r47682);
return r47683;
}
double f(double x) {
double r47684 = 1.0;
double r47685 = x;
double r47686 = r47684 + r47685;
double r47687 = 1.0000001524436755;
bool r47688 = r47686 <= r47687;
double r47689 = r47685 * r47684;
double r47690 = -0.5;
double r47691 = r47684 * r47684;
double r47692 = r47691 / r47685;
double r47693 = r47690 / r47692;
double r47694 = r47685 * r47693;
double r47695 = r47689 + r47694;
double r47696 = log(r47684);
double r47697 = r47695 + r47696;
double r47698 = sqrt(r47686);
double r47699 = log(r47698);
double r47700 = r47699 + r47699;
double r47701 = r47688 ? r47697 : r47700;
return r47701;
}




Bits error versus x
Results
| Original | 39.6 |
|---|---|
| Target | 0.2 |
| Herbie | 0.3 |
if (+ 1.0 x) < 1.0000001524436755Initial program 59.1
Taylor expanded around 0 0.4
Simplified0.4
rmApplied distribute-lft-in0.4
if 1.0000001524436755 < (+ 1.0 x) Initial program 0.2
rmApplied add-sqr-sqrt0.2
Applied log-prod0.2
Final simplification0.3
herbie shell --seed 2019303
(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)))