\frac{1}{\sqrt{x}} - \frac{1}{\sqrt{x + 1}}\frac{\frac{1}{\sqrt{\sqrt{x + 1}} \cdot \left(\sqrt{x} \cdot \sqrt{\sqrt{x + 1}}\right)}}{\sqrt{x} + \sqrt{x + 1}}double f(double x) {
double r6589201 = 1.0;
double r6589202 = x;
double r6589203 = sqrt(r6589202);
double r6589204 = r6589201 / r6589203;
double r6589205 = r6589202 + r6589201;
double r6589206 = sqrt(r6589205);
double r6589207 = r6589201 / r6589206;
double r6589208 = r6589204 - r6589207;
return r6589208;
}
double f(double x) {
double r6589209 = 1.0;
double r6589210 = x;
double r6589211 = r6589210 + r6589209;
double r6589212 = sqrt(r6589211);
double r6589213 = sqrt(r6589212);
double r6589214 = sqrt(r6589210);
double r6589215 = r6589214 * r6589213;
double r6589216 = r6589213 * r6589215;
double r6589217 = r6589209 / r6589216;
double r6589218 = r6589214 + r6589212;
double r6589219 = r6589217 / r6589218;
return r6589219;
}




Bits error versus x
Results
| Original | 19.6 |
|---|---|
| Target | 0.6 |
| Herbie | 0.4 |
Initial program 19.6
rmApplied frac-sub19.6
Simplified19.6
rmApplied flip--19.4
Applied associate-/l/19.4
Simplified0.8
rmApplied associate-/r*0.4
rmApplied add-sqr-sqrt0.4
Applied associate-*r*0.4
Final simplification0.4
herbie shell --seed 2019119
(FPCore (x)
:name "2isqrt (example 3.6)"
:herbie-target
(/ 1 (+ (* (+ x 1) (sqrt x)) (* x (sqrt (+ x 1)))))
(- (/ 1 (sqrt x)) (/ 1 (sqrt (+ x 1)))))