\sqrt{x + 1} - \sqrt{x}\mathsf{expm1}\left(\mathsf{log1p}\left(\frac{1}{\sqrt{x + 1} + \sqrt{x}}\right)\right)double f(double x) {
double r2062061 = x;
double r2062062 = 1.0;
double r2062063 = r2062061 + r2062062;
double r2062064 = sqrt(r2062063);
double r2062065 = sqrt(r2062061);
double r2062066 = r2062064 - r2062065;
return r2062066;
}
double f(double x) {
double r2062067 = 1.0;
double r2062068 = x;
double r2062069 = r2062068 + r2062067;
double r2062070 = sqrt(r2062069);
double r2062071 = sqrt(r2062068);
double r2062072 = r2062070 + r2062071;
double r2062073 = r2062067 / r2062072;
double r2062074 = log1p(r2062073);
double r2062075 = expm1(r2062074);
return r2062075;
}




Bits error versus x
Results
| Original | 30.1 |
|---|---|
| Target | 0.2 |
| Herbie | 0.2 |
Initial program 30.1
rmApplied flip--29.8
Simplified29.4
Simplified29.4
rmApplied expm1-log1p-u29.4
Simplified0.2
Final simplification0.2
herbie shell --seed 2019152 +o rules:numerics
(FPCore (x)
:name "2sqrt (example 3.1)"
:herbie-target
(/ 1 (+ (sqrt (+ x 1)) (sqrt x)))
(- (sqrt (+ x 1)) (sqrt x)))