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 r120459 = 2.0;
double r120460 = atan2(1.0, 0.0);
double r120461 = r120459 * r120460;
double r120462 = 3.0;
double r120463 = r120461 / r120462;
double r120464 = g;
double r120465 = -r120464;
double r120466 = h;
double r120467 = r120465 / r120466;
double r120468 = acos(r120467);
double r120469 = r120468 / r120462;
double r120470 = r120463 + r120469;
double r120471 = cos(r120470);
double r120472 = r120459 * r120471;
return r120472;
}
double f(double g, double h) {
double r120473 = 2.0;
double r120474 = atan2(1.0, 0.0);
double r120475 = r120473 * r120474;
double r120476 = 3.0;
double r120477 = r120475 / r120476;
double r120478 = r120474 / r120476;
double r120479 = r120477 + r120478;
double r120480 = cos(r120479);
double r120481 = g;
double r120482 = h;
double r120483 = r120481 / r120482;
double r120484 = acos(r120483);
double r120485 = r120484 / r120476;
double r120486 = cos(r120485);
double r120487 = r120480 * r120486;
double r120488 = sin(r120479);
double r120489 = sin(r120485);
double r120490 = r120488 * r120489;
double r120491 = r120487 + r120490;
double r120492 = r120473 * r120491;
return r120492;
}



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.1
Final simplification0.1
herbie shell --seed 2020043
(FPCore (g h)
:name "2-ancestry mixing, negative discriminant"
:precision binary64
(* 2 (cos (+ (/ (* 2 PI) 3) (/ (acos (/ (- g) h)) 3)))))