\sqrt{\left|\frac{a \cdot a - b \cdot b}{a \cdot a}\right|}\sqrt{\left|1 - \sqrt[3]{{\left(\frac{b}{a}\right)}^{6}}\right|}double f(double a, double b) {
double r36395 = a;
double r36396 = r36395 * r36395;
double r36397 = b;
double r36398 = r36397 * r36397;
double r36399 = r36396 - r36398;
double r36400 = r36399 / r36396;
double r36401 = fabs(r36400);
double r36402 = sqrt(r36401);
return r36402;
}
double f(double a, double b) {
double r36403 = 1.0;
double r36404 = b;
double r36405 = a;
double r36406 = r36404 / r36405;
double r36407 = 6.0;
double r36408 = pow(r36406, r36407);
double r36409 = cbrt(r36408);
double r36410 = r36403 - r36409;
double r36411 = fabs(r36410);
double r36412 = sqrt(r36411);
return r36412;
}



Bits error versus a



Bits error versus b
Results
Initial program 14.7
Simplified14.7
rmApplied add-cbrt-cube26.9
Applied add-cbrt-cube26.9
Applied cbrt-unprod43.8
Applied add-cbrt-cube43.8
Applied add-cbrt-cube43.8
Applied cbrt-unprod43.8
Applied cbrt-undiv43.8
Simplified0.0
Final simplification0.0
herbie shell --seed 2019325
(FPCore (a b)
:name "Eccentricity of an ellipse"
:precision binary64
:pre (<= 0.0 b a 1)
(sqrt (fabs (/ (- (* a a) (* b b)) (* a a)))))