\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 r709090 = x;
double r709091 = 1.0;
double r709092 = /* ERROR: no posit support in C */;
double r709093 = r709090 + r709092;
double r709094 = sqrt(r709093);
double r709095 = sqrt(r709090);
double r709096 = r709094 - r709095;
return r709096;
}
double f(double x) {
double r709097 = x;
double r709098 = 1.0;
double r709099 = r709098 - r709097;
double r709100 = r709097 + r709099;
double r709101 = r709097 + r709098;
double r709102 = sqrt(r709101);
double r709103 = sqrt(r709097);
double r709104 = r709102 + r709103;
double r709105 = r709100 / r709104;
return r709105;
}



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