\sqrt{x + 1} - \sqrt{x}\frac{\sqrt{\frac{1}{\sqrt{x} + \sqrt{x + 1}}} \cdot \sqrt{1}}{\sqrt{\sqrt{x} + \sqrt{x + 1}}}double f(double x) {
double r122678 = x;
double r122679 = 1.0;
double r122680 = r122678 + r122679;
double r122681 = sqrt(r122680);
double r122682 = sqrt(r122678);
double r122683 = r122681 - r122682;
return r122683;
}
double f(double x) {
double r122684 = 1.0;
double r122685 = x;
double r122686 = sqrt(r122685);
double r122687 = r122685 + r122684;
double r122688 = sqrt(r122687);
double r122689 = r122686 + r122688;
double r122690 = r122684 / r122689;
double r122691 = sqrt(r122690);
double r122692 = sqrt(r122684);
double r122693 = r122691 * r122692;
double r122694 = sqrt(r122689);
double r122695 = r122693 / r122694;
return r122695;
}




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