x + \frac{y \cdot \left(z - x\right)}{t}\begin{array}{l}
\mathbf{if}\;x + \frac{y \cdot \left(z - x\right)}{t} \le -1.5327609267348506 \cdot 10^{261}:\\
\;\;\;\;x + \frac{\frac{y}{t}}{\frac{1}{z - x}}\\
\mathbf{elif}\;x + \frac{y \cdot \left(z - x\right)}{t} \le 2.657445570864299 \cdot 10^{245}:\\
\;\;\;\;x + \frac{y \cdot \left(z - x\right)}{t}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{y}{t} \cdot \left(z - x\right)\\
\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)))) <= -1.5327609267348506e+261)) {
VAR = ((double) (x + ((double) (((double) (y / t)) / ((double) (1.0 / ((double) (z - x))))))));
} else {
double VAR_1;
if ((((double) (x + ((double) (((double) (y * ((double) (z - x)))) / t)))) <= 2.657445570864299e+245)) {
VAR_1 = ((double) (x + ((double) (((double) (y * ((double) (z - x)))) / t))));
} else {
VAR_1 = ((double) (x + ((double) (((double) (y / 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.9 |
|---|---|
| Target | 2.2 |
| Herbie | 1.2 |
if (+ x (/ (* y (- z x)) t)) < -1.5327609267348506e261Initial program 34.3
rmApplied associate-/l*9.1
rmApplied div-inv9.2
Applied associate-/r*2.3
if -1.5327609267348506e261 < (+ x (/ (* y (- z x)) t)) < 2.657445570864299e245Initial program 0.8
if 2.657445570864299e245 < (+ x (/ (* y (- z x)) t)) Initial program 28.5
rmApplied associate-/l*11.0
rmApplied associate-/r/3.0
Final simplification1.2
herbie shell --seed 2020147
(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)) (* (neg z) (/ y t))))
(+ x (/ (* y (- z x)) t)))