\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 r1658141 = x;
double r1658142 = 1.0;
double r1658143 = /* ERROR: no posit support in C */;
double r1658144 = r1658141 + r1658143;
double r1658145 = sqrt(r1658144);
double r1658146 = sqrt(r1658141);
double r1658147 = r1658145 - r1658146;
return r1658147;
}
double f(double x) {
double r1658148 = x;
double r1658149 = 1.0;
double r1658150 = r1658149 - r1658148;
double r1658151 = r1658148 + r1658150;
double r1658152 = r1658148 + r1658149;
double r1658153 = sqrt(r1658152);
double r1658154 = sqrt(r1658148);
double r1658155 = r1658153 + r1658154;
double r1658156 = r1658151 / r1658155;
return r1658156;
}



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