\frac{x}{y} \cdot \left(z - t\right) + t\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 \left(z - t\right)\right) + tdouble f(double x, double y, double z, double t) {
double r316106 = x;
double r316107 = y;
double r316108 = r316106 / r316107;
double r316109 = z;
double r316110 = t;
double r316111 = r316109 - r316110;
double r316112 = r316108 * r316111;
double r316113 = r316112 + r316110;
return r316113;
}
double f(double x, double y, double z, double t) {
double r316114 = x;
double r316115 = cbrt(r316114);
double r316116 = r316115 * r316115;
double r316117 = y;
double r316118 = cbrt(r316117);
double r316119 = r316118 * r316118;
double r316120 = r316116 / r316119;
double r316121 = r316115 / r316118;
double r316122 = z;
double r316123 = t;
double r316124 = r316122 - r316123;
double r316125 = r316121 * r316124;
double r316126 = r316120 * r316125;
double r316127 = r316126 + r316123;
return r316127;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 2.2 |
|---|---|
| Target | 2.4 |
| Herbie | 1.0 |
Initial program 2.2
rmApplied add-cube-cbrt2.7
Applied add-cube-cbrt2.9
Applied times-frac2.9
Applied associate-*l*1.0
Final simplification1.0
herbie shell --seed 2019323 +o rules:numerics
(FPCore (x y z t)
:name "Numeric.Signal.Multichannel:$cget from hsignal-0.2.7.1"
:precision binary64
:herbie-target
(if (< z 2.759456554562692e-282) (+ (* (/ x y) (- z t)) t) (if (< z 2.326994450874436e-110) (+ (* x (/ (- z t) y)) t) (+ (* (/ x y) (- z t)) t)))
(+ (* (/ x y) (- z t)) t))