2 \cdot \tan^{-1} \left(\sqrt{\frac{1 - x}{1 + x}}\right)2 \cdot \tan^{-1} \left(\sqrt{\mathsf{log1p}\left(\mathsf{expm1}\left(\frac{1 - x}{1 + x}\right)\right)}\right)double code(double x) {
return (2.0 * atan(sqrt(((1.0 - x) / (1.0 + x)))));
}
double code(double x) {
return (2.0 * atan(sqrt(log1p(expm1(((1.0 - x) / (1.0 + x)))))));
}



Bits error versus x
Results
Initial program 0.0
rmApplied log1p-expm1-u0.0
Final simplification0.0
herbie shell --seed 2020078 +o rules:numerics
(FPCore (x)
:name "arccos"
:precision binary64
(* 2 (atan (sqrt (/ (- 1 x) (+ 1 x))))))