\cos^{-1} \left(\frac{1 - 5 \cdot \left(v \cdot v\right)}{v \cdot v - 1}\right)\mathsf{expm1}\left(\mathsf{log1p}\left(\cos^{-1} \left(\frac{1}{v + \sqrt{1}} \cdot \frac{1 - 5 \cdot \left(v \cdot v\right)}{v - \sqrt{1}}\right)\right)\right)double code(double v) {
return acos(((1.0 - (5.0 * (v * v))) / ((v * v) - 1.0)));
}
double code(double v) {
return expm1(log1p(acos(((1.0 / (v + sqrt(1.0))) * ((1.0 - (5.0 * (v * v))) / (v - sqrt(1.0)))))));
}



Bits error versus v
Results
Initial program 0.5
rmApplied add-sqr-sqrt0.5
Applied difference-of-squares0.8
Applied *-un-lft-identity0.8
Applied times-frac0.9
rmApplied expm1-log1p-u0.9
Final simplification0.9
herbie shell --seed 2020060 +o rules:numerics
(FPCore (v)
:name "Falkner and Boettcher, Appendix B, 1"
:precision binary64
(acos (/ (- 1 (* 5 (* v v))) (- (* v v) 1))))