x + \frac{y \cdot \left(z - x\right)}{t}\begin{array}{l}
\mathbf{if}\;x + \frac{\left(z - x\right) \cdot y}{t} \le -4.015242760717319860347263184406053785305 \cdot 10^{305}:\\
\;\;\;\;x + \frac{y}{\frac{t}{z - x}}\\
\mathbf{elif}\;x + \frac{\left(z - x\right) \cdot y}{t} \le 2.038097847986346269413799617144504707187 \cdot 10^{304}:\\
\;\;\;\;x + \frac{\left(z - x\right) \cdot y}{t}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{y}{\frac{t}{z - x}}\\
\end{array}double f(double x, double y, double z, double t) {
double r12090413 = x;
double r12090414 = y;
double r12090415 = z;
double r12090416 = r12090415 - r12090413;
double r12090417 = r12090414 * r12090416;
double r12090418 = t;
double r12090419 = r12090417 / r12090418;
double r12090420 = r12090413 + r12090419;
return r12090420;
}
double f(double x, double y, double z, double t) {
double r12090421 = x;
double r12090422 = z;
double r12090423 = r12090422 - r12090421;
double r12090424 = y;
double r12090425 = r12090423 * r12090424;
double r12090426 = t;
double r12090427 = r12090425 / r12090426;
double r12090428 = r12090421 + r12090427;
double r12090429 = -4.01524276071732e+305;
bool r12090430 = r12090428 <= r12090429;
double r12090431 = r12090426 / r12090423;
double r12090432 = r12090424 / r12090431;
double r12090433 = r12090421 + r12090432;
double r12090434 = 2.0380978479863463e+304;
bool r12090435 = r12090428 <= r12090434;
double r12090436 = r12090435 ? r12090428 : r12090433;
double r12090437 = r12090430 ? r12090433 : r12090436;
return r12090437;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 6.7 |
|---|---|
| Target | 2.1 |
| Herbie | 0.9 |
if (+ x (/ (* y (- z x)) t)) < -4.01524276071732e+305 or 2.0380978479863463e+304 < (+ x (/ (* y (- z x)) t)) Initial program 60.1
rmApplied associate-/l*1.3
if -4.01524276071732e+305 < (+ x (/ (* y (- z x)) t)) < 2.0380978479863463e+304Initial program 0.8
Final simplification0.9
herbie shell --seed 2019172 +o rules:numerics
(FPCore (x y z t)
:name "Optimisation.CirclePacking:place from circle-packing-0.1.0.4, D"
:herbie-target
(- x (+ (* x (/ y t)) (* (- z) (/ y t))))
(+ x (/ (* y (- z x)) t)))