\sqrt[3]{\frac{g}{2 \cdot a}}\left(\sqrt[3]{\sqrt[3]{g} \cdot \sqrt[3]{g}} \cdot \sqrt[3]{\frac{1}{2}}\right) \cdot \frac{\sqrt[3]{\sqrt[3]{\sqrt[3]{g} \cdot \sqrt[3]{g}} \cdot \sqrt[3]{\sqrt[3]{g}}}}{\sqrt[3]{a}}double f(double g, double a) {
double r117139 = g;
double r117140 = 2.0;
double r117141 = a;
double r117142 = r117140 * r117141;
double r117143 = r117139 / r117142;
double r117144 = cbrt(r117143);
return r117144;
}
double f(double g, double a) {
double r117145 = g;
double r117146 = cbrt(r117145);
double r117147 = r117146 * r117146;
double r117148 = cbrt(r117147);
double r117149 = 1.0;
double r117150 = 2.0;
double r117151 = r117149 / r117150;
double r117152 = cbrt(r117151);
double r117153 = r117148 * r117152;
double r117154 = cbrt(r117146);
double r117155 = r117148 * r117154;
double r117156 = cbrt(r117155);
double r117157 = a;
double r117158 = cbrt(r117157);
double r117159 = r117156 / r117158;
double r117160 = r117153 * r117159;
return r117160;
}



Bits error versus g



Bits error versus a
Results
Initial program 15.9
rmApplied add-cube-cbrt16.1
Applied times-frac16.1
Applied cbrt-prod5.7
rmApplied cbrt-div1.1
rmApplied div-inv1.1
Applied cbrt-prod1.2
rmApplied add-cube-cbrt1.2
Applied cbrt-prod1.2
Final simplification1.2
herbie shell --seed 2020033 +o rules:numerics
(FPCore (g a)
:name "2-ancestry mixing, zero discriminant"
:precision binary64
(cbrt (/ g (* 2 a))))