\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 r2748098 = x;
double r2748099 = 1.0;
double r2748100 = r2748098 + r2748099;
double r2748101 = sqrt(r2748100);
double r2748102 = sqrt(r2748098);
double r2748103 = r2748101 - r2748102;
return r2748103;
}
double f(double x) {
double r2748104 = 1.0;
double r2748105 = x;
double r2748106 = r2748105 + r2748104;
double r2748107 = sqrt(r2748106);
double r2748108 = sqrt(r2748105);
double r2748109 = r2748107 + r2748108;
double r2748110 = r2748104 / r2748109;
double r2748111 = expm1(r2748110);
double r2748112 = log1p(r2748111);
return r2748112;
}




Bits error versus x
Results
| Original | 29.9 |
|---|---|
| Target | 0.2 |
| Herbie | 0.2 |
Initial program 29.9
rmApplied flip--29.7
Simplified0.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)))