\sqrt[3]{\frac{g}{2 \cdot a}}\sqrt[3]{\frac{1}{2}} \cdot \left(\sqrt[3]{g} \cdot \sqrt[3]{\frac{1}{a}}\right)double f(double g, double a) {
double r119560 = g;
double r119561 = 2.0;
double r119562 = a;
double r119563 = r119561 * r119562;
double r119564 = r119560 / r119563;
double r119565 = cbrt(r119564);
return r119565;
}
double f(double g, double a) {
double r119566 = 1.0;
double r119567 = 2.0;
double r119568 = r119566 / r119567;
double r119569 = cbrt(r119568);
double r119570 = g;
double r119571 = cbrt(r119570);
double r119572 = a;
double r119573 = r119566 / r119572;
double r119574 = cbrt(r119573);
double r119575 = r119571 * r119574;
double r119576 = r119569 * r119575;
return r119576;
}



Bits error versus g



Bits error versus a
Results
Initial program 15.5
rmApplied div-inv15.5
Applied cbrt-prod0.8
rmApplied *-un-lft-identity0.8
Applied times-frac0.8
Applied cbrt-prod0.8
Applied associate-*r*0.8
Simplified0.8
rmApplied associate-*l*0.9
Final simplification0.9
herbie shell --seed 2020047 +o rules:numerics
(FPCore (g a)
:name "2-ancestry mixing, zero discriminant"
:precision binary64
(cbrt (/ g (* 2 a))))