x + \left(y - x\right) \cdot \frac{z}{t}x + \left(\left(\sqrt[3]{\left(y - x\right) \cdot \frac{\sqrt[3]{z} \cdot \sqrt[3]{z}}{\sqrt[3]{t} \cdot \sqrt[3]{t}}} \cdot \sqrt[3]{\left(y - x\right) \cdot \frac{\sqrt[3]{z} \cdot \sqrt[3]{z}}{\sqrt[3]{t} \cdot \sqrt[3]{t}}}\right) \cdot \sqrt[3]{\left(\sqrt[3]{y - x} \cdot \sqrt[3]{y - x}\right) \cdot \left(\sqrt[3]{y - x} \cdot \frac{\sqrt[3]{z} \cdot \sqrt[3]{z}}{\sqrt[3]{t} \cdot \sqrt[3]{t}}\right)}\right) \cdot \frac{\sqrt[3]{z}}{\sqrt[3]{t}}double code(double x, double y, double z, double t) {
return (x + ((y - x) * (z / t)));
}
double code(double x, double y, double z, double t) {
return (x + (((cbrt(((y - x) * ((cbrt(z) * cbrt(z)) / (cbrt(t) * cbrt(t))))) * cbrt(((y - x) * ((cbrt(z) * cbrt(z)) / (cbrt(t) * cbrt(t)))))) * cbrt(((cbrt((y - x)) * cbrt((y - x))) * (cbrt((y - x)) * ((cbrt(z) * cbrt(z)) / (cbrt(t) * cbrt(t))))))) * (cbrt(z) / cbrt(t))));
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 2.0 |
|---|---|
| Target | 2.2 |
| Herbie | 1.0 |
Initial program 2.0
rmApplied add-cube-cbrt2.6
Applied add-cube-cbrt2.7
Applied times-frac2.7
Applied associate-*r*0.9
rmApplied add-cube-cbrt1.0
rmApplied add-cube-cbrt1.0
Applied associate-*l*1.0
Final simplification1.0
herbie shell --seed 2020078 +o rules:numerics
(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))))