\frac{1}{\sqrt{x}} - \frac{1}{\sqrt{x + 1}}\frac{1}{\sqrt{x}} \cdot \frac{1}{\sqrt{x + 1} \cdot \sqrt{x} + \left(x + 1\right)}double f(double x) {
double r109475 = 1.0;
double r109476 = x;
double r109477 = sqrt(r109476);
double r109478 = r109475 / r109477;
double r109479 = r109476 + r109475;
double r109480 = sqrt(r109479);
double r109481 = r109475 / r109480;
double r109482 = r109478 - r109481;
return r109482;
}
double f(double x) {
double r109483 = 1.0;
double r109484 = x;
double r109485 = sqrt(r109484);
double r109486 = r109483 / r109485;
double r109487 = r109484 + r109483;
double r109488 = sqrt(r109487);
double r109489 = r109488 * r109485;
double r109490 = r109489 + r109487;
double r109491 = r109483 / r109490;
double r109492 = r109486 * r109491;
return r109492;
}




Bits error versus x
Results
| Original | 20.2 |
|---|---|
| Target | 0.7 |
| Herbie | 0.3 |
Initial program 20.2
rmApplied frac-sub20.2
Simplified20.2
rmApplied flip--20.0
Simplified19.5
Taylor expanded around 0 0.4
rmApplied times-frac0.4
Simplified0.3
Final simplification0.3
herbie shell --seed 2019291
(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)))))