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 r138654 = 2.0;
double r138655 = atan2(1.0, 0.0);
double r138656 = r138654 * r138655;
double r138657 = 3.0;
double r138658 = r138656 / r138657;
double r138659 = g;
double r138660 = -r138659;
double r138661 = h;
double r138662 = r138660 / r138661;
double r138663 = acos(r138662);
double r138664 = r138663 / r138657;
double r138665 = r138658 + r138664;
double r138666 = cos(r138665);
double r138667 = r138654 * r138666;
return r138667;
}
double f(double g, double h) {
double r138668 = 2.0;
double r138669 = atan2(1.0, 0.0);
double r138670 = r138668 * r138669;
double r138671 = 3.0;
double r138672 = r138670 / r138671;
double r138673 = r138669 / r138671;
double r138674 = r138672 + r138673;
double r138675 = cos(r138674);
double r138676 = g;
double r138677 = h;
double r138678 = r138676 / r138677;
double r138679 = acos(r138678);
double r138680 = r138679 / r138671;
double r138681 = cos(r138680);
double r138682 = r138675 * r138681;
double r138683 = sin(r138674);
double r138684 = sin(r138680);
double r138685 = r138683 * r138684;
double r138686 = r138682 + r138685;
double r138687 = r138668 * r138686;
return r138687;
}



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