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




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 6.8 |
|---|---|
| Target | 2.0 |
| Herbie | 1.9 |
Initial program 6.8
Simplified2.0
rmApplied *-un-lft-identity_binary64_79212.0
Applied add-cube-cbrt_binary64_79562.5
Applied times-frac_binary64_79272.5
Applied add-cube-cbrt_binary64_79562.6
Applied times-frac_binary64_79271.9
Final simplification1.9
herbie shell --seed 2021176
(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)))