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 r166243 = 2.0;
double r166244 = atan2(1.0, 0.0);
double r166245 = r166243 * r166244;
double r166246 = 3.0;
double r166247 = r166245 / r166246;
double r166248 = g;
double r166249 = -r166248;
double r166250 = h;
double r166251 = r166249 / r166250;
double r166252 = acos(r166251);
double r166253 = r166252 / r166246;
double r166254 = r166247 + r166253;
double r166255 = cos(r166254);
double r166256 = r166243 * r166255;
return r166256;
}
double f(double g, double h) {
double r166257 = 2.0;
double r166258 = atan2(1.0, 0.0);
double r166259 = r166257 * r166258;
double r166260 = 3.0;
double r166261 = r166259 / r166260;
double r166262 = r166258 / r166260;
double r166263 = r166261 + r166262;
double r166264 = cos(r166263);
double r166265 = g;
double r166266 = h;
double r166267 = r166265 / r166266;
double r166268 = acos(r166267);
double r166269 = r166268 / r166260;
double r166270 = cos(r166269);
double r166271 = r166264 * r166270;
double r166272 = sin(r166263);
double r166273 = sin(r166269);
double r166274 = r166272 * r166273;
double r166275 = r166271 + r166274;
double r166276 = r166257 * r166275;
return r166276;
}



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