\frac{x - y}{z - y} \cdot t\frac{\sqrt[3]{x - y} \cdot \sqrt[3]{x - y}}{\sqrt[3]{z - y} \cdot \sqrt[3]{z - y}} \cdot \left(\frac{\sqrt[3]{x - y}}{\sqrt[3]{z - y}} \cdot t\right)double f(double x, double y, double z, double t) {
double r450348 = x;
double r450349 = y;
double r450350 = r450348 - r450349;
double r450351 = z;
double r450352 = r450351 - r450349;
double r450353 = r450350 / r450352;
double r450354 = t;
double r450355 = r450353 * r450354;
return r450355;
}
double f(double x, double y, double z, double t) {
double r450356 = x;
double r450357 = y;
double r450358 = r450356 - r450357;
double r450359 = cbrt(r450358);
double r450360 = r450359 * r450359;
double r450361 = z;
double r450362 = r450361 - r450357;
double r450363 = cbrt(r450362);
double r450364 = r450363 * r450363;
double r450365 = r450360 / r450364;
double r450366 = r450359 / r450363;
double r450367 = t;
double r450368 = r450366 * r450367;
double r450369 = r450365 * r450368;
return r450369;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 2.4 |
|---|---|
| Target | 2.3 |
| Herbie | 1.1 |
Initial program 2.4
rmApplied add-cube-cbrt3.4
Applied add-cube-cbrt3.1
Applied times-frac3.1
Applied associate-*l*1.1
Final simplification1.1
herbie shell --seed 2020046 +o rules:numerics
(FPCore (x y z t)
:name "Numeric.Signal.Multichannel:$cput from hsignal-0.2.7.1"
:precision binary64
:herbie-target
(/ t (/ (- z y) (- x y)))
(* (/ (- x y) (- z y)) t))