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



Bits error versus x
Results
Initial program 0.0
Applied expm1-log1p-u_binary640.0
Final simplification0.0
herbie shell --seed 2021275
(FPCore (x)
:name "Hyperbolic secant"
:precision binary64
(/ 2.0 (+ (exp x) (exp (- x)))))