x + \frac{\left(y - x\right) \cdot z}{t}x + \frac{\sqrt[3]{y - x} \cdot \sqrt[3]{y - x}}{\frac{\sqrt[3]{t} \cdot \sqrt[3]{t}}{\sqrt[3]{z} \cdot \sqrt[3]{z}}} \cdot \frac{\sqrt[3]{y - x}}{\frac{\sqrt[3]{t}}{\sqrt[3]{z}}}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((y - x))) / ((cbrt(t) * cbrt(t)) / (cbrt(z) * cbrt(z)))) * (cbrt((y - x)) / (cbrt(t) / cbrt(z)))));
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 6.1 |
|---|---|
| Target | 2.1 |
| Herbie | 1.0 |
Initial program 6.1
rmApplied associate-/l*1.9
rmApplied add-cube-cbrt2.5
Applied add-cube-cbrt2.6
Applied times-frac2.6
Applied add-cube-cbrt2.7
Applied times-frac1.0
Final simplification1.0
herbie shell --seed 2020058
(FPCore (x y z t)
:name "Numeric.Histogram:binBounds from Chart-1.5.3"
:precision binary64
:herbie-target
(if (< x -9.025511195533005e-135) (- x (* (/ z t) (- x y))) (if (< x 4.275032163700715e-250) (+ x (* (/ (- y x) t) z)) (+ x (/ (- y x) (/ t z)))))
(+ x (/ (* (- y x) z) t)))