x + \frac{y \cdot \left(z - t\right)}{a}\begin{array}{l}
\mathbf{if}\;y \cdot \left(z - t\right) \le -2.689959535719968 \cdot 10^{288}:\\
\;\;\;\;x + \frac{y}{a} \cdot \left(z - t\right)\\
\mathbf{elif}\;y \cdot \left(z - t\right) \le 3.12395947034459992 \cdot 10^{186}:\\
\;\;\;\;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 r308089 = x;
double r308090 = y;
double r308091 = z;
double r308092 = t;
double r308093 = r308091 - r308092;
double r308094 = r308090 * r308093;
double r308095 = a;
double r308096 = r308094 / r308095;
double r308097 = r308089 + r308096;
return r308097;
}
double f(double x, double y, double z, double t, double a) {
double r308098 = y;
double r308099 = z;
double r308100 = t;
double r308101 = r308099 - r308100;
double r308102 = r308098 * r308101;
double r308103 = -2.689959535719968e+288;
bool r308104 = r308102 <= r308103;
double r308105 = x;
double r308106 = a;
double r308107 = r308098 / r308106;
double r308108 = r308107 * r308101;
double r308109 = r308105 + r308108;
double r308110 = 3.1239594703446e+186;
bool r308111 = r308102 <= r308110;
double r308112 = r308102 / r308106;
double r308113 = r308105 + r308112;
double r308114 = r308106 / r308101;
double r308115 = r308098 / r308114;
double r308116 = r308105 + r308115;
double r308117 = r308111 ? r308113 : r308116;
double r308118 = r308104 ? r308109 : r308117;
return r308118;
}




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.8 |
| Herbie | 0.5 |
if (* y (- z t)) < -2.689959535719968e+288Initial program 54.0
rmApplied clear-num54.0
rmApplied associate-/r*0.3
rmApplied div-inv0.4
Applied add-cube-cbrt0.4
Applied times-frac0.4
Simplified0.3
Simplified0.2
if -2.689959535719968e+288 < (* y (- z t)) < 3.1239594703446e+186Initial program 0.5
if 3.1239594703446e+186 < (* y (- z t)) Initial program 26.6
rmApplied associate-/l*1.1
Final simplification0.5
herbie shell --seed 2020081
(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)))