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{\pi}{3} + \frac{2 \cdot \pi}{3}\right) \cdot \cos \left(\frac{\cos^{-1} \left(\frac{g}{h}\right)}{3}\right) + \sin \left(\frac{\pi}{3} + \frac{2 \cdot \pi}{3}\right) \cdot \sin \left(\frac{\cos^{-1} \left(\frac{g}{h}\right)}{3}\right)\right)double f(double g, double h) {
double r5886124 = 2.0;
double r5886125 = atan2(1.0, 0.0);
double r5886126 = r5886124 * r5886125;
double r5886127 = 3.0;
double r5886128 = r5886126 / r5886127;
double r5886129 = g;
double r5886130 = -r5886129;
double r5886131 = h;
double r5886132 = r5886130 / r5886131;
double r5886133 = acos(r5886132);
double r5886134 = r5886133 / r5886127;
double r5886135 = r5886128 + r5886134;
double r5886136 = cos(r5886135);
double r5886137 = r5886124 * r5886136;
return r5886137;
}
double f(double g, double h) {
double r5886138 = 2.0;
double r5886139 = atan2(1.0, 0.0);
double r5886140 = 3.0;
double r5886141 = r5886139 / r5886140;
double r5886142 = r5886138 * r5886139;
double r5886143 = r5886142 / r5886140;
double r5886144 = r5886141 + r5886143;
double r5886145 = cos(r5886144);
double r5886146 = g;
double r5886147 = h;
double r5886148 = r5886146 / r5886147;
double r5886149 = acos(r5886148);
double r5886150 = r5886149 / r5886140;
double r5886151 = cos(r5886150);
double r5886152 = r5886145 * r5886151;
double r5886153 = sin(r5886144);
double r5886154 = sin(r5886150);
double r5886155 = r5886153 * r5886154;
double r5886156 = r5886152 + r5886155;
double r5886157 = r5886138 * r5886156;
return r5886157;
}



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