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.3266761392413531e-136)) {
VAR = ((double) fabs(((double) (((double) (((double) (x + 4.0)) / y)) - ((double) (x * ((double) (z / y))))))));
} else {
double VAR_1;
if ((x <= 1.7904176914311191e-62)) {
VAR_1 = ((double) fabs(((double) (((double) (x + ((double) (4.0 - ((double) (x * z)))))) / y))));
} else {
VAR_1 = ((double) fabs(((double) (((double) (4.0 / y)) + ((double) (((double) (x / y)) * ((double) (1.0 - z))))))));
}
VAR = VAR_1;
}
return VAR;
}



Bits error versus x



Bits error versus y



Bits error versus z
Results
if x < -1.3266761392413531e-136Initial program 0.7
Simplified1.5
if -1.3266761392413531e-136 < x < 1.7904176914311191e-62Initial program 2.6
rmApplied associate-*l/0.1
Applied sub-div0.0
Simplified0.0
if 1.7904176914311191e-62 < x Initial program 0.6
Taylor expanded around 0 7.4
Simplified0.6
Final simplification0.7
herbie shell --seed 2020179
(FPCore (x y z)
:name "fabs fraction 1"
:precision binary64
(fabs (- (/ (+ x 4.0) y) (* (/ x y) z))))