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}}}(FPCore (x y z t) :precision binary64 (+ x (/ (* (- y x) z) t)))
(FPCore (x y z t)
:precision binary64
(+
x
(*
(/
(* (cbrt (- y x)) (cbrt (- y x)))
(/ (* (cbrt t) (cbrt t)) (* (cbrt z) (cbrt z))))
(/ (cbrt (- y x)) (/ (cbrt t) (cbrt 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.4 |
|---|---|
| Target | 1.8 |
| Herbie | 1.0 |
Initial program 6.4
rmApplied associate-/l*_binary64_122991.8
rmApplied add-cube-cbrt_binary64_123892.4
Applied add-cube-cbrt_binary64_123892.5
Applied times-frac_binary64_123602.5
Applied add-cube-cbrt_binary64_123892.6
Applied times-frac_binary64_123601.0
Final simplification1.0
herbie shell --seed 2021044
(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)))