double f(double x) {
double r1812252 = x;
double r1812253 = 1.0;
double r1812254 = r1812252 + r1812253;
double r1812255 = sqrt(r1812254);
double r1812256 = sqrt(r1812252);
double r1812257 = r1812255 - r1812256;
return r1812257;
}
double f(double x) {
double r1812258 = x;
double r1812259 = 1.0;
double r1812260 = r1812259 - r1812258;
double r1812261 = r1812258 + r1812260;
double r1812262 = r1812258 + r1812259;
double r1812263 = sqrt(r1812262);
double r1812264 = sqrt(r1812258);
double r1812265 = r1812263 + r1812264;
double r1812266 = r1812261 / r1812265;
return r1812266;
}
\sqrt{x + 1} - \sqrt{x}\frac{x + \left(1 - x\right)}{\sqrt{x + 1} + \sqrt{x}}


Bits error versus x
Initial program 0.8
rmApplied p16-flip--0.6
rmApplied sqrt-sqrd.p160.5
rmApplied sqrt-sqrd.p160.4
rmApplied associate--l+0.3
Final simplification0.3
herbie shell --seed 2019101 +o rules:numerics
(FPCore (x)
:name "2sqrt (example 3.1)"
(-.p16 (sqrt.p16 (+.p16 x (real->posit16 1))) (sqrt.p16 x)))