2 \cdot \cos \left(\frac{2 \cdot \pi}{3} + \frac{\cos^{-1} \left(\frac{-g}{h}\right)}{3}\right)\left(\frac{\sqrt{3} \cdot \sin \left(\frac{\cos^{-1} \left(\frac{g}{h}\right)}{3} - \left(\sqrt[3]{\frac{\pi}{\frac{3}{2}}} \cdot \sqrt[3]{\frac{\pi}{\frac{3}{2}}}\right) \cdot \sqrt[3]{\frac{\pi}{\frac{3}{2}}}\right)}{2} + \cos \left(\frac{\cos^{-1} \left(\frac{g}{h}\right)}{3} - \frac{\pi}{\frac{3}{2}}\right) \cdot \frac{1}{2}\right) \cdot 2double f(double g, double h) {
double r4801535 = 2.0;
double r4801536 = atan2(1.0, 0.0);
double r4801537 = r4801535 * r4801536;
double r4801538 = 3.0;
double r4801539 = r4801537 / r4801538;
double r4801540 = g;
double r4801541 = -r4801540;
double r4801542 = h;
double r4801543 = r4801541 / r4801542;
double r4801544 = acos(r4801543);
double r4801545 = r4801544 / r4801538;
double r4801546 = r4801539 + r4801545;
double r4801547 = cos(r4801546);
double r4801548 = r4801535 * r4801547;
return r4801548;
}
double f(double g, double h) {
double r4801549 = 3.0;
double r4801550 = sqrt(r4801549);
double r4801551 = g;
double r4801552 = h;
double r4801553 = r4801551 / r4801552;
double r4801554 = acos(r4801553);
double r4801555 = r4801554 / r4801549;
double r4801556 = atan2(1.0, 0.0);
double r4801557 = 1.5;
double r4801558 = r4801556 / r4801557;
double r4801559 = cbrt(r4801558);
double r4801560 = r4801559 * r4801559;
double r4801561 = r4801560 * r4801559;
double r4801562 = r4801555 - r4801561;
double r4801563 = sin(r4801562);
double r4801564 = r4801550 * r4801563;
double r4801565 = 2.0;
double r4801566 = r4801564 / r4801565;
double r4801567 = r4801555 - r4801558;
double r4801568 = cos(r4801567);
double r4801569 = 0.5;
double r4801570 = r4801568 * r4801569;
double r4801571 = r4801566 + r4801570;
double r4801572 = r4801571 * r4801565;
return r4801572;
}



Bits error versus g



Bits error versus h
Results
Initial program 1.0
Simplified1.0
rmApplied distribute-frac-neg1.0
Applied acos-neg1.0
Applied div-sub1.0
Applied associate-+l-1.0
Applied cos-diff0.1
Simplified0.1
Simplified0.1
rmApplied add-cube-cbrt0.1
Final simplification0.1
herbie shell --seed 2019152
(FPCore (g h)
:name "2-ancestry mixing, negative discriminant"
(* 2 (cos (+ (/ (* 2 PI) 3) (/ (acos (/ (- g) h)) 3)))))