\log \left(\frac{1}{x} + \frac{\sqrt{1 - x \cdot x}}{x}\right)
\log \left(\frac{1 + \sqrt{1 - x \cdot x}}{x}\right)
(FPCore (x) :precision binary64 (log (+ (/ 1.0 x) (/ (sqrt (- 1.0 (* x x))) x))))
(FPCore (x) :precision binary64 (log (/ (+ 1.0 (sqrt (- 1.0 (* x x)))) x)))
double code(double x) {
return log((1.0 / x) + (sqrt(1.0 - (x * x)) / x));
}
double code(double x) {
return log((1.0 + sqrt(1.0 - (x * x))) / x);
}



Bits error versus x
Results
Initial program 0.0
rmApplied add-exp-log_binary640.0
Simplified0.0
rmApplied rem-exp-log_binary640.0
Final simplification0.0
herbie shell --seed 2021202
(FPCore (x)
:name "Hyperbolic arc-(co)secant"
:precision binary64
(log (+ (/ 1.0 x) (/ (sqrt (- 1.0 (* x x))) x))))