\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 r169260 = 1.0;
double r169261 = x;
double r169262 = sqrt(r169261);
double r169263 = r169260 / r169262;
double r169264 = r169261 + r169260;
double r169265 = sqrt(r169264);
double r169266 = r169260 / r169265;
double r169267 = r169263 - r169266;
return r169267;
}
double f(double x) {
double r169268 = 1.0;
double r169269 = x;
double r169270 = sqrt(r169269);
double r169271 = r169268 / r169270;
double r169272 = r169269 + r169268;
double r169273 = sqrt(r169272);
double r169274 = r169273 * r169270;
double r169275 = r169274 + r169272;
double r169276 = r169268 / r169275;
double r169277 = r169271 * r169276;
return r169277;
}




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