\cos^{-1} \left(\frac{1 - 5 \cdot \left(v \cdot v\right)}{v \cdot v - 1}\right)\frac{\pi}{2} - \sin^{-1} \left(\frac{\sqrt[3]{1 - 5 \cdot \left(v \cdot v\right)} \cdot \sqrt[3]{1 - 5 \cdot \left(v \cdot v\right)}}{\sqrt[3]{v \cdot v - 1} \cdot \sqrt[3]{v \cdot v - 1}} \cdot \frac{\sqrt[3]{1 - 5 \cdot \left(v \cdot v\right)}}{\sqrt[3]{v \cdot v - 1}}\right)(FPCore (v) :precision binary64 (acos (/ (- 1.0 (* 5.0 (* v v))) (- (* v v) 1.0))))
(FPCore (v)
:precision binary64
(-
(/ PI 2.0)
(asin
(*
(/
(* (cbrt (- 1.0 (* 5.0 (* v v)))) (cbrt (- 1.0 (* 5.0 (* v v)))))
(* (cbrt (- (* v v) 1.0)) (cbrt (- (* v v) 1.0))))
(/ (cbrt (- 1.0 (* 5.0 (* v v)))) (cbrt (- (* v v) 1.0)))))))double code(double v) {
return acos((1.0 - (5.0 * (v * v))) / ((v * v) - 1.0));
}
double code(double v) {
return (((double) M_PI) / 2.0) - asin(((cbrt(1.0 - (5.0 * (v * v))) * cbrt(1.0 - (5.0 * (v * v)))) / (cbrt((v * v) - 1.0) * cbrt((v * v) - 1.0))) * (cbrt(1.0 - (5.0 * (v * v))) / cbrt((v * v) - 1.0)));
}



Bits error versus v
Results
Initial program 0.6
rmApplied acos-asin_binary640.6
rmApplied add-cube-cbrt_binary640.6
Applied add-cube-cbrt_binary640.7
Applied times-frac_binary640.7
Final simplification0.7
herbie shell --seed 2020232
(FPCore (v)
:name "Falkner and Boettcher, Appendix B, 1"
:precision binary64
(acos (/ (- 1.0 (* 5.0 (* v v))) (- (* v v) 1.0))))