\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 r117799 = g;
double r117800 = 2.0;
double r117801 = a;
double r117802 = r117800 * r117801;
double r117803 = r117799 / r117802;
double r117804 = cbrt(r117803);
return r117804;
}
double f(double g, double a) {
double r117805 = g;
double r117806 = cbrt(r117805);
double r117807 = r117806 * r117806;
double r117808 = cbrt(r117807);
double r117809 = 1.0;
double r117810 = 2.0;
double r117811 = r117809 / r117810;
double r117812 = cbrt(r117811);
double r117813 = r117808 * r117812;
double r117814 = cbrt(r117806);
double r117815 = r117808 * r117814;
double r117816 = cbrt(r117815);
double r117817 = a;
double r117818 = cbrt(r117817);
double r117819 = r117816 / r117818;
double r117820 = r117813 * r117819;
return r117820;
}



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
(FPCore (g a)
:name "2-ancestry mixing, zero discriminant"
:precision binary64
(cbrt (/ g (* 2 a))))