\frac{1}{\sqrt{x}} - \frac{1}{\sqrt{x + 1}}\frac{\frac{1 \cdot \frac{1}{\sqrt{x + 1} + \sqrt{x}}}{\sqrt{x}}}{\sqrt{x + 1}}double f(double x) {
double r171819 = 1.0;
double r171820 = x;
double r171821 = sqrt(r171820);
double r171822 = r171819 / r171821;
double r171823 = r171820 + r171819;
double r171824 = sqrt(r171823);
double r171825 = r171819 / r171824;
double r171826 = r171822 - r171825;
return r171826;
}
double f(double x) {
double r171827 = 1.0;
double r171828 = x;
double r171829 = r171828 + r171827;
double r171830 = sqrt(r171829);
double r171831 = sqrt(r171828);
double r171832 = r171830 + r171831;
double r171833 = r171827 / r171832;
double r171834 = r171827 * r171833;
double r171835 = r171834 / r171831;
double r171836 = r171835 / r171830;
return r171836;
}




Bits error versus x
Results
| Original | 19.9 |
|---|---|
| Target | 0.7 |
| Herbie | 0.4 |
Initial program 19.9
rmApplied frac-sub19.9
Simplified19.9
rmApplied flip--19.6
Simplified19.2
Taylor expanded around 0 0.4
rmApplied associate-/r*0.4
Final simplification0.4
herbie shell --seed 2020024 +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)))))