\sqrt{x + 1} - \sqrt{x}\frac{\sqrt{\frac{1}{\sqrt{x + 1} + \sqrt{x}}} \cdot \sqrt{1}}{\sqrt{\sqrt{x + 1} + \sqrt{x}}}double f(double x) {
double r103847 = x;
double r103848 = 1.0;
double r103849 = r103847 + r103848;
double r103850 = sqrt(r103849);
double r103851 = sqrt(r103847);
double r103852 = r103850 - r103851;
return r103852;
}
double f(double x) {
double r103853 = 1.0;
double r103854 = x;
double r103855 = r103854 + r103853;
double r103856 = sqrt(r103855);
double r103857 = sqrt(r103854);
double r103858 = r103856 + r103857;
double r103859 = r103853 / r103858;
double r103860 = sqrt(r103859);
double r103861 = sqrt(r103853);
double r103862 = r103860 * r103861;
double r103863 = sqrt(r103858);
double r103864 = r103862 / r103863;
return r103864;
}




Bits error versus x
Results
| Original | 29.8 |
|---|---|
| Target | 0.2 |
| Herbie | 0.2 |
Initial program 29.8
rmApplied flip--29.5
Simplified0.2
rmApplied add-sqr-sqrt0.3
rmApplied sqrt-div0.3
Applied associate-*r/0.2
Final simplification0.2
herbie shell --seed 2019208
(FPCore (x)
:name "2sqrt (example 3.1)"
:precision binary64
:herbie-target
(/ 1 (+ (sqrt (+ x 1)) (sqrt x)))
(- (sqrt (+ x 1)) (sqrt x)))