x + \left(y - x\right) \cdot \frac{z}{t}\begin{array}{l}
\mathbf{if}\;t \leq -6.313002064818998 \cdot 10^{+86}:\\
\;\;\;\;x + \left(\sqrt[3]{y - x} \cdot \sqrt[3]{y - x}\right) \cdot \left(z \cdot \frac{\sqrt[3]{y - x}}{t}\right)\\
\mathbf{elif}\;t \leq 4.633125479398495 \cdot 10^{+59}:\\
\;\;\;\;x + \frac{\left(y - x\right) \cdot z}{t}\\
\mathbf{else}:\\
\;\;\;\;x + \left(y - x\right) \cdot \frac{z}{t}\\
\end{array}(FPCore (x y z t) :precision binary64 (+ x (* (- y x) (/ z t))))
(FPCore (x y z t)
:precision binary64
(if (<= t -6.313002064818998e+86)
(+ x (* (* (cbrt (- y x)) (cbrt (- y x))) (* z (/ (cbrt (- y x)) t))))
(if (<= t 4.633125479398495e+59)
(+ x (/ (* (- y x) z) t))
(+ x (* (- y x) (/ z 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) {
double tmp;
if (t <= -6.313002064818998e+86) {
tmp = x + ((cbrt(y - x) * cbrt(y - x)) * (z * (cbrt(y - x) / t)));
} else if (t <= 4.633125479398495e+59) {
tmp = x + (((y - x) * z) / t);
} else {
tmp = x + ((y - x) * (z / t));
}
return tmp;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 2.3 |
|---|---|
| Target | 2.5 |
| Herbie | 1.8 |
if t < -6.3130020648189983e86Initial program 1.6
rmApplied add-cube-cbrt_binary641.9
Applied associate-*l*_binary641.9
Simplified1.9
rmApplied div-inv_binary641.9
Applied associate-*l*_binary641.1
Simplified1.1
if -6.3130020648189983e86 < t < 4.6331254793984947e59Initial program 3.2
rmApplied associate-*r/_binary642.4
if 4.6331254793984947e59 < t Initial program 1.3
Final simplification1.8
herbie shell --seed 2020253
(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))))