\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 r69428 = a;
double r69429 = r69428 * r69428;
double r69430 = b;
double r69431 = r69430 * r69430;
double r69432 = r69429 - r69431;
double r69433 = r69432 / r69429;
double r69434 = fabs(r69433);
double r69435 = sqrt(r69434);
return r69435;
}
double f(double a, double b) {
double r69436 = 1.0;
double r69437 = a;
double r69438 = b;
double r69439 = r69437 + r69438;
double r69440 = r69437 - r69438;
double r69441 = r69440 / r69437;
double r69442 = r69439 * r69441;
double r69443 = r69437 / r69442;
double r69444 = r69436 / r69443;
double r69445 = fabs(r69444);
double r69446 = sqrt(r69445);
return r69446;
}



Bits error versus a



Bits error versus b
Results
Initial program 14.5
rmApplied difference-of-squares14.5
Applied times-frac0.0
rmApplied associate-*r/0.0
Simplified0.0
rmApplied clear-num0.0
Final simplification0.0
herbie shell --seed 2019354 +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)))))