x - \frac{y \cdot \left(z - t\right)}{a}\begin{array}{l}
\mathbf{if}\;y \cdot \left(z - t\right) \le -1.8448302221257819 \cdot 10^{209}:\\
\;\;\;\;x - \frac{y}{a} \cdot \left(z - t\right)\\
\mathbf{elif}\;y \cdot \left(z - t\right) \le 4.77051970732716704 \cdot 10^{192}:\\
\;\;\;\;\left(x + \frac{t \cdot y}{a}\right) - \frac{z \cdot y}{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 r274840 = x;
double r274841 = y;
double r274842 = z;
double r274843 = t;
double r274844 = r274842 - r274843;
double r274845 = r274841 * r274844;
double r274846 = a;
double r274847 = r274845 / r274846;
double r274848 = r274840 - r274847;
return r274848;
}
double f(double x, double y, double z, double t, double a) {
double r274849 = y;
double r274850 = z;
double r274851 = t;
double r274852 = r274850 - r274851;
double r274853 = r274849 * r274852;
double r274854 = -1.844830222125782e+209;
bool r274855 = r274853 <= r274854;
double r274856 = x;
double r274857 = a;
double r274858 = r274849 / r274857;
double r274859 = r274858 * r274852;
double r274860 = r274856 - r274859;
double r274861 = 4.770519707327167e+192;
bool r274862 = r274853 <= r274861;
double r274863 = r274851 * r274849;
double r274864 = r274863 / r274857;
double r274865 = r274856 + r274864;
double r274866 = r274850 * r274849;
double r274867 = r274866 / r274857;
double r274868 = r274865 - r274867;
double r274869 = r274857 / r274852;
double r274870 = r274849 / r274869;
double r274871 = r274856 - r274870;
double r274872 = r274862 ? r274868 : r274871;
double r274873 = r274855 ? r274860 : r274872;
return r274873;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 6.2 |
|---|---|
| Target | 0.7 |
| Herbie | 0.4 |
if (* y (- z t)) < -1.844830222125782e+209Initial program 30.2
rmApplied associate-/l*0.7
rmApplied associate-/r/0.5
if -1.844830222125782e+209 < (* y (- z t)) < 4.770519707327167e+192Initial program 0.3
Taylor expanded around 0 0.3
if 4.770519707327167e+192 < (* y (- z t)) Initial program 27.1
rmApplied associate-/l*0.9
Final simplification0.4
herbie shell --seed 2020036
(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)))