\frac{1}{\sqrt{x}} - \frac{1}{\sqrt{x + 1}}\frac{\frac{1 \cdot 1}{x + \sqrt{x} \cdot \sqrt{x + 1}}}{\sqrt{x + 1}}double f(double x) {
double r160952 = 1.0;
double r160953 = x;
double r160954 = sqrt(r160953);
double r160955 = r160952 / r160954;
double r160956 = r160953 + r160952;
double r160957 = sqrt(r160956);
double r160958 = r160952 / r160957;
double r160959 = r160955 - r160958;
return r160959;
}
double f(double x) {
double r160960 = 1.0;
double r160961 = r160960 * r160960;
double r160962 = x;
double r160963 = sqrt(r160962);
double r160964 = r160962 + r160960;
double r160965 = sqrt(r160964);
double r160966 = r160963 * r160965;
double r160967 = r160962 + r160966;
double r160968 = r160961 / r160967;
double r160969 = r160968 / r160965;
return r160969;
}




Bits error versus x
Results
| Original | 19.7 |
|---|---|
| Target | 0.7 |
| Herbie | 0.3 |
Initial program 19.7
rmApplied frac-sub19.7
Simplified19.7
rmApplied flip--19.4
Simplified19.0
Taylor expanded around 0 0.4
rmApplied associate-/r*0.4
Simplified0.3
Final simplification0.3
herbie shell --seed 2020043
(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)))))