\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]{\sqrt[3]{y} \cdot \sqrt[3]{y}} \cdot \sqrt[3]{\sqrt[3]{y}}} \cdot \left(z - t\right)\right) + tdouble f(double x, double y, double z, double t) {
double r501694 = x;
double r501695 = y;
double r501696 = r501694 / r501695;
double r501697 = z;
double r501698 = t;
double r501699 = r501697 - r501698;
double r501700 = r501696 * r501699;
double r501701 = r501700 + r501698;
return r501701;
}
double f(double x, double y, double z, double t) {
double r501702 = x;
double r501703 = cbrt(r501702);
double r501704 = r501703 * r501703;
double r501705 = y;
double r501706 = cbrt(r501705);
double r501707 = r501706 * r501706;
double r501708 = r501704 / r501707;
double r501709 = cbrt(r501707);
double r501710 = cbrt(r501706);
double r501711 = r501709 * r501710;
double r501712 = r501703 / r501711;
double r501713 = z;
double r501714 = t;
double r501715 = r501713 - r501714;
double r501716 = r501712 * r501715;
double r501717 = r501708 * r501716;
double r501718 = r501717 + r501714;
return r501718;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 1.8 |
|---|---|
| Target | 2.3 |
| Herbie | 0.9 |
Initial program 1.8
rmApplied add-cube-cbrt2.4
Applied add-cube-cbrt2.5
Applied times-frac2.5
Applied associate-*l*0.9
rmApplied add-cube-cbrt0.9
Applied cbrt-prod0.9
Final simplification0.9
herbie shell --seed 2020025
(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))