\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 r152225 = g;
double r152226 = 2.0;
double r152227 = a;
double r152228 = r152226 * r152227;
double r152229 = r152225 / r152228;
double r152230 = cbrt(r152229);
return r152230;
}
double f(double g, double a) {
double r152231 = g;
double r152232 = cbrt(r152231);
double r152233 = r152232 * r152232;
double r152234 = 2.0;
double r152235 = r152233 / r152234;
double r152236 = cbrt(r152235);
double r152237 = sqrt(r152236);
double r152238 = r152237 * r152237;
double r152239 = cbrt(r152232);
double r152240 = a;
double r152241 = cbrt(r152240);
double r152242 = r152239 / r152241;
double r152243 = r152238 * r152242;
return r152243;
}



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