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{\pi}{3} + \frac{2 \cdot \pi}{3}\right) \cdot \cos \left(\frac{\cos^{-1} \left(\frac{g}{h}\right)}{3}\right) + \sin \left(\frac{\pi}{3} + \frac{2 \cdot \pi}{3}\right) \cdot \sin \left(\frac{\cos^{-1} \left(\frac{g}{h}\right)}{3}\right)\right)double f(double g, double h) {
double r5954067 = 2.0;
double r5954068 = atan2(1.0, 0.0);
double r5954069 = r5954067 * r5954068;
double r5954070 = 3.0;
double r5954071 = r5954069 / r5954070;
double r5954072 = g;
double r5954073 = -r5954072;
double r5954074 = h;
double r5954075 = r5954073 / r5954074;
double r5954076 = acos(r5954075);
double r5954077 = r5954076 / r5954070;
double r5954078 = r5954071 + r5954077;
double r5954079 = cos(r5954078);
double r5954080 = r5954067 * r5954079;
return r5954080;
}
double f(double g, double h) {
double r5954081 = 2.0;
double r5954082 = atan2(1.0, 0.0);
double r5954083 = 3.0;
double r5954084 = r5954082 / r5954083;
double r5954085 = r5954081 * r5954082;
double r5954086 = r5954085 / r5954083;
double r5954087 = r5954084 + r5954086;
double r5954088 = cos(r5954087);
double r5954089 = g;
double r5954090 = h;
double r5954091 = r5954089 / r5954090;
double r5954092 = acos(r5954091);
double r5954093 = r5954092 / r5954083;
double r5954094 = cos(r5954093);
double r5954095 = r5954088 * r5954094;
double r5954096 = sin(r5954087);
double r5954097 = sin(r5954093);
double r5954098 = r5954096 * r5954097;
double r5954099 = r5954095 + r5954098;
double r5954100 = r5954081 * r5954099;
return r5954100;
}



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