e^{x} - 1{x}^{2} \cdot \left(x \cdot \frac{1}{6} + \frac{1}{2}\right) + xdouble f(double x) {
double r166219 = x;
double r166220 = exp(r166219);
double r166221 = 1.0;
double r166222 = r166220 - r166221;
return r166222;
}
double f(double x) {
double r166223 = x;
double r166224 = 2.0;
double r166225 = pow(r166223, r166224);
double r166226 = 0.16666666666666666;
double r166227 = r166223 * r166226;
double r166228 = 0.5;
double r166229 = r166227 + r166228;
double r166230 = r166225 * r166229;
double r166231 = r166230 + r166223;
return r166231;
}




Bits error versus x
Results
| Original | 58.8 |
|---|---|
| Target | 0.4 |
| Herbie | 0.4 |
Initial program 58.8
Taylor expanded around 0 0.4
Simplified0.4
Final simplification0.4
herbie shell --seed 2019318
(FPCore (x)
:name "expm1 (example 3.7)"
:precision binary64
:pre (< -1.7e-4 x)
:herbie-target
(* x (+ (+ 1 (/ x 2)) (/ (* x x) 6)))
(- (exp x) 1))