\frac{1}{\sqrt{x}} - \frac{1}{\sqrt{x + 1}}\frac{1 \cdot \left(\sqrt{\frac{1}{\sqrt{x + 1} + \sqrt{x}}} \cdot \sqrt{\frac{1}{\sqrt{x + 1} + \sqrt{x}}}\right)}{\sqrt{x} \cdot \sqrt{x + 1}}double f(double x) {
double r160788 = 1.0;
double r160789 = x;
double r160790 = sqrt(r160789);
double r160791 = r160788 / r160790;
double r160792 = r160789 + r160788;
double r160793 = sqrt(r160792);
double r160794 = r160788 / r160793;
double r160795 = r160791 - r160794;
return r160795;
}
double f(double x) {
double r160796 = 1.0;
double r160797 = x;
double r160798 = r160797 + r160796;
double r160799 = sqrt(r160798);
double r160800 = sqrt(r160797);
double r160801 = r160799 + r160800;
double r160802 = r160796 / r160801;
double r160803 = sqrt(r160802);
double r160804 = r160803 * r160803;
double r160805 = r160796 * r160804;
double r160806 = r160800 * r160799;
double r160807 = r160805 / r160806;
return r160807;
}




Bits error versus x
Results
| Original | 19.3 |
|---|---|
| Target | 0.7 |
| Herbie | 0.4 |
Initial program 19.3
rmApplied frac-sub19.2
Simplified19.2
rmApplied flip--19.0
Simplified18.6
Taylor expanded around 0 0.4
rmApplied add-sqr-sqrt0.4
Final simplification0.4
herbie shell --seed 2020027 +o rules:numerics
(FPCore (x)
:name "2isqrt (example 3.6)"
:precision binary64
:herbie-target
(/ 1 (+ (* (+ x 1) (sqrt x)) (* x (sqrt (+ x 1)))))
(- (/ 1 (sqrt x)) (/ 1 (sqrt (+ x 1)))))