\frac{1}{\sqrt{x}} - \frac{1}{\sqrt{x + 1}}\frac{1 \cdot \frac{1}{\sqrt{\sqrt{x + 1} + \sqrt{x}} \cdot \sqrt{\sqrt{x + 1} + \sqrt{x}}}}{\sqrt{x} \cdot \sqrt{x + 1}}double f(double x) {
double r160128 = 1.0;
double r160129 = x;
double r160130 = sqrt(r160129);
double r160131 = r160128 / r160130;
double r160132 = r160129 + r160128;
double r160133 = sqrt(r160132);
double r160134 = r160128 / r160133;
double r160135 = r160131 - r160134;
return r160135;
}
double f(double x) {
double r160136 = 1.0;
double r160137 = x;
double r160138 = r160137 + r160136;
double r160139 = sqrt(r160138);
double r160140 = sqrt(r160137);
double r160141 = r160139 + r160140;
double r160142 = sqrt(r160141);
double r160143 = r160142 * r160142;
double r160144 = r160136 / r160143;
double r160145 = r160136 * r160144;
double r160146 = r160140 * r160139;
double r160147 = r160145 / r160146;
return r160147;
}




Bits error versus x
Results
| Original | 20.0 |
|---|---|
| Target | 0.7 |
| Herbie | 0.5 |
Initial program 20.0
rmApplied frac-sub19.9
Simplified19.9
rmApplied flip--19.7
Simplified19.3
Taylor expanded around 0 0.4
rmApplied add-sqr-sqrt0.5
Final simplification0.5
herbie shell --seed 2020057
(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)))))