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 (((x <= -1.3397329542554413e+35) || !(x <= 2.9490020604149927e-58))) {
VAR = ((double) fabs(((double) (((double) (((double) (x + 4.0)) / y)) - ((double) (x * ((double) (z / y))))))));
} else {
VAR = ((double) fabs(((double) (((double) (x + ((double) (4.0 - ((double) (x * z)))))) / y))));
}
return VAR;
}



Bits error versus x



Bits error versus y



Bits error versus z
Results
if x < -1.3397329542554413e35 or 2.9490020604149927e-58 < x Initial program 0.2
Simplified0.3
if -1.3397329542554413e35 < x < 2.9490020604149927e-58Initial program 2.5
rmApplied associate-*l/0.1
Applied sub-div0.1
Simplified0.1
Final simplification0.2
herbie shell --seed 2020181
(FPCore (x y z)
:name "fabs fraction 1"
:precision binary64
(fabs (- (/ (+ x 4.0) y) (* (/ x y) z))))