x + \left(y - x\right) \cdot \frac{z}{t}\begin{array}{l}
\mathbf{if}\;x \leq -0.10068610656125335 \lor \neg \left(x \leq 4.125760857765473 \cdot 10^{-121}\right):\\
\;\;\;\;x + \left(y - x\right) \cdot \frac{z}{t}\\
\mathbf{else}:\\
\;\;\;\;x + \left(\left(y - x\right) \cdot \left(\sqrt[3]{z} \cdot \sqrt[3]{z}\right)\right) \cdot \frac{\sqrt[3]{z}}{t}\\
\end{array}(FPCore (x y z t) :precision binary64 (+ x (* (- y x) (/ z t))))
(FPCore (x y z t) :precision binary64 (if (or (<= x -0.10068610656125335) (not (<= x 4.125760857765473e-121))) (+ x (* (- y x) (/ z t))) (+ x (* (* (- y x) (* (cbrt z) (cbrt z))) (/ (cbrt 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 ((x <= -0.10068610656125335) || !(x <= 4.125760857765473e-121)) {
tmp = x + ((y - x) * (z / t));
} else {
tmp = x + (((y - x) * (cbrt(z) * cbrt(z))) * (cbrt(z) / t));
}
return tmp;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 2.1 |
|---|---|
| Target | 2.4 |
| Herbie | 1.9 |
if x < -0.100687 or 4.12574e-121 < x Initial program 0.5
if -0.100687 < x < 4.12574e-121Initial program 4.3
rmApplied *-un-lft-identity_binary64_25004.3
Applied add-cube-cbrt_binary64_24714.9
Applied times-frac_binary64_24955.0
Applied associate-*r*_binary64_25593.7
Simplified3.7
Final simplification1.9
herbie shell --seed 2020231
(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))))