\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 r519008 = x;
double r519009 = y;
double r519010 = r519008 / r519009;
double r519011 = z;
double r519012 = t;
double r519013 = r519011 - r519012;
double r519014 = r519010 * r519013;
double r519015 = r519014 + r519012;
return r519015;
}
double f(double x, double y, double z, double t) {
double r519016 = x;
double r519017 = cbrt(r519016);
double r519018 = r519017 * r519017;
double r519019 = y;
double r519020 = cbrt(r519019);
double r519021 = r519020 * r519020;
double r519022 = r519018 / r519021;
double r519023 = r519017 / r519020;
double r519024 = z;
double r519025 = t;
double r519026 = r519024 - r519025;
double r519027 = r519023 * r519026;
double r519028 = r519022 * r519027;
double r519029 = r519028 + r519025;
return r519029;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 1.9 |
|---|---|
| Target | 2.1 |
| Herbie | 1.0 |
Initial program 1.9
rmApplied add-cube-cbrt2.4
Applied add-cube-cbrt2.5
Applied times-frac2.5
Applied associate-*l*1.0
Final simplification1.0
herbie shell --seed 2020059 +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))