\frac{e^{x}}{e^{x} - 1}\frac{1}{\frac{\mathsf{expm1}\left(x\right)}{e^{x}}}double f(double x) {
double r6113690 = x;
double r6113691 = exp(r6113690);
double r6113692 = 1.0;
double r6113693 = r6113691 - r6113692;
double r6113694 = r6113691 / r6113693;
return r6113694;
}
double f(double x) {
double r6113695 = 1.0;
double r6113696 = x;
double r6113697 = expm1(r6113696);
double r6113698 = exp(r6113696);
double r6113699 = r6113697 / r6113698;
double r6113700 = r6113695 / r6113699;
return r6113700;
}




Bits error versus x
Results
| Original | 40.1 |
|---|---|
| Target | 39.7 |
| Herbie | 0.3 |
Initial program 40.1
rmApplied expm1-log1p-u40.1
Simplified0.3
rmApplied *-un-lft-identity0.3
Applied associate-/l*0.3
Final simplification0.3
herbie shell --seed 2019125 +o rules:numerics
(FPCore (x)
:name "expq2 (section 3.11)"
:herbie-target
(/ 1 (- 1 (exp (- x))))
(/ (exp x) (- (exp x) 1)))