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 r13264303 = x;
double r13264304 = y;
double r13264305 = z;
double r13264306 = r13264305 - r13264303;
double r13264307 = r13264304 * r13264306;
double r13264308 = t;
double r13264309 = r13264307 / r13264308;
double r13264310 = r13264303 + r13264309;
return r13264310;
}
double f(double x, double y, double z, double t) {
double r13264311 = x;
double r13264312 = z;
double r13264313 = r13264312 - r13264311;
double r13264314 = y;
double r13264315 = r13264313 * r13264314;
double r13264316 = t;
double r13264317 = r13264315 / r13264316;
double r13264318 = r13264311 + r13264317;
double r13264319 = -4.01524276071732e+305;
bool r13264320 = r13264318 <= r13264319;
double r13264321 = r13264316 / r13264313;
double r13264322 = r13264314 / r13264321;
double r13264323 = r13264311 + r13264322;
double r13264324 = 2.0380978479863463e+304;
bool r13264325 = r13264318 <= r13264324;
double r13264326 = r13264325 ? r13264318 : r13264323;
double r13264327 = r13264320 ? r13264323 : r13264326;
return r13264327;
}




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)))