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 r107757 = 2.0;
double r107758 = atan2(1.0, 0.0);
double r107759 = r107757 * r107758;
double r107760 = 3.0;
double r107761 = r107759 / r107760;
double r107762 = g;
double r107763 = -r107762;
double r107764 = h;
double r107765 = r107763 / r107764;
double r107766 = acos(r107765);
double r107767 = r107766 / r107760;
double r107768 = r107761 + r107767;
double r107769 = cos(r107768);
double r107770 = r107757 * r107769;
return r107770;
}
double f(double g, double h) {
double r107771 = 2.0;
double r107772 = atan2(1.0, 0.0);
double r107773 = r107771 * r107772;
double r107774 = 3.0;
double r107775 = r107773 / r107774;
double r107776 = r107772 / r107774;
double r107777 = r107775 + r107776;
double r107778 = cos(r107777);
double r107779 = g;
double r107780 = h;
double r107781 = r107779 / r107780;
double r107782 = acos(r107781);
double r107783 = r107782 / r107774;
double r107784 = cos(r107783);
double r107785 = r107778 * r107784;
double r107786 = sin(r107777);
double r107787 = sin(r107783);
double r107788 = r107786 * r107787;
double r107789 = r107785 + r107788;
double r107790 = r107771 * r107789;
return r107790;
}



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