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:\\
\;\;\;\;x + y \cdot \frac{z - x}{t}\\
\mathbf{elif}\;x + \frac{y \cdot \left(z - x\right)}{t} \leq 6.73935946189453 \cdot 10^{+279}:\\
\;\;\;\;x + \frac{y \cdot \left(z - x\right)}{t}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{1}{\frac{\frac{t}{y}}{z - x}}\\
\end{array}(FPCore (x y z t) :precision binary64 (+ x (/ (* y (- z x)) t)))
(FPCore (x y z t)
:precision binary64
(if (<= (+ x (/ (* y (- z x)) t)) (- INFINITY))
(+ x (* y (/ (- z x) t)))
(if (<= (+ x (/ (* y (- z x)) t)) 6.73935946189453e+279)
(+ x (/ (* y (- z x)) t))
(+ x (/ 1.0 (/ (/ t y) (- z x)))))))double code(double x, double y, double z, double t) {
return ((double) (x + (((double) (y * ((double) (z - x)))) / t)));
}
double code(double x, double y, double z, double t) {
double tmp;
if ((((double) (x + (((double) (y * ((double) (z - x)))) / t))) <= ((double) -(((double) INFINITY))))) {
tmp = ((double) (x + ((double) (y * (((double) (z - x)) / t)))));
} else {
double tmp_1;
if ((((double) (x + (((double) (y * ((double) (z - x)))) / t))) <= 6.73935946189453e+279)) {
tmp_1 = ((double) (x + (((double) (y * ((double) (z - x)))) / t)));
} else {
tmp_1 = ((double) (x + (1.0 / ((t / y) / ((double) (z - x))))));
}
tmp = tmp_1;
}
return tmp;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 6.1 |
|---|---|
| Target | 2.2 |
| Herbie | 0.8 |
if (+.f64 x (/.f64 (*.f64 y (-.f64 z x)) t)) < -inf.0Initial program 64.0
rmApplied *-un-lft-identity_binary6464.0
Applied times-frac_binary640.2
Simplified0.2
if -inf.0 < (+.f64 x (/.f64 (*.f64 y (-.f64 z x)) t)) < 6.7393594618945302e279Initial program 0.7
if 6.7393594618945302e279 < (+.f64 x (/.f64 (*.f64 y (-.f64 z x)) t)) Initial program 40.6
rmApplied clear-num_binary6440.6
rmApplied associate-/r*_binary641.9
Final simplification0.8
herbie shell --seed 2020210
(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)))