\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 r2636761 = x;
double r2636762 = 1.0;
double r2636763 = /* ERROR: no posit support in C */;
double r2636764 = r2636761 + r2636763;
double r2636765 = sqrt(r2636764);
double r2636766 = sqrt(r2636761);
double r2636767 = r2636765 - r2636766;
return r2636767;
}
double f(double x) {
double r2636768 = x;
double r2636769 = 1.0;
double r2636770 = r2636769 - r2636768;
double r2636771 = r2636768 + r2636770;
double r2636772 = r2636768 + r2636769;
double r2636773 = sqrt(r2636772);
double r2636774 = sqrt(r2636768);
double r2636775 = r2636773 + r2636774;
double r2636776 = r2636771 / r2636775;
return r2636776;
}



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