\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 r310076 = x;
double r310077 = 1.0;
double r310078 = r310076 + r310077;
double r310079 = sqrt(r310078);
double r310080 = sqrt(r310076);
double r310081 = r310079 - r310080;
return r310081;
}
double f(double x) {
double r310082 = 1.0;
double r310083 = x;
double r310084 = sqrt(r310083);
double r310085 = r310083 + r310082;
double r310086 = sqrt(r310085);
double r310087 = r310084 + r310086;
double r310088 = r310082 / r310087;
double r310089 = sqrt(r310088);
double r310090 = sqrt(r310082);
double r310091 = r310089 * r310090;
double r310092 = sqrt(r310087);
double r310093 = r310091 / r310092;
return r310093;
}




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