\frac{1}{\sqrt{x}} - \frac{1}{\sqrt{x + 1}}\frac{1}{\sqrt{x}} \cdot \frac{\frac{1}{\sqrt{x + 1} \cdot 1 + 1 \cdot \sqrt{x}}}{\sqrt{x + 1}}double f(double x) {
double r8723467 = 1.0;
double r8723468 = x;
double r8723469 = sqrt(r8723468);
double r8723470 = r8723467 / r8723469;
double r8723471 = r8723468 + r8723467;
double r8723472 = sqrt(r8723471);
double r8723473 = r8723467 / r8723472;
double r8723474 = r8723470 - r8723473;
return r8723474;
}
double f(double x) {
double r8723475 = 1.0;
double r8723476 = x;
double r8723477 = sqrt(r8723476);
double r8723478 = r8723475 / r8723477;
double r8723479 = 1.0;
double r8723480 = r8723476 + r8723479;
double r8723481 = sqrt(r8723480);
double r8723482 = r8723481 * r8723479;
double r8723483 = r8723479 * r8723477;
double r8723484 = r8723482 + r8723483;
double r8723485 = r8723479 / r8723484;
double r8723486 = r8723485 / r8723481;
double r8723487 = r8723478 * r8723486;
return r8723487;
}




Bits error versus x
Results
| Original | 20.0 |
|---|---|
| Target | 0.7 |
| Herbie | 0.4 |
Initial program 20.0
rmApplied frac-sub20.0
rmApplied flip--19.8
Taylor expanded around 0 0.4
rmApplied *-un-lft-identity0.4
Applied *-un-lft-identity0.4
Applied times-frac0.4
Applied times-frac0.4
Simplified0.4
Final simplification0.4
herbie shell --seed 2019171
(FPCore (x)
:name "2isqrt (example 3.6)"
:herbie-target
(/ 1.0 (+ (* (+ x 1.0) (sqrt x)) (* x (sqrt (+ x 1.0)))))
(- (/ 1.0 (sqrt x)) (/ 1.0 (sqrt (+ x 1.0)))))