\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 r438712 = x;
double r438713 = 1.0;
double r438714 = r438712 + r438713;
double r438715 = sqrt(r438714);
double r438716 = sqrt(r438712);
double r438717 = r438715 - r438716;
return r438717;
}
double f(double x) {
double r438718 = 1.0;
double r438719 = 0.0;
double r438720 = r438718 + r438719;
double r438721 = x;
double r438722 = r438721 + r438718;
double r438723 = sqrt(r438722);
double r438724 = sqrt(r438721);
double r438725 = r438723 + r438724;
double r438726 = r438720 / r438725;
double r438727 = sqrt(r438726);
double r438728 = sqrt(r438718);
double r438729 = r438727 * r438728;
double r438730 = sqrt(r438725);
double r438731 = r438729 / r438730;
return r438731;
}




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 "Main:bigenough3 from C"
:precision binary64
:herbie-target
(/ 1 (+ (sqrt (+ x 1)) (sqrt x)))
(- (sqrt (+ x 1)) (sqrt x)))