\sqrt[3]{\frac{g}{2 \cdot a}}\sqrt[3]{g} \cdot \left(\sqrt[3]{\frac{\sqrt[3]{\frac{1}{2}}}{a}} \cdot \sqrt[3]{\sqrt[3]{\frac{1}{2}} \cdot \sqrt[3]{\frac{1}{2}}}\right)double f(double g, double a) {
double r5110833 = g;
double r5110834 = 2.0;
double r5110835 = a;
double r5110836 = r5110834 * r5110835;
double r5110837 = r5110833 / r5110836;
double r5110838 = cbrt(r5110837);
return r5110838;
}
double f(double g, double a) {
double r5110839 = g;
double r5110840 = cbrt(r5110839);
double r5110841 = 0.5;
double r5110842 = cbrt(r5110841);
double r5110843 = a;
double r5110844 = r5110842 / r5110843;
double r5110845 = cbrt(r5110844);
double r5110846 = r5110842 * r5110842;
double r5110847 = cbrt(r5110846);
double r5110848 = r5110845 * r5110847;
double r5110849 = r5110840 * r5110848;
return r5110849;
}



Bits error versus g



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