x + \frac{y \cdot \left(z - x\right)}{t}\begin{array}{l}
\mathbf{if}\;x + \frac{y \cdot \left(z - x\right)}{t} = -\infty \lor \neg \left(x + \frac{y \cdot \left(z - x\right)}{t} \le 1.83014942255142419 \cdot 10^{297}\right):\\
\;\;\;\;x + \frac{y}{\frac{t}{z - x}}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{y \cdot \left(z - x\right)}{t}\\
\end{array}double f(double x, double y, double z, double t) {
double r376982 = x;
double r376983 = y;
double r376984 = z;
double r376985 = r376984 - r376982;
double r376986 = r376983 * r376985;
double r376987 = t;
double r376988 = r376986 / r376987;
double r376989 = r376982 + r376988;
return r376989;
}
double f(double x, double y, double z, double t) {
double r376990 = x;
double r376991 = y;
double r376992 = z;
double r376993 = r376992 - r376990;
double r376994 = r376991 * r376993;
double r376995 = t;
double r376996 = r376994 / r376995;
double r376997 = r376990 + r376996;
double r376998 = -inf.0;
bool r376999 = r376997 <= r376998;
double r377000 = 1.8301494225514242e+297;
bool r377001 = r376997 <= r377000;
double r377002 = !r377001;
bool r377003 = r376999 || r377002;
double r377004 = r376995 / r376993;
double r377005 = r376991 / r377004;
double r377006 = r376990 + r377005;
double r377007 = r377003 ? r377006 : r376997;
return r377007;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 6.2 |
|---|---|
| Target | 2.0 |
| Herbie | 0.7 |
if (+ x (/ (* y (- z x)) t)) < -inf.0 or 1.8301494225514242e+297 < (+ x (/ (* y (- z x)) t)) Initial program 58.1
rmApplied associate-/l*1.7
if -inf.0 < (+ x (/ (* y (- z x)) t)) < 1.8301494225514242e+297Initial program 0.6
Final simplification0.7
herbie shell --seed 2020046
(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)))