\sqrt{x + 1} - \sqrt{x}\frac{\sqrt[3]{1} \cdot \sqrt[3]{1}}{\frac{\sqrt{x + 1} + \sqrt{x}}{\sqrt[3]{1}}}double f(double x) {
double r113509 = x;
double r113510 = 1.0;
double r113511 = r113509 + r113510;
double r113512 = sqrt(r113511);
double r113513 = sqrt(r113509);
double r113514 = r113512 - r113513;
return r113514;
}
double f(double x) {
double r113515 = 1.0;
double r113516 = cbrt(r113515);
double r113517 = r113516 * r113516;
double r113518 = x;
double r113519 = r113518 + r113515;
double r113520 = sqrt(r113519);
double r113521 = sqrt(r113518);
double r113522 = r113520 + r113521;
double r113523 = r113522 / r113516;
double r113524 = r113517 / r113523;
return r113524;
}




Bits error versus x
Results
| Original | 29.3 |
|---|---|
| Target | 0.2 |
| Herbie | 0.2 |
Initial program 29.3
rmApplied flip--29.1
Simplified0.2
rmApplied add-sqr-sqrt0.2
Applied sqrt-prod0.2
rmApplied add-cube-cbrt0.2
Applied associate-/l*0.2
Simplified0.2
Final simplification0.2
herbie shell --seed 2019326
(FPCore (x)
:name "2sqrt (example 3.1)"
:precision binary64
:herbie-target
(/ 1 (+ (sqrt (+ x 1)) (sqrt x)))
(- (sqrt (+ x 1)) (sqrt x)))