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 r154440 = 2.0;
double r154441 = atan2(1.0, 0.0);
double r154442 = r154440 * r154441;
double r154443 = 3.0;
double r154444 = r154442 / r154443;
double r154445 = g;
double r154446 = -r154445;
double r154447 = h;
double r154448 = r154446 / r154447;
double r154449 = acos(r154448);
double r154450 = r154449 / r154443;
double r154451 = r154444 + r154450;
double r154452 = cos(r154451);
double r154453 = r154440 * r154452;
return r154453;
}
double f(double g, double h) {
double r154454 = 2.0;
double r154455 = atan2(1.0, 0.0);
double r154456 = r154454 * r154455;
double r154457 = 3.0;
double r154458 = r154456 / r154457;
double r154459 = r154455 / r154457;
double r154460 = r154458 + r154459;
double r154461 = cos(r154460);
double r154462 = g;
double r154463 = h;
double r154464 = r154462 / r154463;
double r154465 = acos(r154464);
double r154466 = r154465 / r154457;
double r154467 = cos(r154466);
double r154468 = r154461 * r154467;
double r154469 = sin(r154460);
double r154470 = sin(r154466);
double r154471 = r154469 * r154470;
double r154472 = r154468 + r154471;
double r154473 = r154454 * r154472;
return r154473;
}



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