-\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)1 \cdot \frac{2 \cdot \left(-\log \left({\left(\frac{\cosh \left(\pi \cdot \frac{f}{4}\right)}{\sinh \left(\pi \cdot \frac{f}{4}\right)}\right)}^{0.3333333333333333}\right)\right)}{\frac{\pi}{4}} - 1 \cdot \frac{\log \left(\sqrt[3]{\frac{\cosh \left(\pi \cdot \frac{f}{4}\right)}{\sinh \left(\pi \cdot \frac{f}{4}\right)}}\right)}{\frac{\pi}{4}}double code(double f) {
return ((double) -(((double) ((1.0 / (((double) M_PI) / 4.0)) * ((double) log((((double) (((double) exp(((double) ((((double) M_PI) / 4.0) * f)))) + ((double) exp(((double) -(((double) ((((double) M_PI) / 4.0) * f)))))))) / ((double) (((double) exp(((double) ((((double) M_PI) / 4.0) * f)))) - ((double) exp(((double) -(((double) ((((double) M_PI) / 4.0) * f)))))))))))))));
}
double code(double f) {
return ((double) (((double) (1.0 * (((double) (2.0 * ((double) -(((double) log(((double) pow((((double) cosh(((double) (((double) M_PI) * (f / 4.0))))) / ((double) sinh(((double) (((double) M_PI) * (f / 4.0)))))), 0.3333333333333333)))))))) / (((double) M_PI) / 4.0)))) - ((double) (1.0 * (((double) log(((double) cbrt((((double) cosh(((double) (((double) M_PI) * (f / 4.0))))) / ((double) sinh(((double) (((double) M_PI) * (f / 4.0)))))))))) / (((double) M_PI) / 4.0))))));
}



Bits error versus f
Results
Initial program 61.5
rmApplied sinh-undef2.2
Applied associate-/r*2.2
Simplified2.2
rmApplied add-cube-cbrt2.2
Applied log-prod2.2
Applied distribute-lft-in2.2
Simplified2.2
Simplified2.2
rmApplied pow1/32.0
Final simplification2.0
herbie shell --seed 2020196
(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)))))))))