\frac{x}{y} \cdot \left(z - t\right) + t\sqrt{\frac{\sqrt[3]{x} \cdot \sqrt[3]{x}}{\sqrt[3]{y} \cdot \sqrt[3]{y}}} \cdot \left(\left(\left|\frac{\sqrt[3]{x}}{\sqrt[3]{y}}\right| \cdot \left(z - t\right)\right) \cdot \frac{\sqrt[3]{x}}{\sqrt[3]{y}}\right) + tdouble f(double x, double y, double z, double t) {
double r304601 = x;
double r304602 = y;
double r304603 = r304601 / r304602;
double r304604 = z;
double r304605 = t;
double r304606 = r304604 - r304605;
double r304607 = r304603 * r304606;
double r304608 = r304607 + r304605;
return r304608;
}
double f(double x, double y, double z, double t) {
double r304609 = x;
double r304610 = cbrt(r304609);
double r304611 = r304610 * r304610;
double r304612 = y;
double r304613 = cbrt(r304612);
double r304614 = r304613 * r304613;
double r304615 = r304611 / r304614;
double r304616 = sqrt(r304615);
double r304617 = r304610 / r304613;
double r304618 = fabs(r304617);
double r304619 = z;
double r304620 = t;
double r304621 = r304619 - r304620;
double r304622 = r304618 * r304621;
double r304623 = r304622 * r304617;
double r304624 = r304616 * r304623;
double r304625 = r304624 + r304620;
return r304625;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 2.2 |
|---|---|
| Target | 2.5 |
| Herbie | 1.0 |
Initial program 2.2
rmApplied add-cube-cbrt2.7
Applied add-cube-cbrt2.8
Applied times-frac2.8
Applied associate-*l*1.0
rmApplied add-sqr-sqrt1.0
Applied associate-*l*1.0
Simplified1.0
Final simplification1.0
herbie shell --seed 2020046
(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))