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 r102220 = 2.0;
double r102221 = atan2(1.0, 0.0);
double r102222 = r102220 * r102221;
double r102223 = 3.0;
double r102224 = r102222 / r102223;
double r102225 = g;
double r102226 = -r102225;
double r102227 = h;
double r102228 = r102226 / r102227;
double r102229 = acos(r102228);
double r102230 = r102229 / r102223;
double r102231 = r102224 + r102230;
double r102232 = cos(r102231);
double r102233 = r102220 * r102232;
return r102233;
}
double f(double g, double h) {
double r102234 = 2.0;
double r102235 = atan2(1.0, 0.0);
double r102236 = r102234 * r102235;
double r102237 = 3.0;
double r102238 = r102236 / r102237;
double r102239 = r102235 / r102237;
double r102240 = r102238 + r102239;
double r102241 = cos(r102240);
double r102242 = g;
double r102243 = h;
double r102244 = r102242 / r102243;
double r102245 = acos(r102244);
double r102246 = r102245 / r102237;
double r102247 = cos(r102246);
double r102248 = r102241 * r102247;
double r102249 = sin(r102240);
double r102250 = sin(r102246);
double r102251 = r102249 * r102250;
double r102252 = r102248 + r102251;
double r102253 = r102234 * r102252;
return r102253;
}



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