double code(double x, double y, double z) {
return ((double) fabs(((double) (((double) (((double) (x + 4.0)) / y)) - ((double) (((double) (x / y)) * z))))));
}
double code(double x, double y, double z) {
double VAR;
if ((y <= 7.06822395138139e+98)) {
VAR = ((double) fabs(((double) (((double) (4.0 / y)) + ((double) (((double) (x / y)) * ((double) (1.0 - z))))))));
} else {
VAR = ((double) fabs(((double) (((double) (((double) (4.0 + x)) / y)) - ((double) (x * ((double) (z / y))))))));
}
return VAR;
}



Bits error versus x



Bits error versus y



Bits error versus z
Results
if y < 7.0682239513813901e98Initial program 0.8
Taylor expanded around 0 2.7
Simplified0.8
if 7.0682239513813901e98 < y Initial program 2.8
Simplified0.1
Final simplification0.7
herbie shell --seed 2020180
(FPCore (x y z)
:name "fabs fraction 1"
:precision binary64
(fabs (- (/ (+ x 4.0) y) (* (/ x y) z))))