x + \left(y - x\right) \cdot \frac{z}{t}x + \left(\sqrt[3]{y - x} \cdot \sqrt[3]{y - x}\right) \cdot \left(\left(\sqrt[3]{y - x} \cdot \frac{\sqrt[3]{z} \cdot \sqrt[3]{z}}{\sqrt[3]{t} \cdot \sqrt[3]{t}}\right) \cdot \frac{\sqrt[3]{z}}{\sqrt[3]{t}}\right)double f(double x, double y, double z, double t) {
double r567307 = x;
double r567308 = y;
double r567309 = r567308 - r567307;
double r567310 = z;
double r567311 = t;
double r567312 = r567310 / r567311;
double r567313 = r567309 * r567312;
double r567314 = r567307 + r567313;
return r567314;
}
double f(double x, double y, double z, double t) {
double r567315 = x;
double r567316 = y;
double r567317 = r567316 - r567315;
double r567318 = cbrt(r567317);
double r567319 = r567318 * r567318;
double r567320 = z;
double r567321 = cbrt(r567320);
double r567322 = r567321 * r567321;
double r567323 = t;
double r567324 = cbrt(r567323);
double r567325 = r567324 * r567324;
double r567326 = r567322 / r567325;
double r567327 = r567318 * r567326;
double r567328 = r567321 / r567324;
double r567329 = r567327 * r567328;
double r567330 = r567319 * r567329;
double r567331 = r567315 + r567330;
return r567331;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 2.0 |
|---|---|
| Target | 2.1 |
| Herbie | 1.8 |
Initial program 2.0
rmApplied add-cube-cbrt2.5
Applied associate-*l*2.5
rmApplied add-cube-cbrt2.6
Applied add-cube-cbrt2.7
Applied times-frac2.7
Applied associate-*r*1.8
Final simplification1.8
herbie shell --seed 2020018
(FPCore (x y z t)
:name "Graphics.Rendering.Plot.Render.Plot.Axis:tickPosition from plot-0.2.3.4"
:precision binary64
:herbie-target
(if (< (* (- y x) (/ z t)) -1013646692435.887) (+ x (/ (- y x) (/ t z))) (if (< (* (- y x) (/ z t)) -0.0) (+ x (/ (* (- y x) z) t)) (+ x (/ (- y x) (/ t z)))))
(+ x (* (- y x) (/ z t))))