\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 r398658 = x;
double r398659 = y;
double r398660 = r398658 / r398659;
double r398661 = z;
double r398662 = t;
double r398663 = r398661 - r398662;
double r398664 = r398660 * r398663;
double r398665 = r398664 + r398662;
return r398665;
}
double f(double x, double y, double z, double t) {
double r398666 = x;
double r398667 = cbrt(r398666);
double r398668 = r398667 * r398667;
double r398669 = y;
double r398670 = cbrt(r398669);
double r398671 = r398670 * r398670;
double r398672 = r398668 / r398671;
double r398673 = r398667 / r398670;
double r398674 = z;
double r398675 = t;
double r398676 = r398674 - r398675;
double r398677 = r398673 * r398676;
double r398678 = r398672 * r398677;
double r398679 = r398678 + r398675;
return r398679;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 2.2 |
|---|---|
| Target | 2.2 |
| 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
Final simplification1.0
herbie shell --seed 2020002 +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))