\sqrt[3]{\frac{g}{2 \cdot a}}\sqrt[3]{\frac{\sqrt[3]{g} \cdot \sqrt[3]{g}}{2}} \cdot \frac{\sqrt[3]{\sqrt[3]{g}}}{\sqrt[3]{a}}double f(double g, double a) {
double r160941 = g;
double r160942 = 2.0;
double r160943 = a;
double r160944 = r160942 * r160943;
double r160945 = r160941 / r160944;
double r160946 = cbrt(r160945);
return r160946;
}
double f(double g, double a) {
double r160947 = g;
double r160948 = cbrt(r160947);
double r160949 = r160948 * r160948;
double r160950 = 2.0;
double r160951 = r160949 / r160950;
double r160952 = cbrt(r160951);
double r160953 = cbrt(r160948);
double r160954 = a;
double r160955 = cbrt(r160954);
double r160956 = r160953 / r160955;
double r160957 = r160952 * r160956;
return r160957;
}



Bits error versus g



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