\sqrt{\left|\frac{a \cdot a - b \cdot b}{a \cdot a}\right|}\sqrt{\left|\frac{\left(a - b\right) + \frac{a - b}{a} \cdot b}{a}\right|}double f(double a, double b) {
double r76928 = a;
double r76929 = r76928 * r76928;
double r76930 = b;
double r76931 = r76930 * r76930;
double r76932 = r76929 - r76931;
double r76933 = r76932 / r76929;
double r76934 = fabs(r76933);
double r76935 = sqrt(r76934);
return r76935;
}
double f(double a, double b) {
double r76936 = a;
double r76937 = b;
double r76938 = r76936 - r76937;
double r76939 = r76938 / r76936;
double r76940 = r76939 * r76937;
double r76941 = r76938 + r76940;
double r76942 = r76941 / r76936;
double r76943 = fabs(r76942);
double r76944 = sqrt(r76943);
return r76944;
}



Bits error versus a



Bits error versus b
Results
Initial program 14.8
rmApplied difference-of-squares14.8
Applied times-frac0.0
rmApplied associate-*r/0.0
Simplified0.0
rmApplied *-commutative0.0
rmApplied distribute-lft-in0.0
Simplified0.0
Final simplification0.0
herbie shell --seed 2019362
(FPCore (a b)
:name "Eccentricity of an ellipse"
:precision binary64
:pre (<= 0.0 b a 1)
(sqrt (fabs (/ (- (* a a) (* b b)) (* a a)))))