\sqrt{x + 1} - \sqrt{x}\frac{\sqrt{\frac{1 + 0}{\sqrt{x + 1} + \sqrt{x}}} \cdot \sqrt{1}}{\sqrt{\sqrt{x + 1} + \sqrt{x}}}double f(double x) {
double r135443 = x;
double r135444 = 1.0;
double r135445 = r135443 + r135444;
double r135446 = sqrt(r135445);
double r135447 = sqrt(r135443);
double r135448 = r135446 - r135447;
return r135448;
}
double f(double x) {
double r135449 = 1.0;
double r135450 = 0.0;
double r135451 = r135449 + r135450;
double r135452 = x;
double r135453 = r135452 + r135449;
double r135454 = sqrt(r135453);
double r135455 = sqrt(r135452);
double r135456 = r135454 + r135455;
double r135457 = r135451 / r135456;
double r135458 = sqrt(r135457);
double r135459 = sqrt(r135449);
double r135460 = r135458 * r135459;
double r135461 = sqrt(r135456);
double r135462 = r135460 / r135461;
return r135462;
}




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