\sqrt{\left|\frac{a \cdot a - b \cdot b}{a \cdot a}\right|}\log \left(e^{\sqrt{\left|\left(\frac{a + b}{a} \cdot \frac{\sqrt{a - b}}{\sqrt{a}}\right) \cdot \frac{\sqrt{a - b}}{\sqrt{a}}\right|}}\right)double f(double a, double b) {
double r70713 = a;
double r70714 = r70713 * r70713;
double r70715 = b;
double r70716 = r70715 * r70715;
double r70717 = r70714 - r70716;
double r70718 = r70717 / r70714;
double r70719 = fabs(r70718);
double r70720 = sqrt(r70719);
return r70720;
}
double f(double a, double b) {
double r70721 = a;
double r70722 = b;
double r70723 = r70721 + r70722;
double r70724 = r70723 / r70721;
double r70725 = r70721 - r70722;
double r70726 = sqrt(r70725);
double r70727 = sqrt(r70721);
double r70728 = r70726 / r70727;
double r70729 = r70724 * r70728;
double r70730 = r70729 * r70728;
double r70731 = fabs(r70730);
double r70732 = sqrt(r70731);
double r70733 = exp(r70732);
double r70734 = log(r70733);
return r70734;
}



Bits error versus a



Bits error versus b
Results
Initial program 13.7
rmApplied difference-of-squares13.7
Applied times-frac0.0
rmApplied add-log-exp0.0
rmApplied add-sqr-sqrt0.0
Applied add-sqr-sqrt0.0
Applied times-frac0.0
Applied associate-*r*0.0
Final simplification0.0
herbie shell --seed 2020060
(FPCore (a b)
:name "Eccentricity of an ellipse"
:precision binary64
:pre (<= 0.0 b a 1)
(sqrt (fabs (/ (- (* a a) (* b b)) (* a a)))))