\frac{1}{\sqrt{x}} - \frac{1}{\sqrt{x + 1}}\frac{\frac{1}{\sqrt{\sqrt{x + 1}} \cdot \left(\sqrt{x} \cdot \sqrt{\sqrt{x + 1}}\right)}}{\sqrt{x} + \sqrt{x + 1}}double f(double x) {
double r13070610 = 1.0;
double r13070611 = x;
double r13070612 = sqrt(r13070611);
double r13070613 = r13070610 / r13070612;
double r13070614 = r13070611 + r13070610;
double r13070615 = sqrt(r13070614);
double r13070616 = r13070610 / r13070615;
double r13070617 = r13070613 - r13070616;
return r13070617;
}
double f(double x) {
double r13070618 = 1.0;
double r13070619 = x;
double r13070620 = r13070619 + r13070618;
double r13070621 = sqrt(r13070620);
double r13070622 = sqrt(r13070621);
double r13070623 = sqrt(r13070619);
double r13070624 = r13070623 * r13070622;
double r13070625 = r13070622 * r13070624;
double r13070626 = r13070618 / r13070625;
double r13070627 = r13070623 + r13070621;
double r13070628 = r13070626 / r13070627;
return r13070628;
}




Bits error versus x
Results
| Original | 19.6 |
|---|---|
| Target | 0.6 |
| Herbie | 0.4 |
Initial program 19.6
rmApplied frac-sub19.6
Simplified19.6
rmApplied flip--19.4
Applied associate-/l/19.4
Simplified0.8
rmApplied associate-/r*0.4
rmApplied add-sqr-sqrt0.4
Applied associate-*r*0.4
Final simplification0.4
herbie shell --seed 2019119 +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)))))