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 r76916 = 2.0;
double r76917 = atan2(1.0, 0.0);
double r76918 = r76916 * r76917;
double r76919 = 3.0;
double r76920 = r76918 / r76919;
double r76921 = g;
double r76922 = -r76921;
double r76923 = h;
double r76924 = r76922 / r76923;
double r76925 = acos(r76924);
double r76926 = r76925 / r76919;
double r76927 = r76920 + r76926;
double r76928 = cos(r76927);
double r76929 = r76916 * r76928;
return r76929;
}
double f(double g, double h) {
double r76930 = 2.0;
double r76931 = atan2(1.0, 0.0);
double r76932 = r76930 * r76931;
double r76933 = 3.0;
double r76934 = r76932 / r76933;
double r76935 = r76931 / r76933;
double r76936 = r76934 + r76935;
double r76937 = cos(r76936);
double r76938 = g;
double r76939 = h;
double r76940 = r76938 / r76939;
double r76941 = acos(r76940);
double r76942 = r76941 / r76933;
double r76943 = cos(r76942);
double r76944 = r76937 * r76943;
double r76945 = sin(r76936);
double r76946 = sin(r76942);
double r76947 = r76945 * r76946;
double r76948 = r76944 + r76947;
double r76949 = r76930 * r76948;
return r76949;
}



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