double f(double g, double a) {
double r24928556 = g;
double r24928557 = 2.0;
double r24928558 = a;
double r24928559 = r24928557 * r24928558;
double r24928560 = r24928556 / r24928559;
double r24928561 = cbrt(r24928560);
return r24928561;
}
double f(double g, double a) {
double r24928562 = g;
double r24928563 = cbrt(r24928562);
double r24928564 = 0.5;
double r24928565 = cbrt(r24928564);
double r24928566 = a;
double r24928567 = cbrt(r24928566);
double r24928568 = r24928565 / r24928567;
double r24928569 = r24928563 * r24928568;
return r24928569;
}
\sqrt[3]{\frac{g}{2 \cdot a}}\sqrt[3]{g} \cdot \frac{\sqrt[3]{\frac{1}{2}}}{\sqrt[3]{a}}


Bits error versus g



Bits error versus a
Initial program 15.1
rmApplied div-inv15.1
Applied cbrt-prod0.9
rmApplied *-un-lft-identity0.9
Applied times-frac0.8
Applied cbrt-prod0.8
Simplified0.8
rmApplied cbrt-div0.9
Applied associate-*r/0.9
Final simplification0.9
herbie shell --seed 2019102
(FPCore (g a)
:name "2-ancestry mixing, zero discriminant"
(cbrt (/ g (* 2 a))))