\frac{1}{\sqrt{x}} - \frac{1}{\sqrt{x + 1}}\frac{\frac{1}{x + \sqrt{x + 1} \cdot \sqrt{x}}}{\sqrt{x + 1}}double f(double x) {
double r5393795 = 1.0;
double r5393796 = x;
double r5393797 = sqrt(r5393796);
double r5393798 = r5393795 / r5393797;
double r5393799 = r5393796 + r5393795;
double r5393800 = sqrt(r5393799);
double r5393801 = r5393795 / r5393800;
double r5393802 = r5393798 - r5393801;
return r5393802;
}
double f(double x) {
double r5393803 = 1.0;
double r5393804 = x;
double r5393805 = r5393804 + r5393803;
double r5393806 = sqrt(r5393805);
double r5393807 = sqrt(r5393804);
double r5393808 = r5393806 * r5393807;
double r5393809 = r5393804 + r5393808;
double r5393810 = r5393803 / r5393809;
double r5393811 = r5393810 / r5393806;
return r5393811;
}




Bits error versus x
Results
| Original | 19.7 |
|---|---|
| Target | 0.6 |
| Herbie | 0.3 |
Initial program 19.7
rmApplied frac-sub19.7
rmApplied flip--19.4
Simplified18.9
Simplified18.9
rmApplied associate-/r*18.9
Simplified0.3
Final simplification0.3
herbie shell --seed 2019163
(FPCore (x)
:name "2isqrt (example 3.6)"
:herbie-target
(/ 1 (+ (* (+ x 1) (sqrt x)) (* x (sqrt (+ x 1)))))
(- (/ 1 (sqrt x)) (/ 1 (sqrt (+ x 1)))))