2 \cdot \cos \left(\frac{2 \cdot \pi}{3} + \frac{\cos^{-1} \left(\frac{-g}{h}\right)}{3}\right)2 \cdot \left(\cos \left(\frac{2 \cdot \pi}{3} + \frac{\pi}{3}\right) \cdot \cos \left(\frac{\cos^{-1} \left(\frac{g}{h}\right)}{3}\right) + \sin \left(\frac{2 \cdot \pi}{3} + \frac{\pi}{3}\right) \cdot \sin \left(\frac{\cos^{-1} \left(\frac{g}{h}\right)}{3}\right)\right)double f(double g, double h) {
double r162816 = 2.0;
double r162817 = atan2(1.0, 0.0);
double r162818 = r162816 * r162817;
double r162819 = 3.0;
double r162820 = r162818 / r162819;
double r162821 = g;
double r162822 = -r162821;
double r162823 = h;
double r162824 = r162822 / r162823;
double r162825 = acos(r162824);
double r162826 = r162825 / r162819;
double r162827 = r162820 + r162826;
double r162828 = cos(r162827);
double r162829 = r162816 * r162828;
return r162829;
}
double f(double g, double h) {
double r162830 = 2.0;
double r162831 = atan2(1.0, 0.0);
double r162832 = r162830 * r162831;
double r162833 = 3.0;
double r162834 = r162832 / r162833;
double r162835 = r162831 / r162833;
double r162836 = r162834 + r162835;
double r162837 = cos(r162836);
double r162838 = g;
double r162839 = h;
double r162840 = r162838 / r162839;
double r162841 = acos(r162840);
double r162842 = r162841 / r162833;
double r162843 = cos(r162842);
double r162844 = r162837 * r162843;
double r162845 = sin(r162836);
double r162846 = sin(r162842);
double r162847 = r162845 * r162846;
double r162848 = r162844 + r162847;
double r162849 = r162830 * r162848;
return r162849;
}



Bits error versus g



Bits error versus h
Results
Initial program 1.0
rmApplied distribute-frac-neg1.0
Applied acos-neg1.0
Applied div-sub1.0
Applied associate-+r-1.0
Applied cos-diff0.0
Final simplification0.0
herbie shell --seed 2019351
(FPCore (g h)
:name "2-ancestry mixing, negative discriminant"
:precision binary64
(* 2 (cos (+ (/ (* 2 PI) 3) (/ (acos (/ (- g) h)) 3)))))