x + \frac{\left(y - x\right) \cdot z}{t}\begin{array}{l}
\mathbf{if}\;x \leq -5.816891157054731 \cdot 10^{-255}:\\
\;\;\;\;x + \frac{y - x}{\frac{t}{z}}\\
\mathbf{elif}\;x \leq 5.688424274737943 \cdot 10^{-196}:\\
\;\;\;\;x + \frac{y - x}{\sqrt[3]{t} \cdot \sqrt[3]{t}} \cdot \frac{z}{\sqrt[3]{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 (<= x -5.816891157054731e-255)
(+ x (/ (- y x) (/ t z)))
(if (<= x 5.688424274737943e-196)
(+ x (* (/ (- y x) (* (cbrt t) (cbrt t))) (/ z (cbrt 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 (x <= -5.816891157054731e-255) {
tmp = x + ((y - x) / (t / z));
} else if (x <= 5.688424274737943e-196) {
tmp = x + (((y - x) / (cbrt(t) * cbrt(t))) * (z / cbrt(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 | 6.3 |
|---|---|
| Target | 1.9 |
| Herbie | 1.8 |
if x < -5.81689115705473114e-255Initial program 6.0
rmApplied associate-/l*_binary641.5
if -5.81689115705473114e-255 < x < 5.68842427473794347e-196Initial program 5.9
rmApplied add-cube-cbrt_binary646.7
Applied times-frac_binary644.8
if 5.68842427473794347e-196 < x Initial program 6.8
rmApplied *-un-lft-identity_binary646.8
Applied times-frac_binary641.1
Simplified1.1
Final simplification1.8
herbie shell --seed 2020233
(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)))