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 r131810 = 2.0;
double r131811 = atan2(1.0, 0.0);
double r131812 = r131810 * r131811;
double r131813 = 3.0;
double r131814 = r131812 / r131813;
double r131815 = g;
double r131816 = -r131815;
double r131817 = h;
double r131818 = r131816 / r131817;
double r131819 = acos(r131818);
double r131820 = r131819 / r131813;
double r131821 = r131814 + r131820;
double r131822 = cos(r131821);
double r131823 = r131810 * r131822;
return r131823;
}
double f(double g, double h) {
double r131824 = 2.0;
double r131825 = atan2(1.0, 0.0);
double r131826 = r131824 * r131825;
double r131827 = 3.0;
double r131828 = r131826 / r131827;
double r131829 = r131825 / r131827;
double r131830 = r131828 + r131829;
double r131831 = cos(r131830);
double r131832 = g;
double r131833 = h;
double r131834 = r131832 / r131833;
double r131835 = acos(r131834);
double r131836 = r131835 / r131827;
double r131837 = cos(r131836);
double r131838 = r131831 * r131837;
double r131839 = sin(r131830);
double r131840 = sin(r131836);
double r131841 = r131839 * r131840;
double r131842 = r131838 + r131841;
double r131843 = r131824 * r131842;
return r131843;
}



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 2019304
(FPCore (g h)
:name "2-ancestry mixing, negative discriminant"
:precision binary64
(* 2 (cos (+ (/ (* 2 PI) 3) (/ (acos (/ (- g) h)) 3)))))