\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 r135746 = g;
double r135747 = 2.0;
double r135748 = a;
double r135749 = r135747 * r135748;
double r135750 = r135746 / r135749;
double r135751 = cbrt(r135750);
return r135751;
}
double f(double g, double a) {
double r135752 = g;
double r135753 = cbrt(r135752);
double r135754 = r135753 * r135753;
double r135755 = 2.0;
double r135756 = r135754 / r135755;
double r135757 = cbrt(r135756);
double r135758 = sqrt(r135757);
double r135759 = r135758 * r135758;
double r135760 = cbrt(r135753);
double r135761 = a;
double r135762 = cbrt(r135761);
double r135763 = r135760 / r135762;
double r135764 = r135759 * r135763;
return r135764;
}



Bits error versus g



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