\frac{1}{\sqrt{x}} - \frac{1}{\sqrt{x + 1}}\frac{\frac{1}{x}}{\frac{1}{\sqrt{1 + x}} + \frac{1}{\sqrt{x}}} \cdot \frac{1 \cdot 1}{1 + x}double f(double x) {
double r6251504 = 1.0;
double r6251505 = x;
double r6251506 = sqrt(r6251505);
double r6251507 = r6251504 / r6251506;
double r6251508 = r6251505 + r6251504;
double r6251509 = sqrt(r6251508);
double r6251510 = r6251504 / r6251509;
double r6251511 = r6251507 - r6251510;
return r6251511;
}
double f(double x) {
double r6251512 = 1.0;
double r6251513 = x;
double r6251514 = r6251512 / r6251513;
double r6251515 = r6251512 + r6251513;
double r6251516 = sqrt(r6251515);
double r6251517 = r6251512 / r6251516;
double r6251518 = sqrt(r6251513);
double r6251519 = r6251512 / r6251518;
double r6251520 = r6251517 + r6251519;
double r6251521 = r6251514 / r6251520;
double r6251522 = r6251512 * r6251512;
double r6251523 = r6251522 / r6251515;
double r6251524 = r6251521 * r6251523;
return r6251524;
}




Bits error versus x
Results
| Original | 20.0 |
|---|---|
| Target | 0.7 |
| Herbie | 0.4 |
Initial program 20.0
rmApplied flip--20.0
rmApplied frac-times25.2
Applied frac-times20.1
Applied frac-sub19.9
Simplified19.5
Simplified19.5
Taylor expanded around 0 5.9
rmApplied *-un-lft-identity5.9
Applied times-frac5.5
Applied times-frac0.4
Simplified0.4
Final simplification0.4
herbie shell --seed 2019171 +o rules:numerics
(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)))))