\frac{x}{1 + \sqrt{x + 1}}\begin{array}{l}
\mathbf{if}\;x \le 28329089.9676262401044368743896484375:\\
\;\;\;\;\frac{x}{{1}^{3} + {\left(\sqrt{x + 1}\right)}^{3}} \cdot \left(1 \cdot 1 + \left(\sqrt{x + 1} \cdot \sqrt{x + 1} - 1 \cdot \sqrt{x + 1}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\sqrt{x} \cdot \frac{\sqrt{x}}{1 + \sqrt{x + 1}}\\
\end{array}double f(double x) {
double r149486 = x;
double r149487 = 1.0;
double r149488 = r149486 + r149487;
double r149489 = sqrt(r149488);
double r149490 = r149487 + r149489;
double r149491 = r149486 / r149490;
return r149491;
}
double f(double x) {
double r149492 = x;
double r149493 = 28329089.96762624;
bool r149494 = r149492 <= r149493;
double r149495 = 1.0;
double r149496 = 3.0;
double r149497 = pow(r149495, r149496);
double r149498 = r149492 + r149495;
double r149499 = sqrt(r149498);
double r149500 = pow(r149499, r149496);
double r149501 = r149497 + r149500;
double r149502 = r149492 / r149501;
double r149503 = r149495 * r149495;
double r149504 = r149499 * r149499;
double r149505 = r149495 * r149499;
double r149506 = r149504 - r149505;
double r149507 = r149503 + r149506;
double r149508 = r149502 * r149507;
double r149509 = sqrt(r149492);
double r149510 = r149495 + r149499;
double r149511 = r149509 / r149510;
double r149512 = r149509 * r149511;
double r149513 = r149494 ? r149508 : r149512;
return r149513;
}



Bits error versus x
Results
if x < 28329089.96762624Initial program 0.0
rmApplied flip3-+0.0
Applied associate-/r/0.0
if 28329089.96762624 < x Initial program 0.5
rmApplied *-un-lft-identity0.5
Applied add-sqr-sqrt0.1
Applied times-frac0.0
Simplified0.0
Final simplification0.0
herbie shell --seed 2019350
(FPCore (x)
:name "Numeric.Log:$clog1p from log-domain-0.10.2.1, B"
:precision binary64
(/ x (+ 1 (sqrt (+ x 1)))))