\sqrt{x + 1} - \sqrt{x}\mathsf{log1p}\left(\mathsf{expm1}\left(\frac{1}{\sqrt{x + 1} + \sqrt{x}}\right)\right)double f(double x) {
double r1636046 = x;
double r1636047 = 1.0;
double r1636048 = r1636046 + r1636047;
double r1636049 = sqrt(r1636048);
double r1636050 = sqrt(r1636046);
double r1636051 = r1636049 - r1636050;
return r1636051;
}
double f(double x) {
double r1636052 = 1.0;
double r1636053 = x;
double r1636054 = r1636053 + r1636052;
double r1636055 = sqrt(r1636054);
double r1636056 = sqrt(r1636053);
double r1636057 = r1636055 + r1636056;
double r1636058 = r1636052 / r1636057;
double r1636059 = expm1(r1636058);
double r1636060 = log1p(r1636059);
return r1636060;
}




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