\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]{g}}}{\sqrt[3]{a}}double f(double g, double a) {
double r150183 = g;
double r150184 = 2.0;
double r150185 = a;
double r150186 = r150184 * r150185;
double r150187 = r150183 / r150186;
double r150188 = cbrt(r150187);
return r150188;
}
double f(double g, double a) {
double r150189 = g;
double r150190 = cbrt(r150189);
double r150191 = r150190 * r150190;
double r150192 = cbrt(r150191);
double r150193 = 1.0;
double r150194 = 2.0;
double r150195 = r150193 / r150194;
double r150196 = cbrt(r150195);
double r150197 = r150192 * r150196;
double r150198 = cbrt(r150190);
double r150199 = a;
double r150200 = cbrt(r150199);
double r150201 = r150198 / r150200;
double r150202 = r150197 * r150201;
return r150202;
}



Bits error versus g



Bits error versus a
Results
Initial program 15.1
rmApplied add-cube-cbrt15.3
Applied times-frac15.3
Applied cbrt-prod5.6
rmApplied cbrt-div1.2
rmApplied div-inv1.2
Applied cbrt-prod1.2
Final simplification1.2
herbie shell --seed 2020056 +o rules:numerics
(FPCore (g a)
:name "2-ancestry mixing, zero discriminant"
:precision binary64
(cbrt (/ g (* 2 a))))