\sqrt[3]{\frac{g}{2 \cdot a}}\sqrt[3]{\sqrt[3]{g} \cdot \sqrt[3]{g}} \cdot \left(\sqrt[3]{\frac{1}{2}} \cdot \frac{\sqrt[3]{\sqrt[3]{g}}}{\sqrt[3]{a}}\right)double f(double g, double a) {
double r109338 = g;
double r109339 = 2.0;
double r109340 = a;
double r109341 = r109339 * r109340;
double r109342 = r109338 / r109341;
double r109343 = cbrt(r109342);
return r109343;
}
double f(double g, double a) {
double r109344 = g;
double r109345 = cbrt(r109344);
double r109346 = r109345 * r109345;
double r109347 = cbrt(r109346);
double r109348 = 1.0;
double r109349 = 2.0;
double r109350 = r109348 / r109349;
double r109351 = cbrt(r109350);
double r109352 = cbrt(r109345);
double r109353 = a;
double r109354 = cbrt(r109353);
double r109355 = r109352 / r109354;
double r109356 = r109351 * r109355;
double r109357 = r109347 * r109356;
return r109357;
}



Bits error versus g



Bits error versus a
Results
Initial program 15.5
rmApplied add-cube-cbrt15.7
Applied times-frac15.6
Applied cbrt-prod5.7
rmApplied cbrt-div1.2
rmApplied div-inv1.2
Applied cbrt-prod1.2
rmApplied *-un-lft-identity1.2
Final simplification1.2
herbie shell --seed 2019304
(FPCore (g a)
:name "2-ancestry mixing, zero discriminant"
:precision binary64
(cbrt (/ g (* 2 a))))