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 r143185 = 2.0;
double r143186 = atan2(1.0, 0.0);
double r143187 = r143185 * r143186;
double r143188 = 3.0;
double r143189 = r143187 / r143188;
double r143190 = g;
double r143191 = -r143190;
double r143192 = h;
double r143193 = r143191 / r143192;
double r143194 = acos(r143193);
double r143195 = r143194 / r143188;
double r143196 = r143189 + r143195;
double r143197 = cos(r143196);
double r143198 = r143185 * r143197;
return r143198;
}
double f(double g, double h) {
double r143199 = 2.0;
double r143200 = atan2(1.0, 0.0);
double r143201 = r143199 * r143200;
double r143202 = 3.0;
double r143203 = r143201 / r143202;
double r143204 = r143200 / r143202;
double r143205 = r143203 + r143204;
double r143206 = cos(r143205);
double r143207 = g;
double r143208 = h;
double r143209 = r143207 / r143208;
double r143210 = acos(r143209);
double r143211 = r143210 / r143202;
double r143212 = cos(r143211);
double r143213 = r143206 * r143212;
double r143214 = sin(r143205);
double r143215 = sin(r143211);
double r143216 = r143214 * r143215;
double r143217 = r143213 + r143216;
double r143218 = r143199 * r143217;
return r143218;
}



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