2 \cdot \cos \left(\frac{2 \cdot \pi}{3} + \frac{\cos^{-1} \left(\frac{-g}{h}\right)}{3}\right)2 \cdot \sqrt[3]{{\left(\sqrt[3]{\cos \left(\frac{\frac{3}{2 \cdot \pi} \cdot \frac{\cos^{-1} \left(\frac{-g}{h}\right)}{\sqrt[3]{3} \cdot \sqrt[3]{3}} + \sqrt[3]{3}}{\frac{3}{2 \cdot \pi} \cdot \sqrt[3]{3}}\right) \cdot \cos \left(\frac{\frac{3}{2 \cdot \pi} \cdot \frac{\cos^{-1} \left(\frac{-g}{h}\right)}{\sqrt[3]{3} \cdot \sqrt[3]{3}} + \sqrt[3]{3}}{\frac{3}{2 \cdot \pi} \cdot \sqrt[3]{3}}\right)} \cdot \sqrt[3]{\cos \left(\frac{\frac{3}{2 \cdot \pi} \cdot \frac{\cos^{-1} \left(\frac{-g}{h}\right)}{\sqrt[3]{3} \cdot \sqrt[3]{3}} + \sqrt[3]{3}}{\frac{3}{2 \cdot \pi} \cdot \sqrt[3]{3}}\right)}\right)}^{3}}double code(double g, double h) {
return ((double) (2.0 * ((double) cos(((double) (((double) (((double) (2.0 * ((double) M_PI))) / 3.0)) + ((double) (((double) acos(((double) (((double) -(g)) / h)))) / 3.0))))))));
}
double code(double g, double h) {
return ((double) (2.0 * ((double) cbrt(((double) pow(((double) (((double) cbrt(((double) (((double) cos(((double) (((double) (((double) (((double) (3.0 / ((double) (2.0 * ((double) M_PI))))) * ((double) (((double) acos(((double) (((double) -(g)) / h)))) / ((double) (((double) cbrt(3.0)) * ((double) cbrt(3.0)))))))) + ((double) cbrt(3.0)))) / ((double) (((double) (3.0 / ((double) (2.0 * ((double) M_PI))))) * ((double) cbrt(3.0)))))))) * ((double) cos(((double) (((double) (((double) (((double) (3.0 / ((double) (2.0 * ((double) M_PI))))) * ((double) (((double) acos(((double) (((double) -(g)) / h)))) / ((double) (((double) cbrt(3.0)) * ((double) cbrt(3.0)))))))) + ((double) cbrt(3.0)))) / ((double) (((double) (3.0 / ((double) (2.0 * ((double) M_PI))))) * ((double) cbrt(3.0)))))))))))) * ((double) cbrt(((double) cos(((double) (((double) (((double) (((double) (3.0 / ((double) (2.0 * ((double) M_PI))))) * ((double) (((double) acos(((double) (((double) -(g)) / h)))) / ((double) (((double) cbrt(3.0)) * ((double) cbrt(3.0)))))))) + ((double) cbrt(3.0)))) / ((double) (((double) (3.0 / ((double) (2.0 * ((double) M_PI))))) * ((double) cbrt(3.0)))))))))))), 3.0))))));
}



Bits error versus g



Bits error versus h
Results
Initial program 1.0
rmApplied add-cube-cbrt1.0
Applied associate-/r*1.0
Applied clear-num1.0
Applied frac-add1.0
Simplified1.0
rmApplied add-cbrt-cube1.0
Simplified1.0
rmApplied add-cube-cbrt1.0
Simplified0.1
Final simplification0.1
herbie shell --seed 2020114
(FPCore (g h)
:name "2-ancestry mixing, negative discriminant"
:precision binary64
(* 2 (cos (+ (/ (* 2 PI) 3) (/ (acos (/ (- g) h)) 3)))))