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 r104977 = 2.0;
double r104978 = atan2(1.0, 0.0);
double r104979 = r104977 * r104978;
double r104980 = 3.0;
double r104981 = r104979 / r104980;
double r104982 = g;
double r104983 = -r104982;
double r104984 = h;
double r104985 = r104983 / r104984;
double r104986 = acos(r104985);
double r104987 = r104986 / r104980;
double r104988 = r104981 + r104987;
double r104989 = cos(r104988);
double r104990 = r104977 * r104989;
return r104990;
}
double f(double g, double h) {
double r104991 = 2.0;
double r104992 = atan2(1.0, 0.0);
double r104993 = r104991 * r104992;
double r104994 = 3.0;
double r104995 = r104993 / r104994;
double r104996 = r104992 / r104994;
double r104997 = r104995 + r104996;
double r104998 = cos(r104997);
double r104999 = g;
double r105000 = h;
double r105001 = r104999 / r105000;
double r105002 = acos(r105001);
double r105003 = r105002 / r104994;
double r105004 = cos(r105003);
double r105005 = r104998 * r105004;
double r105006 = sin(r104997);
double r105007 = sin(r105003);
double r105008 = r105006 * r105007;
double r105009 = r105005 + r105008;
double r105010 = r104991 * r105009;
return r105010;
}



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