\frac{1}{\sqrt{x}} - \frac{1}{\sqrt{x + 1}}\frac{1}{\sqrt{\sqrt{x + 1} + \sqrt{x}}} \cdot \frac{\frac{1}{\sqrt{x + 1} \cdot \sqrt{x}}}{\sqrt{\sqrt{x + 1} + \sqrt{x}}}double f(double x) {
double r11252751 = 1.0;
double r11252752 = x;
double r11252753 = sqrt(r11252752);
double r11252754 = r11252751 / r11252753;
double r11252755 = r11252752 + r11252751;
double r11252756 = sqrt(r11252755);
double r11252757 = r11252751 / r11252756;
double r11252758 = r11252754 - r11252757;
return r11252758;
}
double f(double x) {
double r11252759 = 1.0;
double r11252760 = x;
double r11252761 = r11252760 + r11252759;
double r11252762 = sqrt(r11252761);
double r11252763 = sqrt(r11252760);
double r11252764 = r11252762 + r11252763;
double r11252765 = sqrt(r11252764);
double r11252766 = r11252759 / r11252765;
double r11252767 = r11252762 * r11252763;
double r11252768 = r11252759 / r11252767;
double r11252769 = r11252768 / r11252765;
double r11252770 = r11252766 * r11252769;
return r11252770;
}




Bits error versus x
Results
| Original | 20.3 |
|---|---|
| Target | 0.6 |
| Herbie | 0.5 |
Initial program 20.3
rmApplied frac-sub20.3
Simplified20.3
rmApplied flip--20.0
Applied associate-/l/20.0
Simplified0.8
rmApplied associate-/r*0.4
rmApplied add-sqr-sqrt0.5
Applied div-inv0.5
Applied times-frac0.5
Final simplification0.5
herbie shell --seed 2019124 +o rules:numerics
(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)))))