double f(double x) {
double r1600502 = x;
double r1600503 = 1.0;
double r1600504 = r1600502 + r1600503;
double r1600505 = sqrt(r1600504);
double r1600506 = sqrt(r1600502);
double r1600507 = r1600505 - r1600506;
return r1600507;
}
double f(double x) {
double r1600508 = x;
double r1600509 = 1.0;
double r1600510 = r1600509 - r1600508;
double r1600511 = r1600508 + r1600510;
double r1600512 = r1600508 + r1600509;
double r1600513 = sqrt(r1600512);
double r1600514 = sqrt(r1600508);
double r1600515 = r1600513 + r1600514;
double r1600516 = r1600511 / r1600515;
return r1600516;
}
\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
(FPCore (x)
:name "2sqrt (example 3.1)"
(-.p16 (sqrt.p16 (+.p16 x (real->posit16 1))) (sqrt.p16 x)))