\frac{1}{\sqrt{x}} - \frac{1}{\sqrt{x + 1}}\frac{\frac{1}{x + \sqrt{x + 1} \cdot \sqrt{x}}}{\sqrt{x + 1}}double f(double x) {
double r4915252 = 1.0;
double r4915253 = x;
double r4915254 = sqrt(r4915253);
double r4915255 = r4915252 / r4915254;
double r4915256 = r4915253 + r4915252;
double r4915257 = sqrt(r4915256);
double r4915258 = r4915252 / r4915257;
double r4915259 = r4915255 - r4915258;
return r4915259;
}
double f(double x) {
double r4915260 = 1.0;
double r4915261 = x;
double r4915262 = r4915261 + r4915260;
double r4915263 = sqrt(r4915262);
double r4915264 = sqrt(r4915261);
double r4915265 = r4915263 * r4915264;
double r4915266 = r4915261 + r4915265;
double r4915267 = r4915260 / r4915266;
double r4915268 = r4915267 / r4915263;
return r4915268;
}




Bits error versus x
Results
| Original | 19.7 |
|---|---|
| Target | 0.6 |
| Herbie | 0.3 |
Initial program 19.7
rmApplied frac-sub19.7
rmApplied flip--19.4
Simplified18.9
Simplified18.9
rmApplied associate-/r*18.9
Simplified0.3
Final simplification0.3
herbie shell --seed 2019163
(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)))))