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 r144589 = 2.0;
double r144590 = atan2(1.0, 0.0);
double r144591 = r144589 * r144590;
double r144592 = 3.0;
double r144593 = r144591 / r144592;
double r144594 = g;
double r144595 = -r144594;
double r144596 = h;
double r144597 = r144595 / r144596;
double r144598 = acos(r144597);
double r144599 = r144598 / r144592;
double r144600 = r144593 + r144599;
double r144601 = cos(r144600);
double r144602 = r144589 * r144601;
return r144602;
}
double f(double g, double h) {
double r144603 = 2.0;
double r144604 = atan2(1.0, 0.0);
double r144605 = r144603 * r144604;
double r144606 = 3.0;
double r144607 = r144605 / r144606;
double r144608 = r144604 / r144606;
double r144609 = r144607 + r144608;
double r144610 = cos(r144609);
double r144611 = g;
double r144612 = h;
double r144613 = r144611 / r144612;
double r144614 = acos(r144613);
double r144615 = r144614 / r144606;
double r144616 = cos(r144615);
double r144617 = r144610 * r144616;
double r144618 = sin(r144609);
double r144619 = sin(r144615);
double r144620 = r144618 * r144619;
double r144621 = r144617 + r144620;
double r144622 = r144603 * r144621;
return r144622;
}



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