x + \frac{y \cdot \left(z - x\right)}{t}\begin{array}{l}
\mathbf{if}\;x + \frac{y \cdot \left(z - x\right)}{t} \leq -\infty \lor \neg \left(x + \frac{y \cdot \left(z - x\right)}{t} \leq 5.935478454692787 \cdot 10^{+299}\right):\\
\;\;\;\;x + y \cdot \frac{z - x}{t}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{y \cdot \left(z - x\right)}{t}\\
\end{array}(FPCore (x y z t) :precision binary64 (+ x (/ (* y (- z x)) t)))
(FPCore (x y z t)
:precision binary64
(if (or (<= (+ x (/ (* y (- z x)) t)) (- INFINITY))
(not (<= (+ x (/ (* y (- z x)) t)) 5.935478454692787e+299)))
(+ x (* y (/ (- z x) t)))
(+ x (/ (* y (- z x)) t))))double code(double x, double y, double z, double t) {
return x + ((y * (z - x)) / t);
}
double code(double x, double y, double z, double t) {
double tmp;
if (((x + ((y * (z - x)) / t)) <= -((double) INFINITY)) || !((x + ((y * (z - x)) / t)) <= 5.935478454692787e+299)) {
tmp = x + (y * ((z - x) / t));
} else {
tmp = x + ((y * (z - x)) / t);
}
return tmp;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 6.8 |
|---|---|
| Target | 2.0 |
| Herbie | 1.0 |
if (+.f64 x (/.f64 (*.f64 y (-.f64 z x)) t)) < -inf.0 or 5.93547845469278713e299 < (+.f64 x (/.f64 (*.f64 y (-.f64 z x)) t)) Initial program 59.9
rmApplied *-un-lft-identity_binary64_1030859.9
Applied times-frac_binary64_103141.8
if -inf.0 < (+.f64 x (/.f64 (*.f64 y (-.f64 z x)) t)) < 5.93547845469278713e299Initial program 0.9
Final simplification1.0
herbie shell --seed 2021093
(FPCore (x y z t)
:name "Optimisation.CirclePacking:place from circle-packing-0.1.0.4, D"
:precision binary64
:herbie-target
(- x (+ (* x (/ y t)) (* (- z) (/ y t))))
(+ x (/ (* y (- z x)) t)))