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 r140342 = 2.0;
double r140343 = atan2(1.0, 0.0);
double r140344 = r140342 * r140343;
double r140345 = 3.0;
double r140346 = r140344 / r140345;
double r140347 = g;
double r140348 = -r140347;
double r140349 = h;
double r140350 = r140348 / r140349;
double r140351 = acos(r140350);
double r140352 = r140351 / r140345;
double r140353 = r140346 + r140352;
double r140354 = cos(r140353);
double r140355 = r140342 * r140354;
return r140355;
}
double f(double g, double h) {
double r140356 = 2.0;
double r140357 = atan2(1.0, 0.0);
double r140358 = r140356 * r140357;
double r140359 = 3.0;
double r140360 = r140358 / r140359;
double r140361 = r140357 / r140359;
double r140362 = r140360 + r140361;
double r140363 = cos(r140362);
double r140364 = g;
double r140365 = h;
double r140366 = r140364 / r140365;
double r140367 = acos(r140366);
double r140368 = r140367 / r140359;
double r140369 = cos(r140368);
double r140370 = r140363 * r140369;
double r140371 = sin(r140362);
double r140372 = sin(r140368);
double r140373 = r140371 * r140372;
double r140374 = r140370 + r140373;
double r140375 = r140356 * r140374;
return r140375;
}



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