\left|\frac{x + 4}{y} - \frac{x}{y} \cdot z\right|
\left|\left(\frac{4}{y} + \frac{x}{y}\right) - \frac{x}{y} \cdot z\right|
(FPCore (x y z) :precision binary64 (fabs (- (/ (+ x 4.0) y) (* (/ x y) z))))
(FPCore (x y z) :precision binary64 (fabs (- (+ (/ 4.0 y) (/ x y)) (* (/ x y) z))))
double code(double x, double y, double z) {
return fabs(((x + 4.0) / y) - ((x / y) * z));
}
double code(double x, double y, double z) {
return fabs(((4.0 / y) + (x / y)) - ((x / y) * z));
}



Bits error versus x



Bits error versus y



Bits error versus z
Results
Initial program 1.5
Taylor expanded in x around 0 1.5
Taylor expanded in y around 0 1.5
Simplified1.5
Final simplification1.5
herbie shell --seed 2022061
(FPCore (x y z)
:name "fabs fraction 1"
:precision binary64
(fabs (- (/ (+ x 4.0) y) (* (/ x y) z))))