\frac{1}{\sqrt{x}} - \frac{1}{\sqrt{x + 1}}\frac{1}{\sqrt{x + 1} + \sqrt{x}} \cdot \frac{1}{\sqrt{x + 1} \cdot \sqrt{x}}double f(double x) {
double r2972623 = 1.0;
double r2972624 = x;
double r2972625 = sqrt(r2972624);
double r2972626 = r2972623 / r2972625;
double r2972627 = r2972624 + r2972623;
double r2972628 = sqrt(r2972627);
double r2972629 = r2972623 / r2972628;
double r2972630 = r2972626 - r2972629;
return r2972630;
}
double f(double x) {
double r2972631 = 1.0;
double r2972632 = x;
double r2972633 = r2972632 + r2972631;
double r2972634 = sqrt(r2972633);
double r2972635 = sqrt(r2972632);
double r2972636 = r2972634 + r2972635;
double r2972637 = r2972631 / r2972636;
double r2972638 = r2972634 * r2972635;
double r2972639 = r2972631 / r2972638;
double r2972640 = r2972637 * r2972639;
return r2972640;
}




Bits error versus x
Results
| Original | 19.7 |
|---|---|
| Target | 0.6 |
| Herbie | 0.4 |
Initial program 19.7
rmApplied frac-sub19.7
Simplified19.7
rmApplied flip--19.4
Simplified18.9
Taylor expanded around 0 0.4
rmApplied div-inv0.4
Final simplification0.4
herbie shell --seed 2019163 +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)))))