double code(double x) {
return ((double) log(((double) (((double) (1.0 / x)) + ((double) (((double) sqrt(((double) (1.0 - ((double) (x * x)))))) / x))))));
}
double code(double x) {
return ((double) (((double) log(((double) sqrt(((double) (((double) (1.0 + ((double) sqrt(((double) (1.0 - ((double) (x * x)))))))) / x)))))) + ((double) log(((double) sqrt(((double) (((double) (1.0 + ((double) sqrt(((double) (1.0 - ((double) (x * x)))))))) / x))))))));
}



Bits error versus x
Results
Initial program 0.1
Simplified0.1
rmApplied add-sqr-sqrt0.1
Applied log-prod0.1
Final simplification0.1
herbie shell --seed 2020191
(FPCore (x)
:name "Hyperbolic arc-(co)secant"
:precision binary64
(log (+ (/ 1.0 x) (/ (sqrt (- 1.0 (* x x))) x))))