\sqrt{x + 1} - \sqrt{x}\frac{1}{\sqrt{\sqrt{x + 1}} \cdot \sqrt{\sqrt{x + 1}} + \sqrt{x}}double f(double x) {
double r503622 = x;
double r503623 = 1.0;
double r503624 = r503622 + r503623;
double r503625 = sqrt(r503624);
double r503626 = sqrt(r503622);
double r503627 = r503625 - r503626;
return r503627;
}
double f(double x) {
double r503628 = 1.0;
double r503629 = x;
double r503630 = r503629 + r503628;
double r503631 = sqrt(r503630);
double r503632 = sqrt(r503631);
double r503633 = r503632 * r503632;
double r503634 = sqrt(r503629);
double r503635 = r503633 + r503634;
double r503636 = r503628 / r503635;
return r503636;
}




Bits error versus x
Results
| Original | 29.6 |
|---|---|
| Target | 0.2 |
| Herbie | 0.3 |
Initial program 29.6
rmApplied flip--29.4
Simplified0.2
rmApplied add-sqr-sqrt0.2
Applied sqrt-prod0.3
Final simplification0.3
herbie shell --seed 2020059
(FPCore (x)
:name "Main:bigenough3 from C"
:precision binary64
:herbie-target
(/ 1 (+ (sqrt (+ x 1)) (sqrt x)))
(- (sqrt (+ x 1)) (sqrt x)))