\sqrt{x + 1} - \sqrt{x}\frac{\sqrt{1 + 0}}{\frac{\sqrt{x + 1} + \sqrt{x}}{\sqrt{1}}}double f(double x) {
double r110012 = x;
double r110013 = 1.0;
double r110014 = r110012 + r110013;
double r110015 = sqrt(r110014);
double r110016 = sqrt(r110012);
double r110017 = r110015 - r110016;
return r110017;
}
double f(double x) {
double r110018 = 1.0;
double r110019 = 0.0;
double r110020 = r110018 + r110019;
double r110021 = sqrt(r110020);
double r110022 = x;
double r110023 = r110022 + r110018;
double r110024 = sqrt(r110023);
double r110025 = sqrt(r110022);
double r110026 = r110024 + r110025;
double r110027 = sqrt(r110018);
double r110028 = r110026 / r110027;
double r110029 = r110021 / r110028;
return r110029;
}




Bits error versus x
Results
| Original | 29.9 |
|---|---|
| Target | 0.2 |
| Herbie | 0.2 |
Initial program 29.9
rmApplied flip--29.7
Simplified0.2
rmApplied add-sqr-sqrt0.2
Applied associate-/l*0.2
Simplified0.2
Final simplification0.2
herbie shell --seed 2019354 +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)))