\sqrt[3]{\frac{g}{2 \cdot a}}\frac{1}{\frac{\sqrt[3]{2 \cdot a}}{\sqrt[3]{g}}}double f(double g, double a) {
double r116413 = g;
double r116414 = 2.0;
double r116415 = a;
double r116416 = r116414 * r116415;
double r116417 = r116413 / r116416;
double r116418 = cbrt(r116417);
return r116418;
}
double f(double g, double a) {
double r116419 = 1.0;
double r116420 = 2.0;
double r116421 = a;
double r116422 = r116420 * r116421;
double r116423 = cbrt(r116422);
double r116424 = g;
double r116425 = cbrt(r116424);
double r116426 = r116423 / r116425;
double r116427 = r116419 / r116426;
return r116427;
}



Bits error versus g



Bits error versus a
Results
Initial program 15.9
rmApplied cbrt-div0.8
rmApplied clear-num0.9
Final simplification0.9
herbie shell --seed 2020046
(FPCore (g a)
:name "2-ancestry mixing, zero discriminant"
:precision binary64
(cbrt (/ g (* 2 a))))