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 r279156 = 2.0;
double r279157 = atan2(1.0, 0.0);
double r279158 = r279156 * r279157;
double r279159 = 3.0;
double r279160 = r279158 / r279159;
double r279161 = g;
double r279162 = -r279161;
double r279163 = h;
double r279164 = r279162 / r279163;
double r279165 = acos(r279164);
double r279166 = r279165 / r279159;
double r279167 = r279160 + r279166;
double r279168 = cos(r279167);
double r279169 = r279156 * r279168;
return r279169;
}
double f(double g, double h) {
double r279170 = 2.0;
double r279171 = atan2(1.0, 0.0);
double r279172 = r279170 * r279171;
double r279173 = 3.0;
double r279174 = r279172 / r279173;
double r279175 = r279171 / r279173;
double r279176 = r279174 + r279175;
double r279177 = cos(r279176);
double r279178 = g;
double r279179 = h;
double r279180 = r279178 / r279179;
double r279181 = acos(r279180);
double r279182 = r279181 / r279173;
double r279183 = cos(r279182);
double r279184 = r279177 * r279183;
double r279185 = sin(r279176);
double r279186 = sin(r279182);
double r279187 = r279185 * r279186;
double r279188 = r279184 + r279187;
double r279189 = r279170 * r279188;
return r279189;
}



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