-\frac{1}{\frac{\pi}{4}} \cdot \log \left(\frac{e^{\frac{\pi}{4} \cdot f} + e^{-\frac{\pi}{4} \cdot f}}{e^{\frac{\pi}{4} \cdot f} - e^{-\frac{\pi}{4} \cdot f}}\right)-\frac{\log \left(\left(\sqrt{\frac{e^{\frac{\pi}{4} \cdot f} + e^{-\frac{\pi}{4} \cdot f}}{\sinh \left(\frac{\pi}{4} \cdot f\right)}} \cdot \sqrt{e^{\log \left(\frac{e^{\frac{\pi}{4} \cdot f} + e^{-\frac{\pi}{4} \cdot f}}{\sinh \left(\frac{\pi}{4} \cdot f\right)}\right)}}\right) \cdot 0.5\right)}{\frac{\pi}{4}}(FPCore (f)
:precision binary64
(-
(*
(/ 1.0 (/ PI 4.0))
(log
(/
(+ (exp (* (/ PI 4.0) f)) (exp (- (* (/ PI 4.0) f))))
(- (exp (* (/ PI 4.0) f)) (exp (- (* (/ PI 4.0) f)))))))))(FPCore (f)
:precision binary64
(-
(/
(log
(*
(*
(sqrt
(/
(+ (exp (* (/ PI 4.0) f)) (exp (- (* (/ PI 4.0) f))))
(sinh (* (/ PI 4.0) f))))
(sqrt
(exp
(log
(/
(+ (exp (* (/ PI 4.0) f)) (exp (- (* (/ PI 4.0) f))))
(sinh (* (/ PI 4.0) f)))))))
0.5))
(/ PI 4.0))))double code(double f) {
return -((1.0 / (((double) M_PI) / 4.0)) * log((exp((((double) M_PI) / 4.0) * f) + exp(-((((double) M_PI) / 4.0) * f))) / (exp((((double) M_PI) / 4.0) * f) - exp(-((((double) M_PI) / 4.0) * f)))));
}
double code(double f) {
return -(log((sqrt((exp((((double) M_PI) / 4.0) * f) + exp(-((((double) M_PI) / 4.0) * f))) / sinh((((double) M_PI) / 4.0) * f)) * sqrt(exp(log((exp((((double) M_PI) / 4.0) * f) + exp(-((((double) M_PI) / 4.0) * f))) / sinh((((double) M_PI) / 4.0) * f))))) * 0.5) / (((double) M_PI) / 4.0));
}





Bits error versus f
Results
| Alternative 1 | |
|---|---|
| Accuracy | 2.0 |
| Cost | 3392 |
| Alternative 2 | |
|---|---|
| Accuracy | 2.0 |
| Cost | 4992 |
Initial program 61.7
rmApplied sinh-undef_binary642.0
Applied *-un-lft-identity_binary642.0
Applied times-frac_binary642.1
Simplified2.1
rmApplied associate-*l/_binary642.0
Simplified2.0
rmApplied add-sqr-sqrt_binary642.0
rmApplied add-exp-log_binary642.0
Applied add-exp-log_binary642.0
Applied div-exp_binary642.0
Simplified2.0
Final simplification2.0
herbie shell --seed 2020322
(FPCore (f)
:name "VandenBroeck and Keller, Equation (20)"
:precision binary64
(- (* (/ 1.0 (/ PI 4.0)) (log (/ (+ (exp (* (/ PI 4.0) f)) (exp (- (* (/ PI 4.0) f)))) (- (exp (* (/ PI 4.0) f)) (exp (- (* (/ PI 4.0) f)))))))))