\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]{{\left(\sqrt[3]{x}\right)}^{3}}}{\sqrt[3]{y}} \cdot \left(z - t\right)\right) + tdouble f(double x, double y, double z, double t) {
double r630352 = x;
double r630353 = y;
double r630354 = r630352 / r630353;
double r630355 = z;
double r630356 = t;
double r630357 = r630355 - r630356;
double r630358 = r630354 * r630357;
double r630359 = r630358 + r630356;
return r630359;
}
double f(double x, double y, double z, double t) {
double r630360 = x;
double r630361 = cbrt(r630360);
double r630362 = r630361 * r630361;
double r630363 = y;
double r630364 = cbrt(r630363);
double r630365 = r630364 * r630364;
double r630366 = r630362 / r630365;
double r630367 = 3.0;
double r630368 = pow(r630361, r630367);
double r630369 = cbrt(r630368);
double r630370 = r630369 / r630364;
double r630371 = z;
double r630372 = t;
double r630373 = r630371 - r630372;
double r630374 = r630370 * r630373;
double r630375 = r630366 * r630374;
double r630376 = r630375 + r630372;
return r630376;
}




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.5
Applied add-cube-cbrt2.6
Applied times-frac2.6
Applied associate-*l*0.9
rmApplied add-cbrt-cube0.9
Simplified0.9
Final simplification0.9
herbie shell --seed 2020043 +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))