\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 r333020 = x;
double r333021 = y;
double r333022 = r333020 / r333021;
double r333023 = z;
double r333024 = t;
double r333025 = r333023 - r333024;
double r333026 = r333022 * r333025;
double r333027 = r333026 + r333024;
return r333027;
}
double f(double x, double y, double z, double t) {
double r333028 = x;
double r333029 = cbrt(r333028);
double r333030 = r333029 * r333029;
double r333031 = y;
double r333032 = cbrt(r333031);
double r333033 = r333032 * r333032;
double r333034 = r333030 / r333033;
double r333035 = r333029 / r333032;
double r333036 = z;
double r333037 = t;
double r333038 = r333036 - r333037;
double r333039 = r333035 * r333038;
double r333040 = r333034 * r333039;
double r333041 = r333040 + r333037;
return r333041;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 2.1 |
|---|---|
| Target | 2.4 |
| Herbie | 0.9 |
Initial program 2.1
rmApplied add-cube-cbrt2.6
Applied add-cube-cbrt2.8
Applied times-frac2.8
Applied associate-*l*0.9
Final simplification0.9
herbie shell --seed 2019208 +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.7594565545626922e-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))