\sqrt[3]{\frac{g}{2 \cdot a}}\left(\sqrt{\sqrt[3]{\frac{\sqrt[3]{g} \cdot \sqrt[3]{g}}{2}}} \cdot \sqrt{\sqrt[3]{\frac{\sqrt[3]{g} \cdot \sqrt[3]{g}}{2}}}\right) \cdot \frac{\sqrt[3]{\sqrt[3]{g}}}{\sqrt[3]{a}}double f(double g, double a) {
double r172968 = g;
double r172969 = 2.0;
double r172970 = a;
double r172971 = r172969 * r172970;
double r172972 = r172968 / r172971;
double r172973 = cbrt(r172972);
return r172973;
}
double f(double g, double a) {
double r172974 = g;
double r172975 = cbrt(r172974);
double r172976 = r172975 * r172975;
double r172977 = 2.0;
double r172978 = r172976 / r172977;
double r172979 = cbrt(r172978);
double r172980 = sqrt(r172979);
double r172981 = r172980 * r172980;
double r172982 = cbrt(r172975);
double r172983 = a;
double r172984 = cbrt(r172983);
double r172985 = r172982 / r172984;
double r172986 = r172981 * r172985;
return r172986;
}



Bits error versus g



Bits error versus a
Results
Initial program 15.4
rmApplied add-cube-cbrt15.6
Applied times-frac15.6
Applied cbrt-prod5.6
rmApplied cbrt-div1.2
rmApplied add-sqr-sqrt1.2
Final simplification1.2
herbie shell --seed 2020025
(FPCore (g a)
:name "2-ancestry mixing, zero discriminant"
:precision binary64
(cbrt (/ g (* 2 a))))