\sqrt[3]{\frac{g}{2 \cdot a}}\left(\sqrt[3]{g} \cdot \sqrt[3]{\frac{1}{a}}\right) \cdot \sqrt[3]{0.5}double f(double g, double a) {
double r5293861 = g;
double r5293862 = 2.0;
double r5293863 = a;
double r5293864 = r5293862 * r5293863;
double r5293865 = r5293861 / r5293864;
double r5293866 = cbrt(r5293865);
return r5293866;
}
double f(double g, double a) {
double r5293867 = g;
double r5293868 = cbrt(r5293867);
double r5293869 = 1.0;
double r5293870 = a;
double r5293871 = r5293869 / r5293870;
double r5293872 = cbrt(r5293871);
double r5293873 = r5293868 * r5293872;
double r5293874 = 0.5;
double r5293875 = cbrt(r5293874);
double r5293876 = r5293873 * r5293875;
return r5293876;
}



Bits error versus g



Bits error versus a
Results
Initial program 16.0
rmApplied cbrt-div0.8
rmApplied div-inv0.9
rmApplied add-cbrt-cube0.9
Applied cbrt-undiv0.9
Simplified0.8
Taylor expanded around 0 49.1
Simplified0.8
Final simplification0.8
herbie shell --seed 2019168 +o rules:numerics
(FPCore (g a)
:name "2-ancestry mixing, zero discriminant"
(cbrt (/ g (* 2.0 a))))