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




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 2.1 |
|---|---|
| Target | 2.3 |
| Herbie | 1.7 |
Initial program 2.1
Taylor expanded around 0 6.5
Simplified6.2
rmApplied add-cube-cbrt_binary646.6
Applied add-cube-cbrt_binary646.8
Applied times-frac_binary646.8
Applied associate-*r*_binary641.7
Final simplification1.7
herbie shell --seed 2021118
(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.8867) (+ 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))))