x - \frac{y \cdot \left(z - t\right)}{a}\begin{array}{l}
\mathbf{if}\;a \le -1.598281406670949442324391521520754455565 \cdot 10^{-10}:\\
\;\;\;\;x + \frac{y}{\frac{a}{t - z}}\\
\mathbf{elif}\;a \le 2.399671998072701683137690037832409821452 \cdot 10^{-142}:\\
\;\;\;\;\left(x + \frac{t \cdot y}{a}\right) - \frac{y \cdot z}{a}\\
\mathbf{elif}\;a \le 2.181728582825251587184966617016023093761 \cdot 10^{52}:\\
\;\;\;\;x + \frac{\frac{y}{a}}{\frac{1}{t - z}}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{y}{\frac{a}{t - z}}\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r826717 = x;
double r826718 = y;
double r826719 = z;
double r826720 = t;
double r826721 = r826719 - r826720;
double r826722 = r826718 * r826721;
double r826723 = a;
double r826724 = r826722 / r826723;
double r826725 = r826717 - r826724;
return r826725;
}
double f(double x, double y, double z, double t, double a) {
double r826726 = a;
double r826727 = -1.5982814066709494e-10;
bool r826728 = r826726 <= r826727;
double r826729 = x;
double r826730 = y;
double r826731 = t;
double r826732 = z;
double r826733 = r826731 - r826732;
double r826734 = r826726 / r826733;
double r826735 = r826730 / r826734;
double r826736 = r826729 + r826735;
double r826737 = 2.3996719980727017e-142;
bool r826738 = r826726 <= r826737;
double r826739 = r826731 * r826730;
double r826740 = r826739 / r826726;
double r826741 = r826729 + r826740;
double r826742 = r826730 * r826732;
double r826743 = r826742 / r826726;
double r826744 = r826741 - r826743;
double r826745 = 2.1817285828252516e+52;
bool r826746 = r826726 <= r826745;
double r826747 = r826730 / r826726;
double r826748 = 1.0;
double r826749 = r826748 / r826733;
double r826750 = r826747 / r826749;
double r826751 = r826729 + r826750;
double r826752 = r826746 ? r826751 : r826736;
double r826753 = r826738 ? r826744 : r826752;
double r826754 = r826728 ? r826736 : r826753;
return r826754;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 6.0 |
|---|---|
| Target | 0.7 |
| Herbie | 0.8 |
if a < -1.5982814066709494e-10 or 2.1817285828252516e+52 < a Initial program 9.6
Simplified1.7
rmApplied *-un-lft-identity1.7
Applied *-un-lft-identity1.7
Applied times-frac1.7
Applied associate-*l*1.7
Simplified0.5
if -1.5982814066709494e-10 < a < 2.3996719980727017e-142Initial program 1.1
Simplified4.4
rmApplied sub-neg4.4
Applied distribute-lft-in4.4
Simplified4.4
rmApplied *-un-lft-identity4.4
Applied *-un-lft-identity4.4
Applied distribute-lft-out4.4
Simplified1.1
if 2.3996719980727017e-142 < a < 2.1817285828252516e+52Initial program 1.1
Simplified1.2
rmApplied *-un-lft-identity1.2
Applied *-un-lft-identity1.2
Applied times-frac1.2
Applied associate-*l*1.2
Simplified5.2
rmApplied div-inv5.3
Applied associate-/r*1.3
Final simplification0.8
herbie shell --seed 2019179
(FPCore (x y z t a)
:name "Optimisation.CirclePacking:place from circle-packing-0.1.0.4, F"
:herbie-target
(if (< y -1.0761266216389975e-10) (- x (/ 1.0 (/ (/ a (- z t)) y))) (if (< y 2.894426862792089e-49) (- x (/ (* y (- z t)) a)) (- x (/ y (/ a (- z t))))))
(- x (/ (* y (- z t)) a)))