\frac{1}{\sqrt{x}} - \frac{1}{\sqrt{x + 1}}\frac{1 \cdot 1}{\sqrt{x + 1} \cdot x + \left(x + 1\right) \cdot \sqrt{x}}double f(double x) {
double r126180 = 1.0;
double r126181 = x;
double r126182 = sqrt(r126181);
double r126183 = r126180 / r126182;
double r126184 = r126181 + r126180;
double r126185 = sqrt(r126184);
double r126186 = r126180 / r126185;
double r126187 = r126183 - r126186;
return r126187;
}
double f(double x) {
double r126188 = 1.0;
double r126189 = r126188 * r126188;
double r126190 = x;
double r126191 = r126190 + r126188;
double r126192 = sqrt(r126191);
double r126193 = r126192 * r126190;
double r126194 = sqrt(r126190);
double r126195 = r126191 * r126194;
double r126196 = r126193 + r126195;
double r126197 = r126189 / r126196;
return r126197;
}




Bits error versus x
Results
| Original | 19.6 |
|---|---|
| Target | 0.7 |
| Herbie | 0.7 |
Initial program 19.6
rmApplied frac-sub19.6
Simplified19.6
rmApplied flip--19.4
Applied associate-*r/19.4
Applied associate-/l/19.4
Simplified19.4
Taylor expanded around 0 0.7
rmApplied *-un-lft-identity0.7
Final simplification0.7
herbie shell --seed 2020046 +o rules:numerics
(FPCore (x)
:name "2isqrt (example 3.6)"
:precision binary64
:herbie-target
(/ 1 (+ (* (+ x 1) (sqrt x)) (* x (sqrt (+ x 1)))))
(- (/ 1 (sqrt x)) (/ 1 (sqrt (+ x 1)))))