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




Bits error versus x
Results
| Original | 58.5 |
|---|---|
| Target | 0.5 |
| Herbie | 0.5 |
Initial program 58.5
Taylor expanded around 0 0.5
Simplified0.5
rmApplied distribute-lft-in0.5
Applied associate-+l+0.5
Simplified0.5
Final simplification0.5
herbie shell --seed 2020130
(FPCore (x)
:name "expm1 (example 3.7)"
:precision binary64
:pre (< -0.00017 x)
:herbie-target
(* x (+ (+ 1.0 (/ x 2.0)) (/ (* x x) 6.0)))
(- (exp x) 1.0))