e^{a \cdot x} - 1(e^{a \cdot x} - 1)^*double f(double a, double x) {
double r5057027 = a;
double r5057028 = x;
double r5057029 = r5057027 * r5057028;
double r5057030 = exp(r5057029);
double r5057031 = 1.0;
double r5057032 = r5057030 - r5057031;
return r5057032;
}
double f(double a, double x) {
double r5057033 = a;
double r5057034 = x;
double r5057035 = r5057033 * r5057034;
double r5057036 = expm1(r5057035);
return r5057036;
}




Bits error versus a




Bits error versus x
Results
| Original | 28.8 |
|---|---|
| Target | 0.2 |
| Herbie | 0.0 |
Initial program 28.8
Simplified0.0
Final simplification0.0
herbie shell --seed 2019119 +o rules:numerics
(FPCore (a x)
:name "expax (section 3.5)"
:herbie-target
(if (< (fabs (* a x)) 1/10) (* (* a x) (+ 1 (+ (/ (* a x) 2) (/ (pow (* a x) 2) 6)))) (- (exp (* a x)) 1))
(- (exp (* a x)) 1))