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 r22107 = 2.0;
double r22108 = atan2(1.0, 0.0);
double r22109 = r22107 * r22108;
double r22110 = 3.0;
double r22111 = r22109 / r22110;
double r22112 = g;
double r22113 = -r22112;
double r22114 = h;
double r22115 = r22113 / r22114;
double r22116 = acos(r22115);
double r22117 = r22116 / r22110;
double r22118 = r22111 + r22117;
double r22119 = cos(r22118);
double r22120 = r22107 * r22119;
return r22120;
}
double f(double g, double h) {
double r22121 = 2.0;
double r22122 = atan2(1.0, 0.0);
double r22123 = r22121 * r22122;
double r22124 = 3.0;
double r22125 = r22123 / r22124;
double r22126 = r22122 / r22124;
double r22127 = r22125 + r22126;
double r22128 = cos(r22127);
double r22129 = g;
double r22130 = h;
double r22131 = r22129 / r22130;
double r22132 = acos(r22131);
double r22133 = r22132 / r22124;
double r22134 = cos(r22133);
double r22135 = r22128 * r22134;
double r22136 = sin(r22127);
double r22137 = sin(r22133);
double r22138 = r22136 * r22137;
double r22139 = r22135 + r22138;
double r22140 = r22121 * r22139;
return r22140;
}



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