\frac{1}{\sqrt{x}} - \frac{1}{\sqrt{x + 1}}\frac{\frac{\sqrt{1 + \left(x - x\right)}}{\sqrt{\sqrt{x} + \sqrt{x + 1}}}}{\frac{\sqrt{x} \cdot \sqrt{x + 1}}{\frac{\sqrt{1 + \left(x - x\right)}}{\sqrt{\sqrt{x} + \sqrt{x + 1}}}}}double f(double x) {
double r5037947 = 1.0;
double r5037948 = x;
double r5037949 = sqrt(r5037948);
double r5037950 = r5037947 / r5037949;
double r5037951 = r5037948 + r5037947;
double r5037952 = sqrt(r5037951);
double r5037953 = r5037947 / r5037952;
double r5037954 = r5037950 - r5037953;
return r5037954;
}
double f(double x) {
double r5037955 = 1.0;
double r5037956 = x;
double r5037957 = r5037956 - r5037956;
double r5037958 = r5037955 + r5037957;
double r5037959 = sqrt(r5037958);
double r5037960 = sqrt(r5037956);
double r5037961 = r5037956 + r5037955;
double r5037962 = sqrt(r5037961);
double r5037963 = r5037960 + r5037962;
double r5037964 = sqrt(r5037963);
double r5037965 = r5037959 / r5037964;
double r5037966 = r5037960 * r5037962;
double r5037967 = r5037966 / r5037965;
double r5037968 = r5037965 / r5037967;
return r5037968;
}




Bits error versus x
Results
| Original | 19.3 |
|---|---|
| Target | 0.6 |
| Herbie | 0.5 |
Initial program 19.3
rmApplied frac-sub19.3
Simplified19.3
rmApplied flip--19.1
Simplified0.4
rmApplied add-sqr-sqrt0.4
Applied add-sqr-sqrt0.4
Applied times-frac0.5
Applied associate-/l*0.5
Final simplification0.5
herbie shell --seed 2019130 +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)))))