\sqrt{x + 1} - \sqrt{x}\frac{1}{\sqrt{\sqrt{x + 1}} \cdot \sqrt{\sqrt{x + 1}} + \sqrt{x}}double f(double x) {
double r1762325 = x;
double r1762326 = 1.0;
double r1762327 = r1762325 + r1762326;
double r1762328 = sqrt(r1762327);
double r1762329 = sqrt(r1762325);
double r1762330 = r1762328 - r1762329;
return r1762330;
}
double f(double x) {
double r1762331 = 1.0;
double r1762332 = x;
double r1762333 = r1762332 + r1762331;
double r1762334 = sqrt(r1762333);
double r1762335 = sqrt(r1762334);
double r1762336 = r1762335 * r1762335;
double r1762337 = sqrt(r1762332);
double r1762338 = r1762336 + r1762337;
double r1762339 = r1762331 / r1762338;
return r1762339;
}




Bits error versus x
Results
| Original | 29.7 |
|---|---|
| Target | 0.2 |
| Herbie | 0.3 |
Initial program 29.7
rmApplied flip--29.5
Simplified29.0
Taylor expanded around 0 0.2
rmApplied add-sqr-sqrt0.2
Applied sqrt-prod0.3
Final simplification0.3
herbie shell --seed 2019154 +o rules:numerics
(FPCore (x)
:name "2sqrt (example 3.1)"
:herbie-target
(/ 1 (+ (sqrt (+ x 1)) (sqrt x)))
(- (sqrt (+ x 1)) (sqrt x)))