\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 r264664 = x;
double r264665 = y;
double r264666 = r264664 / r264665;
double r264667 = z;
double r264668 = t;
double r264669 = r264667 - r264668;
double r264670 = r264666 * r264669;
double r264671 = r264670 + r264668;
return r264671;
}
double f(double x, double y, double z, double t) {
double r264672 = x;
double r264673 = cbrt(r264672);
double r264674 = r264673 * r264673;
double r264675 = y;
double r264676 = cbrt(r264675);
double r264677 = r264676 * r264676;
double r264678 = r264674 / r264677;
double r264679 = r264673 / r264676;
double r264680 = z;
double r264681 = t;
double r264682 = r264680 - r264681;
double r264683 = r264679 * r264682;
double r264684 = r264678 * r264683;
double r264685 = r264684 + r264681;
return r264685;
}




Bits error versus x




Bits error versus y




Bits error versus z




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