x + \frac{\left(y - x\right) \cdot z}{t}x + \left(\left(y - x\right) \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}}(FPCore (x y z t) :precision binary64 (+ x (/ (* (- y x) z) t)))
(FPCore (x y z t) :precision binary64 (+ x (* (* (- y x) (/ (* (cbrt z) (cbrt z)) (* (cbrt t) (cbrt t)))) (/ (cbrt z) (cbrt 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 + (((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 | 6.3 |
|---|---|
| Target | 2.2 |
| Herbie | 1.0 |
Initial program 6.3
rmApplied *-un-lft-identity_binary64_126956.3
Applied times-frac_binary64_127012.1
Simplified2.1
rmApplied add-cube-cbrt_binary64_127302.6
Applied add-cube-cbrt_binary64_127302.8
Applied times-frac_binary64_127012.8
Applied associate-*r*_binary64_126351.0
Final simplification1.0
herbie shell --seed 2021098
(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)))