\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 r383077 = x;
double r383078 = y;
double r383079 = r383077 / r383078;
double r383080 = z;
double r383081 = t;
double r383082 = r383080 - r383081;
double r383083 = r383079 * r383082;
double r383084 = r383083 + r383081;
return r383084;
}
double f(double x, double y, double z, double t) {
double r383085 = x;
double r383086 = cbrt(r383085);
double r383087 = r383086 * r383086;
double r383088 = y;
double r383089 = cbrt(r383088);
double r383090 = r383089 * r383089;
double r383091 = r383087 / r383090;
double r383092 = r383086 / r383089;
double r383093 = z;
double r383094 = t;
double r383095 = r383093 - r383094;
double r383096 = r383092 * r383095;
double r383097 = r383091 * r383096;
double r383098 = r383097 + r383094;
return r383098;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 2.1 |
|---|---|
| Target | 2.5 |
| 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 2019354 +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))