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 r98770 = 2.0;
double r98771 = atan2(1.0, 0.0);
double r98772 = r98770 * r98771;
double r98773 = 3.0;
double r98774 = r98772 / r98773;
double r98775 = g;
double r98776 = -r98775;
double r98777 = h;
double r98778 = r98776 / r98777;
double r98779 = acos(r98778);
double r98780 = r98779 / r98773;
double r98781 = r98774 + r98780;
double r98782 = cos(r98781);
double r98783 = r98770 * r98782;
return r98783;
}
double f(double g, double h) {
double r98784 = 2.0;
double r98785 = atan2(1.0, 0.0);
double r98786 = r98784 * r98785;
double r98787 = 3.0;
double r98788 = r98786 / r98787;
double r98789 = r98785 / r98787;
double r98790 = r98788 + r98789;
double r98791 = cos(r98790);
double r98792 = g;
double r98793 = h;
double r98794 = r98792 / r98793;
double r98795 = acos(r98794);
double r98796 = r98795 / r98787;
double r98797 = cos(r98796);
double r98798 = r98791 * r98797;
double r98799 = sin(r98790);
double r98800 = sin(r98796);
double r98801 = r98799 * r98800;
double r98802 = r98798 + r98801;
double r98803 = r98784 * r98802;
return r98803;
}



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