\sqrt[3]{\frac{g}{2 \cdot a}}\sqrt[3]{\frac{\sqrt[3]{\frac{1}{2}}}{a}} \cdot \left(\sqrt[3]{g} \cdot \sqrt[3]{\sqrt[3]{\frac{1}{2}} \cdot \sqrt[3]{\frac{1}{2}}}\right)double f(double g, double a) {
double r14655943 = g;
double r14655944 = 2.0;
double r14655945 = a;
double r14655946 = r14655944 * r14655945;
double r14655947 = r14655943 / r14655946;
double r14655948 = cbrt(r14655947);
return r14655948;
}
double f(double g, double a) {
double r14655949 = 0.5;
double r14655950 = cbrt(r14655949);
double r14655951 = a;
double r14655952 = r14655950 / r14655951;
double r14655953 = cbrt(r14655952);
double r14655954 = g;
double r14655955 = cbrt(r14655954);
double r14655956 = r14655950 * r14655950;
double r14655957 = cbrt(r14655956);
double r14655958 = r14655955 * r14655957;
double r14655959 = r14655953 * r14655958;
return r14655959;
}



Bits error versus g



Bits error versus a
Results
Initial program 14.9
rmApplied cbrt-div0.9
rmApplied div-inv0.9
rmApplied add-cbrt-cube0.9
Applied cbrt-undiv0.9
Simplified0.8
rmApplied *-un-lft-identity0.8
Applied add-cube-cbrt0.8
Applied times-frac0.8
Applied cbrt-prod0.9
Applied associate-*r*0.9
Final simplification0.9
herbie shell --seed 2019112 +o rules:numerics
(FPCore (g a)
:name "2-ancestry mixing, zero discriminant"
(cbrt (/ g (* 2 a))))