{ "cells": [ { "cell_type": "code", "execution_count": 106, "metadata": { "collapsed": true }, "outputs": [], "source": [ "from pathlib import Path\n", "import pandas as pd\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "import scipy\n", "from glob import glob" ] }, { "cell_type": "markdown", "source": [ "Prior code from mva:" ], "metadata": { "collapsed": false } }, { "cell_type": "code", "execution_count": 112, "outputs": [], "source": [ "def read_xy(path: str):\n", " df = pd.read_csv(path, skipinitialspace=True)\n", " df['N'] = df.index + 1\n", " df['r'] = (df.x ** 2 + df.y ** 2) ** 0.5\n", " df['cr'] = df.r.cummax()\n", " df['fd'] = np.log(df.N) / np.log(df.cr)\n", "\n", " return df\n", "\n", "\n", "def read_load_dir(load_dir: str):\n", " paths = glob(f'{load_dir}/*.csv')\n", " return [read_xy(path) for path in paths]\n", "\n", "\n", "# We expect the format of sp data to be:\n", "# ROOT/{probability}/{RUN_ID}.csv\n", "\n", "def read_sp_load_dir(load_dir: str):\n", " probability = float(Path(load_dir).name)\n", " return probability, read_load_dir(load_dir)\n", "\n", "\n", "def read_sp_full(root: str):\n", " if not Path(root).exists():\n", " raise Exception(\"Root does not exist\")\n", "\n", " raw_data = [read_sp_load_dir(load_dir) for load_dir in glob(f'{root}/*')]\n", " fd_data = [fd_stats(dfs) for (p, dfs) in raw_data]\n", " probabilities = [p for (p, dfs) in raw_data]\n", " fd_means, fd_stds = list(zip(*fd_data))\n", " out = pd.DataFrame(zip(probabilities, fd_means, fd_stds), columns=['p', 'fd_mean', 'fd_std'])\n", " out.sort_values(by='p', inplace=True)\n", "\n", " return out" ], "metadata": { "collapsed": false } }, { "cell_type": "code", "execution_count": null, "outputs": [], "source": [ "def convergent_tail_index(series, tol):\n", " diffs = np.abs(np.ediff1d(series))\n", " for i in range(0, len(diffs)):\n", " if np.max(diffs[i:]) <= tol:\n", " return i\n", "\n", " # No convergence found\n", " return None\n", "\n", "\n", "def mean_of_tail(series, tol=0.05):\n", " tail_index = convergent_tail_index(series, tol)\n", " if tail_index is None:\n", " raise Exception(\"No convergence found.\")\n", "\n", " return np.mean(series[tail_index:])\n", "\n", "\n", "def fd_stats(dfs):\n", " fds = [mean_of_tail(df.fd, 0.001) for df in dfs]\n", " fds_clean = [f for f in fds if f < np.inf]\n", " return np.mean(fds_clean), np.std(fds_clean)" ], "metadata": { "collapsed": false } }, { "cell_type": "code", "execution_count": 113, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/Users/joshuacoles/Library/Caches/pypoetry/virtualenvs/data-analysis-B4Au_hWl-py3.10/lib/python3.10/site-packages/numpy/lib/arraysetops.py:89: RuntimeWarning: invalid value encountered in subtract\n", " return ary[1:] - ary[:-1]\n" ] } ], "source": [ "c_sticking_probability = read_sp_full(\"../data/stick-probability\")" ], "metadata": { "collapsed": false } }, { "cell_type": "code", "execution_count": 121, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/Users/joshuacoles/Library/Caches/pypoetry/virtualenvs/data-analysis-B4Au_hWl-py3.10/lib/python3.10/site-packages/numpy/lib/arraysetops.py:89: RuntimeWarning: invalid value encountered in subtract\n", " return ary[1:] - ary[:-1]\n" ] } ], "source": [ "rust_sticking_probability = read_sp_full(\"../data/rust-sticking-probability\")" ], "metadata": { "collapsed": false } }, { "cell_type": "code", "execution_count": 151, "outputs": [ { "data": { "text/plain": "" }, "execution_count": 151, "metadata": {}, "output_type": "execute_result" }, { "data": { "text/plain": "
", "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiMAAAGdCAYAAADAAnMpAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy88F64QAAAACXBIWXMAAA9hAAAPYQGoP6dpAABXc0lEQVR4nO3dd3wUdf7H8dfspncCBBJI6L33KggiiMjJWe9AAUVswYZnQc/6U1EPPfVEsFBsCCpFRSyI0nsJAgk9kAChk942u/P7YyUSaQkkmZT38/GYx7FTPzPusW++853vGKZpmoiIiIhYxGZ1ASIiIlK5KYyIiIiIpRRGRERExFIKIyIiImIphRERERGxlMKIiIiIWEphRERERCylMCIiIiKW8rC6gMJwuVwcOnSIwMBADMOwuhwREREpBNM0SUtLIyIiApvt/O0f5SKMHDp0iMjISKvLEBERkUuQmJhI7dq1z7u8XISRwMBAwH0yQUFBFlcjIiIihZGamkpkZGT+7/j5lIswcvrWTFBQkMKIiIhIOXOxLhbqwCoiIiKWUhgRERERSymMiIiIiKXKRZ8RERG5NKZpkpeXh9PptLoUqYDsdjseHh6XPeyGwoiISAWVm5tLUlISmZmZVpciFZifnx/h4eF4eXld8j4URkREKiCXy0V8fDx2u52IiAi8vLw0aKQUK9M0yc3N5dixY8THx9OoUaMLDmx2IQojIiIVUG5uLi6Xi8jISPz8/KwuRyooX19fPD092b9/P7m5ufj4+FzSftSBVUSkArvUf6mKFFZxfMf0LRURERFLKYyIiIiIpRRGRERExFJFCiPjx4+nU6dOBAYGEhYWxpAhQ9ixY8cFt9m2bRs33ngjdevWxTAM3nrrrcupV0REKoHDhw/zwAMPUL9+fby9vYmMjGTw4MEsWrTI6tKkBBQpjCxZsoTo6GhWr17NwoULcTgc9O/fn4yMjPNuk5mZSf369Xn11VepWbPmZRdcrLZ8Dd8+CInrrK5ERET+sG/fPjp06MCvv/7Kf/7zH7Zs2cKPP/5Inz59iI6Otro8KQFFerT3xx9/LPB5+vTphIWFsWHDBnr16nXObTp16kSnTp0AePLJJy+xzBKyfT5smwvVGkFkJ6urEREpUaZpkuUo/ZFYfT3tRRrj5P7778cwDNauXYu/v3/+/BYtWnDnnXeWRIliscsaZyQlJQWA0NDQYinmtJycHHJycvI/p6amFuv+T4vNq0Vz4PjeTVTrXiKHEBEpM7IcTpo/+1OpHzf2xQH4eRXu5+bkyZP8+OOPvPzyywWCyGkhISHFXJ2UBZfcgdXlcvHwww/To0cPWrZsWZw1MX78eIKDg/OnyMjIYt3/actTwwAwj8SWyP5FRKRodu/ejWmaNG3a1OpSpBRdcstIdHQ0W7duZfny5cVZDwDjxo1j7Nix+Z9TU1NLJJB4RbSAJAhO3wsuJ9jsxX4MEZGywtfTTuyLAyw5bmGZplmClUhZdUlhZMyYMcyfP5+lS5dSu3bt4q4Jb29vvL29i32/f1UtqgnZ6z3xIQdO7YOqDUr8mCIiVjEMo9C3S6zSqFEjDMNg+/btVpcipahIt2lM02TMmDHMnTuXX3/9lXr16pVUXaWiYY1gdpm1ADCP6laNiIjVQkNDGTBgABMnTjznk5rJycmlX5SUuCKFkejoaD777DNmzJhBYGAghw8f5vDhw2RlZeWvM3z4cMaNG5f/OTc3l5iYGGJiYsjNzeXgwYPExMSwe/fu4juLS1Svmj87Tfftn4zELRZXIyIiABMnTsTpdNK5c2dmz57Nrl27iIuL45133qFbt25WlycloEhhZNKkSaSkpHDllVcSHh6eP82aNSt/nYSEBJKSkvI/Hzp0iHbt2tGuXTuSkpKYMGEC7dq146677iq+s7hE3h52jvnWByDr4FaLqxEREYD69euzceNG+vTpw6OPPkrLli25+uqrWbRoEZMmTbK6PCkBRbp5WJiORYsXLy7wuW7dumW6Q1J2lSZwBDyO6/6kiEhZER4ezrvvvsu7775rdSlSCir9u2m8ItyPJQdl7IO8XGuLERERqYQqfRipUas+qaYvdpxwwvp+LCIiIpVNpQ8jDWsE5ndiRU/UiIiIlLpKH0YahAWw0+UeKyX7kDqxioiIlLZKH0YCvD047OMeLyXrgMKIiIhIaav0YQT+eKIG8DihJ2pERERKm8II4BneAoCAzAOQe/aIfyIiIlJyFEaAWrWiOGYGYWDCsR1WlyMiIlKpKIwADcMC2Ok6/URNnLXFiIiIVDIKI/wRRkz3EzWOw9ssrkZEpHIbOXIkhmFgGAaenp7Uq1ePxx9/nOzs7GLZf926dXnrrbeKZV9SPMr2u6RLSai/Fwc964IJ2Qe34ml1QSIildw111zDtGnTcDgcbNiwgREjRmAYBq+99prVpUkJUMvIH7JC/3iiRu+oERGxnLe3NzVr1iQyMpIhQ4bQr18/Fi5cCJy7ZaNt27Y8//zzgPs9as8//zxRUVF4e3sTERHBgw8+CMCVV17J/v37eeSRR/JbX8R6ahn5g1d4CzgBvtlHIOsU+FaxuiQRkeJlmuDILP3jevrBZfzob926lZUrV1KnTp1CrT979mz++9//MnPmTFq0aMHhw4fZvHkzAHPmzKFNmzbcfffdjB49+pJrkuKlMPKHqPAaHNhSjdrGcTi6Hep0s7okEZHi5ciEVyJK/7hPHQIv/yJtMn/+fAICAsjLyyMnJwebzVboN/gmJCRQs2ZN+vXrh6enJ1FRUXTu3BmA0NBQ7HY7gYGB1KxZs8inIiVDt2n+0PCMYeH1jhoREWv16dOHmJgY1qxZw4gRI7jjjju48cYbC7XtzTffTFZWFvXr12f06NHMnTuXvLy8Eq5YLodaRv7QMCyAb8xI+hKD80gsdqsLEhEpbp5+7lYKK45bRP7+/jRs2BCAqVOn0qZNG6ZMmcKoUaOw2WyYpllgfYfDkf/nyMhIduzYwS+//MLChQu5//77+c9//sOSJUvw9NQjCmWRwsgfagb5sN/uvh+Zc2grRf+/johIGWcYRb5dUhbYbDaeeuopxo4dy9ChQ6levTpJSUn5y1NTU4mPjy+wja+vL4MHD2bw4MFER0fTtGlTtmzZQvv27fHy8sLpdJb2acgF6DbNHwzDIOv0O2qOx7k7eomISJlw8803Y7fbmThxIn379uXTTz9l2bJlbNmyhREjRmC3/9mePX36dKZMmcLWrVvZu3cvn332Gb6+vvkdYOvWrcvSpUs5ePAgx48ft+qU5AxqGTmDd82m5J2y4ZWbAmlJEGRBRy8RETmLh4cHY8aM4fXXX2fXrl3Ex8dz3XXXERwczP/93/8VaBkJCQnh1VdfZezYsTidTlq1asV3331H1apVAXjxxRe55557aNCgATk5OWfd8pHSZ5jl4L9CamoqwcHBpKSkEBQUVGLHmbR4D1f9OpjGtoMw9Cto3L/EjiUiUpKys7OJj4+nXr16+Pj4WF2OVGAX+q4V9vdbt2nO0CgsgFjzj+fYD/9ubTEiIiKVhMLIGRqGBRDrcocR1+GtFlcjIiJSOSiMnCEy1I/dtroAOA9ttrYYERGRSkJh5Ax2m0F2aHMAPJLjITfD4opEREQqPoWRvwgLj+SoGYKBCUc0EquIiEhJUxj5i0Y1AvP7jagTq4iISMlTGPmLAk/UHFEnVhERkZKmMPIXjWoEEueKAsBM2mJxNSIiIhWfwshfRIX6sdtWDwDzyFZw6f0FIlK5ZToyafVxK1p93IpMR6bV5UgFpDDyF3abgb1aQ7JML2x5WXAy/uIbiYhIuTd9+nRCQkKsLqPELV68GMMwSE5OLvQ2devW5a233iqxmhRGzqFBzWB2mJHuD+rEKiJSqkaOHIlhGBiGgZeXFw0bNuTFF18kLy+vRI976623snPnzhI9xsWcPvd77733rGXR0dEYhsHIkSNLv7ASpjByDo3PfKJGnVhFpJJ6L+Y9Jm+efM5lkzdP5r2Y90rs2Ndccw1JSUns2rWLRx99lOeff57//Oc/51w3Nze3WI7p6+tLWFhYsezrckRGRjJz5kyysrLy52VnZzNjxgyioqIsrKzkKIycQ8MC76hRJ1YRqZxsho2JMROZsmVKgfmTN09mYsxEbEbJ/YR4e3tTs2ZN6tSpw3333Ue/fv349ttvAXfrwZAhQ3j55ZeJiIigSZMmACQmJnLLLbcQEhJCaGgo119/Pfv27QPg559/xsfH56xbEw899BB9+/YFzn2bZtKkSTRo0AAvLy+aNGnCp59+mr9s3759GIZBTExM/rzk5GQMw2Dx4sUAnDp1imHDhlG9enV8fX1p1KgR06ZNu+C5t2/fnsjISObMmZM/b86cOURFRdGuXbsC6+bk5PDggw8SFhaGj48PPXv2ZN26dQXWWbBgAY0bN8bX15c+ffrkX5MzLV++nCuuuAJfX18iIyN58MEHycgovYE/FUbO4cyWEVPvqBGRCsI0TTIdmYWehjcfzt2t7uaDLR/k7+P9ze8zMWYid7e6m+HNhxdqP8XxcnhfX98CLSCLFi1ix44dLFy4kPnz5+NwOBgwYACBgYEsW7aMFStWEBAQwDXXXENubi5XXXUVISEhzJ49O38fTqeTWbNmMWzYsHMec+7cuTz00EM8+uijbN26lXvuuYc77riD3377rdB1P/PMM8TGxvLDDz8QFxfHpEmTqFat2kW3u/POOwuElqlTp3LHHXectd7jjz/O7Nmz+fjjj9m4cSMNGzZkwIABnDx5EnAHtBtuuIHBgwcTExPDXXfdxZNPPllgH3v27OGaa67hxhtv5Pfff2fWrFksX76cMWPGFPo8L5dHqR2pHIkK9SPeXgeXaWBLOwQZx8H/4l8eEZGyLCsviy4zulzWPqZumwrAB1s+KBBSLmTN0DX4efpd0vFM02TRokX89NNPPPDAA/nz/f39+eijj/Dy8gLgs88+w+Vy8dFHH2EYBgDTpk0jJCSExYsX079/f/7xj38wY8YMRo0aBbgDTXJyMjfeeOM5jz1hwgRGjhzJ/fffD8DYsWNZvXo1EyZMoE+fPoWqPyEhgXbt2tGxY0fA3RG0MG677TbGjRvH/v37AVixYgUzZ87Mb3EByMjIYNKkSUyfPp2BAwcC8OGHH7Jw4UKmTJnCY489lt+y88YbbwDQpEkTtmzZwmuvvZa/n/HjxzNs2DAefvhhABo1asQ777xD7969mTRpEj4+PoWq+XIojJyD3WZQs3p19p8Io55xxH2rpkHhvngiInL55s+fT0BAAA6HA5fLxdChQ3n++efzl7dq1So/iABs3ryZ3bt3ExgYWGA/2dnZ7NmzB4Bhw4bRtWtXDh06REREBJ9//jmDBg067xM0cXFx3H333QXm9ejRg7fffrvQ53Hfffdx4403snHjRvr378+QIUPo3r37RberXr06gwYNYvr06ZimyaBBg85qUdmzZw8Oh4MePXrkz/P09KRz587ExcXln0OXLgUDaLdu3Qp83rx5M7///juff/55/jzTNHG5XMTHx9OsWbNCn++lUhg5j8Y1Aog9Xod6HHF3YlUYEZFyztfDlzVD1xR5u/c3v5/fIgJwd6u7GdVqVJGOW1R9+vRh0qRJeHl5ERERgYdHwZ8rf3//Ap/T09Pp0KFDgR/U06pXrw5Ap06daNCgATNnzuS+++5j7ty5TJ8+vci1nWazuXs6nHkbyuFwFFhn4MCB7N+/nwULFrBw4UKuuuoqoqOjmTBhwkX3f+edd+bfKpk4ceIl13kx6enp3HPPPTz44INnLSutDrMKI+fRKCyAuC11GGRfq06sIlIhGIZR5NslkzdPPiuIfLDlAzztntzb5uzHT4uLv78/DRs2LPT67du3Z9asWYSFhREUFHTe9YYNG8bnn39O7dq1sdlsDBo06LzrNmvWjBUrVjBixIj8eStWrKB5c/fb3U+HnKSkpPyOpWd2Zj2tevXqjBgxghEjRnDFFVfw2GOPFSqMnO7vYhgGAwYMOGv56Y61K1asoE4ddz9Hh8PBunXr8m+5NGvWLL/j72mrV68u8Ll9+/bExsYW6XoXN3VgPY9GNQLPeKJGnVhFpPI5/dTM3a3+vFUxqtUoottGMzFm4nkf+7XCsGHDqFatGtdffz3Lli0jPj6exYsX8+CDD3LgwIEC623cuJGXX36Zm266CW9v7/Pu87HHHmP69OlMmjSJXbt28eabbzJnzhz+9a9/Ae5OtV27duXVV18lLi6OJUuW8O9//7vAPp599lm++eYbdu/ezbZt25g/f36hb3vY7Xbi4uKIjY3Fbreftdzf35/77ruPxx57jB9//JHY2FhGjx5NZmZmfr+Ye++9l127dvHYY4+xY8cOZsyYcVZr0BNPPMHKlSsZM2YMMTEx7Nq1i2+++aZUO7AqjJxHgSdqju+AvByLKxIRKV0u00V02+izbsnc2+ZeottG4zJdFlV2Nj8/P5YuXUpUVBQ33HADzZo1Y9SoUWRnZxdoKWnYsCGdO3fm999/P+9TNKcNGTKEt99+mwkTJtCiRQvef/99pk2bxpVXXpm/ztSpU8nLy6NDhw48/PDDvPTSSwX24eXlxbhx42jdujW9evXCbrczc+bMQp9XUFDQBVt6Xn31VW688UZuv/122rdvz+7du/npp5+oUqUK4L7NMnv2bObNm0ebNm2YPHkyr7zySoF9tG7dmiVLlrBz506uuOIK2rVrx7PPPktERESh67xchlkcz1yVsNTUVIKDg0lJSbngf5Ti5HSZNHv2B9bYR1PFSId7lkJ4m1I5tojI5crOziY+Pp569epd9tMQmY7M/KdwLufJGKmYLvRdK+zvt1pGzsNuM2hQ/YyRWNVvREREpEQojFxAozNHYj240dpiREQs4ufpx5YRW9gyYotaRaREKIxcQOMaAax3NXZ/SFhlbTEiIiIVlMLIBTSqEcg6V1P3h6OxkHHC2oJEREQqIIWRC2gUFsBJgthl1nbPSFhpbUEiIiIVkMLIBUSF+uHlYWO184/Wkf0KIyIiIsVNYeQCPOw26lfzZ43rjwFq9i23tiAREZEKSGHkIhrXCGTN6X4jh7dAVrKl9YiIiFQ0CiMX0aFOFY5RhQO2CMCEhNUX3UZEpELJyADDcE8ZGVZXIxWQwshFDGlbCx9PG8tym7hn7F9hbUEiIpXc4sWLMQyD5ORkq0uRYqIwchHBfp4Mbh3xZ78RhRERqWyczj//vHRpwc/FzDCMC07PP/98iR1brKMwUgjDutZh7R/9RsxDMZCTZm1BIiKlZc4caN78z8/XXgt167rnl4CkpKT86a233iIoKKjAvNNvzLVCbm6uZceu6BRGCqFN7WBCa9Un0VUdw3RC4hqrSxIRKXlz5sBNN8HBgwXnHzzonl8CgaRmzZr5U3BwMIZhFJgXEBCQv+6GDRvo2LEjfn5+dO/enR07dhTY1zfffEP79u3x8fGhfv36vPDCC+Tl5eUvT0hI4PrrrycgIICgoCBuueUWjhw5kr/8+eefp23btnz00Uf5L4H75JNPqFq1Kjk5Bd/kPmTIEG6//fZivx6VhcJIIRiGwW1d6rDGdN+qMfdpvBERKYdM090BtTBTaio8+KB7m3PtB+Chh9zrXWxfJfRy+Keffpo33niD9evX4+HhwZ133pm/bNmyZQwfPpyHHnqI2NhY3n//faZPn87LL78MgMvl4vrrr+fkyZMsWbKEhQsXsnfvXm699dYCx9i9ezezZ89mzpw5xMTEcPPNN+N0Ovn222/z1zl69Cjff/99geNLEZnlQEpKigmYKSkpltWQkeMw//3s46b5XJCZ/G4fy+oQESmMrKwsMzY21szKyvpzZnq6abqjQelO6emXdA7Tpk0zg4ODz5r/22+/mYD5yy+/5M/7/vvvTSD/fK+66irzlVdeKbDdp59+aoaHh5umaZo///yzabfbzYSEhPzl27ZtMwFz7dq1pmma5nPPPWd6enqaR48eLbCf++67zxw4cGD+5zfeeMOsX7++6XK5Luk8y7tzftf+UNjfb7WMFJKflwdVW/QBwP/YZsjNtLgiEZHKrXXr1vl/Dg8PB9ytFACbN2/mxRdfJCAgIH8aPXo0SUlJZGZmEhcXR2RkJJGRkfn7aN68OSEhIcTFxeXPq1OnDtWrVy9w3NGjR/Pzzz9z8I/bV9OnT2fkyJEYhlFi51rReVhdQHlyba/uJG0LJdw4yckdKwhtdbXVJYmIFJ6fH6SnF27dpUvdnVUvZsEC6NXr4sctAZ6envl/Ph0EXC4XAOnp6bzwwgvccMMNZ23n4+NT6GP4+/ufNa9du3a0adOGTz75hP79+7Nt2za+//77opYvZ1AYKYLGNYNY5tuG8Ozf2LH2J7opjIhIeWIYcI4f13Pq3x9q13Z3Vj1Xnw/DcC/v3x/s9uKtsxi0b9+eHTt20LBhw3Mub9asGYmJiSQmJua3jsTGxpKcnEzzM58eOo+77rqLt956i4MHD9KvX78CLSxSdLpNU0RBTXsD4HVgJXlOl8XViIiUELsd3n7b/ee/3n44/fmtt8pkEAF49tln+eSTT3jhhRfYtm0bcXFxzJw5k3//+98A9OvXj1atWjFs2DA2btzI2rVrGT58OL1796Zjx44X3f/QoUM5cOAAH374oTquFgOFkSJq1vUaAFq4drJ+z2GLqxERKUE33ABffw0REQXn167tnn+OWyBlxYABA5g/fz4///wznTp1omvXrvz3v/+lTp06gPu2zjfffEOVKlXo1asX/fr1o379+syaNatQ+w8ODubGG28kICCAIUOGlOCZVA6GaZbQM1fFKDU1leDgYFJSUggKCrK2GNMk5eUGBOed4LNGb3PbsJHW1iMicg7Z2dnEx8fnj49xWVJTITjY/ecFC8rsrZnSdtVVV9GiRQveeecdq0ux1IW+a4X9/S5Sy8j48ePp1KkTgYGBhIWFMWTIkLMGmTmXr776iqZNm+Lj40OrVq1YsGBBUQ5bthgGqbWvBMA7fiHlIMuJiFyeM4NHr16VPoicOnWKuXPnsnjxYqKjo60up0IoUhhZsmQJ0dHRrF69moULF+JwOOjfvz8ZF3iL48qVK/nnP//JqFGj2LRpE0OGDGHIkCFs3br1sou3SrX2gwHomLuO3UcL2TNdRKS88vf/c9SQwnaArcDatWvHyJEjee2112jSpInV5VQIl3Wb5tixY4SFhbFkyRJ6nefRrltvvZWMjAzmz5+fP69r1660bduWyZMnF+o4Zeo2DUBOGnnj6+JBHp93nsuwa/taXZGISAHFeptG5AJK/TbNX6WkpAAQGhp63nVWrVpFv379CswbMGAAq1atOu82OTk5pKamFpjKFO9AjlXtAEBO7A8WFyMiIlK+XXIYcblcPPzww/To0YOWLVued73Dhw9To0aNAvNq1KjB4cPnfxJl/PjxBAcH509l8flt/5buwYAap67kWFrORdYWERGR87nkMBIdHc3WrVuZOXNmcdYDwLhx40hJScmfEhMTi/0Ylyuo1XUAdDbiWLY13uJqRETOTZ3spaQVx3fsksLImDFjmD9/Pr/99hu1a9e+4Lo1a9Ys8EpmgCNHjlCzZs3zbuPt7U1QUFCBqcyp2oBkn9p4GU6SYn60uhoRkQJOD5Weman3aEnJOv0dO3N4/qIq0nDwpmnywAMP5D/SVK9evYtu061bNxYtWsTDDz+cP2/hwoV069atyMWWKYaBs+HVsHUa1ZOWkO14AB/Pyv24m4iUHXa7nZCQkPwXx/n5+elFblKsTNMkMzOTo0ePEhISgv0yHvkuUhiJjo5mxowZfPPNNwQGBub3+wgODsbX1xeA4cOHU6tWLcaPHw/AQw89RO/evXnjjTcYNGgQM2fOZP369XzwwQeXXHRZEdp2MGydRi9jEyt2HeOq5udv7RERKW2nW6BPBxKRkhASEnLBux2FUaQwMmnSJACuvPLKAvOnTZvGyJEjAUhISMBm+/PuT/fu3ZkxYwb//ve/eeqpp2jUqBHz5s27YKfX8sKo25Mcmy81Xaf4ctNyrmp+k9UliYjkMwyD8PBwwsLCcDgcVpcjFZCnp+dltYicVuTbNBezePHis+bdfPPN3HzzzUU5VPng4U1aRA+8D/yCx55fcLluxGZTM6iIlC12u71YfjBESopelHeZglsPAqBL3np+P5hicTUiIiLlj8LIZfJsMgCAdsZulsVst7gaERGR8kdh5HIF1yIluCk2w+TE5gXkOV1WVyQiIlKuKIwUg9OjsfbMWcpvO45ZXI2IiEj5ojBSDDzaDcPEoJ99E0uWLba6HBERkXJFYaQ4VGtIZkP38PAdD3xM4kmNeCgiIlJYCiPFxP+qxwAYbFvJD8vO/0ZiERERKUhhpLiEt+FYjSuwGyahMZNxqCOriIhIoSiMFKMq1zwJwGDXryzdsMXiakRERMoHhZFi5FGvJwcD2+Bt5JG19B2ryxERESkXFEaKmXdfd9+RK9Pmsy8x0eJqREREyj6FkWJWre11JHg2IMDIZv8Pb1ldjoiISJmnMFLcDIOUjmMAaHNoJtkZel+NiIjIhSiMlIDmVw3nADUJIZ3YXz61uhwREZEyTWGkBNg9PNhX5wYAfLbNsrgaERGRsk1hpIQ0uGoULtOgee7vHIqPs7ocERGRMkthpISERzVkq297APYv+sjiakRERMouhZES5Gj1TwDqHvgGp9NpcTUiIiJlk8JICWrZ95+k4Uc4x9iy4nuryxERESmTFEZKkLdvADuq9QcgZ90nFlcjIiJSNimMlLCqPe8AoHXqEk6eOG5xNSIiImWPwkgJq9emN4n22vgauWz9Ra0jIiIif6UwUtIMg2MNbgKgys6vME3T4oJERETKFoWRUtDo6lE4TYNWzlhit8VYXY6IiEiZojBSCgKrR7EzoDMAhxZPsbgaERGRskVhpJR4dRoOQNtj37IsNtHiakRERMoOhZFS0uCKW0n2qkl1I4UVX7/NifQcq0sSEREpExRGSovdE/8rHwFgaN48npodo86sIiIiKIyUKs+Ow8nzqUqU7RjeO75l1jrdrhEREVEYKU1efnh0vw+A+zy+5YXvthF/PMPiokRERKylMFLaOt2F6RVAM1sCXZwbeHjmJhxOl9VViYiIWEZhpLT5VsHo6B4i/gGv+Ww+kMLcTQctLkpERMQ6CiNW6BoNdi86EEcHYwefrNqnzqwiIlJpKYxYISgc2vwTgGjP79h6MJVNicnW1iQiImIRhRGr9HgIMOhr20hTI4FPVu6zuiIRERFLKIxYpWoDaDEEgIc9ZrNgy2GOpWkgNBERqXwURqzU+0nA4Br7Ohq7djNrXYLVFYmIiJQ6hRErhTWF1rcAMNbjaz5fk0CeHvMVEZFKRmHEar2fwDTs9LXHEJ76O7/EHbG6IhERkVKlMGK1qg0w2rqfrBnr8RWfrNpvcUEiIiKlS2GkLOj1OKbNk572bbjil7HrSJrVFYmIiJQahZGyoEodjPbDAXfryKer9llbj4iISClSGCkrev0Ll82LzrYdJG1cQFq2w+qKRERESoXCSFkRFIHRaRQA0cxi7sYDFhckIiJSOhRGyhDjirHk2Xxoa9tD7PJ5el+NiIhUCgojZUlAGM4OIwG4KX0Gq/Yct7YeERGRUqAwUsZ493oEh+FFR9tO1v46z+pyRERESpzCSFkTWJP0FkMB6HpgKkkpWRYXJCIiUrIURsqgKlc/Rh4edLXFsmThN1aXIyIiUqIURsqi4NocrHcjAHW3vUdunt5XIyIiFZfCSBkVcd1T5GGnq7mZ1Ut/tLocERGREqMwUkZ5Vq3L9hqDAPBf/WaBZdkOpx77FRGRCsPD6gLk/MIHPU3elPl0yF3HA69PZrWjISlZDnLzXHSsU4Uv7+mGzWZYXaaIiMhlUctIGVY1qikbQvoD8ETGf8hLO5bff2T9/lOsiT9pZXkiIiLFQmGkjGt553tkBtShtnGcZQ0+Y/ljV3BLx9oAfLk+0eLqRERELp/CSBnnH1wVv9tngqc/AQeXU3vjBIZ2qQPAgi1JpGTphXoiIlK+KYyUBzWaw/Xvuv+84m3apPxKkxqB5OS5+G7zIWtrExERuUwKI+VFyxug+4MAGN9Ec3dT98isulUjIiLlncJIeXLVc1D/SnBkMmTnkwTZc/j9QApxSalWVyYiInLJFEbKE7sH3DQNgmphP7WX/1afD8CsdWodERGR8kthpLzxC4XB7wDQN3kOHY3tzIs5SE6e0+LCRERELk2Rw8jSpUsZPHgwERERGIbBvHnzLrrNxIkTadasGb6+vjRp0oRPPvnkUmqV0xr1g3a3YWDyps+HZGemszD2iNVViYiIXJIih5GMjAzatGnDxIkTC7X+pEmTGDduHM8//zzbtm3jhRdeIDo6mu+++67IxcoZ+r8MgeFEmUk86vGVbtWIiEi5VeTh4AcOHMjAgQMLvf6nn37KPffcw6233gpA/fr1WbduHa+99hqDBw8u6uHlNN8QGPw2zLiFUfYf+HFPZw4mt6ZWiK/VlYmIiBRJifcZycnJwcfHp8A8X19f1q5di8OhAbsuS+MB0GYoNsPkdY/3+WLFDqsrEhERKbISDyMDBgzgo48+YsOGDZimyfr16/noo49wOBwcP378nNvk5OSQmppaYJLzuOYVcnzCaGBLImjNGySezLS6IhERkSIp8TDyzDPPMHDgQLp27YqnpyfXX389I0aMcB/cdu7Djx8/nuDg4PwpMjKypMssv3yr4DXkLQDuML7nk29+sLYeERGRIirxMOLr68vUqVPJzMxk3759JCQkULduXQIDA6levfo5txk3bhwpKSn5U2KiOmdeiNF0EGl1B+BpOOkf/yrr48/d4iQiIlIWldo4I56entSuXRu73c7MmTO57rrrztsy4u3tTVBQUIFJLizw72+SY/Olk20nq2a/jctlWl2SiIhIoRQ5jKSnpxMTE0NMTAwA8fHxxMTEkJCQALhbNYYPH56//s6dO/nss8/YtWsXa9eu5R//+Adbt27llVdeKZ4zELfg2jiueBKA29Km8v3qLRYXJCIiUjhFDiPr16+nXbt2tGvXDoCxY8fSrl07nn32WQCSkpLygwmA0+nkjTfeoE2bNlx99dVkZ2ezcuVK6tatWzxnIPkCeo3heEBjqhjpGL88Q0ZOntUliYiIXJRhmmaZb89PTU0lODiYlJQU3bK5iNz9a/CYNgAbJu/VeRvq9ST+WAb7TmRwNC2HR/s34W9tIqwuU0REKoHC/n7r3TQVjFedLhyo/w8Abtr3HAG/PEnuppkk7d/B/hMZvPbDdpzqTyIiImWIWkYqIDPrFMn/7UaV3KQC849ShSdzR/GPYaPp36KmRdWJiEhloZaRSszwrUKVsWvglk+gazTU6gg2T8I4xRuek5mzXJ1bRUSk7Cjyu2mknPAJhubXuyeA3Exy3+9LlRNx9EycxO6jXWkYFmhtjSIiIqhlpPLw8sNr8BsADLX/yi+LfrK4IBERETeFkcqkbg+O1f0bNsOk6/ZXScvKsboiERERhZHKptoNr5GJD22NXWyeP9nqckRERBRGKhsjKIK4xvcB0Dz2DcysUxZXJCIilZ3CSCXUZMjj7DVrEWqmcGje81aXIyIilZzCSCUU4OfHsob/AiB8xyew+FVwOiyuSkREKiuFkUqqx4Bb+DSvHzZcsHg8x9+6gsO7NlpdloiIVEIaZ6SSahgWwKcdn+fBtU150XMa1dLiyPnsaqb4DmNZ9X+S6zJwOF3kOk2q+Hky4eY2VAvwtrpsERGpgDQcfCWXeDKT5Zu20njtv+mQswaAOFckE/JuYZGrPWAAcFOH2ky4uY2FlYqISHlT2N9vhRFxM00y1n6C5y//xsuRCsCp0LZsbPQgo5b4ADAvugdtI0MsLFJERMoTvZtGisYw8O8yAq9HNkPPR8DDlyonY7hqzZ38UO1tAsnk+W+34dIbf0VEpJgpjEhBfqHQ73l4KAY63QU2D5qlr+EZ7y+ISUxmXsxBqysUEZEKRmFEzi2wJgx6A4Z/A8AtxiI6Gtt59YftpOfkWVyciIhUJAojcmF1e0L74QD8x2cap9IymPjbbouLEhGRikRhRC6u3wvgV416ZiJ32+czZVk8+45nWF2ViIhUEAojcnF+oXDNeAAe8pxLuOsQT8z+nQzdrhERkWKgMCKF0+pmqN8HLxyM95rKmvgT3Dx5FYdTsq2uTEREyjmFESkcw4Dr3gQPH7obW7nP71fikpIZMnEF2w6lWF2diIiUYwojUnih9aHXYwA84ZrCCt9/cVPGF4yZPJ9ftx+xuDgRESmvNAKrFI3TAQufg02fQo57pFanabDE1Ybv/f9OTuQVtKodQqvawbSNDMHPS68/EhGprDQcvJSs3EyI/QbXxk+wJazMn/27qx6T8v7GT65OBPt5M6ZvI27rGoW3h93CYkVExAoKI1J6ju8me+UkPDd/ht3p7tC6nwhezb2ZH1xdqF3Fl8cGNGFw6whsNsPiYkVEpLQojEjpyzgOa96Hte9DtrtT66v2e5ic0RuAFhFBTBzanrrV/K2sUkRESolelCelz78a9H0aHtkGne8G4Enn+0xtHUuAtwfbDqXy2NebKQf5V0RESpHCiBQ/70AY+Dp0uQ+AvjtfZunVB/H1tLNu3ynm/55kcYEiIlKWKIxIyTAM96itne8GTEJ/GcvbzeIAGL8gjqxcp7X1iYhImaHnLqXkGIa7hcTlhPVTuHrnC3zt15J1GXVZNGc31w28DoIi3OuJiEilpQ6sUvJcLljwKKyfevayao2h9xPQ4gawqaFORKQi0dM0UvYcicU8sI5Fv/xAeEYsTW2J2HG5l9VoCX2fgcYD1FIiIlJBKIxImbX1YAqD312Ov5nJj91iqR37Uf5orkR2gf4vQ2Qna4sUEZHLpkd7pcxqWSuYf3SKJB0/7tnfl6z7NkKPh8DDFxLXwJSr4buHIPPkObc/mpaNw+kq5apFRKSkKIyIJR7t34RAH/fYI7d/sYuUns/Ag5ug7TDAhA3T4d2OEDMDzmi8m7PxAF1fWcSYGRs1XomISAWhMCKWqBbgzfQ7OhPk48H6/af45werOW4LhSHvwR0/QPVmkHkC5t0H066Fw1tYtecET8z+HZcJP207wk/bDlt9GiIiUgwURsQyHepUYebd3agW4EVsUiq3TF7FoeQsqNMd7l0GV78Inn6QsBLz/V4kfHov/s5UagR5A/DCd7Gk5+RZfBYiInK5FEbEUs0jgvjynm7UCvFl7/EMbp68iuW7juM0PNz9SKLXktN0CIbp4lZ+ZpnvoyzutYt6VbxJSsnmrYU7rT4FERG5THqaRsqEQ8lZ3DZlDXuPZQBQM8iH69tFcF2rCF74bhv2hJW87PMpDc19AGQG1ueJE9fyA934dkwvmkfoeyEiUtbo0V4pd06k5/Dmwp18t/kQqdkFb78E+ngw997ONEycDb++DFnuJ222uyL5JmQ4jz30L2x2NfSJiJQlCiNSbuXkOflt+1HmbjrIr9uPYmAw7Y5O9GhYzb1CdiqseR/Xynew/TE+yYngFlQd9bV7eHkRESkTFEakQkjJcpCT5yQs0OfshVmn2PTlyzTa+ykBRjZH/Juwb/DXtGlQCx9Pe+kXKyIiBWjQM6kQgn09zx1EAHyr0Oq213kk9F2Om0HUyNhByucjafv8j9zy/iqW7jxWusWKiMglURiRcs3DbuO/9/6dDd3exWF40t++gYeNGayNP8m9n23gRHqO1SWKiMhFKIxIuRfg7cGAa67H8++TALjXYz6PhK4iM9fJ5CV7LK5OREQuxsPqAkSKTeub4cRuWPIqD2RNItmezN7VtTnRNIuqNaLAryrYlL9FRMoadWCVisU0YfZdsPXrs5cZNvAJ/mMKcf9vne7Q/UHw8iv1UkVEKrrC/n6rZUQqFsOA6ydClTok715L0sF9hBkpVDVSwXRB1in3dFr8Eoj5HK6dAI0HWFe3iEglppYRqdCGfbSaFbtPMLRjOK8MiICsZMhOITv9BOs2b6VDwlT8spLcKzcbDNe8BsG18rfffjiVHIeLNpEhltQvIlKe6dFeEWDs1U0AmLXxMPtyAiGsKYsy6nLVtz7cvrkFHU69zP6md4Fhh7jvYGJn2PAxmCYbE04x+H/LuX7iCl7/cTtOV+Fye1q2g/eX7GFjwqmLrywiImoZkYrvjmlr+W3HMfo1q4HdBj9tOwKAj6eNbIeLIB8Pfh5ajZrLxkHiGgCymwzh2r03sjftz8HTejWuzjv/aEuIn9d5j7XrSBr3fLaBvccyMAy4q2c9Hu3fRIOwiUilpJYRkT+cbh35Je4IP207gofN4N7eDVgzrh9tIkNIzc5j9E9Z5Az/Hvq9gGnY8dkxj6k5/2Jg1cO8fmNrfDxtLN15jL+9u4K4pNRzHmfBliSun7iCvccyCPTxwDThw2XxDHpnGZsTk0vxjEVEyhe1jEil8OAXm/h28yE61qnCS39vSdOa7u/RweQsBr2zjORMB7d1jeKlIa34cMZMBu54mtrGcUybJ0a7YZzMzGPdzkTM3Ew8bCb7Im/A2egamtQMpHGNQD5euY/3l+4FoFv9qvxvaDs2Jybz5JwtHEvLwW4zuP/KBjzcrzF2m2HlpRARKTV6N43IGRxOFzuPpNGsZhC2v4SBxTuOcsf0dZgmDGkbwbyYQwSRzsIGX1Hj4MLz7nNa3gDG5w0lF8/8eff0qs9jA5rg8ccbhE9l5PLst9v4bvMhwH3b5t/XNS+BMxQRKXsURkSK4M2fd/DOr7vzPz90VSMe6dcIts6GpM3g5Q+efrg8/Tm0ayO1d30GwE5bA+7Nvp9jXpG8ekNrBrUOP+f+v1yXyOOzfwfgpSEtua1rnZI/KRERiymMiBSB02UyYupalu8+zlVNw/hweMezWlAK2PkTzL0Xsk5iegXg6joGuyMdTu1zT+lHoOVN0P8lsLuH8/nfol28sXAndpvBtJGd6NW4eqmcm4iIVRRGRIoo2+Fk5Z7j9GhYDW+PQjz9knIQ5oyG/SvOv07ja+CmqeDlj2maPPrVZuZsPEigtwez7+9O4xqBxXcCIiJljMKISGlw5sGaSZC4FkKioEpd95R5Ar57CPKyoVYHGPol+FcjJ8/J7VPWsjb+JLVCfJkX3YPqgd5Wn4WISIlQGBGxWsIa+OJW9/DzofXhttkQWp9TGbn8/b0V7DuRSbuoEGbe3bVwLTEiIuWMxhkRsVpUFxi10N1icnIvTOkPR7dTxd+LqSM7EeTjwaaEZP49dyvl4N8EIiIlRmFEpCRVawSjfoGarSHjGHx2I6Qeon71AN4d2h6bAV9tOMD0lfusrlRExDJFDiNLly5l8ODBREREYBgG8+bNu+g2n3/+OW3atMHPz4/w8HDuvPNOTpw4cSn1ipQ/gTVg+DdQtRGkHoDPboKsZHo1rs5T1zYD4KXv41ix+7jFhYqIWKPIYSQjI4M2bdowceLEQq2/YsUKhg8fzqhRo9i2bRtfffUVa9euZfTo0UUuVqTc8gt19xkJqAFHt8HMYeDIZlTPetzQrhaGy8Gszz8g+cdXIENBXUQql8vqwGoYBnPnzmXIkCHnXWfChAlMmjSJPXv25M/73//+x2uvvcaBAwcKdRx1YJUKI+l3mHYt5KZB8+uh9xPkbfyctLWfU8VMBsDhF4bHjR9iNLjS0lJFRC5XmenA2q1bNxITE1mwYAGmaXLkyBG+/vprrr322vNuk5OTQ2pqaoFJpEIIbw3/+BxsnhD7DUzqjseaiVQxkzlBMPtdYXhmHsX8dAgrJo1h3e7DOF3q3CoiFZtHSR+gR48efP7559x6661kZ2eTl5fH4MGDL3ibZ/z48bzwwgslXZqINer3hr9Phtl3gc3uHhit7VBOBnXjvV/j6LbrDW4xFtHjyKds/mQFg5zR5IY2pE6oH3Wq+tMwLIC/t6uFv/fF/+/rcpnsOprO2vgTeHnYuKVjJIahF/WJSNlS4rdpYmNj6devH4888ggDBgwgKSmJxx57jE6dOjFlypRzbpOTk0NOTk7+59TUVCIjI3WbRiqWk3vBOwj8qxWYne1wsv3Xz2m85in8XGm4TINlrlbMcF7FIlc78vCgde1gpt/RmVB/r7N2m5GTx6x1iazae4J1+06SnOnIX/bMdc0Z1bNeiZ+aiAiU0qBnhQkjt99+O9nZ2Xz11Vf585YvX84VV1zBoUOHCA8/94vFzqQ+I1IppRzAnD8WY9dP+bPSPKsxM683k7KuJjQsgk9HdSY82Dd/eUxiMg/N3MT+E5n583w97TSqEcDvB1LwsBnMuqcrHeqEluqpiEjlVGb6jGRmZmKzFTyM3e4ebVIDPYlcQHBtjGFfwoOboMfD4F+dQMdxRpuzWeIzlt4nZvGP95YSfzwDp8vk3V93ceOklew/kUlEsA9PDmzK3Pu78/vz/fkmugfXtQ4nz2UyZsYmTqTnXPTwIiKlpcgtI+np6eze7X7Vert27XjzzTfp06cPoaGhREVFMW7cOA4ePMgnn3wCwPTp0xk9ejTvvPNO/m2ahx9+GJvNxpo1awp1TLWMiAB5ubDje1j2BhzeAkC8qwbveowgofqVrNufDMCg1uG8MqQVwX6eBTZPz8njb+8u58CxZO6IPMwTbXKw1e8NEW1L+UREpLIosds0ixcvpk+fPmfNHzFiBNOnT2fkyJHs27ePxYsX5y/73//+x+TJk4mPjyckJIS+ffvy2muvUatWrWI9GZFKweWEmBm4fnkBW+YxAOJcUcQaDWjQsjNt2nfDqN7EvW5uxp/T0VjSt/2EsW8p/sYfLSOGHfqMg55j3Z1pRUSKkV6UJ1LR5aSRvfgN7KvfxdN0XHz9Mxwzg9ln1qSTbYd7Rt0r4O/vQ/C5/4GQletk/f6TrNh9gqOp2Tx+TVNqBvtc7hmISAWnMCJSWaQdhoTVcDTWPR2JdT+pY9jAKwC8/MHLDwLDoUFfaNiPp1a6mLHuADfYlvGi53QCjGwy7UFs6/gy+8P6kpGTR3pOHqnZDjYlJLMp4RQO559/VbSpHcyse7rh46nWFBE5P4URkcrM5XSHkfOMKZKT52TK8ngWxh4hOXE7b3m+SxvbXgD+47iFic7rgYLbhgf70K1BVX7dfpTkTAe3dKzNaze21rglInJeCiMiUijJmbms2HmIwGUv0+vElwAsDh7CgloP4efjTaMaAfRoUI06Vf0wDINlu44xYupaXCa8NKQlt3WtY/EZiEhZpTAiIkW3ejL8+CRgQvMhcMMH4OENzjzY9ROsnwYH13PUI5ylydXYQ23+PuBqGkfVgrwscGSBIxNME+r0gKCLjyMkIhWXwoiIXJqts2HOPeByQL1eENUdNn4CaYeKvq9aHaHZddB0MFRrWPy1ikiZpjAiIpduz28w6zbITf9znl9VaDsMmv0N0g6RmxTLqtXLqZmzD397Hj6+AQQGBuHt6+/e7tCmgvts2A/+ORPsBcc/EZGKS2FERC7PoU3uFpKAMOgwEpoNdt+yOcO+4xncNHklx9Nz8+e1jwrh7+1q8bcGNoL3L4Tt8yF+KbjyYNAb0OmuUj4REbGKwoiIlIq0bAc/bTvCvE0HWbHnOKf/RvGy2+jfoga3doqkx8m52H54DPzD3MPbewcUat8ZOXn4edn/fGIn/aj78eW6V4DNTk6eky/WJNC4ZiDdG1S78M5EpNQpjIhIqTuSms13mw/x9YYDbD+clj+/TrAH3xhjCck+AFeOgyufPHvjvYvh1H5odzumYTBleTyv/rCdtpEh/PfWtkR6psKHV0HqAajakJROD3HH+rpsPOA+zrAuUTw9qBl+Xh6ldLYicjEKIyJiGdM02XYolS/XJzJv00FSs/MYZFvNRK93yDZ8+LX/T/Tp0Apfrz8GTYv7Dr4cDqYLV9vbeJF7mL46MX9/Yd55/FzlNUKStxU4zn5XGB/yd2bl9sCBB3Wr+vHGLW3pUKdKaZ6uiJyHwoiIlAnZDic/bTvM1+sT+VfCfbSx7eXjvKuZYB/NdW3CGVnrII1/HoHh/PNNwl/m9eaJvNE8cFUTVu48zOjDLzDAvp50ezCLu3xE3LI53GH7nmpGKgCZQQ14OGsUP6fVxWbA3b0acFOH2jSo7q9B2UQspDAiImXOsS2/UH32jeRhp1/O63jj4CuvFwkyMtlbtRdzcjrzSNob2A2TA3X+Tu0RU3D9/Cy21e+SY3ryz9yn2Wg2BmBwsyAm1N+E9+p3IOMYJgbLqvyde5OuIxP3e3Oq+nvRuV4oneuFEh7sQ67TxJHnwuF0YbMZXNm4OmFBeseOSElRGBGRsumzm2D3QpJr9sA4voPgvOOsdTXh9txx5ODFP3zXMZ53MEwn1OoABzcAsKfX2wxfG8XB5Cwe6deYB/o2xGYzIPMk/PwMxHwGQJZfBJO9RrL3RCY1XMeobbgnP3IwAMNw/5XnNG0cJwSv0NrUrdeIRo2a4FGrHYREWnVlRCochRERKZsOb4XJPQH3Xz2uak35tv0UZm1NJyfPyZu3tKXu4Z9g9l1gOt3bXPkUXPkE2Q4nJzNyiQjxPXu/uxfB/IchOeGSS3MadsxOd+Fx5ZPgF3rxDfJy3E/3VG8GnmphEfkrhRERKbvm3gebZ0BwJIz6GYIizl5n21z47iFoeZN7fJLC9P3ISYffXoEdC8C/OoRE/TFFgk/wGSsa4MrjeNJ+9u7dRdrR/dRwHqalbR8Aed4hePR9GjrecfYgbaf2we5fYNcv7vFTHBlQ/0q4bS7YbJd4QUQqJoURESm7spJh3UfQ6iaoUvf867mcYLOXeDkOp4sFW5L4+dsZPJg3nSa2AwCYoQ0xQiIh69SfU07quXcy4BXoFl3itYqUJwojIiJFdDIjl2fnxhASN4OxHl8RaqSftY4TG9vszVjkaMUvjta0s+3iJc9p5BleZN+5iIDI1gXW33c8gzXxJ2gREUzLWsFn7U+kIivs77dGBxIR+UOovxf/G9aJ736vzfXzetMldxVO00YyASSbASQTwBGzSv7TOoYB++z16evcRF97DAc+uo3vu3xGz6a1WLbrOD/HHmbnkfT8dUf1qMe/BjTBx7PkW3tEyhO1jIiInMPR1Gy+2nCAbIcTm2HgYTOw2QwCvD2ICvUjqqoftUJ88bAZ/Lx2M91//hshZiqT8gbzWt4/8/fjYTNoXCOQ2CT37Z361fz5z81tNDCbVAq6TSMiUopcsd9i+/J2XBiMMl7At2FPrmkSQp8aOQTmHOL3fUnMWLWPtKxcPAwXXZo3pOeAm4mqVrj39IiURwojIiKlbV60e7wT72D3o77pRy64+o/OTrxX5XF6tajD1c1r0KpWsHvsFJEKQmFERKS0Zae6x1BJ3v/nPK8ACKkD3oFg2MBm52Smg8CjG/DEQYyrPqNz/8UxQmgXFcKUEZ0I9fOEhFWw4wf3o8Ve/uAVyAmHJ87gKMKa9zr7kWORMkhhRETECqlJcHA9BNVyP7bsW+XcY6TsX4Vr5lBsWSc54RHGHTmPscURzm0hsTxT5We8ktaf9xA5HgF4Nbkao/E10PBq8K8KQEqmg7cW7eRoWg6P9GtEw7DAEjpJqSgyHZl0mdEFgDVD1+Dn6Ves+1cYEREp607sgc9vhpN7cHkGkJAXTF3zIAAuuze2ljeQ6xnIhp2JnDh1kkCyaGXbW/CRY8OG2eFOfgi/l2d/2M/xdPcLBz1sBnf2rMcDfRsS6KNWFDm3shJG9GiviIhVqjaAu36BWbdh27+CuqSTjj8f513FN7bBPFC/J2/9spM9xzLwsBk8ObApCXaT7xbM5wo2MsBzM43NeIz1H9HC9S31HPcSXL0Ddar68+v2o3ywdC9zNx3kyWuaUreaH/uOZ7LvRAbxxzM4mpaD02WS5zJxulyYJgxqHc69vRqo34qUOrWMiIhYLS8HVrwDXv6caHwLI2bEsfXgnyO91gzyYeKwdnSo435fztaDKTzwxSbij2fQ3baV1z0/oLZxHBMDZ5d78ej3HL/tSeOF77ax70RmgUNVJYXHPGbR1RZHMv6cNIM4SRDHzSDizXCqNerEg/8YjLfPRf6FbJqwYRps+Rp6PgKNrsbhdPHS/Fj2Hs+gVogvtUJ8iQjxJaqqH+2jqmBXyClzykrLiMKIiEgZk5bt4K6P17Mm/iRXNKrGW7e2pWqAd4F10nPyeGbeVuZuOsiVdbx5J/RrguK+cC+sUhc6jSan5a1M2ZjC1OXxeNvgDr9lDE2bip8z7YLHz8MDI6wZ9siO7vfzhLcpuELWKfj2AYj77o8ZBuZVz/LE4b58ueHgOffZvUFV3hvWnhA/r3Muz81z4eWhd/uUNoWRIlAYEZHKJs/pYseRNJrWDLpgi8Lx9Byq+nthGAbs/Am+fRDSD7sX2r2hxRBoci2sehcOrHPPr9ka+jzlbt3IPA4ZxyHjGKf2bcZI2kzIX4fBb9AXejwE9XpDwmr3G5VTD4DNE+pdAXt+BWC+swtP5N3DfVe3Js9lcig5i0PJ2WzYf4osh5N61fyZMqIj9av/ObZKcmYu/124k8/XJNCtQVVeGtKSOlX9i/NSygUojBSBwoiISCHlpMOWr2D9VDj8e8FlXoHQ92noNBrs5+4yuD0phXHTFhCWvoMh3usYwCpsptO9sHozOL4TTCdUqQc3TYVa7Vnz5QTabXsFL8PJqcDGVBn8EqQlwcm9cGIPWScP8sOJGnyf1YKtXm14Y1gPutYP5Yu1CbyxcCfJmY784/t42ni4X2Pu6h6FB073eC0l7FByFj9sPUz/5jWIDC3eH+OyLjP5GH5Vwtx/PnUUv5Dqxbp/hRERkcrMNOHQRncoif0OGvSBa8ZDUMRFNz2cks0d09cRl5RKbeMoY3x+4kbjVzxd7id1aHULDHoDfIL4adth7vtsA+3ZzscB7+LvOHnBfeeYHqwxmxPn045l6RFsd0VRtUZtovs2ZPaaXXjuW0J/23oGeG4iwMOJ/Z9fQP3ehT/v7FRIPeR+u3J2inuye0Gj/mcFm6xcJx8s3cukJbvJdrjw87LzxDVNub1rnUrTiVdhpAgURkRESldWrpNZ6xL4cFk8B5OzCCWVoV5L8ahWnw1+V2C327AbBst3Hycnz8U/O0fyylWhGPPHwondEFrfPVVtAH5VIWEVrp0/Y0tJOOtYpn8YRmh9zMO/YzgKdrjNs/viMeIbiOpy1nY7j6Sxeu8Jbu4QiW96Aqx8BzZ9Ds6cs08osgv84wvwr4ppmny/JYnxC7ZzMDkLgGoB3vmPRXeuF8rrN7ambrWKf7tIYaQIFEZERKyR53Tx/ZYk3l+yN/9lf3/Vr1kNJt/WHg/7RTqgmibmsR1s+vUr/I9toqFrH/ZTe4EzfoaCI8msfw3vH25EhwOf0su+hVyPALzu/B4i2uav9k3MQR7/+nfqOPfzVNCP9M5dinH6dpJPMPiEgE+Qe2j+w1sgJ4WswDp8Vn8CcxN8888lItiHpwY1Y2DLcD5bvZ/XftxOZq6Thp7HuLdnHW64+soK3UqiMFIECiMiItYyTZP1+09xKDkLl2mS5zRxukyCfD3p16zGpT8Jk5sBR+PcrSlhzdydaw0Dl8vkxbnrGbh5DF1s28n2DMFn9I84qzXl7e83kLR6FkNsK+hh35a/q6yo3vj2eQzq9gTDIM/p4puYQyxftYKxx54m0jjGSTOA0bmPss2jGff1bsjdverj62XP38ehfTvY9eXT9Mz4BRc2xoe/zZjbbyXU/+yngBJOZHI4NZtOdau4OxCXQwojRaAwIiJS+ZimyX++XU//DffQ1raHDM+qbPNsSeuMlfgY7k6vJgaLbV15I2sQh/2bMm1kZ5qGBzJ340He/W03CSfdt32qkcInvm/Q3NyN0+ZFVs8nCYhqC8GREFwLHFmw7A1Y9xE4c/NriHNFMtp7Am/8sxNd6ruH3d99NI13f93Nt5sP4TLhhna1eOnvLfHzKn/jiCqMFIHCiIhI5WSaJpN+WE+f1XfSzPZnf5O0gPoEdh4GrW7miL0GI6e5O9z6edmp4ueV3xck1N+LUT3rMahVOHWCwJg9GnZ8f/aBDLv7KSGAuldAj4fIm303HtknedNxE++6buC+Kxuw/0Qm329J4vQvp80AlwkNwwJ4b1h7GtcoX+8DUhgpAoUREZHKbfrP6whZ9ixZXlXp9Ld7adi6R4EXEKZlO7j3sw2s2H0CcHdIvadXfYZ1jSrYYuFywqqJsG8ZpByA5ETI/WMQuPC20O85qN/Hve8tX8PsUeThwbU5L7PTjMzfTf/mNXigbyOyHE4e+GIjR1Jz8PW0839DWnJTh9qlcUmKhcJIESiMiIjI3mPp1Ajywd/73LdDcvNcfLB0D8G+ntzcMRIfT/s51ztLdop7VNngKLCd0ffFNOGLf8LOHzgZ0orrMp6lbd2qjOnTiOYRf/4WHU/P4ZFZMSzbdRyAgS1r8ujVjWiYusZ92+donPvR6la3EO/fmhlrE4mq6s+wzlHn7By7+2g6T8z+HZdpck+vBvRvXuOs9bIdThbFHcVlmlzXOvyS+6wojBSBwoiIiFgi9RBM7OIet6T/y9B9zDlXc7lM3vt1B58s2sAQ2zKGeSyijnH0rPUOmaF85+zGbGcvQuq04T83ty4w4uzXGw7wzLytZDmc+fOa1gxkTN+GDGwZzu6j6cxcl8DcTQfzB4sb1Cqc125qTcB5QtqFKIwUgcKIiIhYZsN0+O4h8PCFPuMg/SgkJ7hv82Qccz8R5MiEvOwCm6WYfsxx9eZAUHuapCznGvtagoys/OUrnC2YYVxLt4FDGdI+ime/2cqcje53+3RvUJV2USF8snI/aTl5AO5h/zOO0dSWQBMjgXbeBznq8OM1xy1EhoXy/u0daHDGUPtnMU330P3rp4JXANTpRmZgI/ya9AQg+7u5+AwcDPZCtigVgsKIiIhIcTBN+Hiwu59JYYS3IbHBUF7Y15xfdrv7o3jaDW5qXY0HovYRsf8bzB0LMEwXAPtdYXxJf1Y4GrOHWtzdrw3392mI3WaQeuwgq37+Eteun+lobqO6kXLW4dYYbRiR9TAe3v5MuLkNVzSqRuKpTBJPZpF4MhNfLzs3hO7De+l4SFj554ZxDlw/ZGNLOyMG1K4Nb78NN9xwyZfrTAojIiIixSU5EX54HOyefzwOHAkhkRBQA7z83ZOnP3j5gadv/mbr950k7nAaVzerQc1gnwL7M9d+SO7aaXjn/WUwueAoCGsK6UcgaXOBRSYGRtUGUKOFe4TbNR+AI4NYr5bckvow6RR8t05HYzsPe8ym5+nxWOze7jcxr03AfH4mAAV6m5zue/L118USSBRGREREyrrcTI6v/JSc3+cQnhOPLePI2euEt4GGV0PDq9xP/HidETgS18JnN0FOCgf9m3PtiYdJxY/rfTZzj8d3NMvb7j6MaWeWsw87G9/DA3/rQXCzxngdPsQ5u70ahruFJD7+sm/ZKIyIiIiUN5kn4dh29xM4XgHup3ACwi68zaEY+PTvkHWSvKpNMEwn9pO73cvsXuS2vJX38obwzsYcXCZccXArn3725MVr+e03uPLKyzqdwv5+l7/h4kRERCoqv1Co0909FVZEWxj5PXxyPR4ndrjn+QRDx1HQ5V68AmvwMHBVtxSemruFKrEnCrffpKQiFn/pFEZERETKuxrN4c4f4bdXIKIddBgB3gVHg21VO5h50T2ICTwK3xVin+HhJVPrOeg2jYiISGXidELdunDwIJwrAljQZ+QSX7MoIiIi5ZLd7n58FzD/OnLr6c9vvVWs441cjMKIiIhIZXPDDfD115h/vRVTu3axPdZbFAojIiIildENN5D9+4b8j9nfzXXfminlIAIKIyIiIpXXGbdiXD17lOqtmTMpjIiIiIilFEZERETEUgojIiIiYimFEREREbGUwoiIiIhYSmFERERELKUwIiIiIpZSGBERERFL6a29IiIilZW/P62mtwRgjb+/ZWWoZUREREQspTAiIiIillIYEREREUsVOYwsXbqUwYMHExERgWEYzJs374Lrjxw5EsMwzppatGhxqTWLiIjIZXgv5j0mb558zmWTN0/mvZj3SrWeIoeRjIwM2rRpw8SJEwu1/ttvv01SUlL+lJiYSGhoKDfffHORixUREZHLZzNsTIyZyJQtUwrMn7x5MhNjJmIzSvfGSZGfphk4cCADBw4s9PrBwcEEBwfnf543bx6nTp3ijjvuKOqhRUREpBjc2+ZeACbG/NmwMGXLFD7Y8gHRbaPzl5eWUn+0d8qUKfTr1486deqcd52cnBxycnLyP6emppZGaSIiIpXGvW3uxeF08MGWDwAsCyJQyh1YDx06xA8//MBdd911wfXGjx+f36ISHBxMZGRkKVUoIiJSeYxqNSr/z542T0uCCJRyGPn4448JCQlhyJAhF1xv3LhxpKSk5E+JiYmlU6CIiEglcmafEYfLcd5OrSWt1G7TmKbJ1KlTuf322/Hy8rrgut7e3nh7e5dSZSIiIpXP5M2T82/RANzd6u78PiSl3UJSai0jS5YsYffu3YwaNeriK4uIiEiJOf3UzN2t7s6fN6rVKKLbRjMxZmKpt5AUuWUkPT2d3bt353+Oj48nJiaG0NBQoqKiGDduHAcPHuSTTz4psN2UKVPo0qULLVu2vPyqRURE5JK5TBfRbaMZ3nx4gdaR0y0iLtNVqvUUOYysX7+ePn365H8eO3YsACNGjGD69OkkJSWRkJBQYJuUlBRmz57N22+/fZnlioiIyOW6v+39AGQ6Ms9aZkUn1iKHkSuvvBLTNM+7fPr06WfNCw4OJjPz7BMWERER0btpRERExFIKIyIiImIphRERERGxlMKIiIiIWEphRERERCylMCIiIiKWUhgRERERSymMiIiIiKUURkRERMRSCiMiIiJiKYURERERsZTCiIiIiFhKYUREREQspTAiIiIillIYEREREUspjIiIiIilFEZERETEUgojIiIiYimFEREREbGUwoiIiIhYSmFERERELKUwIiIiIpZSGBERERFLKYyIiIiIpRRGRERExFIKIyIiImIphRERERGxlMKIiIiIWEphRERERCylMCIiIiKWUhgRERERSymMiIiIiKUURkRERMRSCiMiIiJiKYURERERsZTCiIiIiFjKw+oCRERExBp+nn5sGbHF6jLUMiIiIiLWUhgRERERSymMiIiIiKUURkRERMRSCiMiIiJiKYURERERsZTCiIiIiFhKYUREREQspTAiIiIillIYEREREUspjIiIiIilFEZERETEUgojIiIiYimFEREREbGUwoiIiIhYysPqAgrDNE0AUlNTLa5ERERECuv07/bp3/HzKRdhJC0tDYDIyEiLKxEREZGiSktLIzg4+LzLDfNicaUMcLlcHDp0iMDAQAzDuOT9pKamEhkZSWJiIkFBQcVYofyVrnXp0bUuPbrWpUfXuvSU5LU2TZO0tDQiIiKw2c7fM6RctIzYbDZq165dbPsLCgrSl7uU6FqXHl3r0qNrXXp0rUtPSV3rC7WInKYOrCIiImIphRERERGxVKUKI97e3jz33HN4e3tbXUqFp2tdenStS4+udenRtS49ZeFal4sOrCIiIlJxVaqWERERESl7FEZERETEUgojIiIiYimFEREREbFUhQsjEydOpG7duvj4+NClSxfWrl17wfW/+uormjZtio+PD61atWLBggWlVGn5V5Rr/eGHH3LFFVdQpUoVqlSpQr9+/S7630b+VNTv9WkzZ87EMAyGDBlSsgVWIEW91snJyURHRxMeHo63tzeNGzfW3yOFVNRr/dZbb9GkSRN8fX2JjIzkkUceITs7u5SqLb+WLl3K4MGDiYiIwDAM5s2bd9FtFi9eTPv27fH29qZhw4ZMnz69ZIs0K5CZM2eaXl5e5tSpU81t27aZo0ePNkNCQswjR46cc/0VK1aYdrvdfP31183Y2Fjz3//+t+np6Wlu2bKllCsvf4p6rYcOHWpOnDjR3LRpkxkXF2eOHDnSDA4ONg8cOFDKlZc/Rb3Wp8XHx5u1atUyr7jiCvP6668vnWLLuaJe65ycHLNjx47mtddeay5fvtyMj483Fy9ebMbExJRy5eVPUa/1559/bnp7e5uff/65GR8fb/70009meHi4+cgjj5Ry5eXPggULzKefftqcM2eOCZhz58694Pp79+41/fz8zLFjx5qxsbHm//73P9Nut5s//vhjidVYocJI586dzejo6PzPTqfTjIiIMMePH3/O9W+55RZz0KBBBeZ16dLFvOeee0q0zoqgqNf6r/Ly8szAwEDz448/LqkSK4xLudZ5eXlm9+7dzY8++sgcMWKEwkghFfVaT5o0yaxfv76Zm5tbWiVWGEW91tHR0Wbfvn0LzBs7dqzZo0ePEq2zoilMGHn88cfNFi1aFJh36623mgMGDCixuirMbZrc3Fw2bNhAv3798ufZbDb69evHqlWrzrnNqlWrCqwPMGDAgPOuL26Xcq3/KjMzE4fDQWhoaEmVWSFc6rV+8cUXCQsLY9SoUaVRZoVwKdf622+/pVu3bkRHR1OjRg1atmzJK6+8gtPpLK2yy6VLudbdu3dnw4YN+bdy9u7dy4IFC7j22mtLpebKxIrfxnLxorzCOH78OE6nkxo1ahSYX6NGDbZv337ObQ4fPnzO9Q8fPlxidVYEl3Kt/+qJJ54gIiLirC+8FHQp13r58uVMmTKFmJiYUqiw4riUa713715+/fVXhg0bxoIFC9i9ezf3338/DoeD5557rjTKLpcu5VoPHTqU48eP07NnT0zTJC8vj3vvvZennnqqNEquVM7325iamkpWVha+vr7FfswK0zIi5cerr77KzJkzmTt3Lj4+PlaXU6GkpaVx++238+GHH1KtWjWry6nwXC4XYWFhfPDBB3To0IFbb72Vp59+msmTJ1tdWoWzePFiXnnlFd577z02btzInDlz+P777/m///s/q0uTYlBhWkaqVauG3W7nyJEjBeYfOXKEmjVrnnObmjVrFml9cbuUa33ahAkTePXVV/nll19o3bp1SZZZIRT1Wu/Zs4d9+/YxePDg/HkulwsADw8PduzYQYMGDUq26HLqUr7X4eHheHp6Yrfb8+c1a9aMw4cPk5ubi5eXV4nWXF5dyrV+5plnuP3227nrrrsAaNWqFRkZGdx99908/fTT2Gz6t3VxOd9vY1BQUIm0ikAFahnx8vKiQ4cOLFq0KH+ey+Vi0aJFdOvW7ZzbdOvWrcD6AAsXLjzv+uJ2Kdca4PXXX+f//u//+PHHH+nYsWNplFruFfVaN23alC1bthATE5M//e1vf6NPnz7ExMQQGRlZmuWXK5fyve7Rowe7d+/OD3wAO3fuJDw8XEHkAi7lWmdmZp4VOE6HQFOvWCtWlvw2lljXWAvMnDnT9Pb2NqdPn27Gxsaad999txkSEmIePnzYNE3TvP32280nn3wyf/0VK1aYHh4e5oQJE8y4uDjzueee06O9hVTUa/3qq6+aXl5e5tdff20mJSXlT2lpaVadQrlR1Gv9V3qapvCKeq0TEhLMwMBAc8yYMeaOHTvM+fPnm2FhYeZLL71k1SmUG0W91s8995wZGBhofvHFF+bevXvNn3/+2WzQoIF5yy23WHUK5UZaWpq5adMmc9OmTSZgvvnmm+amTZvM/fv3m6Zpmk8++aR5++23569/+tHexx57zIyLizMnTpyoR3uL6n//+58ZFRVlenl5mZ07dzZXr16dv6x3797miBEjCqz/5Zdfmo0bNza9vLzMFi1amN9//30pV1x+FeVa16lTxwTOmp577rnSL7wcKur3+kwKI0VT1Gu9cuVKs0uXLqa3t7dZv3598+WXXzbz8vJKueryqSjX2uFwmM8//7zZoEED08fHx4yMjDTvv/9+89SpU6VfeDnz22+/nfPv39PXd8SIEWbv3r3P2qZt27aml5eXWb9+fXPatGklWqNhmmrfEhEREetUmD4jIiIiUj4pjIiIiIilFEZERETEUgojIiIiYimFEREREbGUwoiIiIhYSmFERERELKUwIiIiIpZSGBERERFLKYyIiIiIpRRGRERExFIKIyIiImKp/wcA6Nq0OOzbbgAAAABJRU5ErkJggg==" }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.plot(c_sticking_probability.p, c_sticking_probability.fd_mean)\n", "plt.plot(rust_sticking_probability.p, rust_sticking_probability.fd_mean)\n", "plt.errorbar([1], [1.70], marker='x', yerr=[0.06])\n", "plt.errorbar([1], [1.7311567309540636], marker='o', color='red', yerr=[0.020648406595153657])\n", "plt.legend([\"C\", \"Rust\", \"Previous Model\", \"Theory\"])" ], "metadata": { "collapsed": false } }, { "cell_type": "code", "execution_count": 125, "outputs": [ { "data": { "text/plain": "" }, "execution_count": 125, "metadata": {}, "output_type": "execute_result" }, { "data": { "text/plain": "
", "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiMAAAGdCAYAAADAAnMpAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy88F64QAAAACXBIWXMAAA9hAAAPYQGoP6dpAABUiklEQVR4nO3de3xTVbo38N/eadILbVMKlLbSUkCxFBQVqAKiwAEpKsh7HGGGEXDEG6BY9eiAo7QwM4KDOuCIqKMcBEVAVERAEGWQo4AFhBlqKYpQirYFAXujtEmz1/tHmrS5tUmae37fzyfndO+snb2yZejDWs96liSEECAiIiLyE9nfHSAiIqLwxmCEiIiI/IrBCBEREfkVgxEiIiLyKwYjRERE5FcMRoiIiMivGIwQERGRXzEYISIiIr+K8HcHnKEoCsrKyhAXFwdJkvzdHSIiInKCEAI1NTVITU2FLDse/wiKYKSsrAxpaWn+7gYRERG54fTp0+jWrZvD94MiGImLiwNg/DLx8fF+7g0RERE5o7q6Gmlpaebf444ERTBimpqJj49nMEJERBRk2kqxYAIrERER+RWDESIiIvIrBiNERETkV0GRM0JERORtBoMBer3e390IKiqVChEREe0uu8FghIiIwl5tbS1++uknCCH83ZWgExMTg5SUFGg0Grc/g8EIERGFNYPBgJ9++gkxMTHo0qULi2s6SQgBnU6HX375BSdPnsQVV1zRamGz1jAYISKisKbX6yGEQJcuXRAdHe3v7gSV6OhoqNVqnDp1CjqdDlFRUW59DhNYiYiI0HYtDLLP3dEQi8/wQD+IiIjCXp2uERlztiBjzhbU6Rr93Z2gwmCEiIiI/IrBCBERkQcYlOaVOAUnL1gcU+sYjBAREbXTtsJyjHrpS/PxPf+7Hzc+vxPbCsu9fu+Kigo88sgj6NmzJyIjI5GWloZx48bhiy++8Pq9PYXBCBERUTtsKyzHjHe+xZnqBovzFVX1mPHOt14NSEpKSjBgwADs3LkTixcvxpEjR7Bt2zaMGDECs2bN8tp9PS1sl/bW6RqRNW87AKBowRjEaML2URARkZsMisD8T4pgb0JGAJAAzP+kCKOzkqGSPb9aZ+bMmZAkCQUFBejQoYP5fN++fXHvvfd6/H7eEr4jI7qLKImajJKoyYDuor97Q0REQajg5AWUV9U7fF8AKK+qR8HJCx6/94ULF7Bt2zbMmjXLIhAxSUhI8Pg9vSV8gxEiIqJ2OlvjOBBxp50rjh8/DiEEMjMzPf7ZvsZghIiIyE1Jcc5VHHW2nStCaR8dBiNERERuyu6RiBRtFBxlg0gAUrRRyO6R6PF7X3HFFZAkCcXFxR7/bF9jMEJEROQmlSwhb1wWANgEJKbjvHFZXkleTUxMxJgxY7Bs2TJcvGib+1hZWenxe3oLgxEiIqJ2yOmXguV3X4ek+EiL88naKCy/+zrk9Evx2r2XLVsGg8GA7OxsfPDBB/jhhx9w9OhRvPzyyxg8eLDX7utpXM9KRETUTjn9UjD08s64Kv8zAMDKPwzCsCu6eGVEpKWePXvi22+/xV//+lc88cQTKC8vR5cuXTBgwAAsX77cq/f2JAYjREREHtAy8Mjukej1QMQkJSUFr7zyCl555RWf3M8bGIwQERF5QIwmAiWLbvN3N4ISc0aIiIjIrxiMEBERkV8xGCEiIiK/YjBCREREfsVghIiIiPyKwQgRERH5FYMRIiIiT9BdBPK1xpfOtjw7OcZghIiIiPyKwQgREZEnKIbmn0/tsTz2gnvuuQeSJEGSJKjVavTo0QNPPfUU6uvrPfL5GRkZWLJkiUc+qy0MRoiIiNqraBOwLLv5+N3fAEv6Gc97UU5ODsrLy3HixAn8/e9/x+uvv468vDyv3tMbGIwQERG1R9EmYP1UoKbc8nx1ufG8FwOSyMhIJCcnIy0tDRMmTMCoUaOwY8cOAPZHNq655hrk5+cDAIQQyM/PR3p6OiIjI5GamorZs2cDAIYPH45Tp07hscceM4++eBODESIiIncpBmDbHwEIO282nds2x+tTNgBQWFiIPXv2QKPRONX+gw8+MI+m/PDDD9i4cSOuuuoqAMCHH36Ibt26YcGCBSgvL0d5eXkbn9Y+4btRXos/GHLpXiBzNCCr/NghIiIKOqf2ANVlrTQQQPXPxnY9hnn89ps3b0ZsbCwaGxvR0NAAWZad3r23tLQUycnJGDVqFNRqNdLT05GdbZxqSkxMhEqlQlxcHJKTkz3eb2vhOTJStAlRbwwxH0atn+STuT0iIgoxtWc8285FI0aMwOHDh/HNN99g2rRp+MMf/oA777zTqWvvuusuXLp0CT179sT999+Pjz76CI2NjV7pZ1tcCkYWLlyIQYMGIS4uDklJSZgwYQKOHTvW6jXfffcd7rzzTmRkZECSJJ9l5jrUNLcn1VoOOQkfzO0REVGIie3q2XYu6tChAy6//HL0798fK1aswDfffIO33noLACDLMoSwnD7S6/Xmn9PS0nDs2DG8+uqriI6OxsyZM3HTTTdZtPEVl4KRL7/8ErNmzcK+ffuwY8cO6PV63HLLLbh40XFxl7q6OvTs2ROLFi3yyVBPq5rm9gQErFNxJAjj7J6P5vaIiCgEdB8CxKcCNr9VTCQg/jJjOy+TZRlPP/00nnnmGVy6dAldunSxyPWorq7GyZMnLa6Jjo7GuHHj8PLLL2PXrl3Yu3cvjhw5AgDQaDQwGHzz+9ClYGTbtm2455570LdvX/Tv3x8rV65EaWkpDh486PCaQYMGYfHixfjtb3+LyMjIdne4XZrm9hz/kWkxt0dERNQWWQXkPN90YPvPXABAziKf5STeddddUKlUWLZsGUaOHInVq1fj//7v/3DkyBFMmzYNKlVzP1auXIm33noLhYWFOHHiBN555x1ER0eje/fuAIyrcXbv3o2ff/4Z586d82q/25UzUlVVBcCY6BIMlJoK59qtvJ2lfImIyDlZ44GJq4A4q9H/+FTj+azxPutKREQEHn74Yfztb3/DnDlzcPPNN+P222/HbbfdhgkTJqBXr17mtgkJCfjnP/+JoUOH4uqrr8bnn3+OTz75BJ06dQIALFiwACUlJejVqxe6dOni1X5LwnpCyUmKomD8+PGorKzEV1995dQ1GRkZyM3NRW5ubqvtGhoa0NDQYD6urq5GWloaqqqqEB8f7053AQDffb0FfXdMbrudIR19n90HaDq4fS8iIgoO9fX1OHnyJHr06IGoqKh2fFA1sCjN+PPvNwC9RobFKs3Wnl91dTW0Wm2bv7/dHhmZNWsWCgsLsXbtWnc/wqGFCxdCq9WaX2lpaR753OMxV6FMJEJxEH4pAigTnXAcqR65HxERhZGWgUf3IWERiHiKW8HIww8/jM2bN+Nf//oXunXr5uk+Ye7cuaiqqjK/Tp8+7ZHPTYrvgPn6qQBgE5CYjufrpyAJVajT+Wd5ExERBSlNByC/yvjiyLpLXApGhBB4+OGH8dFHH2Hnzp3o0aOHVzoVGRmJ+Ph4i5cnZPdIxH/ibsJMfS4qYJnnUoFOmKnPxX+UnsiWiz1yPyIiImqbSxVYZ82ahTVr1uDjjz9GXFwcKiqMCaFarRbR0dEAgKlTp+Kyyy7DwoULAQA6nQ5FRUXmn3/++WccPnwYsbGxuPzyyz35XdqkkiXkjcvCjHfqsaNhIAbJxUhCJc4iAfuVTCiQsVz9d6gkt9JoiIiIyA0ujYwsX74cVVVVGD58OFJSUsyvdevWmduUlpZarGsuKyvDtddei2uvvRbl5eV44YUXcO211+K+++7z3LdwQU6/FCy/+zp0jovCPiULm5Qh2KdkIUkbgyW/yUSOar9f+kVERBSuXBoZcWbhza5duyyOMzIynLrOl3L6pWBoejQKF+fgLBKgnfQ6hmWlo6GuGtjs794REZE/BNrvqmDhiecWthvlqWQJg1VHAQB1GQlQyd7dHpmIiAKTqRCYTqczpxyQ8+rq6gAAarXa7c8I22CEiIgIMBYKi4mJwS+//AK1Wg1ZDs89ZF0lhEBdXR3Onj2LhIQEi+qurmIwQkREYU2SJKSkpODkyZM4deqUv7sTdBISEtq99xyDESIiCnsajQZXXHEFdDqdv7sSVNRqdbtGREwYjBAREcG46227ysGT2zgxRkRERH7FYISIiIj8isEIERER+RWDESIiIvIrBiNERETkVwxGiIiIyK8YjBAREZFfMRghIiIiv2IwQkRERH7FYISIiIj8isFIS4rB/KNcutd4rLsI5GuNL91FP3aOiIgoNHFvGpOiTYja+pT5MGr9JCA+FRi1wI+dIiIiCn3hG4xoOiCjfg0A4NiJncCHf4AEYdmmuhz48H4/dI6IiCh8hP00jQwF6h1PAxCQbN4VthcQERGRR4V9MJItF0OuKWulBQMSIiIibwr7YCQJlf7uAhERUVgL+2DkLBL83QUiIqKwFrbBSIwmAiWLbsO7Cx7DGXSC4mA2RhFAhUiAgbM1REREXhG2wYhJwakqzNNNAQCbgMR0nKe/BwVKHx/3jIiIKDyEfTBytqYe25VszNDnogKJFu9VoBNm6HOxXcnGWSSgTtfop14SERGFrvCtM9IkKS4KALBdycaOhoHIlouRhEqcRQIKlEwoTfEaE12JiIi8I+yDkeweiUjRRqGiqh4KZOxTsizelwAk4xyy5WI0+KeLREREIS3sp2lUsoS8ccYAxLromek4T70aKokZrERERN4Q9sEIAOT0S8Hyu69DUpzG4nyyNgpLfpOJHNV+P/WMiIgo9IX9NI1JTr8UDE2PRuHiHJxFArSTXsewrHQ01FUDm/3dOyIiotDFYKQFlSxhsOooAKAuIwEq2Xa3GiIiIvIsTtMQERGRXzEYISIiIr9iMEJERER+xWCEiIiI/IrBiIvqdI3ImLMFGXO2sDw8ERGRB3A1jat0F1ESNRkAUKcrBTRaP3eIiIgouHFkhIiIiPyKwQgRERH5FYMRIiIi8isGI0RERORXDEbaohjMP8qley2OiYiIqP24mqY1RZsQtfUp82HU+klQYlP82CEiIqLQw2DEAdX3nwKbHoIEYXFeqq3wU4+IiIhCE6dpHFDvzAcgYL1vr0VwwikbIiKidmMw4oBcW952m58KfNATIiKi0MZgpB2ki2f93QUiIqKgx2CkJU0HZNSvwW91zzjVXHRIAsD9aoiIiNqDwYgdBUomlLhUOxkjRqIpbUTplu3DXhEREYUmBiN2KJBxoM8fIYSAIqzeE4AAcMjQC5BVfukfERFRKGEw4sAjh7thhj4XFUi0OF+BTpipfxQz9Y/CYIpUmnbyLYmaDOgu+qG3REREwYt1Rhw4U92A7cjGjoaByJaLkYRKnEWCcQqnKYY7WFqF4f0S/NtRIiKiIMdgpIUYTQRKFt2Gjw//jEfXHgZgnLLZp2TZbf9Lrd6HvSMiIgpNnKaxIykuyql2XWLVXu4JERFR6GMwYkd2j0SkaKMcrKUBJChIwTkMSNf6tF9EREShiMGIHSpZQt4449SMbTl44//NU6+GSnYUrhAREZGzGIw4kNMvBcvvvg5JcRqL813jNViuXoIc1X4/9YyIiCi0MIG1FTn9UjA0PRqFi3NwFgnQTnodA9K1iHuJgQgREZGnMBhpg0qWMFh1FABQl5Hg384QERGFIAYjLjIoAnsNfYwjJSWVGJYV7+8uERERBTUGIy7YUXwOi7b9gO6GO41F0Nasx59i++OpW3rhDn93joiIKEgxGHHSNsMgfPrhKnygXoVUzQXz+bKGRCz4cCoiVYOY1EpEROQGrqZxgkFI2GW4Gq+qlyAZFyzeS8YFvKpegl2Gq2EQTUt9dReBfK3xxb1qiIiIWsVgxAkFypWYrf4IAGBdWsR0PFu9EQXKlT7uGRERUfBjMOIEA2SkShdsAhETWQJSpfMw8HESERG5jL89nZAkVXq0HRERETVzKRhZuHAhBg0ahLi4OCQlJWHChAk4duxYm9e9//77yMzMRFRUFK666ips3brV7Q77Qy+pzKPtiIiIqJlLwciXX36JWbNmYd++fdixYwf0ej1uueUWXLzoOElzz549+N3vfofp06fj0KFDmDBhAiZMmIDCwsJ2d95XVBJwSaihCPvvKwK4JDRQNU3j1Okaze+1/JmIiIhsSUIIB79i2/bLL78gKSkJX375JW666Sa7bSZNmoSLFy9i8+bN5nM33HADrrnmGrz22mtO3ae6uhparRZVVVWIj/dtkbG62irEvJAOABCQICAsIjgFgAQJEoyPse5/SgHAfE3d/5QiJpa7+xIRUfhx9vd3u3JGqqqqAACJiYkO2+zduxejRo2yODdmzBjs3bvX4TUNDQ2orq62eAUC3fjXgNgUi3NS/GXG80REROQWt4MRRVGQm5uLoUOHol+/fg7bVVRUoGvXrhbnunbtioqKCofXLFy4EFqt1vxKS0tzt5seZeg9FvUP7DEf109cByn3CAy9x/qxV0RERMHN7WBk1qxZKCwsxNq1az3ZHwDA3LlzUVVVZX6dPn3a4/dwm6wy/6ikD7Y4JiIiIte5VQ7+4YcfxubNm7F7925069at1bbJyck4c+aMxbkzZ84gOTnZ4TWRkZGIjIx0p2sBq2XuCfNIiIiImrk0MiKEwMMPP4yPPvoIO3fuRI8ePdq8ZvDgwfjiiy8szu3YsQODBw92radEREQUklwaGZk1axbWrFmDjz/+GHFxcea8D61Wi+joaADA1KlTcdlll2HhwoUAgEcffRQ333wzXnzxRdx2221Yu3YtDhw4gDfeeMPDX4WIiIiCkUsjI8uXL0dVVRWGDx+OlJQU82vdunXmNqWlpSgvLzcfDxkyBGvWrMEbb7yB/v37Y8OGDdi4cWOrSa/BxqAI7DX0wceGwSgoqYTBUUESIiIisuHSyIgzJUl27dplc+6uu+7CXXfd5cqtAoemAzLq1wAAijQdbHbh3VZYjryPC3FG/6zxxDuF6BqnwXzDIOSo9vu6t0REREHHrQRWMtpRfA65G4ohQcENcjGSUImzSMD+mkzMQC6WYwnsl4IjIiIiEwYjbjIICc9tP4Fb5ALkqVchVbpgfq9MJGKBfgrm66fgMwdTNnW6RmTN2w4AKFowBjEa/qcgIqLwxF173VSgZOKai19huXoJknHB4r1kXMCr6qW4Wj6Bg6VVfuohERFRcGAw4qaz0CJPvQoAIEuW75mO89Sr8UtNvY97RkREFFw4N9CGGE0EShbdZj6u0xn//+Uos5iasSZLQCrO44pLRwBc6eVeEhERBS+OjLipj1zqVLsrYy+23YiIiCiMMRhxk/XUjCNSbNe2GxEREYUxBiPtoMSmQMB+VGIqyaJ0y/Zhj4iIiIIPg5F20I/MhwTYBCQCEiTTKe7qS0RE1CoGI+1g6D0WmLgKItZyB2IR53hHYiIiIrLEYKS9ssaj/oE95sP6ietw8Q+78Z0hHQeUK/B9wQ4YGhv92EEiIqLAxqW9ntBiKuY/P5ai+/qZ6Ks6bzyxezrO7H4KpwbMgUX2iO4iSqImAwDqdKWARuu7/hIREQUQjox40CFDLww88CS6iPMW57uI8xh44EkcMvTyU8+IiIgCF4MRDzEIIFU2BiGOKrKmyOdhaHvjYyIiorDCYMRDipV0dJUqHdYfkSUgWapEsZLu244REREFOAYjHnJJimx3uzpdIzLmbEHGnC2o0zHplYiIwgODERfFaCLs/hwtGpy63tl2RERE4YKraVyl6QDkV9mczpRLcUYkoAvsT9UoAjiLBGTKpWA4QkRE1IwjIx6ikoAypRMAY+DRkum4XOkEVWt72jQt9y2JmgzouMEeERGFBwYjHnSt6kccGLgYv0idLM6flTrhwMDFuFb1o596RkREFLgYjHhYvxGTEJO731yB9fBNb6HLM9+j34hJ/u4aERFRQGLOiBeoIiLQV1UKAKjLHg1VBB8zERGRIxwZ8RFDYyP3qyEiIrKD/2T3gUPb30bq3vlW+9V0wnFDL+aREBFR2GMw4mWF/1qHgQeeNB60WEnTRZxHF/k8Dhl64Ur/dI2IiCggcJrGiwwC6H5wEYA29qvhlA0REYUxBiNeVKykoyvOt7lfzY/f/st+A91FIF9rfLHuCBERhSgGI17k7H419ZUVTn8m968hIqJQw2DEi5zdhyYqIdnLPSEiIgpcDEa8KFMuxRl0sikPb6IIoFx0xIXO18HgqBEREVGIYzDiRSoJODVgDgDH+9Xk66fh3jVHcePzO7Gj+JxFm5bTMJySISKiUMVgxMv6jZiEfw952Wa/mgp0wgx9LrYr2cbjqnrkbijGNsMgf3STiIjIb1hnxAeuHTMNNdffjiMvjsC7hlEoQQoKlEwoLWJBAWMZkvn6KRgtH/BbX4mIiHyNwYiPqCIiUIsYrFX+CwAgQ8ENchGSUImzSDAHJ+XojAIlE/0dfZDuIkqiJgMA6nSlgEbrmy9ARETkJQxGfOgsEgAAY+QC5KlXIVW6YH6vTCRivn4qtivZ5nZEREThgDkjPpSESoyRC7BcvQTJuGDxXjIuYLl6CcbIBUhCpX86SERE5AccGfGAmFgtkF9l/BlAXW2V3XbZ8lH0lMsB2C8PrwhgvnoluqASzlUoISIiCn4MRnxIJQFdWxn1kCUgmaMiREQUZjhNQ0RERH7FYISIiIj8isGIF8RoIuz+bCJgfxtfAQkNIgIfG25AQUklS8QTEVFYYM6InwhIkCAsjgGB2fqHjVVZ3ylE1zgN5hsGIUe13+Hn1OkakTVvOwCgaMEYu8EPERFRIOPIiB/oxr8GEWu5U2+5SMQM3WxUIRbj5T24QS7CuZp6zNDntl4ivqkIWknUZEB30cs9JyIi8jz+M9oPDL3HwnD5LYh5qQcUAcyKyAcu/Yp56tU2hdAW6Kdgvn4KPuOUDRERhSiOjHiDpoOx7kh+lfFne2QVAOAbpQ+US5VYpl5qtxDaq+qluFo+gYOl9muX2FOna0TGnC3ImLOFu/0SEVHAYzDiZ2ehRZ56FQD7hdAAIE+9Gr/U1Pu4Z0RERL7BYMTPLkcZUqULNoGIiSwBqdJ5XHHpiG87RkRE5CMMRvysj1zqVLsrY5mcSkREoYnBiJ85GhGxJsV29W5HiIiI/ITBSABQYlMcF0JrWkSjdMv2YY+IiIh8h8FIANCPzG8qeWYZkAhIkEynmlbfuKuutgrI1wL5Woe7ChMREfkDg5EAYOg9Fpi4yqYQmohLdnAFERFR6GAwEiiyxqP+gT3mw/qJ63DxD7vxnSEdB5Qr8H3BDhgaG2FQBPYa+uBjw2DH+9ewKisREQURVmANJC2mYv7zYym6r5+JvqrzxhO7p6Ni91P4G+7Bh/pnjeec3L+GiIgokHFkJAAdMvTCwANPoos4b3E+SZzHC+JFjJELzOfO1uja3r+GiIgogDEYCTAGAaTKxiCktYqsMhQAaNr3V2C+for9KRsiIqIAx2AkwBQr6egqVbZZkTVbLjafE5BRjs4u7V9DREQUKBiMBJhLUqRT7ZJQaXPul1q9h3tDRETkfQxGAky0aHCq3eXyT7hBLjJP1wBAl1i10/fhzr5ERBQouJrGR2I0ERY/OwoAMuVSnBEJ6AL7UzVCAJIEzI7YiNnYiDKRiAX6Kfi30gsD0rXOd6hp+S8A1OlKAY0WdbpGZM3bDgAoWjDGos9ERETewt82vqLpAOS3yOnQ2c/vUElAmaETusiVUIRlEquwk5+ajAt4Vb0U/1Z6QSVP8GiXGZwQEZEvcJomAF2r+hEHBi7GL1InyzckNJeHbyJLgNR0DRRD+27MYmlEROQHDEYCVL8RkxCTux/fGdJxTLkMABxspdccoMg/FThoQUREFLg47h7AVBER6Ksqdbq9dPGsF3tDRETkHRwZCSFHa2JY+IyIiIIOR0b8RdMBGfVrAABFmg5tNldiUyDXVsBUc9XiPQFUIRbPf/4jSvftwFO39MIdnuijnRU3REREnubyyMju3bsxbtw4pKamQpIkbNy4sc1rli1bhj59+iA6OhpXXnklVq1a5U5fw5p+ZD4AQFhljoimFTcdpVq8p3kO7zc8iE8/XMW9aoiIKGi4HIxcvHgR/fv3x7Jly5xqv3z5csydOxf5+fn47rvvMH/+fMyaNQuffPKJy50NJTGaCJQsug0li25zasmsofdYYOIqiNjkVtsZl/ouwS7D1TAIRymvREREgcPlaZqxY8di7NixTrdfvXo1HnzwQUyaNAkA0LNnT+zfvx/PP/88xo0b5+rtw1vWeNR3uxEnF9+EVPlXJKDW7lJfRQCz1RtRoFyJ/v7pKRERkdO8nsDa0NCAqKgoi3PR0dEoKCiAXm9/L5WGhgZUV1dbvKiJrMKviENHyTYQMTdp2kzPwPxkIiIKAl7/bTVmzBi8+eabOHjwIIQQOHDgAN58803o9XqcO3fO7jULFy6EVqs1v9LS0rzdzaCSJFV6tB0REZE/eT0YefbZZzF27FjccMMNUKvVuOOOOzBt2jTjzWX7t587dy6qqqrMr9OnT3u7m0Gll1Tm0XZERET+5PVgJDo6GitWrEBdXR1KSkpQWlqKjIwMxMXFoUuXLnaviYyMRHx8vMWLmqkk4JJQw1FJEUUANSIam5UbUFBSCYMiYFAE9hr64GPDYPM5t+guAvla40t3kbv/EhFRu/mszoharUa3bt0AAGvXrsXtt9/ucGSE2hYt6SEgQYGwiCgVAJCAN/S34pRIwdk16/F4ZD8okFGpf9bY6J1CdI3TYL5hEHJU+9vXEdYiISKidnI5GKmtrcXx48fNxydPnsThw4eRmJiI9PR0zJ07Fz///LO5lsj333+PgoICXH/99fj111/x0ksvobCwEG+//bbnvkWY0o1/Deqd+UBtuflcpYgFADyh/sB8rkxJxAL93aiU45GESpxFAvbXZGIGcrEcS3CTC/es0zUipsXPRERE7eVyMHLgwAGMGDHCfPz4448DAKZNm4aVK1eivLwcpaXN+6kYDAa8+OKLOHbsGNRqNUaMGIE9e/YgIyOj/b0Pc4beY2G4/BbEvNQDigD+qZqI+w3rbdol4wKWq1+2WH1TJhKxQD8F8/VT8BlLyBMRkR+5HIwMHz4cQjj+5bVy5UqL4z59+uDQoUMudywcxcRqgfwq488A6mqr2r5IVgEAvlEyMQ6fG0/ZqT1i/Z/MWBxtKWboc3GwtArD+yW0s/fN6nSNyJq3HQBQtGCMU0XdiIgofDFpI0QYICNVumATiJjYK44GAHnq1filpt6znWnKIymJmmxMeCUiImoFg5EQ4U5NEVNxtCsuHfF8h1rgihsiImoNg5EQ0Z6aIlfGcvSCiIj8h8FIiFBJwKWoJIe1R1ojxXb1fIeIiIicxGAkhMi3/BmSJBlrjbTgKD4xJbUq3bK92S3mkBARUasYjIQQQ++xkCauAmJTbN6zDkgEmpNa5Z8KAMXg9f4RERHZw2Ak1GSNR/0De8yH9RPXQbrrbQjrAEVq/k8ftX4SxJJ+UH3/qa96aVNWnoiIwhcLQASwlvU5XKrV0VR7BACU9MFArBb1aTch5qUeAIyjIkIoaLnaV1SXQbPpoXb2mIiIyHUcGQkXTQHKJaGGELb/4WUAQgjUCQ0OKJfj+4Id0On0ntlcj4iIqBUcGQky1qMlrtTtMAjjBnuOyBIQAx0GSseB3dNR/mUiVhqmYruSDbxTiBTtccwZnYE72vUNjKz3uInReOBDiYgoKHFkJJBpOhjLw+dXGX9upx9Fqkvtu+IClquXYIxcAACoqKpH7oZibDMMandfiIiITBiMhJGzIsGl9qaS8YvUb2KwXAjTouH5+ikwCAd159vDiaRWVnMlIgo9DEbCiAoKykSiS4XRZAnoKNXiPc1z+CpyNm6RC1COzihQMr3XUSIiCisMRsJItnwML+v/HwC4Vak1ucW0zVkkeLZz9tgbKWEBNSKikMNgJIyoJIHhqv9gpj4XFUh0+fqWO/0m4VcP9w4W0y6cgiEiCh9cTRNmclT70XDH67hz22B0rzuCrriAeerV6Igac7DRGlkCUnEeXeXzaPB+d9tUV1uFmBfSjT//TyliYrV+7hEREbmKIyPBxgMrbEZndsZns2/Ao6oPMFJ1CF9d8UcArk3dHBJXeL32CEdKiIjCA0dGwpRKljBYdRQAUDdhHf79dTxS985HV5x36voXGidh3zuF6BqnwXzDIOSo9nuzu0REFMI4MkIAgGvHTENM7n58Z0iDTkQ43OlXEcA5EY+uuIAb5CKcq6nHDH0ua48QEZHbODJCZqqICPRVnQYACEjGnX1bhCWKACQAnaVqLNW8CgAoE4lYoJ+C+fop+CxQysXrLgLPNRV4e7rMIwXjiIjIezgyEuTc3kyvDbrxr0HEJlucs5ffmowLeFW9FNMjtmLnZ5/A0Oj/3A7mmhARBReOjAQ7U0Krhxl6j4Xh8lsQ81IPHFAuRw/pjN0VN6bj+yK2AYXbcKYwDyXXzYFi6IOzSIC2pBID0rnChYiIHGMwQo417fQbLXToJNc4dUkXcR5dDj6JGYZc8wZ7THIlIqLWcJqG2tRHLnW6rWmk5K/qFbhD/opJrkRE1CaOjFCbnCmGZt2+MyyTXOfrp2COfjpeOlmJm/vGe6GXjtXpGpE1bzsAoGjBGI/m1hARUftxZIScosSmQNhNYW2bcU+bpbheLsa97xbixud3YkfxOQ/3sBX29rOx3vfGiR2DiYjIOxiMkFP0I/MhAW4FJC33tJGhoKKqHrkbijltQ0REABiMhCZNB2TUr0FG/RpzjY32LgE29B4LTFxls9zXWbIEpErncY9qGyQoAID5+ikwCPdGW3yhTteIjDlbkDFni0tLhN29jogoXDEYIedljUf9A3sAACVKVwjh2n42ADBP/Q6+ipyNMfI+dJfP4keRArl0L6AYvNBhx1iLhIgocDCTj1zTtNw3Qz6DggGL0f3gIqf3szEx5pC8DMk0KLJ+EhCfCtXwPA931vO4SzARkedxZITc1m/EpKb9bNJxQLkcOk1Hp3JK7K7OqS6HZtNDnu9kAOC0DRFR6zgyEi6sK7XqPFO11bifjbEOSUPO65A2PQQByWJPG3skm4DEsr1cuhdKt2yP9NFtTatwAKBOVwpoOApCROQNHBkhj2lvkmvL+CRq/SREvTHE9h6KwF5DH3xsGIyCkkoY/Lw5H0c9iIjajyMj5FlZ41Hf7UbEvNSj3R8l1VZYHG8rLEfex4U4o3/WeCIQSs3bGz2xOcddg4mIWsNgJATFaCJQsug2/3WgKckVgFNTNo60vO7br7Zi1q5YGKwG887W6DADuViOJbjJvd4CMI5wxLT8WdOODyMiIpdwmoa8zjqpVbgRm9y47yH8X+RsjJELrD7b+H/n66f4fcqGiIjcw2CEvEo3/jW7OSTuBCTGJcFL7AQkMsrRGQdLPZOUS0REvsVpGvIqQ++xMFx+C2Je6oEDyhV4oXESOqIGz6pXIxUXXPosWTIWWfuregWi9DqcQSIKlEwoTTH1L7V6b3yF9rOTV8LN+4iImvFvwHDVVDIeAIq8nWDZlEOiFxHYp2QBALY3DEK2XIyuuIB56tXoiBqndge2vyPwVGxXstElVu21r0BERN7DYIR8JlsuRtc4Dc7W6KBANgcm9XoNlquXQBEOCqK1wjR1s84wHNnSQ4DSzQs9JyIib2LOCPmMShJ4ekxPAJY1RbYr2Zihz0UFEl3+TFkyvn4XsQsx7/8WYkk/qL7/1EM99h3WKyGicMZgJEyZlv+WLLrNp/kKozM7Y/nd1yEpznLt7HYlGzc2vIzf6Z7GryLW5Q34TER1GTSbHsQxw2UBUxiNiIhax2CEHGoZpHgyYMnpl4LPHxmE99R/xsOqD83nFcjYq/TDHP19xmM3YggZxpU6cXI9HtPPwj3vFGLUP/Zjm2GQh3rvJU1JriVRkwHdRfM55GuNL9M5IqIQxGCEzGJitcb9a/KrvL4brUqWMFh1FI9FfICucRqPTdsAxmmbVOk8suViAE2F0fS5gRWQ2As+iIjCFIMRcsy0uV5+lddKmreWR3Jjw8v4re4ZPKqbifMizuWRkiRUAmBhNCKiQMdghFzijakbR3kkphU3Hys34mn9dOM5F2KJy+WfcINchAg04nq5GIPk7/Hjgc8BxeCRfvuas0muTIYlomDDpb0UEHL6pWBoejQKF+dghzIAKwy3QgLMu9NsV7IxU5+LeepVbRZLEwKQJGB2xEbMxkYYhASV1PRJuwEcfhaqm581t5dL90Lplu2V70VERG3jyAgFDFMeyTz1O1j6m0ybkZJ/x92E/bdvh0EA7+mHQwjbkRJ7ZeZl6436qsug+WSG+TBq/SREvTHEU1/DY1qOatTpGv2eZ8IRFyLyFo6MUEAandkZI3t3QuHiHJxFArSTXsewrHQ01FVDtQXoKF/EDAcjJZJV4TTrY3uk2gqbcwZFYK+hj/H+JZUYlhXfnq/kNyw9T0SBjn8rUcAyjZQAQF1GAlQtyrPmqPaj4Y7Xcee2wehedwRDpELMVm90KvAAAOtmEoR5ekc+9TW2KVcjb9NRnNE3Tee8U4gU7XHMGZ2BO9r4bOsgZkC6F1cm6S4Cz6Uaf366zPlEY3evIyLyAgYjFLSaR0/yoZbaP21gCmSi3v8drhaJuEY/FdvRnEtSUVWP3A3FiFQPQo5qv93P2FZYjryPCy2CmK5xGsw3OL6mPep0jYhp+bMGthvzATYb9dlgcEJEfsScEQpqptGTgfIPHv1c0543Y+QC8zlT5sl8/RQYhO0QzI7ic5jxzrc4U6OzOB+QdU6IiAIIgxFyjRO1R7xVubUtSmwKbCdg3GOaEfqregXukL/CDXIRZCiQoKC7fBafGgbh+4IdMDQaR2QMQsJz209Yp8oCMAYxAgJz9NOx72Rglqe3SZYFE1aJyHc4TUMhQz8yH5GbHoKABKlFWGDKBXF07IgsAZ1RjaWaVwEAF0QsACBRqjU22F2AM7s74cvGbBSLdJxp0Dn6KAAyKhGPe991PvfEm6ynd4iI/IkjI9Q+mg7IqF+DjPo1fs8zMPQeC0xcBRGbbHFesfPH3N4S4LZ0RC06otbiXBdxHreqCvCD6Ob055hyTzhtQ0RkxJERahfT7r8BI2s86rvdiJOLb8brhvE4iwQcUHpjoPw9klCJs0hAR9TgWfXqNounWbM3miJLxloni9Rvolofg2+ULABAtlxsvl+BkmkREAkYJ5Pm66dgtHygHV/Wx6wTY5kIS0QewmCEPM+UV9IkBs0FunySQyKrcByXYZPSXMhsX1OQYLK9YRCul4vwqvplaFELuR2pJrJkHDV5T/Oc7VQOgDKRiPn6qdiuNK/MEQDK0RkFSib6u39rv6urrULMC+nGn/+nlDVMiMgtnKahkGTaJM8RBTL2Kv0wR3+f8dhDOaX2pnLsrcwxOYsEz9yYiCiIMRihkJQtFyMF51tdWxOHWmxXsjFDn4sKJHrkvpJkO51jGnXJU6+GDMXivR+Uy1BQEpgrbDxWfl53EcjXGl8ufA5X8xCFDwYjFJJUkkCeehUAe9VWja+Jqt0AjJvw3djwMn6rewaP6mbivIjz2EiJiSwBqdJ53KPaZhGQvGL4b9zzTiFG/WN/WCW0MtAgopYYjFDIylHtxxI7G+4la6Ow5DeZGCUfNJ9TIGOfkoWPlRvxtH668ZydTfjcWYXT0jz1O/gqcrbNlE2oFEazV6/EKRw9IQprzDYj77NKaDWdy6hfAwAoMq24sD7ngZ1pW9twL/IT41ROBTpZFCszTd3kWW3Cp4tMwCVdIxKsckJcZcohmaHPNSe1GlfYCMzXT8FngThl0w52S9a7gyt1iEIWR0Yo5JlKxt+h2ovsFhvutTaV85mSjWENL6PAcCUAoH7iOkTOOQFV7n/wnSENOhFht9qqMxzlkAjIKEdnHCytcnAluYOjJ0SBj8EIhbXWpnJe+k0WslXHAABK+mBAVkEVEYG+qtPQSI0wjmVYhjECcCpIMeWQZMvFxmMouEEuwnh5D3BqL6AYPPDtwoBVki0DD6LgxGka8gt3i6VZ73vjiV84rU3lYLPj63TjX4N6Zz6k2nLzOSk6EUIIoP5Xp+6dIxdgkFSMyRE7kSI1TQntB8Sx+VANz2vP1wo5HpvuIaKA4/LIyO7duzFu3DikpqZCkiRs3LixzWveffdd9O/fHzExMUhJScG9996L8+fPu9NfCic+LDXvaCqnNYbeY1H/wB7zcf3EdcCTx3Fp5rdO3/eeiM/whHoDkq2qwYrqMmg2PYhjhsvwsWFw4C7/dZenlg27cy97ybJuJtASkWe4HIxcvHgR/fv3x7Jly5xq//XXX2Pq1KmYPn06vvvuO7z//vsoKCjA/fff73JniQKOrDL/aJrKaXnOehrHfN4qrrCpTdLUppNcg52Ga/HamvUYvXQvlur/n8vBiUER2GvoE9BBjdurcAL8XkTkHJenacaOHYuxY8c63X7v3r3IyMjA7NmzAQA9evTAgw8+iOeff97VWxMFJZtdhOHersFljYmYL6bi7/q7gHcc7/5rCj7OIgE/7D6F9w+dwRn9s8Y3W7ku4FjvhePONfb2z/EkrvAh8givJ7AOHjwYp0+fxtatWyGEwJkzZ7BhwwbceuutDq9paGhAdXW1xYsoGOnGv2azi7C72+CYlgQ/E7EKN8hFOFtVZ7P777bCcoz6x378Tv8sHtU/gld2n8aZGp3F53DXYFt2R0usp27sTOVwlIXIM7wejAwdOhTvvvsuJk2aBI1Gg+TkZGi12laneRYuXAitVmt+paWlebubFABMSa0li24LmQ3XrPNK2kOWjK/7IrZhreYv+L+m4mnz9VNgEBJ2FJ/DjHe+tQk+rJnGaEzXERH5m9eDkaKiIjz66KOYN28eDh48iG3btqGkpAQPPfSQw2vmzp2Lqqoq8+v06dPe7iYFILvBiamAWn5V8AyJt8gh8aRkXMCr6iW4Wj6BF/V3In/rcadrn0hQ0F0+ix9FCuRSLiW2x1OjHnW1VeYRlbpa1pAhssfr//xcuHAhhg4diieffBIAcPXVV6NDhw4YNmwY/vKXvyAlJcXmmsjISERGRnq7a0Q+p8SmQK6tgHPVSFonS8aS9X9Vr8Cf9Xfjyvr/oACZUKz+jSFDQbZcjCRU4iwS0BE1eFa9GqmmpcTrJ0GJtf3fIVny1NLiOl0jsuZtBwAULRgTMqOARO3h9f8V1NXVISLC8jYqlfFfiqK9G30QBRn9yHxEbnrINqlVOJfUas0myVUkYoH+blQiHkmoRHepwrKGCezvryPVVrR5r5aJsdqSSgzLine9wyHOqUDD10m2REHA5WCktrYWx48fNx+fPHkShw8fRmJiItLT0zF37lz8/PPPWLXKWGZ73LhxuP/++7F8+XKMGTMG5eXlyM3NRXZ2NlJTUz33TYiaeKMwmqcYeo8FJq6C2PqURbG0M1Ii3tWPQKnoinnq1eiIGjhR6sSGMcn1ZYvAxm7wYfXZEgSEAC5Bg4KT5zGsb3zorMoJAhwtoXDn8p/4AwcOYMSIEebjxx9/HAAwbdo0rFy5EuXl5SgtbV6Gd88996CmpgavvPIKnnjiCSQkJGDkyJFc2kueYW/DvUCXNR713W5EzEs9ABiLpXVIHYQhL96Gy1GG8iHzkbj3CSgQLid1yVLbNUwcjcBIEhADHXavfRGPa26HAhmVpuBjt23elmlVTqR6EHJU+13saYjiqAeRW1wORoYPH97q9MrKlSttzj3yyCN45JFHXL0VUeiyKpamAjBYdRQAUDd0HaS0LhBbnwJajJ44y53pnpbmqd/BfcpWzNdPxXZkO2xn3GnYuCpntHygfTclS6xfQmGGY4EU1Jza48be6In1uUArAW41emIM/y3zTLzJVNNkhj4X2xVjQGKdCFugGJNly9EZBUom+vukZ+HB2WTZutoqxLyQbvz5f0oRE8uRGApODEaIAlWL0RPd+NdtNuXz6q2bVurkqVdjR8NAjJYPIE+9qnkFDiyTZdVSI+TSvVBSBzR/RuleIHO0T/ob9JyoNutM4OFu7glzVsjf+CeOwpL1iEpd63XC/M7QeywMl99iHilp0HSEuuFXt5JcAedW78gSkIrz+If6ZdwqF9iMydgky66fBCE1Z7lErZ8ExKdCdfOzzZ/JACUw2QRDEZwmIp/yetEzomBlvSrH71oWT8tZBEmSoFg1cWYSx9V98m5TFUBqqv5q0R17wYyw6lF1GTSfzDAfRq2fBCzpB9X3n7rWCXKOL3dDJvIgBiNEQcjQeyykiasAO8XKrGMN6+MKdMJM/Wws0N/d7n7YLhF2QnU5NJscV2Amz6rTNSJjzhZkzNmCOl2j5yrC2tmrx5n7E9nDYIRCj59LxvtsRCVrvMW+N/UT10G6620IqwCl5UZ99RPX4T+/2Y1vO9yElYYclIlEl0dKXGUboDh3Q1Odk48Ng1FQUglDY/MvMlMJe5s23v4y1DonAxQiawEw9kwUAIKxXglgs0QYsVrUp91kUcNE6ZZtPlbSByMnVouh3WNRuDgHPymdkaK6AAW+/ZeJzaqgFnvjyKV7sc3QD3mbjpqLrI1ZswJ9NKsR19Qmav0kXIpOxl90d2Ndi0JsXeM0mG9g3RMLTiTHeordVUBBUHuFCbz+xydO5IhphCXYWAcodqhkyVzXpKFppY47NU3a64ByBaK+3oCsI4vN56LWT8LVIhHXNNU5GSMXYLl6SXNhkyaRdRVYiBeQpLoTp0QKziIB+2syMQO5WI4luMnn3ybEWdU+8dRePabrLYIBNLSZQMsAIrTwvx4RnKxX4uA6ez+72safrFfqOMvd/XRa+tLQH4/tedxmLsdU52Smfjbmqd8BYD+BVgjgCfUH5nNlIhHz9VMwRz8dL52sxM19uX+Op1gHHwHJulgcwFVBQSLw/mYkCmAxsVrzaEkMEDrz4i1HU1zYWdg6IGnr2HwPAVQgEZMjdhpvb92dpjonf1GvRGep2uH9rT/bGMQsxQx9Lu59Nx4p2uN4alR3JLfY4O+abvH4Nzf8a5UnRz2cYjWVUwf49v7kd0xgJfI0PyfQtpd+ZD4AQFgNV7QMTb664Q3M1D+KciRatFHs/JVivXuEKcf0vcYRSJEuOKyVIktoNRBxdA0ALFK/icFyIc5U1eGJD45iqeFOfGEYgNfWrEf287stjm9atAM7jp5p/oym5FjyDo+t5nH3/s6s7nEiEZerhDyLIyMU8tydgglXjnYWFnEpkGqMx9fdmIPabjfhzo1D0L3uiLlE/AGlNwbK35uPO6IGz6pXIxXNlVsr0AkL9L9HluydZEpZAjqiFu9pnsMFEQsASJRqm7+fkKCSmiOkCw2xwGaYp4qi1k+CiE+Fanhe2zezSrxVujney4d8z+cjPOQ2BiNE7RGsSa5tsbOzcMtVOQCQ0y8FQ9OjUbg4HzuUAdhkGAIJwD4ly+KjtjcMstjTpiNqME+92qK0vLd0RK3NOdlq+sleG1Fd1nYtlKJNiNr6lPkwav0kKHbqvpBnWJfDB+DeKh17q3tszvFXo69xmobIDwI9qRWA7aqclhVgm5hW5cxTv4Olv8lEUpztPz0VyNinZGGTMgRa1GKZeimS4f1ABDDmlNgUZrNzbH1OBiCEQIMw/reRS/cCjc17BkTsWQKsn2qzV5BUW+GhnpM7Wk6XtGfqxFOf0577h9sUUID+LUhEwWZ0ZmeM7N0JhYtzcBYJ+GHw37Bs92kAxnwTGQry1KsAOCgl7yJPrOZpjSwBkTD+Ioiy2ndHs+cl65XGAIz1U8z9aprCMRVmMyXMDkgPvDobRP7GYITI20J1KseOlvVL6m7qjn7pXZD3cSHO1OiQLRd7bGpmS2M2xqoKIIRnAhunKIpF9OHotqYASf6pANuU/sbvz8Js4cl6qbGmA+ujOMCnQERe05xXkgO15Lnh5tXKLfhEGYJ56lUWybHe5OoozI4DRZhdZLBZIH22RsfCbCHKOq/Fm4FGqAU1wd17olARrOXondBytMQZ50UcOqLG7oiHsT5JJxQomeiqjcH+Ub/HHVsGAbAp0Op37xYbzNNTLRN4jSuOirHLcDWSCnag/00T/N1V8iU3y+NbBx/BUGbfFQxGiAKUdZKru4lsgZQs21pBNVOg8Wf93XhVsxQCksUeNqbjn5ROWHH31RiWlY6Gumpgi6mF1BSQtLzG9E4ze4XZAM/lnygCqEIsJCgYK39jXNrcYnrKYmnx7rdwZvdTODVgDrgoOLB5LRiwM5XjqJ1H9hgK0Cq1XE1DRD7juKCaBAnAT0onTJo8Hcpdqyx2GwYAEZcMSQKyVd8jOyMBKquhE93412yukaITIaISLDshW/619yti8Sti3f1Kln1symHpKBnrnLyqXooUq2kk66XFXcR5DDzwJDY3Zre6QzH5WVMwUBI12enKy86syrHbxo17Ofrstlbl+HvlkAmDESJ/CPIqre4yF1RrK9DoewfqH9hjfr9+4jrU37/H+uNsPtv6Gjx5HPUzD1mcu/To9xbHBb/5Brep38Jvdc/gUd1MnBdx5iqxbbGuLmtPW0uLTTHVdaof8Zh+Ft5bswJ1SwaZ349aPwlY0g9S8WZ8Z0jHAeUKfF+wwyJgMTGt3DEHNc5+ESI/4zQNUZAIpOmWdnGioBoAp3YftmGvNkorn6OkD0ZOrBZDu8eicHE+KkRHPCc/gMXKi1CsVupYT++Yfs//Xf/f+EPEZ0hAbZuBh8NuS0AqzmOWaiMei9hgkwAjqssQtXkm+pq+yu7pxumd654yT+98+9U2PHlAi/IWK3dStMcxZ3QG7nCuGxTkgjmvJIj/RiMKHWFXst6dQMOLWibZyre9gpkfGmxW6iiQoYJiPq5AJ8zXT0EVYvG49KFH+vFQxCeQYGc0BU2jMC3OdxHnkXTwj+ZzN+57AB+IRMyXp2K7YgxRKqrqkbuhGJHq5qXE1nVPhmXFsxZKqPBUXokfMBghChZO1isJmREUPxmd2RmRkx/CnRsHW+y7863ojeuk5n13Tsf2x//c3hNZm8Z77N4dpAaH79mb3rGeJjLuWrwEM/S52K5km7NTntbfi0tCjR93n8IH35ahu+FO4/dYsx6PR/aDAhmVrtZCsdqXB5mj7QY6FFis9+sJFPybiihQhVGxtEDTct+ds0iAdtLruKZbPI6+uNB8bFrNE7Pl53bfz91qsvYCFEUYdy2u1sfgGyULCmRcgBaPNT6MMV99gA/Uq5CqaR7xKVMSsUB/NyrleHOgtb8ms9VaKKrvP4V6Z775OGr9JFyKTsZfdHdjndU00VOjuiOZoy7UBgYjREHM3VGQmFitOdCJAfyylXugs6gmm5EAABbH1qt5rJciN59ve2mxJ8vat9y1uEwkYr7eOG0zRi7AcvUSm/bG0ZSXLfpQJjpibeNI7DFk4aZTe4A+t1hcY9xE0PK7RtZVYCFeQKWca54mKq+qx2MfHAPQ+qgLp4mIwQhRMOPoSUCxrY1iS5Esc0+8ub+OadpmSaMxydZePoq96Z4U/IrH1R8YD95/GyI+FarheS1aWC/Obh6VyVOvxo6GgVCaFmtaF32zHnXZVljOkvnEYIQo5NgLUBiweJ1u/GtQ78y32MlXir8MDTc/g8hPZgAwrhzSJ1+H0r//F6IlHXrKzu3y6+7oiSnQeFzdeoJtW6uARHUZNJsebH6/lful4jxyIzZgj9IPHVFjU/StTCRigX4K5uun4Kmic3j8w2KboO1cTT1WyrcgRTqPK5vyUbyJIzP+x2CEiMgDDL3HwnD5LRZLlqMyR8NQVwt8YmyjpA+GCkBflesrHdyd3vHEyIsM29U8rZkdsRGzsdFuHZZkXMCr6qVYYcjBpi3nIKE3RIuSV2PkAuSpVzUHMOsn2RmZccAqqVZJHWBxbC+o4chMYGAwQkR2c088VY4+rNirc9KKtsrjVyEWM/WzkYCLxhEGq2quns43aY2793FU5O2+iG24D9tQFtmUQIt4jJIOYnrEpzZPwzgy85DFOZuVO9IuRH36R/P7UesnQUiyxTGsgpodxeeQu8F2ZIabGfoegxEiIj/Rj8xH5KaHbHJNTAXV5ujvw16lHwBge8Mg/Jd8EPfJW/BP9e+g1lXaDVACiTMBjL0EWpt8FACKEKiHBtGSDjuKz+G57SfMoxlj1qzAzZolth8uFMvj6nJoNj2EY4bLUCS6Y8HW43bzeiQoyJaLuJmhDzEYISKnEmGtV+BQ+5nL4299yiLXxFRQzbQqRQIgIONO1Ve4XvU9Km67A7kbivFZwyBky0V4Vf0ytKi1u9NxW/y927G9BFpH7aKhwz7DlcjdUAwJCm6Qi9EVFzBPvdq8L1BLtt9LQACIk+vxRMNMKHZK6ttMEznYzJA1VTyLwQgR2efMSh1NB2TUrwEAFDXtsWN9TG2wKo//1Q1vGMu6NzTnPyRrozBndAZyNhtzGEZndsbyu69D3seF2FvTD3P092G5egkUuLfhmL8DElemgEpEMh5WfYjJETuRIrk+KiQDSJWak2wLlEzzyh9Hy5+7iPPocuBJHFJ64VrVjzYjMyy9334MRojIY+yWtbcTsJCVFrkl192Yg8+GCBQuzrEpsIbNzZc0F2Yztisa8hKyjiwGWoywtKX+rvcQBR2UT+dYjMz4MhfFVb+N+BKAc6MprTEl2RpX99yNKsRikfpNh8ufFQGkyOextXGg3TwTe6X3yXkMRojIbT7dUyeMghrrgmvWBdbsthu6DvWD/595hEU35HFo9vy9adSjZe2T5vwUpftQIFaL+rSbEPNSD/xb6YnPDdcZRx1a5KI4Ck4UAVQgERKArrhgd5rIU6Mu1n3wVLBkL2fFHlkCklGJ95XhdvNMTN9zvn4KRssHPNO5MMJghIiCFlf8WGkxwtI4JBeabtfa5KOIuGRINeV2r+svn0DJ+A347213mffl6S5V4LGIDTY5GaYk2/n6qQBgnCayk7fhCd4cqXE2Z8UkDvUO3xMAytEZf2+8EwNLKjGgW2zzfRzs32Ovpokz+SihlrPCYISIKFRZ5aPUT1wHpVu2+die0ZmdMbJ3JxQuzkeF6IgFmlz8UN/NZhdj6yTbGfpcY+Jny9U98ZdB11T0rUTpiu7SGQhYBizOBBruBiLeqMVyFglttnnF8N8Ys2YF+mhWI67pnGn/nr82TEZJi40KT8VchbwWNU3s1T2xzkcJxZwVBiNE5FU+ncohW9a1T5zQcvpHvvUV5G5oxI6GgRjUoqx7y8RPAPhMycbnDQPxnvrPyFYdsyn6liGfQcGAxeh+cBG64nzzzSTbqRwFAASwS7kaI1X/cfurK7Asvd8eximpTihQMgHYlrk3nc+Wi5vrpVh9MeP+PS9B0jSfK9MnYoHBOLrU4KDuSUVVPR7dUIwS1f9DFTpgxYZiyE2ricxl9qsynctZsS4M1y3bcVsfYjBCRH5nb9mwzRSMvQutV/zo2i577+7mguGq5cqdfTVZ5vOycb2xmWnFT/bmYwDsF33rN2ISDMPuxHcvDsclKRIRw+fgqiQVpG1zLaaSEJeCHypllIhkAM4HI4oALiAef9bfjTNIxAGlNwbK32OIVIjZ6o1ufPvmzwWA9xqH43Z5H7pLFTareS4I45RMolRrPmdv/x7rKSFjRdol+B/9g/iXg7onAsbgZ6/oiyRU4hE7q4mMibhTLXJWnCkMp3RIbu6fqUptG8X6vIH/SySioODMCAtroXiH9cod7aTXcU23eBx9cWyrK37sUUVEmMvh12WPhipWi7r0m22mkq58qQcuGOJQJhKR7CA5tiVTwPAn/b3mqSMA2KdkoQCZ+E3Ebqc+x55KGAONJ0ybB8I2qOiIWjjD0UqdZ9Xvory+E76BMeBrOepib48fR0HNDH0u/t54J7D7FN4/dKbNwnDSxeb9kaKaSu9LOc8DWeOd+j6ewmCEiAKTt3YktvO5HC1pm/UKHwBOrfhxioOppGz5GJ6OeAILG1+wSY61zgexzmFpSYGM+fqpeE2zxGZKyNFqHyGAS1DjtcZxyI2w3Wiwrc0FXSFLxmDmPc1zdkdY7CXYOgpq8tSrcWPDUii7T5unclwpDCeqy4D1UyFNXOXTgIT/qyOi0ObMkmB3Ax/rz9ZddLeXZIdKEhgydjJmfthok0BbjkSs1Y/Af6kOoWbIHEzbFWORw9JMYLrqUwz73VNQpOsgffpHy52VYRuQCEiQJIFIocfv1V8CXlolZI+jERZngh3TrsnZcjG0qLWsJOskc+n9T55EdOZtPpuyYTBCRD7ny5EIJtAGt9GZnRE5+SHcuXGwebnxWSSgtMNVmKe8gP7yCdTdmINl3eqMq1BqdOZrk+M1yL/0PHJU+1GX8RxUsXegLm2Y5c7Koh7Cuuhb0/JnlQRjsq0PC8DZCzpcHXV5QPUJRsj/tpt/4gxZAqIvVcBQ8jVUPX2zVSCDESLyPW9NwXgR81H8pzlnJd+cozIgXYu4l/bbaZPjsA0A2ymhFkXfAOeWP7vLV5VtR6r+DaD9MdSPJ35EbwYjRER+FEYVX4OBvZwVd9rY5cbyZ3d5MyDx9GefFQno7bmPaxWDESIiO5xavWPTxnb6iVVig5sSmwK5tgJwYtLDpl5K0yUrDDn4XBnYvCoGrm/w1+a9PRiImGqqqDKGeuYDneDOBo9ERERhQT8yH4AxqbUlm9AkOhGI6mhx6qzUCTP0ufhL41TsU7JwOO5m7L99OwBAJ1TmYMUZ1k2tj9tTpbYlU59eVk9Hdq8u7n2oGxiMEBEROWDoPRaYuAoiNtnifMvj+onrgCeP49LMby3Odcjdj3tU27FU/Q+svLsfvvrjSIzu0xUAoJEMkCTJpj6sgDNjMLDpj7MUAZwT8XhUNxMv6n+DCiRavF+BTpipz8XwCfe2b7m2izhNQ0TkKV5MzGUCrR+1scePudqsVe6JCq3XYtGNfw3qnflAy+qz0YnG4Yr6X5vPtdjjx979nWWvMNwywwSLAmunY/vj2buuQk6/FJc+u70YjBAR+RITY4OTF5JcDb3HwnD5LZZLjTNHo662yuacaY8fe/d3Nq+lAp2wwKownAIZ+5Q+mK76FP/9uwcwLCvdpyMiJgxGiIh8yKm6J24GLM7Ub2G12QBjHeTYGWFpq/CYfmQ+Ijc9ZCzW1iIgaZlQWz9xHf5j6IdDm44CrdVi8UMgAjAYISLyLiembpzaFDAIMNDxD3Ney9anrIq3pUCqMR4r6YORE6vF0O6xbddi8QP+aSEi8jeb3YedXP5rM4LSYNPEJkDwVF4Lp5sCSxt5LSZu12LxMgYjREREocCHxds8jUt7iYjCnPXoiSenW5jHQs7gf3UiolDRjt2HLaeJuPsw+RaDESKiUBZgmxIGQ70UjtT4Hp8yEVGAcWr5rwvtXGYvgHE3YdWNYMjXwQCDD//jUyciosDSjukm64CJK36CA4MRIiJqk/UojN3dhz2Vs+LkNb4MNLj7sncxGCEiIpd5bYrIg+z10akpmQDLswkHDEaIiCjoBEMwRM5jMEJERIGPFV9DGoMRIiIKHw5WCnlixQ/zStzHYISIiAKeM9MynLoJXgxGiIiIvMHO1JKnRk9CrTaKy3vT7N69G+PGjUNqaiokScLGjRtbbX/PPfdAkiSbV9++fd3tMxERkW+ZpnLyq4w/Wx/bYRqpKVl0W3PA4OR19n721PfIqF9jDJICKPfG5WDk4sWL6N+/P5YtW+ZU+6VLl6K8vNz8On36NBITE3HXXXe53FkiIqKA5USgYc3ZwCPUNxx0ubdjx47F2LFjnW6v1Wqh1WrNxxs3bsSvv/6KP/zhD67emoiIKDzZS7JtI+k2GPYBMvF56PTWW29h1KhR6N69u8M2DQ0NaGhoMB9XV1f7omtERES+5ecCa4ESsPg0GCkrK8Onn36KNWvWtNpu4cKFmD9/vo96RUREFEC8GKAE6oojl3NG2uPtt99GQkICJkyY0Gq7uXPnoqqqyvw6ffq0bzpIREQUrNzIWQkUPhsZEUJgxYoVmDJlCjQaTattIyMjERkZ6aOeERERkT/5bGTkyy+/xPHjxzF9+nRf3ZKIiIiCgMsjI7W1tTh+/Lj5+OTJkzh8+DASExORnp6OuXPn4ueff8aqVassrnvrrbdw/fXXo1+/fu3vNREREYUMl4ORAwcOYMSIEebjxx9/HAAwbdo0rFy5EuXl5SgtLbW4pqqqCh988AGWLl3azu4SERFRqJGEEMLfnWhLdXU1tFotqqqqEB8f7+/uEBERkROc/f3t09U0RERERNYYjBAREZFfMRghIiIiv2IwQkRERH7FYISIiIj8isEIERER+RWDESIiIvIrBiNERETkVwxGiIiIyK8YjBAREZFfubw3jT+YKtZXV1f7uSdERETkLNPv7bZ2ngmKYKSmpgYAkJaW5ueeEBERkatqamqg1Wodvh8UG+UpioKysjLExcVBkiS3P6e6uhppaWk4ffo0N9zzMj5r3+Gz9h0+a9/hs/Ydbz5rIQRqamqQmpoKWXacGRIUIyOyLKNbt24e+7z4+Hj+4fYRPmvf4bP2HT5r3+Gz9h1vPevWRkRMmMBKREREfsVghIiIiPwqrIKRyMhI5OXlITIy0t9dCXl81r7DZ+07fNa+w2ftO4HwrIMigZWIiIhCV1iNjBAREVHgYTBCREREfsVghIiIiPyKwQgRERH5VcgFI8uWLUNGRgaioqJw/fXXo6CgoNX277//PjIzMxEVFYWrrroKW7du9VFPg58rz/qf//wnhg0bho4dO6Jjx44YNWpUm/9tqJmrf65N1q5dC0mSMGHCBO92MIS4+qwrKysxa9YspKSkIDIyEr179+bfI05y9VkvWbIEV155JaKjo5GWlobHHnsM9fX1Pupt8Nq9ezfGjRuH1NRUSJKEjRs3tnnNrl27cN111yEyMhKXX345Vq5c6d1OihCydu1aodFoxIoVK8R3330n7r//fpGQkCDOnDljt/3XX38tVCqV+Nvf/iaKiorEM888I9RqtThy5IiPex58XH3WkydPFsuWLROHDh0SR48eFffcc4/QarXip59+8nHPg4+rz9rk5MmT4rLLLhPDhg0Td9xxh286G+RcfdYNDQ1i4MCB4tZbbxVfffWVOHnypNi1a5c4fPiwj3sefFx91u+++66IjIwU7777rjh58qTYvn27SElJEY899piPex58tm7dKv70pz+JDz/8UAAQH330UavtT5w4IWJiYsTjjz8uioqKxD/+8Q+hUqnEtm3bvNbHkApGsrOzxaxZs8zHBoNBpKamioULF9ptP3HiRHHbbbdZnLv++uvFgw8+6NV+hgJXn7W1xsZGERcXJ95++21vdTFkuPOsGxsbxZAhQ8Sbb74ppk2bxmDESa4+6+XLl4uePXsKnU7nqy6GDFef9axZs8TIkSMtzj3++ONi6NChXu1nqHEmGHnqqadE3759Lc5NmjRJjBkzxmv9CplpGp1Oh4MHD2LUqFHmc7IsY9SoUdi7d6/da/bu3WvRHgDGjBnjsD0ZufOsrdXV1UGv1yMxMdFb3QwJ7j7rBQsWICkpCdOnT/dFN0OCO89606ZNGDx4MGbNmoWuXbuiX79+eO6552AwGHzV7aDkzrMeMmQIDh48aJ7KOXHiBLZu3Ypbb73VJ30OJ/743RgUG+U549y5czAYDOjatavF+a5du6K4uNjuNRUVFXbbV1RUeK2focCdZ23tj3/8I1JTU23+wJMld571V199hbfeeguHDx/2QQ9DhzvP+sSJE9i5cyd+//vfY+vWrTh+/DhmzpwJvV6PvLw8X3Q7KLnzrCdPnoxz587hxhtvhBACjY2NeOihh/D000/7osthxdHvxurqaly6dAnR0dEev2fIjIxQ8Fi0aBHWrl2Ljz76CFFRUf7uTkipqanBlClT8M9//hOdO3f2d3dCnqIoSEpKwhtvvIEBAwZg0qRJ+NOf/oTXXnvN310LObt27cJzzz2HV199Fd9++y0+/PBDbNmyBX/+85/93TXygJAZGencuTNUKhXOnDljcf7MmTNITk62e01ycrJL7cnInWdt8sILL2DRokX4/PPPcfXVV3uzmyHB1Wf9448/oqSkBOPGjTOfUxQFABAREYFjx46hV69e3u10kHLnz3VKSgrUajVUKpX5XJ8+fVBRUQGdTgeNRuPVPgcrd571s88+iylTpuC+++4DAFx11VW4ePEiHnjgAfzpT3+CLPPf1p7i6HdjfHy8V0ZFgBAaGdFoNBgwYAC++OIL8zlFUfDFF19g8ODBdq8ZPHiwRXsA2LFjh8P2ZOTOswaAv/3tb/jzn/+Mbdu2YeDAgb7oatBz9VlnZmbiyJEjOHz4sPk1fvx4jBgxAocPH0ZaWpovux9U3PlzPXToUBw/ftwc8AHA999/j5SUFAYirXDnWdfV1dkEHKYgUHCLNY/yy+9Gr6XG+sHatWtFZGSkWLlypSgqKhIPPPCASEhIEBUVFUIIIaZMmSLmzJljbv/111+LiIgI8cILL4ijR4+KvLw8Lu11kqvPetGiRUKj0YgNGzaI8vJy86umpsZfXyFouPqsrXE1jfNcfdalpaUiLi5OPPzww+LYsWNi8+bNIikpSfzlL3/x11cIGq4+67y8PBEXFyfee+89ceLECfHZZ5+JXr16iYkTJ/rrKwSNmpoacejQIXHo0CEBQLz00kvi0KFD4tSpU0IIIebMmSOmTJlibm9a2vvkk0+Ko0ePimXLlnFpr6v+8Y9/iPT0dKHRaER2drbYt2+f+b2bb75ZTJs2zaL9+vXrRe/evYVGoxF9+/YVW7Zs8XGPg5crz7p79+4CgM0rLy/P9x0PQq7+uW6JwYhrXH3We/bsEddff72IjIwUPXv2FH/9619FY2Ojj3sdnFx51nq9XuTn54tevXqJqKgokZaWJmbOnCl+/fVX33c8yPzrX/+y+/ev6flOmzZN3HzzzTbXXHPNNUKj0YiePXuK//3f//VqHyUhOL5FRERE/hMyOSNEREQUnBiMEBERkV8xGCEiIiK/YjBCREREfsVghIiIiPyKwQgRERH5FYMRIiIi8isGI0RERORXDEaIiIjIrxiMEBERkV8xGCEiIiK/YjBCREREfvX/Af6jSdLJxXMoAAAAAElFTkSuQmCC" }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.errorbar(c_sticking_probability.p, c_sticking_probability.fd_mean, c_sticking_probability.fd_std, fmt=\"o\")\n", "plt.errorbar(rust_sticking_probability.p, rust_sticking_probability.fd_mean, rust_sticking_probability.fd_std, fmt=\"o\")\n", "plt.legend([\"C\", \"Rust\"])" ], "metadata": { "collapsed": false } }, { "cell_type": "markdown", "source": [ "Hmmm even with a different RNG I still get a dip around 0.6, I think theres something hiding there" ], "metadata": { "collapsed": false } } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.6" } }, "nbformat": 4, "nbformat_minor": 0 }