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 r118640 = 2.0;
double r118641 = atan2(1.0, 0.0);
double r118642 = r118640 * r118641;
double r118643 = 3.0;
double r118644 = r118642 / r118643;
double r118645 = g;
double r118646 = -r118645;
double r118647 = h;
double r118648 = r118646 / r118647;
double r118649 = acos(r118648);
double r118650 = r118649 / r118643;
double r118651 = r118644 + r118650;
double r118652 = cos(r118651);
double r118653 = r118640 * r118652;
return r118653;
}
double f(double g, double h) {
double r118654 = 2.0;
double r118655 = atan2(1.0, 0.0);
double r118656 = r118654 * r118655;
double r118657 = 3.0;
double r118658 = r118656 / r118657;
double r118659 = r118655 / r118657;
double r118660 = r118658 + r118659;
double r118661 = cos(r118660);
double r118662 = g;
double r118663 = h;
double r118664 = r118662 / r118663;
double r118665 = acos(r118664);
double r118666 = r118665 / r118657;
double r118667 = cos(r118666);
double r118668 = r118661 * r118667;
double r118669 = sin(r118660);
double r118670 = sin(r118666);
double r118671 = r118669 * r118670;
double r118672 = r118668 + r118671;
double r118673 = r118654 * r118672;
return r118673;
}



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