\frac{2}{e^{x} + e^{-x}}
\mathsf{log1p}\left(\mathsf{expm1}\left(\frac{1}{\cosh x}\right)\right)
(FPCore (x) :precision binary64 (/ 2.0 (+ (exp x) (exp (- x)))))
(FPCore (x) :precision binary64 (log1p (expm1 (/ 1.0 (cosh x)))))
double code(double x) {
return 2.0 / (exp(x) + exp(-x));
}
double code(double x) {
return log1p(expm1(1.0 / cosh(x)));
}



Bits error versus x
Results
Initial program 0.0
Applied cosh-undef_binary640.0
Applied add-sqr-sqrt_binary640.5
Applied times-frac_binary640.5
Applied log1p-expm1-u_binary640.5
Simplified0.0
Final simplification0.0
herbie shell --seed 2021332
(FPCore (x)
:name "Hyperbolic secant"
:precision binary64
(/ 2.0 (+ (exp x) (exp (- x)))))