\frac{1}{\sqrt{x}} - \frac{1}{\sqrt{x + 1}}\frac{\sqrt{\frac{1}{\sqrt{x} \cdot 1 + \sqrt{x + 1} \cdot 1}}}{\sqrt{x}} \cdot \frac{\sqrt{\frac{1}{\sqrt{x} \cdot 1 + \sqrt{x + 1} \cdot 1}}}{\sqrt{x + 1}}double f(double x) {
double r6332459 = 1.0;
double r6332460 = x;
double r6332461 = sqrt(r6332460);
double r6332462 = r6332459 / r6332461;
double r6332463 = r6332460 + r6332459;
double r6332464 = sqrt(r6332463);
double r6332465 = r6332459 / r6332464;
double r6332466 = r6332462 - r6332465;
return r6332466;
}
double f(double x) {
double r6332467 = 1.0;
double r6332468 = x;
double r6332469 = sqrt(r6332468);
double r6332470 = r6332469 * r6332467;
double r6332471 = r6332468 + r6332467;
double r6332472 = sqrt(r6332471);
double r6332473 = r6332472 * r6332467;
double r6332474 = r6332470 + r6332473;
double r6332475 = r6332467 / r6332474;
double r6332476 = sqrt(r6332475);
double r6332477 = r6332476 / r6332469;
double r6332478 = r6332476 / r6332472;
double r6332479 = r6332477 * r6332478;
return r6332479;
}




Bits error versus x
Results
| Original | 20.0 |
|---|---|
| Target | 0.7 |
| Herbie | 0.4 |
Initial program 20.0
rmApplied frac-sub20.0
rmApplied flip--19.8
Taylor expanded around 0 0.4
rmApplied add-sqr-sqrt0.4
Applied times-frac0.4
Final simplification0.4
herbie shell --seed 2019171
(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)))))