\frac{x}{y} \cdot \left(z - t\right) + t\left(\frac{\sqrt[3]{x}}{\sqrt[3]{y}} \cdot \left(z - t\right)\right) \cdot \frac{\sqrt[3]{x} \cdot \sqrt[3]{x}}{\sqrt[3]{y} \cdot \sqrt[3]{y}} + tdouble f(double x, double y, double z, double t) {
double r20281249 = x;
double r20281250 = y;
double r20281251 = r20281249 / r20281250;
double r20281252 = z;
double r20281253 = t;
double r20281254 = r20281252 - r20281253;
double r20281255 = r20281251 * r20281254;
double r20281256 = r20281255 + r20281253;
return r20281256;
}
double f(double x, double y, double z, double t) {
double r20281257 = x;
double r20281258 = cbrt(r20281257);
double r20281259 = y;
double r20281260 = cbrt(r20281259);
double r20281261 = r20281258 / r20281260;
double r20281262 = z;
double r20281263 = t;
double r20281264 = r20281262 - r20281263;
double r20281265 = r20281261 * r20281264;
double r20281266 = r20281258 * r20281258;
double r20281267 = r20281260 * r20281260;
double r20281268 = r20281266 / r20281267;
double r20281269 = r20281265 * r20281268;
double r20281270 = r20281269 + r20281263;
return r20281270;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 2.1 |
|---|---|
| Target | 2.2 |
| Herbie | 0.9 |
Initial program 2.1
rmApplied add-cube-cbrt2.6
Applied add-cube-cbrt2.7
Applied times-frac2.7
Applied associate-*l*0.9
Final simplification0.9
herbie shell --seed 2019163
(FPCore (x y z t)
:name "Numeric.Signal.Multichannel:$cget from hsignal-0.2.7.1"
: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))