x - \frac{y \cdot \left(z - t\right)}{a}\begin{array}{l}
\mathbf{if}\;y \cdot \left(z - t\right) \le -5.275596583469645 \cdot 10^{271}:\\
\;\;\;\;x - \frac{1}{\frac{\frac{a}{y}}{z - t}}\\
\mathbf{elif}\;y \cdot \left(z - t\right) \le 1.72919691515906818 \cdot 10^{262}:\\
\;\;\;\;x - \frac{1}{a \cdot \frac{1}{y \cdot \left(z - t\right)}}\\
\mathbf{else}:\\
\;\;\;\;x - \frac{y}{\frac{a}{z - t}}\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r286690 = x;
double r286691 = y;
double r286692 = z;
double r286693 = t;
double r286694 = r286692 - r286693;
double r286695 = r286691 * r286694;
double r286696 = a;
double r286697 = r286695 / r286696;
double r286698 = r286690 - r286697;
return r286698;
}
double f(double x, double y, double z, double t, double a) {
double r286699 = y;
double r286700 = z;
double r286701 = t;
double r286702 = r286700 - r286701;
double r286703 = r286699 * r286702;
double r286704 = -5.275596583469645e+271;
bool r286705 = r286703 <= r286704;
double r286706 = x;
double r286707 = 1.0;
double r286708 = a;
double r286709 = r286708 / r286699;
double r286710 = r286709 / r286702;
double r286711 = r286707 / r286710;
double r286712 = r286706 - r286711;
double r286713 = 1.7291969151590682e+262;
bool r286714 = r286703 <= r286713;
double r286715 = r286707 / r286703;
double r286716 = r286708 * r286715;
double r286717 = r286707 / r286716;
double r286718 = r286706 - r286717;
double r286719 = r286708 / r286702;
double r286720 = r286699 / r286719;
double r286721 = r286706 - r286720;
double r286722 = r286714 ? r286718 : r286721;
double r286723 = r286705 ? r286712 : r286722;
return r286723;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 6.4 |
|---|---|
| Target | 0.7 |
| Herbie | 0.5 |
if (* y (- z t)) < -5.275596583469645e+271Initial program 49.3
rmApplied clear-num49.3
rmApplied associate-/r*0.3
if -5.275596583469645e+271 < (* y (- z t)) < 1.7291969151590682e+262Initial program 0.4
rmApplied clear-num0.5
rmApplied div-inv0.5
if 1.7291969151590682e+262 < (* y (- z t)) Initial program 44.4
rmApplied associate-/l*0.2
Final simplification0.5
herbie shell --seed 2020060
(FPCore (x y z t a)
:name "Optimisation.CirclePacking:place from circle-packing-0.1.0.4, F"
: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)))