double code(double g, double a) {
return cbrt(g / (2.0 * a));
}
double code(double g, double a) {
return cbrt(g) / cbrt(2.0 * a);
}



Bits error versus g



Bits error versus a
Results
Initial program 15.6
rmApplied cbrt-div0.9
Final simplification0.9
herbie shell --seed 2020322
(FPCore (g a)
:name "2-ancestry mixing, zero discriminant"
:precision binary64
(cbrt (/ g (* 2.0 a))))