\frac{e^{x}}{e^{x} - 1}\frac{\frac{1}{\mathsf{expm1}\left(x\right)}}{\frac{1}{e^{x}}}double f(double x) {
double r6472580 = x;
double r6472581 = exp(r6472580);
double r6472582 = 1.0;
double r6472583 = r6472581 - r6472582;
double r6472584 = r6472581 / r6472583;
return r6472584;
}
double f(double x) {
double r6472585 = 1.0;
double r6472586 = x;
double r6472587 = expm1(r6472586);
double r6472588 = r6472585 / r6472587;
double r6472589 = exp(r6472586);
double r6472590 = r6472585 / r6472589;
double r6472591 = r6472588 / r6472590;
return r6472591;
}




Bits error versus x
Results
| Original | 39.8 |
|---|---|
| Target | 39.3 |
| Herbie | 0.4 |
Initial program 39.8
rmApplied expm1-def0.4
rmApplied *-un-lft-identity0.4
Applied associate-/l*0.4
rmApplied div-inv0.4
Applied associate-/r*0.4
Final simplification0.4
herbie shell --seed 2019124 +o rules:numerics
(FPCore (x)
:name "expq2 (section 3.11)"
:herbie-target
(/ 1 (- 1 (exp (- x))))
(/ (exp x) (- (exp x) 1)))