e^{x} - 1{x}^{2} \cdot \left(x \cdot \frac{1}{6} + \frac{1}{2}\right) + xdouble code(double x) {
return (exp(x) - 1.0);
}
double code(double x) {
return ((pow(x, 2.0) * ((x * 0.16666666666666666) + 0.5)) + x);
}




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