\sqrt{\left|\frac{a \cdot a - b \cdot b}{a \cdot a}\right|}\sqrt{\left|\frac{1}{\frac{a}{\left(a + b\right) \cdot \frac{a - b}{a}}}\right|}double f(double a, double b) {
double r88690 = a;
double r88691 = r88690 * r88690;
double r88692 = b;
double r88693 = r88692 * r88692;
double r88694 = r88691 - r88693;
double r88695 = r88694 / r88691;
double r88696 = fabs(r88695);
double r88697 = sqrt(r88696);
return r88697;
}
double f(double a, double b) {
double r88698 = 1.0;
double r88699 = a;
double r88700 = b;
double r88701 = r88699 + r88700;
double r88702 = r88699 - r88700;
double r88703 = r88702 / r88699;
double r88704 = r88701 * r88703;
double r88705 = r88699 / r88704;
double r88706 = r88698 / r88705;
double r88707 = fabs(r88706);
double r88708 = sqrt(r88707);
return r88708;
}



Bits error versus a



Bits error versus b
Results
Initial program 14.6
rmApplied difference-of-squares14.6
Applied times-frac0.0
rmApplied associate-*r/0.0
Simplified0.0
rmApplied clear-num0.0
Final simplification0.0
herbie shell --seed 2020036 +o rules:numerics
(FPCore (a b)
:name "Eccentricity of an ellipse"
:precision binary64
:pre (<= 0.0 b a 1)
(sqrt (fabs (/ (- (* a a) (* b b)) (* a a)))))