\left(\sqrt{\left(\frac{x}{\left(1\right)}\right)}\right) - \left(\sqrt{x}\right)\frac{x + \left(1 - x\right)}{\sqrt{x + 1} + \sqrt{x}}double f(double x) {
double r3166925 = x;
double r3166926 = 1.0;
double r3166927 = /* ERROR: no posit support in C */;
double r3166928 = r3166925 + r3166927;
double r3166929 = sqrt(r3166928);
double r3166930 = sqrt(r3166925);
double r3166931 = r3166929 - r3166930;
return r3166931;
}
double f(double x) {
double r3166932 = x;
double r3166933 = 1.0;
double r3166934 = r3166933 - r3166932;
double r3166935 = r3166932 + r3166934;
double r3166936 = r3166932 + r3166933;
double r3166937 = sqrt(r3166936);
double r3166938 = sqrt(r3166932);
double r3166939 = r3166937 + r3166938;
double r3166940 = r3166935 / r3166939;
return r3166940;
}



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 2019107
(FPCore (x)
:name "2sqrt (example 3.1)"
(-.p16 (sqrt.p16 (+.p16 x (real->posit16 1))) (sqrt.p16 x)))