x + \frac{y \cdot \left(z - t\right)}{a}\begin{array}{l}
\mathbf{if}\;y \cdot \left(z - t\right) = -\infty:\\
\;\;\;\;x + \frac{y}{a} \cdot \left(z - t\right)\\
\mathbf{elif}\;y \cdot \left(z - t\right) \le 1.219869957807508 \cdot 10^{198}:\\
\;\;\;\;x + \frac{y \cdot \left(z - t\right)}{a}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{y}{\frac{a}{z - t}}\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r323460 = x;
double r323461 = y;
double r323462 = z;
double r323463 = t;
double r323464 = r323462 - r323463;
double r323465 = r323461 * r323464;
double r323466 = a;
double r323467 = r323465 / r323466;
double r323468 = r323460 + r323467;
return r323468;
}
double f(double x, double y, double z, double t, double a) {
double r323469 = y;
double r323470 = z;
double r323471 = t;
double r323472 = r323470 - r323471;
double r323473 = r323469 * r323472;
double r323474 = -inf.0;
bool r323475 = r323473 <= r323474;
double r323476 = x;
double r323477 = a;
double r323478 = r323469 / r323477;
double r323479 = r323478 * r323472;
double r323480 = r323476 + r323479;
double r323481 = 1.219869957807508e+198;
bool r323482 = r323473 <= r323481;
double r323483 = r323473 / r323477;
double r323484 = r323476 + r323483;
double r323485 = r323477 / r323472;
double r323486 = r323469 / r323485;
double r323487 = r323476 + r323486;
double r323488 = r323482 ? r323484 : r323487;
double r323489 = r323475 ? r323480 : r323488;
return r323489;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 6.5 |
|---|---|
| Target | 0.8 |
| Herbie | 0.4 |
if (* y (- z t)) < -inf.0Initial program 64.0
rmApplied associate-/l*0.2
rmApplied associate-/r/0.2
if -inf.0 < (* y (- z t)) < 1.219869957807508e+198Initial program 0.4
rmApplied *-un-lft-identity0.4
Applied times-frac6.8
Simplified6.8
rmApplied associate-*r/0.4
if 1.219869957807508e+198 < (* y (- z t)) Initial program 28.6
rmApplied associate-/l*0.9
Final simplification0.4
herbie shell --seed 2020056
(FPCore (x y z t a)
:name "Optimisation.CirclePacking:place from circle-packing-0.1.0.4, E"
:precision binary64
:herbie-target
(if (< y -1.0761266216389975e-10) (+ x (/ 1 (/ (/ a (- z t)) y))) (if (< y 2.894426862792089e-49) (+ x (/ (* y (- z t)) a)) (+ x (/ y (/ a (- z t))))))
(+ x (/ (* y (- z t)) a)))