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 r149949 = 2.0;
double r149950 = atan2(1.0, 0.0);
double r149951 = r149949 * r149950;
double r149952 = 3.0;
double r149953 = r149951 / r149952;
double r149954 = g;
double r149955 = -r149954;
double r149956 = h;
double r149957 = r149955 / r149956;
double r149958 = acos(r149957);
double r149959 = r149958 / r149952;
double r149960 = r149953 + r149959;
double r149961 = cos(r149960);
double r149962 = r149949 * r149961;
return r149962;
}
double f(double g, double h) {
double r149963 = 2.0;
double r149964 = atan2(1.0, 0.0);
double r149965 = r149963 * r149964;
double r149966 = 3.0;
double r149967 = r149965 / r149966;
double r149968 = r149964 / r149966;
double r149969 = r149967 + r149968;
double r149970 = cos(r149969);
double r149971 = g;
double r149972 = h;
double r149973 = r149971 / r149972;
double r149974 = acos(r149973);
double r149975 = r149974 / r149966;
double r149976 = cos(r149975);
double r149977 = r149970 * r149976;
double r149978 = sin(r149969);
double r149979 = sin(r149975);
double r149980 = r149978 * r149979;
double r149981 = r149977 + r149980;
double r149982 = r149963 * r149981;
return r149982;
}



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