\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 r365784 = x;
double r365785 = y;
double r365786 = r365784 / r365785;
double r365787 = z;
double r365788 = t;
double r365789 = r365787 - r365788;
double r365790 = r365786 * r365789;
double r365791 = r365790 + r365788;
return r365791;
}
double f(double x, double y, double z, double t) {
double r365792 = x;
double r365793 = cbrt(r365792);
double r365794 = r365793 * r365793;
double r365795 = y;
double r365796 = cbrt(r365795);
double r365797 = r365796 * r365796;
double r365798 = r365794 / r365797;
double r365799 = r365793 / r365796;
double r365800 = z;
double r365801 = t;
double r365802 = r365800 - r365801;
double r365803 = r365799 * r365802;
double r365804 = r365798 * r365803;
double r365805 = r365804 + r365801;
return r365805;
}




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))