x + \frac{y \cdot \left(z - t\right)}{a}\begin{array}{l}
\mathbf{if}\;y \cdot \left(z - t\right) \le -3.7422788828719246 \cdot 10^{161}:\\
\;\;\;\;x + y \cdot \frac{z - t}{a}\\
\mathbf{elif}\;y \cdot \left(z - t\right) \le 2.424193312752178 \cdot 10^{192}:\\
\;\;\;\;x + \frac{1}{\frac{a}{y \cdot \left(z - t\right)}}\\
\mathbf{else}:\\
\;\;\;\;\left(x + z \cdot \frac{y}{a}\right) + \frac{y}{a} \cdot \left(-t\right)\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r584018 = x;
double r584019 = y;
double r584020 = z;
double r584021 = t;
double r584022 = r584020 - r584021;
double r584023 = r584019 * r584022;
double r584024 = a;
double r584025 = r584023 / r584024;
double r584026 = r584018 + r584025;
return r584026;
}
double f(double x, double y, double z, double t, double a) {
double r584027 = y;
double r584028 = z;
double r584029 = t;
double r584030 = r584028 - r584029;
double r584031 = r584027 * r584030;
double r584032 = -3.7422788828719246e+161;
bool r584033 = r584031 <= r584032;
double r584034 = x;
double r584035 = a;
double r584036 = r584030 / r584035;
double r584037 = r584027 * r584036;
double r584038 = r584034 + r584037;
double r584039 = 2.424193312752178e+192;
bool r584040 = r584031 <= r584039;
double r584041 = 1.0;
double r584042 = r584035 / r584031;
double r584043 = r584041 / r584042;
double r584044 = r584034 + r584043;
double r584045 = r584027 / r584035;
double r584046 = r584028 * r584045;
double r584047 = r584034 + r584046;
double r584048 = -r584029;
double r584049 = r584045 * r584048;
double r584050 = r584047 + r584049;
double r584051 = r584040 ? r584044 : r584050;
double r584052 = r584033 ? r584038 : r584051;
return r584052;
}




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.6 |
| Herbie | 0.7 |
if (* y (- z t)) < -3.7422788828719246e+161Initial program 21.2
rmApplied *-un-lft-identity21.2
Applied times-frac1.7
Simplified1.7
if -3.7422788828719246e+161 < (* y (- z t)) < 2.424193312752178e+192Initial program 0.4
rmApplied clear-num0.4
if 2.424193312752178e+192 < (* y (- z t)) Initial program 26.2
rmApplied associate-/l*1.1
rmApplied associate-/r/1.0
rmApplied sub-neg1.0
Applied distribute-lft-in1.0
Applied associate-+r+1.0
Simplified1.0
Final simplification0.7
herbie shell --seed 2020042
(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)))