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 r194517 = 2.0;
double r194518 = atan2(1.0, 0.0);
double r194519 = r194517 * r194518;
double r194520 = 3.0;
double r194521 = r194519 / r194520;
double r194522 = g;
double r194523 = -r194522;
double r194524 = h;
double r194525 = r194523 / r194524;
double r194526 = acos(r194525);
double r194527 = r194526 / r194520;
double r194528 = r194521 + r194527;
double r194529 = cos(r194528);
double r194530 = r194517 * r194529;
return r194530;
}
double f(double g, double h) {
double r194531 = 2.0;
double r194532 = atan2(1.0, 0.0);
double r194533 = r194531 * r194532;
double r194534 = 3.0;
double r194535 = r194533 / r194534;
double r194536 = r194532 / r194534;
double r194537 = r194535 + r194536;
double r194538 = cos(r194537);
double r194539 = g;
double r194540 = h;
double r194541 = r194539 / r194540;
double r194542 = acos(r194541);
double r194543 = r194542 / r194534;
double r194544 = cos(r194543);
double r194545 = r194538 * r194544;
double r194546 = sin(r194537);
double r194547 = sin(r194543);
double r194548 = r194546 * r194547;
double r194549 = r194545 + r194548;
double r194550 = r194531 * r194549;
return r194550;
}



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