\frac{e^{x}}{e^{x} - 1}\frac{0.5}{\sinh x} + \frac{0.5}{\sinh x} \cdot e^{x}(FPCore (x) :precision binary64 (/ (exp x) (- (exp x) 1.0)))
(FPCore (x) :precision binary64 (+ (/ 0.5 (sinh x)) (* (/ 0.5 (sinh x)) (exp x))))
double code(double x) {
return exp(x) / (exp(x) - 1.0);
}
double code(double x) {
return (0.5 / sinh(x)) + ((0.5 / sinh(x)) * exp(x));
}




Bits error versus x
Results
| Original | 41.0 |
|---|---|
| Target | 40.6 |
| Herbie | 0.5 |
Initial program 41.0
rmApplied flip--_binary6441.1
Applied associate-/r/_binary6441.1
Simplified41.0
rmApplied sinh-undef_binary640.5
Applied associate-/r*_binary640.5
Simplified0.5
rmApplied distribute-rgt-in_binary640.5
Simplified0.5
Simplified0.5
Final simplification0.5
herbie shell --seed 2020288
(FPCore (x)
:name "expq2 (section 3.11)"
:precision binary64
:herbie-target
(/ 1.0 (- 1.0 (exp (- x))))
(/ (exp x) (- (exp x) 1.0)))