\frac{1}{\sqrt{x}} - \frac{1}{\sqrt{x + 1}}\frac{\frac{1}{\sqrt{x + 1}}}{\sqrt{x} \cdot \left(\sqrt{x + 1} + \sqrt{x}\right)}double f(double x) {
double r13871371 = 1.0;
double r13871372 = x;
double r13871373 = sqrt(r13871372);
double r13871374 = r13871371 / r13871373;
double r13871375 = r13871372 + r13871371;
double r13871376 = sqrt(r13871375);
double r13871377 = r13871371 / r13871376;
double r13871378 = r13871374 - r13871377;
return r13871378;
}
double f(double x) {
double r13871379 = 1.0;
double r13871380 = x;
double r13871381 = r13871380 + r13871379;
double r13871382 = sqrt(r13871381);
double r13871383 = r13871379 / r13871382;
double r13871384 = sqrt(r13871380);
double r13871385 = r13871382 + r13871384;
double r13871386 = r13871384 * r13871385;
double r13871387 = r13871383 / r13871386;
return r13871387;
}




Bits error versus x
Results
| Original | 19.7 |
|---|---|
| Target | 0.6 |
| Herbie | 0.4 |
Initial program 19.7
rmApplied frac-sub19.7
Simplified19.7
rmApplied flip--19.5
Simplified0.4
rmApplied div-inv0.4
Applied times-frac0.4
rmApplied *-un-lft-identity0.4
Applied associate-*l*0.4
Simplified0.4
Final simplification0.4
herbie shell --seed 2019128
(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)))))