\left|\frac{x + 4}{y} - \frac{x}{y} \cdot z\right|\left|\mathsf{fma}\left(4, \frac{1}{y}, \frac{x}{y}\right) - \frac{\sqrt[3]{x} \cdot \sqrt[3]{x}}{\sqrt[3]{y} \cdot \sqrt[3]{y}} \cdot \left(\frac{\sqrt[3]{x}}{\sqrt[3]{y}} \cdot z\right)\right|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((fma(4.0, (1.0 / y), (x / y)) - (((cbrt(x) * cbrt(x)) / (cbrt(y) * cbrt(y))) * ((cbrt(x) / cbrt(y)) * z))));
}



Bits error versus x



Bits error versus y



Bits error versus z
Results
Initial program 1.5
Taylor expanded around 0 1.6
Simplified1.6
rmApplied add-cube-cbrt1.8
Applied add-cube-cbrt1.9
Applied times-frac1.9
Applied associate-*l*0.6
Final simplification0.6
herbie shell --seed 2020057 +o rules:numerics
(FPCore (x y z)
:name "fabs fraction 1"
:precision binary64
(fabs (- (/ (+ x 4) y) (* (/ x y) z))))