\cos^{-1} \left(\frac{1 - 5 \cdot \left(v \cdot v\right)}{v \cdot v - 1}\right)\sqrt{\cos^{-1} \left(\frac{\mathsf{fma}\left(-5 \cdot v, v, 1\right)}{\mathsf{fma}\left(v, v, -1\right)}\right)} \cdot \left(\sqrt{\sqrt[3]{\sqrt{\cos^{-1} \left(\frac{\mathsf{fma}\left(-5 \cdot v, v, 1\right)}{\mathsf{fma}\left(v, v, -1\right)}\right)} \cdot \cos^{-1} \left(\frac{\mathsf{fma}\left(-5 \cdot v, v, 1\right)}{\mathsf{fma}\left(v, v, -1\right)}\right)}} \cdot \sqrt{\sqrt[3]{\sqrt{\cos^{-1} \left(\frac{\mathsf{fma}\left(-5 \cdot v, v, 1\right)}{\mathsf{fma}\left(v, v, -1\right)}\right)} \cdot \cos^{-1} \left(\frac{\mathsf{fma}\left(-5 \cdot v, v, 1\right)}{\mathsf{fma}\left(v, v, -1\right)}\right)}}\right)double f(double v) {
double r3800421 = 1.0;
double r3800422 = 5.0;
double r3800423 = v;
double r3800424 = r3800423 * r3800423;
double r3800425 = r3800422 * r3800424;
double r3800426 = r3800421 - r3800425;
double r3800427 = r3800424 - r3800421;
double r3800428 = r3800426 / r3800427;
double r3800429 = acos(r3800428);
return r3800429;
}
double f(double v) {
double r3800430 = -5.0;
double r3800431 = v;
double r3800432 = r3800430 * r3800431;
double r3800433 = 1.0;
double r3800434 = fma(r3800432, r3800431, r3800433);
double r3800435 = -1.0;
double r3800436 = fma(r3800431, r3800431, r3800435);
double r3800437 = r3800434 / r3800436;
double r3800438 = acos(r3800437);
double r3800439 = sqrt(r3800438);
double r3800440 = r3800439 * r3800438;
double r3800441 = cbrt(r3800440);
double r3800442 = sqrt(r3800441);
double r3800443 = r3800442 * r3800442;
double r3800444 = r3800439 * r3800443;
return r3800444;
}



Bits error versus v
Initial program 0.6
Simplified0.6
rmApplied add-sqr-sqrt1.5
rmApplied add-cbrt-cube1.5
Simplified1.5
rmApplied add-sqr-sqrt0.6
Final simplification0.6
herbie shell --seed 2019154 +o rules:numerics
(FPCore (v)
:name "Falkner and Boettcher, Appendix B, 1"
(acos (/ (- 1 (* 5 (* v v))) (- (* v v) 1))))