\frac{x}{y} \cdot \left(z - t\right) + t\left(x \cdot \frac{\sqrt[3]{z - t} \cdot \sqrt[3]{z - t}}{\sqrt[3]{y} \cdot \sqrt[3]{y}}\right) \cdot \frac{\sqrt[3]{z - t}}{\sqrt[3]{y}} + tdouble f(double x, double y, double z, double t) {
double r422518 = x;
double r422519 = y;
double r422520 = r422518 / r422519;
double r422521 = z;
double r422522 = t;
double r422523 = r422521 - r422522;
double r422524 = r422520 * r422523;
double r422525 = r422524 + r422522;
return r422525;
}
double f(double x, double y, double z, double t) {
double r422526 = x;
double r422527 = z;
double r422528 = t;
double r422529 = r422527 - r422528;
double r422530 = cbrt(r422529);
double r422531 = r422530 * r422530;
double r422532 = y;
double r422533 = cbrt(r422532);
double r422534 = r422533 * r422533;
double r422535 = r422531 / r422534;
double r422536 = r422526 * r422535;
double r422537 = r422530 / r422533;
double r422538 = r422536 * r422537;
double r422539 = r422538 + r422528;
return r422539;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 2.1 |
|---|---|
| Target | 2.3 |
| Herbie | 1.6 |
Initial program 2.1
rmApplied div-inv2.1
Applied associate-*l*6.6
Simplified6.5
rmApplied add-cube-cbrt7.0
Applied add-cube-cbrt7.1
Applied times-frac7.1
Applied associate-*r*1.6
Final simplification1.6
herbie shell --seed 2020024
(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))