\sqrt{x + 1} - \sqrt{x}\frac{1}{\mathsf{fma}\left(\sqrt{\left|\sqrt[3]{x}\right|} \cdot \sqrt{\sqrt{\sqrt[3]{x}}}, \sqrt{\sqrt{x}}, \sqrt{x + 1}\right)}double f(double x) {
double r187450 = x;
double r187451 = 1.0;
double r187452 = r187450 + r187451;
double r187453 = sqrt(r187452);
double r187454 = sqrt(r187450);
double r187455 = r187453 - r187454;
return r187455;
}
double f(double x) {
double r187456 = 1.0;
double r187457 = x;
double r187458 = cbrt(r187457);
double r187459 = fabs(r187458);
double r187460 = sqrt(r187459);
double r187461 = sqrt(r187458);
double r187462 = sqrt(r187461);
double r187463 = r187460 * r187462;
double r187464 = sqrt(r187457);
double r187465 = sqrt(r187464);
double r187466 = r187457 + r187456;
double r187467 = sqrt(r187466);
double r187468 = fma(r187463, r187465, r187467);
double r187469 = r187456 / r187468;
return r187469;
}




Bits error versus x
| Original | 30.2 |
|---|---|
| Target | 0.2 |
| Herbie | 0.2 |
Initial program 30.2
rmApplied flip--29.9
Simplified0.2
Simplified0.2
rmApplied add-sqr-sqrt0.2
Applied sqrt-prod0.3
Applied fma-def0.2
rmApplied add-cube-cbrt0.2
Applied sqrt-prod0.2
Applied sqrt-prod0.2
Simplified0.2
Final simplification0.2
herbie shell --seed 2020047 +o rules:numerics
(FPCore (x)
:name "2sqrt (example 3.1)"
:precision binary64
:herbie-target
(/ 1 (+ (sqrt (+ x 1)) (sqrt x)))
(- (sqrt (+ x 1)) (sqrt x)))