x + \frac{y \cdot \left(z - x\right)}{t}\begin{array}{l}
\mathbf{if}\;x + \frac{y \cdot \left(z - x\right)}{t} = -inf.0:\\
\;\;\;\;x + \frac{y}{t} \cdot \left(z - x\right)\\
\mathbf{elif}\;x + \frac{y \cdot \left(z - x\right)}{t} \le 1.7094698165966556 \cdot 10^{306}:\\
\;\;\;\;x + \frac{1}{\frac{t}{y \cdot \left(z - x\right)}}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{y}{\frac{t}{z - x}}\\
\end{array}double code(double x, double y, double z, double t) {
return ((double) (x + ((double) (((double) (y * ((double) (z - x)))) / t))));
}
double code(double x, double y, double z, double t) {
double VAR;
if ((((double) (x + ((double) (((double) (y * ((double) (z - x)))) / t)))) <= -inf.0)) {
VAR = ((double) (x + ((double) (((double) (y / t)) * ((double) (z - x))))));
} else {
double VAR_1;
if ((((double) (x + ((double) (((double) (y * ((double) (z - x)))) / t)))) <= 1.7094698165966556e+306)) {
VAR_1 = ((double) (x + ((double) (1.0 / ((double) (t / ((double) (y * ((double) (z - x))))))))));
} else {
VAR_1 = ((double) (x + ((double) (y / ((double) (t / ((double) (z - x))))))));
}
VAR = VAR_1;
}
return VAR;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 6.6 |
|---|---|
| Target | 2.1 |
| Herbie | 0.8 |
if (+ x (/ (* y (- z x)) t)) < -inf.0Initial program 64.0
rmApplied associate-/l*0.2
rmApplied associate-/r/0.2
if -inf.0 < (+ x (/ (* y (- z x)) t)) < 1.7094698165966556e+306Initial program 0.8
rmApplied clear-num0.8
if 1.7094698165966556e+306 < (+ x (/ (* y (- z x)) t)) Initial program 61.7
rmApplied associate-/l*0.8
Final simplification0.8
herbie shell --seed 2020128
(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)))