x + \left(y - x\right) \cdot \frac{z}{t}\begin{array}{l}
\mathbf{if}\;t \leq -3.947928725274646 \cdot 10^{+86}:\\
\;\;\;\;x + \frac{y - x}{\sqrt[3]{t} \cdot \sqrt[3]{t}} \cdot \frac{z}{\sqrt[3]{t}}\\
\mathbf{elif}\;t \leq 9.835726741021531 \cdot 10^{+41}:\\
\;\;\;\;x + \left(\left(y - x\right) \cdot z\right) \cdot \frac{1}{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 -3.947928725274646e+86)
(+ x (* (/ (- y x) (* (cbrt t) (cbrt t))) (/ z (cbrt t))))
(if (<= t 9.835726741021531e+41)
(+ x (* (* (- y x) z) (/ 1.0 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 <= -3.947928725274646e+86) {
tmp = x + (((y - x) / (cbrt(t) * cbrt(t))) * (z / cbrt(t)));
} else if (t <= 9.835726741021531e+41) {
tmp = x + (((y - x) * z) * (1.0 / 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 | 1.9 |
|---|---|
| Target | 2.2 |
| Herbie | 1.5 |
if t < -3.94792872527464621e86Initial program 1.0
rmApplied add-cube-cbrt_binary641.3
Applied *-un-lft-identity_binary641.3
Applied times-frac_binary641.3
Applied associate-*r*_binary640.9
Simplified0.9
if -3.94792872527464621e86 < t < 9.835726741021531e41Initial program 3.0
rmApplied div-inv_binary643.1
Applied associate-*r*_binary642.2
if 9.835726741021531e41 < t Initial program 0.9
Final simplification1.5
herbie shell --seed 2020224
(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))))