\sqrt{x + 1} - \sqrt{x}\sqrt{\frac{1 \cdot 1}{\left(\sqrt{x} + \sqrt{x + 1}\right) \cdot \left(\sqrt{x} + \sqrt{x + 1}\right)}}double f(double x) {
double r96985 = x;
double r96986 = 1.0;
double r96987 = r96985 + r96986;
double r96988 = sqrt(r96987);
double r96989 = sqrt(r96985);
double r96990 = r96988 - r96989;
return r96990;
}
double f(double x) {
double r96991 = 1.0;
double r96992 = r96991 * r96991;
double r96993 = x;
double r96994 = sqrt(r96993);
double r96995 = r96993 + r96991;
double r96996 = sqrt(r96995);
double r96997 = r96994 + r96996;
double r96998 = r96997 * r96997;
double r96999 = r96992 / r96998;
double r97000 = sqrt(r96999);
return r97000;
}




Bits error versus x
Results
| Original | 29.5 |
|---|---|
| Target | 0.2 |
| Herbie | 0.3 |
Initial program 29.5
rmApplied flip--29.4
Simplified0.2
Simplified0.2
rmApplied add-sqr-sqrt0.3
Simplified0.3
Simplified0.3
rmApplied sqrt-unprod0.2
Simplified0.3
Final simplification0.3
herbie shell --seed 2020046 +o rules:numerics
(FPCore (x)
:name "2sqrt (example 3.1)"
:precision binary64
:herbie-target
(/ 1 (+ (sqrt (+ x 1)) (sqrt x)))
(- (sqrt (+ x 1)) (sqrt x)))