x + \frac{\left(y - x\right) \cdot z}{t}\begin{array}{l}
\mathbf{if}\;x + \frac{\left(y - x\right) \cdot z}{t} \leq -3.3826999295111967 \cdot 10^{+304}:\\
\;\;\;\;x + \frac{y - x}{\frac{t}{z}}\\
\mathbf{elif}\;x + \frac{\left(y - x\right) \cdot z}{t} \leq 3.227055028057096 \cdot 10^{+294}:\\
\;\;\;\;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 (<= (+ x (/ (* (- y x) z) t)) -3.3826999295111967e+304)
(+ x (/ (- y x) (/ t z)))
(if (<= (+ x (/ (* (- y x) z) t)) 3.227055028057096e+294)
(+ 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 ((x + (((y - x) * z) / t)) <= -3.3826999295111967e+304) {
tmp = x + ((y - x) / (t / z));
} else if ((x + (((y - x) * z) / t)) <= 3.227055028057096e+294) {
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 | 6.9 |
|---|---|
| Target | 2.1 |
| Herbie | 0.9 |
if (+.f64 x (/.f64 (*.f64 (-.f64 y x) z) t)) < -3.38269992951119667e304Initial program 58.7
rmApplied associate-/l*_binary64_144720.5
if -3.38269992951119667e304 < (+.f64 x (/.f64 (*.f64 (-.f64 y x) z) t)) < 3.22705502805709617e294Initial program 0.9
if 3.22705502805709617e294 < (+.f64 x (/.f64 (*.f64 (-.f64 y x) z) t)) Initial program 50.2
rmApplied *-un-lft-identity_binary64_1452750.2
Applied times-frac_binary64_145331.3
Simplified1.3
Final simplification0.9
herbie shell --seed 2020280
(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)))