\left|\frac{x + 4}{y} - \frac{x}{y} \cdot z\right|\left|\frac{x + 4}{y} - \frac{x}{\sqrt[3]{y} \cdot \sqrt[3]{y}} \cdot \frac{z}{\sqrt[3]{y}}\right|double f(double x, double y, double z) {
double r28774 = x;
double r28775 = 4.0;
double r28776 = r28774 + r28775;
double r28777 = y;
double r28778 = r28776 / r28777;
double r28779 = r28774 / r28777;
double r28780 = z;
double r28781 = r28779 * r28780;
double r28782 = r28778 - r28781;
double r28783 = fabs(r28782);
return r28783;
}
double f(double x, double y, double z) {
double r28784 = x;
double r28785 = 4.0;
double r28786 = r28784 + r28785;
double r28787 = y;
double r28788 = r28786 / r28787;
double r28789 = cbrt(r28787);
double r28790 = r28789 * r28789;
double r28791 = r28784 / r28790;
double r28792 = z;
double r28793 = r28792 / r28789;
double r28794 = r28791 * r28793;
double r28795 = r28788 - r28794;
double r28796 = fabs(r28795);
return r28796;
}



Bits error versus x



Bits error versus y



Bits error versus z
Results
Initial program 1.6
rmApplied div-inv1.6
Applied associate-*l*3.4
Simplified3.4
rmApplied add-cube-cbrt3.7
Applied *-un-lft-identity3.7
Applied times-frac3.7
Applied associate-*r*2.0
Simplified2.0
Final simplification2.0
herbie shell --seed 2020034 +o rules:numerics
(FPCore (x y z)
:name "fabs fraction 1"
:precision binary64
(fabs (- (/ (+ x 4) y) (* (/ x y) z))))