Skip to content

Instantly share code, notes, and snippets.

@maedoc
Created October 12, 2023 11:58
Show Gist options
  • Save maedoc/add7c3206f81d59105753a04f7c1fcf4 to your computer and use it in GitHub Desktop.
Save maedoc/add7c3206f81d59105753a04f7c1fcf4 to your computer and use it in GitHub Desktop.
comparing point vs extended geometry for source sensor forward model
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "markdown",
"id": "e855bcfe-15e1-4b59-8628-d6f684b00ff4",
"metadata": {},
"source": [
"# gain matrix for less trivial geometry\n",
"\n",
"the typical gain matrix we used for non-oriented sEEG is \n",
"\n",
"$$ G_{ij} = 1/|xyz_{s_i} - xyz_{c_j}|^2 $$\n",
"\n",
"for source j and sensor i. this assume point geometry, but that's less reasonable for neural fields. as a simple approach, we can write the integral with both non-point geometries, considering just distances along the vector from source to contact,"
]
},
{
"cell_type": "markdown",
"id": "a8ffd7f4-1a67-459b-8e92-ffb5510fcd64",
"metadata": {},
"source": [
"$$ G_{ij} = \\int_{-D_s}^{D_s} \\int_{-D_c}^{D_c} ds dc / |c - (d_{ij} - s)|^2 $$\n",
"\n",
"where $D_s$ is the diameter of the sensor and $D_c$ is the diameter of the contact and $$ d_{ij} = |xyz_{s_i} - xyz_{c_j}| $$"
]
},
{
"cell_type": "code",
"execution_count": 1,
"id": "5305343b-eca8-48a9-b4e2-28aa1b992740",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"%pylab is deprecated, use %matplotlib inline and import the required libraries.\n",
"Populating the interactive namespace from numpy and matplotlib\n"
]
}
],
"source": [
"%pylab inline\n",
"import numpy as np"
]
},
{
"cell_type": "code",
"execution_count": 49,
"id": "1350e3c1-53ee-4bc2-b157-35977e20644e",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"0.1016131519696175"
]
},
"execution_count": 49,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"def f(d, c, s, nug=1e-3, att=2):\n",
" # att 2 is parabolic attentuation, for nfm use 1\n",
" k = np.abs(c - (d - s))**att\n",
" return np.reciprocal(nug + k)\n",
"\n",
"def G(d, c, s, nug=1e-3, att=2):\n",
" f_ = f(d, c, s, nug=nug, att=att)\n",
" return np.trapz(np.trapz(f_, c[:,0], axis=-2), s[0], axis=-1)\n",
"\n",
"dij = 25.0 # mm\n",
"Ds = 3.0\n",
"Dc = 5.0\n",
"ng = 64\n",
"\n",
"c, s = np.mgrid[-Dc:Dc:1j*ng, -Ds:Ds:1j*(ng+1)]\n",
"G(dij, c, s)"
]
},
{
"cell_type": "code",
"execution_count": 62,
"id": "c588167e-368e-404b-9db2-7cec432edcee",
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAxcAAAGGCAYAAAAJhXzzAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy81sbWrAAAACXBIWXMAAA9hAAAPYQGoP6dpAAC5tElEQVR4nOzdd1hTdxcH8O9NGGGDgIBsBSxYxIW7iqtUrVato/Wtq4i2YrWuOurscLTW1VLrqlpbrW21aqt17z1x4QAEEUUUgbBX8nv/wFwJSSBAQhI4n+fxkdzce3MSwj05+S2OMcZACCGEEEIIIdUk0HUAhBBCCCGEkNqBigtCCCGEEEKIRlBxQQghhBBCCNEIKi4IIYQQQgghGkHFBSGEEEIIIUQjqLgghBBCCCGEaAQVF4QQQgghhBCNoOKCEEIIIYQQohFUXBBCCCGEEEI0gooLotT8+fPBcRxSU1N1HUqVaCP+kJAQhISE8LcTEhLAcRw2bdqksccglTNy5Eh4eXnpOgxCiJ7y8vLC22+/rbHzKbvuy/IN0R2O4zB//nxdh0FeouKCEAOSm5uL+fPn4/jx4wr37du3r8YurtHR0Zg/fz4SEhK0/lhPnjzB/PnzERUVpfXHIoQQQ3D27FnMnz8fGRkZCvctXLgQu3btqpE4tm7dihUrVtTIY9VkjiPVQ8UFIVXk6emJvLw8DBs2rMYeMzc3FwsWLFBZXCxYsKBG4oiOjsaCBQtqrLhYsGCB0uJi3bp1uHfvntZjIIQQVWbPno28vLwafcyzZ89iwYIFda64UJXj8vLyMHv27BqJg1SMigtSq+Tk5NTYY3EcB5FIBKFQWGOPSeQZGxvD1NRU12EQQmoAY6zGP8Srw8jICCKRSNdh1GkikQhGRka6DoO8RMUFKVdGRgZGjhwJW1tb2NjYYNSoUcjNzeXv79y5M4KCgpQe27hxY4SGhgJ41U916dKlWL58OTw9PWFmZobOnTvj1q1bCsfevXsXAwcORL169SASidCqVSvs2bNHbp9NmzaB4zicOHEC48aNQ/369eHm5ia3T2pqKgYPHgxra2vY29tj4sSJyM/Pl9unuLgYX375JRo1agRTU1N4eXlh1qxZKCgoKPe1UTXm4u7duxg8eDAcHR1hZmaGxo0b4/PPPy/3XIWFhZg7dy5atmwJGxsbWFhY4I033sCxY8fkHs/R0REAsGDBAnAcx/czHTlyJCIjIwGA3166D7BUKsWKFSvQpEkTiEQiODk5YezYsUhPT5eLQ9Y/+fTp02jdujVEIhEaNmyIX375Re51HzRoEACgS5cu/GPJWlNU9X318vLCyJEj+dtpaWmYOnUqAgMDYWlpCWtra/Ts2RPXr1/n9zl+/DiCg4MBAKNGjeIfS/aaKxtzkZOTgylTpsDd3R2mpqZo3Lgxli5dCsaY3H4cx2H8+PHYtWsXXn/9dZiamqJJkybYv3+/ql8TIUQF2bgD2fWvvGvuxo0b0bVrV9SvXx+mpqYICAjA6tWrFc4pux4dOHAArVq1gpmZGdasWVOpc8gcPHgQzZo1g0gkQkBAAHbu3Kmwz4MHDzBo0CDUq1cP5ubmaNu2Lfbu3av2cy/r119/RevWrWFubg47Ozt06tQJBw8eLPdcN27cwMiRI9GwYUOIRCI4Ozvjww8/xIsXL+Qeb9q0aQAAb29v/rooy0k5OTnYvHkzv730dffx48f48MMP4eTkxF/zfv75Z7kYjh8/Do7j8Mcff+Drr7+Gm5sbRCIRunXrhtjYWH6/kJAQ7N27Fw8fPuQfS3Y9luXnsq3bsnOXbn0/deoUBg0aBA8PD5iamsLd3R2TJk2SKyQrynHK8s61a9fQs2dPWFtbw9LSEt26dcP58+fl9pHFeebMGUyePBmOjo6wsLBA//798fz583J/V0Q1KvNIuQYPHgxvb28sWrQIV69exfr161G/fn0sWbIEADBs2DCEh4fj1q1beP311/njLl26hPv37ys0U/7yyy/IyspCREQE8vPzsXLlSnTt2hU3b96Ek5MTAOD27dvo0KEDXF1dMWPGDFhYWOCPP/5Av379sGPHDvTv31/unOPGjYOjoyPmzp2r0HIxePBgeHl5YdGiRTh//jxWrVqF9PR0uQ/Lo0ePxubNmzFw4EBMmTIFFy5cwKJFi3Dnzh38/ffflXq9bty4gTfeeAPGxsYYM2YMvLy8EBcXh3/++Qdff/21yuMyMzOxfv16vP/++wgPD0dWVhY2bNiA0NBQXLx4Ec2aNYOjoyNWr16Njz/+GP3798eAAQMAAE2bNkVOTg6ePHmCQ4cOYcuWLQrnHzt2LDZt2oRRo0ZhwoQJiI+Pxw8//IBr167hzJkzMDY25veNjY3FwIEDERYWhhEjRuDnn3/GyJEj0bJlSzRp0gSdOnXChAkTsGrVKsyaNQv+/v4AwP+vrgcPHmDXrl0YNGgQvL29kZKSgjVr1qBz586Ijo5GgwYN4O/vjy+++AJz587FmDFj8MYbbwAA2rdvr/ScjDH07dsXx44dQ1hYGJo1a4YDBw5g2rRpePz4MZYvXy63/+nTp7Fz506MGzcOVlZWWLVqFd59910kJibC3t6+Us+HEKLeNXf16tVo0qQJ+vbtCyMjI/zzzz8YN24cpFIpIiIi5M537949vP/++xg7dizCw8PRuHHjSp8jJiYGQ4YMwUcffYQRI0Zg48aNGDRoEPbv348ePXoAAFJSUtC+fXvk5uZiwoQJsLe3x+bNm9G3b1/89ddfCnmnIgsWLMD8+fPRvn17fPHFFzAxMcGFCxdw9OhRvPnmmyqPO3ToEB48eIBRo0bB2dkZt2/fxtq1a3H79m2cP38eHMdhwIABuH//PrZt24bly5fDwcEBAODo6IgtW7Zg9OjRaN26NcaMGQMAaNSoEf8c27Zty3+x4ujoiP/++w9hYWHIzMzEp59+KhfL4sWLIRAIMHXqVIjFYnzzzTf43//+hwsXLgAAPv/8c4jFYiQlJfHXVktLy0q9TgDw559/Ijc3Fx9//DHs7e1x8eJFfP/990hKSsKff/4JoCSHlZfjyrp9+zbeeOMNWFtb47PPPoOxsTHWrFmDkJAQnDhxAm3atJHb/5NPPoGdnR3mzZuHhIQErFixAuPHj8f27dsr/XwIAEaIEvPmzWMA2Icffii3vX///sze3p6/nZGRwUQiEZs+fbrcfhMmTGAWFhYsOzubMcZYfHw8A8DMzMxYUlISv9+FCxcYADZp0iR+W7du3VhgYCDLz8/nt0mlUta+fXvm6+vLb9u4cSMDwDp27MiKi4uVxt+3b1+57ePGjWMA2PXr1xljjEVFRTEAbPTo0XL7TZ06lQFgR48e5bd17tyZde7cmb8te04bN27kt3Xq1IlZWVmxhw8fyp1PKpWy8hQXF7OCggK5benp6czJyUnud/D8+XMGgM2bN0/hHBEREUzZn/SpU6cYAPbbb7/Jbd+/f7/Cdk9PTwaAnTx5kt/27NkzZmpqyqZMmcJv+/PPPxkAduzYMYXHUxWfp6cnGzFiBH87Pz+fSSQSuX3i4+OZqakp++KLL/htly5dUnidZUaMGME8PT3527t27WIA2FdffSW338CBAxnHcSw2NlYuThMTE7lt169fZwDY999/r/BYhBDV1L3mMsZYbm6uwvGhoaGsYcOGcttk16P9+/cr7F/Zc+zYsYPfJhaLmYuLC2vevDm/7dNPP2UA2KlTp/htWVlZzNvbm3l5efHXKmXXfdlzl4mJiWECgYD1799f4RpXUS5Q9ry2bdumcF3+9ttvGQAWHx+vsL+FhYXctVYmLCyMubi4sNTUVLnt7733HrOxseEf+9ixYwwA8/f3l8tLK1euZADYzZs3+W29e/eWuwbLyPJz2fhk5y6dO5Q950WLFjGO4+Ryqaocx5hi3unXrx8zMTFhcXFx/LYnT54wKysr1qlTJ4U4u3fvLve7mTRpEhMKhSwjI0Pp45HyUbcoUq6PPvpI7vYbb7yBFy9eIDMzEwBgY2ODd955B9u2beO7nUgkEmzfvh39+vWDhYWF3PH9+vWDq6srf7t169Zo06YN9u3bB6Ckq8zRo0cxePBgZGVlITU1FampqXjx4gVCQ0MRExODx48fy50zPDxc5biHst9gffLJJwDAP57s/8mTJ8vtN2XKFABQq0lc5vnz5zh58iQ+/PBDeHh4yN1X0TSFQqEQJiYmAEq6MKWlpaG4uBitWrXC1atX1Y5BmT///BM2Njbo0aMH/3qmpqaiZcuWsLS0lOt6BQABAQF8CwFQ8m1Y48aN8eDBg2rFUZapqSkEgpJLkEQiwYsXL2BpaYnGjRtX+Tnv27cPQqEQEyZMkNs+ZcoUMMbw33//yW3v3r07/60eUNIKZG1trfHnSkhdUdE1FwDMzMz4n8ViMVJTU9G5c2c8ePAAYrFY7nhvb2++e21plTlHgwYN5FoerK2tMXz4cFy7dg1Pnz7l42vdujU6duzI72dpaYkxY8YgISEB0dHRar8Gu3btglQqxdy5c/lrnExFuaD088rPz0dqairatm0LANXKBYwx7NixA3369AFjTC4XhIaGQiwWK5x/1KhRfF4CwOcFTV8fSz/nnJwcpKamon379mCM4dq1a5U+n0QiwcGDB9GvXz80bNiQ3+7i4oKhQ4fi9OnT/GcYmTFjxsj9bt544w1IJBI8fPiwCs+IUHFBylX2Q7KdnR0AyPXVHz58OBITE3Hq1CkAwOHDh5GSkqJ0FiVfX1+FbX5+fny/zNjYWDDGMGfOHDg6Osr9mzdvHgDg2bNncsd7e3urjL/s4zVq1AgCgYB/vIcPH0IgEMDHx0duP2dnZ9ja2lbqwiK74JbuHlYZmzdvRtOmTSESiWBvbw9HR0fs3btXIVFWVkxMDMRiMerXr6/wmmZnZyu8nmV/50DJ773s+IzqkkqlWL58OXx9fWFqagoHBwc4Ojrixo0bVX7ODx8+RIMGDWBlZSW3XdZlq+zvs6aeKyF1RUXXXAA4c+YMunfvDgsLC9ja2sLR0RGzZs0CAKXFhTKVOYePj4/Ch3o/Pz8AkMsFsi5Xpam6dpQnLi4OAoEAAQEBah8jk5aWhokTJ8LJyQlmZmZwdHTkX4Pq5ILnz58jIyMDa9euVcgDo0aNAqCYW9XJ/5qQmJiIkSNHol69erC0tISjoyM6d+4MoGrP+fnz58jNzVX5+5RKpXj06JHc9pp6rnUFjbkg5VLVIsBKDY4NDQ2Fk5MTfv31V3Tq1Am//vornJ2d0b1790o/nlQqBQBMnTpV6bdVABQKgdLfelRE1bdGul4A6ddff8XIkSPRr18/TJs2DfXr14dQKMSiRYsQFxdXrXNLpVLUr18fv/32m9L7ZYPEZdT5nVeFRCKRu71w4ULMmTMHH374Ib788kvUq1cPAoEAn376Kf8+0DZtPVdCSImy19a4uDh069YNr732GpYtWwZ3d3eYmJhg3759WL58ucLfvrLre2XPYUgGDx6Ms2fPYtq0aWjWrBksLS0hlUrx1ltvVet5yY794IMPMGLECKX7NG3aVO52da6PqnJq2TwgkUjQo0cPpKWlYfr06XjttddgYWGBx48fY+TIkZQLDBQVF6TahEIhhg4dik2bNmHJkiXYtWuXyq5KMTExCtvu37/PzzAha8I0NjauUnGi7PFKf/MVGxsLqVTKP56npyekUiliYmLkBiSnpKQgIyMDnp6eaj+WLHZls19V5K+//kLDhg2xc+dOuYuyrLVGprwiSNV9jRo1wuHDh9GhQ4dKFWLlKS8OOzs7hbnXCwsLkZycLLftr7/+QpcuXbBhwwa57RkZGfwAxYoeqyxPT08cPnwYWVlZcq0Xd+/e5e8nhGhPRdfcf/75BwUFBdizZ4/ct8Vlu2eWp7LnkLWIl76W3L9/HwDkcoGyNXOqcu1o1KgRpFIpoqOj0axZM7WPS09Px5EjR7BgwQLMnTuX364sb1Y2Fzg6OsLKygoSiUQjubWiOGTf/JfNBWVbgG7evIn79+9j8+bNGD58OL/90KFDaj9WWY6OjjA3N1f5+xQIBHB3d1frXKRqqFsU0Yhhw4YhPT0dY8eORXZ2Nj744AOl++3atUtuzMTFixdx4cIF9OzZEwBQv359hISEYM2aNQofRgFUemo42dR1Mt9//z0A8I/Xq1cvAFBYBGjZsmUAgN69e6v9WI6OjujUqRN+/vlnJCYmyt1X0bcfskKs9H4XLlzAuXPn5PYzNzcHoHjBBsCPbyl73+DBgyGRSPDll18qHFNcXKz0XBVR9VhASWI9efKk3La1a9cqfGMlFAoVXpc///xTYUxNeY9VVq9evSCRSPDDDz/IbV++fDk4juN/74QQ7ajomqvsWicWi7Fx40a1H6Oy53jy5InczH+ZmZn45Zdf0KxZMzg7OwMouXZcvHhR7pqbk5ODtWvXwsvLq1JdnPr16weBQIAvvvhC4Zv38nKBsucFKOYnoPzrooWFhcJ2oVCId999Fzt27FD6BVhVp121sLBQ2nVJNpatdC6QSCRYu3atQlyA/HNmjGHlypVKHwuoOBcIhUK8+eab2L17t1x3vJSUFGzduhUdO3aEtbV1+U+MVAu1XBCNaN68OV5//XX8+eef8Pf3R4sWLZTu5+Pjg44dO+Ljjz9GQUEBVqxYAXt7e3z22Wf8PpGRkejYsSMCAwMRHh6Ohg0bIiUlBefOnUNSUpLcOggViY+PR9++ffHWW2/h3Llz+PXXXzF06FB+bY6goCCMGDECa9euRUZGBjp37oyLFy9i8+bN6NevH7p06VKp12HVqlXo2LEjWrRogTFjxsDb2xsJCQnYu3ev0hWmZd5++23s3LkT/fv3R+/evREfH4+ffvoJAQEByM7O5vczMzNDQEAAtm/fDj8/P9SrVw+vv/46Xn/9dbRs2RIAMGHCBISGhkIoFOK9995D586dMXbsWCxatAhRUVF48803YWxsjJiYGPz5559YuXIlBg4cWKnn2axZMwiFQixZsgRisRimpqb8nPOjR4/GRx99hHfffRc9evTA9evXceDAAbnWCNlz/uKLLzBq1Ci0b98eN2/exG+//SY3AA8oSVK2trb46aefYGVlBQsLC7Rp00ZpX+w+ffqgS5cu+Pzzz5GQkICgoCAcPHgQu3fvxqeffio3eJsQonkVXXPffPNNmJiYoE+fPvyXUevWrUP9+vWVfqGkTGXP4efnh7CwMFy6dAlOTk74+eefkZKSIleMzJgxA9u2bUPPnj0xYcIE1KtXD5s3b0Z8fDx27NihMDC7PD4+Pvj888/x5Zdf4o033sCAAQNgamqKS5cuoUGDBli0aJHS46ytrdGpUyd88803KCoqgqurKw4ePIj4+HiFfWXX+88//xzvvfcejI2N0adPH1hYWKBly5Y4fPgwli1bhgYNGsDb2xtt2rTB4sWLcezYMbRp0wbh4eEICAhAWloarl69isOHDyMtLU3t51g6ju3bt2Py5MkIDg6GpaUl+vTpgyZNmqBt27aYOXMm0tLSUK9ePfz+++8oLi6WO/61115Do0aNMHXqVDx+/BjW1tbYsWOH0rEOqnKcMl999RUOHTqEjh07Yty4cTAyMsKaNWtQUFCAb775ptLPk1RSjc5NRQyGbGq958+fy21XNb0cY4x98803DABbuHChwn2y6fu+/fZb9t133zF3d3dmamrK3njjDbkpCmXi4uLY8OHDmbOzMzM2Nmaurq7s7bffZn/99ZdCLJcuXVIZf3R0NBs4cCCzsrJidnZ2bPz48SwvL09u36KiIrZgwQLm7e3NjI2Nmbu7O5s5c6bcVLiMqTcVLWOM3bp1i/Xv35/Z2toykUjEGjduzObMmaMQY2lSqZQtXLiQeXp6MlNTU9a8eXP277//Kky1yhhjZ8+eZS1btmQmJiZy0+8VFxezTz75hDk6OjKO4xSm7Fu7di1r2bIlMzMzY1ZWViwwMJB99tln7MmTJ/w+np6erHfv3grxlX3ujDG2bt061rBhQyYUCuWmFpRIJGz69OnMwcGBmZubs9DQUBYbG6t0KtopU6YwFxcXZmZmxjp06MDOnTun9LF2797NAgICmJGRkdxrruz1ycrKYpMmTWINGjRgxsbGzNfXl3377bcKU0ACYBEREQrPtWychJCKVeaau2fPHta0aVMmEomYl5cXW7JkCfv5558Vcouq61FVznHgwAHWtGlTZmpqyl577TX2559/KpwzLi6ODRw4kL92t27dmv37779y+6gzFa3Mzz//zJo3b85MTU2ZnZ0d69y5Mzt06FC5r2NSUhKfP2xsbNigQYPYkydPlE7x/eWXXzJXV1cmEAjknvfdu3dZp06dmJmZGQMgdz1LSUlhERERzN3dnRkbGzNnZ2fWrVs3tnbtWn4f2XSxZV8jZc89OzubDR06lNna2jIActfjuLg41r17d2ZqasqcnJzYrFmz2KFDhxSmoo2Ojmbdu3dnlpaWzMHBgYWHh/PTgpd+rPJynLLX5+rVqyw0NJRZWloyc3Nz1qVLF3b27Fm5fVR9jlA2ZS5RH8cYjVYhmrFy5UpMmjQJCQkJCjMvJCQkwNvbG99++y2mTp2qowgJIYRow/z587FgwQI8f/5coZWSEFK30JgLohGMMWzYsAGdO3dWOr0nIYQQQgip/Wp9cZGRkYFWrVqhWbNmeP3117Fu3Tpdh1Sr5OTkYNu2bRg7dixu3ryJSZMm6TokQghRC+UHQgjRvFo/oNvKygonT56Eubk5cnJy8Prrr2PAgAGwt7fXdWi1wvPnzzF06FDY2tpi1qxZ6Nu3r65DIoQQtVB+IIQQzatTYy7S0tLQokULXL58mfqEEkII4VF+IIQQzdD7blEnT55Enz590KBBA3Ach127dinsExkZCS8vL4hEIrRp0wYXL16Uuz8jIwNBQUFwc3PDtGnTKHEQQkgtQPmBEEL0j94XFzk5OQgKClJYmEdGNr/yvHnzcPXqVQQFBSE0NBTPnj3j97G1tcX169cRHx+PrVu3IiUlpabCJ4QQoiWUHwghRP8YVLcojuPw999/o1+/fvy2Nm3aIDg4mF+RVyqVwt3dHZ988glmzJihcI5x48aha9euKhcNKygoQEFBAX9bKpUiLS0N9vb2ai89TwghhoQxhqysLDRo0KBSi4XpE8oPhBCiHZXOETpbYaMKALC///6bv11QUMCEQqHcNsYYGz58OOvbty9jjLGnT5+yzMxMxhhjGRkZrEmTJuzGjRsqH0O2GA79o3/0j/7VtX+PHj3S+HW7pgCUH+gf/aN/9E+b/9TNEQY9W1RqaiokEgmcnJzktjs5OeHu3bsAgIcPH2LMmDFgjIExhk8++QSBgYEqzzlz5kxMnjyZvy0Wi+Hh4YH4+HhYWVlVGNPjx4/RrFkzsFINQgKBANeuXYOrq2tln6LOFBUV4dixY+jSpQuMjY0Vd8jJgbGnZ8m+Dx8CFhY1HKFyOYU58FxVElfcx3GwsbDRcUSVV+Frr+cMOX5Djh2oevxZWVnw9vZW6xpnKCg/aI9a7zPKEVpjyNcpQ44dqLvxVzZHGHRxoY7WrVsjKipK7f1NTU1hamqqsL1evXqwtrau8Hh7e3usWLECEydOBAAIhUKsWbMGTZs2VTsGfVBUVARzc3PY29srfwOKRK9+trfXm8QhKhQBL0OrZ18Ptha2Oo2nKip87fWcIcdvyLEDVY9ftm9d69qji/ywbt06hIeHgzEGjuOwdu3a2pcfAMoRWmTI1ylDjh2ou/FXNkcYZufalxwcHCAUChUG4KWkpMDZ2bla546MjERAQACCg4MrfWx4eDj/87Vr1xAWFlatWAghhFSOvuaHsLAwjB8/HgDwv//9j/IDIaTWMejiwsTEBC1btsSRI0f4bVKpFEeOHEG7du2qde6IiAhER0fj0qVLlT629DdbtBgTIYTUPH3NDwD44sZQB88TQkh59L5bVHZ2NmJjY/nb8fHxiIqKQr169eDh4YHJkydjxIgRaNWqFVq3bo0VK1YgJycHo0aN0lnMAoEAxsbGKCoqkptZhBBCiOYYYn4AAIuXXYRycnJ0GgchhGiD3hcXly9fRpcuXfjbssF0I0aMwKZNmzBkyBA8f/4cc+fOxdOnT9GsWTPs379fYRBfTTM1NUVRURHy8/N1GgdRn0QiQVFRka7DQFFREYyMjJCfnw+JRKLrcCrNkOM35NiB8uM3NjaGUCjUUWTaYaj5wdzcHACQl5en0zhI5VCOqD5Djh2ovfFrOj/ofXEREhIiN7OGMuPHj+f7sGpKZGQkIiMjq/zmEYlEyM7OppYLA8AYw9OnT5GRkaHrUACUxOPs7IxHjx4Z5ABbQ47fkGMHKo7f1tYWzs7OBvnclDHU/GBmZgYAyM3N1WRYREsoR2iOIccO1O74NZkf9L640JWIiAhEREQgMzMTNjaVn6pONu6Cigv9J0sa9evXh7m5uc4vGFKpFNnZ2bC0tDTIPtmGHL8hxw6ojp8xhtzcXH5lahcXF12FWCtUNz/IWi6ouDAMlCM0x5BjB2pn/NrID1RcaIno5TR8VFzoN4lEwicNfRl8L5VKUVhYCJFIZLAXL0ON35BjB8qPX/Zt+bNnz1C/fv1a10XKkFBxYTgoR2iWIccO1N74NZ0fDO+VqSHVmWoQeNVyQWMu9Jus/6ws2RNSm8ne5/rQb9yQVTc/yAZ0U3Gh/yhHkLpCk/mBigsVqjvVIHWLMiy6buYmpCbQ+1wzqpsfaMyF4aG/HVLbafI9TsWFlsi6RVHLBSGEkNJyJCVdDnJyXhUXyeI8nI1LRbI4T+W2mtyHEEKqisZcaAl1iyKk6jZt2oRPP/2Un51l/vz52LVrF6KionQaFyHVtf1SIqZtuQEAyMrJwfZLiQCAmTtvQsoAAQcsGhCosK1/c1f8fe1xjeyzaEAgOvk5Ij41B94OFnCxMavpl4kQlSg/6D8qLrRE1nJRWFio40gIMXxTp07FJ598ouswAJQ0Hf/999/o16+frkMhBiZZnIeZO2+CGZXkB0gl+OyPa+CEr1KxlAHTd9yUO07KgB1XH9fYPtN33AQHgOFVsTEk2APJ4jzEp+bAzca0Kk+fEK2g/KB/qLhQobrzmJuYmACglgtSe0gkEnAcp5MZMiwtLWFpaVnjj6tNRUVFMDY21nUYpAqqmh/iU3NKWgeMX304Z0X54IT6996WrR4iKzYORqfg6N1nYC9bNwZ7c2guzkeSWEytG4Tyg4YZen6gMRcqVHfAHk1FS7Ttr7/+QmBgIMzMzGBvb4/u3bsjJycHQMl0c1988QXc3NxgamrKr0wsc/z4cXAcJ7coVFRUFDiOQ0JCAoCSpmdbW1vs2bMHAQEBMDU1RWJiIgoKCjB9+nS4u7vD1NQUPj4+2LBhA3+eW7duYeDAgbC2toaTkxOGDRuG1NTUcp/Lpk2b4OHhAXNzc/Tv3x8vXryQu3/+/Plo1qwZf/vSpUvo0aMHHBwcYGNjg86dO+Pq1atyx3AchzVr1uDtt9+Gubk5/P39ce7cOcTGxiIkJAQWFhZo37494uLi5I7bt28fWrVqBZFIhIYNG2LBggUoLi4GAHh5eQEA+vfvD47j+NsAsHv3brRo0ULpcbJ4Vq9ejb59+8LCwgJff/11ua8J0V9VzQ/eDhYQcAAERgD3co75ogKUHUfJAQrbylJ7nwpiUmcfADhyp6SwAEoKjt8fCND5u5MYuu4COiw+ynfv4sduZNIXa7pE+UE7+WH37t3o3LkzzM3NKT+Ug4oLLaExFwaMMSAnRzf/KlhtWCY5ORnvv/8+PvzwQ9y5cwfHjx/HgAED+NWKV65cie+++w5Lly7FjRs3EBoair59+yImJqZSL0Vubi6WLFmC9evX4/bt26hfvz6GDx+Obdu2YdWqVbhz5w7WrFnDf2uUkZGB7t27o2nTprh48SL279+PlJQUDB48WOVjXLhwAWFhYRg/fjyioqLQpUsXfPXVV+XGlZWVhREjRuD06dM4f/48fH190atXL2RlZcnt9+WXX2L48OGIiorCa6+9hqFDh2Ls2LGYOXMmLl++DMaY3OrNp06dwkcffYRPPvkE0dHRWLNmDTZt2sRf6GUfJjdu3Ijk5GT+9qlTpzB8+HBMnDhR6XEy8+fPR//+/XHz5k18+OGHav4WSG3hYmOGRQMCIRRw4F62Xkzt5oXFAwIhfFkpCDkOi98NVNj2bgvXqu3zbuX2EUC9YgPg5IqNmTtvYuXh++iw+GhJwbHqPLY37VHNV0wPUX7g1cX8MHLkSHz00Ue4desW5YfyMFIusVjMADCxWFyp48LDwxkANnfuXC1Fpl2FhYVs165drLCwUPkO2dmMlVzqSn7WE9kF2QzzwTAfLD07vcL98/LyWHR0NMvLyyt1klLPrab/ZWcziUTC0tPTmUQiURn3lStXGACWkJCg9P4GDRqwr7/+Wm5bcHAwGzduHGOMsWPHjjEALD391Wt07do1BoDFx8czxhjbuHEjA8CioqL4fe7du8cAsEOHDil93C+//JL16NFDLv5Hjx4xAOzevXtKj3n//fdZr1695LYNGTKE2djY8LfnzZvHgoKClB7PGGMSiYRZWVmxf/75h98GgM2ePZu/fe7cOQaAbdiwgd+2bds2JhKJ+NvdunVjc+bMkXvtt2zZwlxcXOTO+/fff8s9frdu3djChQvltik77tNPP1X5HDShoveO0vc7q/p1rq6r6uv28Hkms7SxZQDY9evXGWOMPcnIZWdjU9mTjFx+v7Lbamqf3y8+ZA1n7GWe0/9lDWfsZT8dj2XeM/5lntMr96/htN0sytmHnbn1SO6xdKnaOULH+YGxiv/OKT/I02R++Prrr+XiN6T8wFj57x1V+YGxyl/raMyFllDLBdGmoKAgdOvWDYGBgQgNDcWbb76JgQMHws7ODpmZmXjy5Ak6dOggd0yHDh1w/fr1Sj2OiYkJmjZtyt+OioqCUChE586dle5//fp1HD9+HG5ubgr3xcXFwc/PT2H7nTt30L9/f7lt7dq1k2umLyslJQWzZ8/G8ePH8ezZM0gkEuTm5iIxMVFuv9KxOzk5AQACAwPltuXn5yMzMxPW1ta4fv06zpw5g2XLlvH7SCQS5OfnIzc3V+VCWrLjSn8Tpey4Vq1aqXxOpO5wsRHBQmSKbPGrtS5cbMwUxi2U3VZT+wwJ9kAnP0ckpObCy8EcLjZmsDU3xqydtyBh7OVgb4aK2jgkAiHeGb4c2HJd6cBwGquhHZQftJsfFi5cyO9D+UE5Ki5UqO6AbllxQbNFGSBzcyA7W3ePrUbTt1AoxKFDh3D27FkcPHgQ33//PT7//HNcuHAB9vb2FR4vG3THSj2WslU5zczM5BbWkS3+pUp2djbefvttzJ49G5aWlnKD+1xcXCqMS10jRozAixcvsHLlSnh6esLU1BTt2rVT+HsrPSBO9jyUbZNKpXz8M2bMwPvvv68wMFE2jkqZ7OxsLFiwAAMGDFC4r/RxspWZiWGrbn4AXuUIWT94fVO2ACldcLjamODHncfwR7wQUlbSveqztxpjyf67kJa9fMn+xl4ODD9y5xkO30mRm/Z2SLBHDT6zatJ1flAD5Qft5Yf58+ejR48eCvFTfpBHxYUKERERiIiIQGZmJmxsbCp9PC2iZ8A4DtDlH7ma/Wo5jkOHDh3QoUMHzJ07F56envj7778xefJkNGjQAGfOnJH7BunMmTNo3bo1AMDR0RFASd9cOzs7AFBrjvDAwEBIpVKcOHEC3bt3V7i/RYsW2LFjBzw8PFCvXj21Zg7x9/fHhQsX5LadP3++3GPOnDmDH3/8Eb169QIAPHr0qMJBgepo0aIFYmNj4ePjozJ2Y2NjhQ+VLVq0wL179+Dj41PtGIj+q25+AF4VF4a0Sres4CgqKkI7J4ZxAzrhsbiwnNYNRQejU/ifZWM1OvmVXI8MojVD1/lBTZQftJMf7t27h7Fjx8La2lpp/JQfSlBxoSWyxEGzRRFtuHDhAo4cOYI333wT9evXx4ULF/D8+XP4+/sDAKZNm4Z58+ahUaNGaNasGTZu3IioqCj89ttvAAAfHx+4u7tj/vz5+Prrr3H//n189913FT6ul5cXRowYgQ8//BCrVq1CUFAQHj58iGfPnmHw4MGIiIjAunXrMHr0aMycORMODg6IjY3F77//jvXr10MoFCqcc8KECejQoQOWLl2Kd955BwcOHCi3yRsAfH19sWXLFrRq1QqZmZmYNm1ahd+aqWP27Nno27cvGjVqhEGDBkEgEOD69eu4desWP4jQy8sLR44cQYcOHWBqago7OzvMnTsXb7/9Njw8PDBw4EClxxFSmiEWF2W52Ijg4WDF3y7dumEuKUD/9VcgreADpJQBb686jbScQrl1NWgRv6qj/KCd/CC7zjs5OWHo0KEwMjKi/KACzRalJdRyQbTJ2toaJ0+eRK9eveDn54fZs2fju+++Q8+ePQGUXJAnT56MKVOmIDAwEPv378eePXvg6+sLoOTblW3btuHu3bto2rQplixZovZFbvXq1Rg4cCDGjRuH1157DeHh4XzXjgYNGuDUqVOQSCR46623EBgYiE8//RS2trYqv6Vq27Yt1q1bh5UrVyIoKAgHDx7E7Nmzy41hw4YNSE9PR4sWLTBs2DBMmDAB9evXV/flUyk0NBS///47Dh06hODgYLRt2xbLly+Hp6cnv893332HQ4cOwd3dHc2bN+eP+/fff3Hw4EGVxxFSWm0oLpRxsTFDu0b2CHK1waID30MoLfkWV8hxmNnztZKpeMt48bKwAF51n2ovm3Wq1DS3RD2UH7SXH/bs2YOjR4+iTZs2lB/KwTGmZh+MOkrW7C0Wi2Ftba32catWrcLEiRPRv39/7Ny5U4sRakdRURH27duHXr16KV/IJScHkC1ak52tN83EOYU5sFxUElf61HTYWtiWu39+fj7i4+Ph7e1dbp/JmiSVSvkBZLpYkKi6DDl+Q44dqDh+Ve/3ql7n6rqqvm5FRUUICQnB2bNn8cMPPyAiIkKLUWpehfkB4HNEspU9Es5FwcvNHi42Zth+KZHvOiXkOAxt44Et5x+W+3gCDjgzoyuA6nedohyhW4YcO1C74y/vvV7Zax11i9ISWkSPEEKIKvo+oFtTXLJewMXLDrBQPhMVAPx24aHiQPBSpAwY9NM5PE7Pk+s6ZVADwQmpQwyv7KohkZGRCAgIQHBwcJWONzExAUDFBSGE1DbVzQ9A7e0WpQ5Z1ynZAPFFAypexC/pZWEBvBoIHpWYVrIauDivpkInhKiBWi5U0NRsUVRcEEJI7VJXZ4vSlrKtGSfvP5frOjWwpRu2X34kd4yUAf1+PAcAtIYGIXqGigstoUX0CCGEqELFhbzyFvEDgD+vPFLZdUrKgBk7buJGkhjbLiYa7hoahNQS1C1KS6jlghBCiCpUXJSvoq5TZTEAv11I5AsQKQNm7byF64/SqesUITWMWi60hKaiJYQQokpdGdCtKXJraJgI0P/Hs+UOAgcACWPoF3lWbhD420EVr1BNCKkearnQElniKLvcPCGEEEItF5XHr6HhbifXklHeGhqlB4HP2nkLT8X5EDJ7mEoCkZJJPQsI0QZqudASGnNBCCFEFSouqqfsuAwXGzPYmhvzA8E5vCosZCSM4Zv9cXDN3wgOAvRcdZnGZRCiBVRcaIkscdCYC0IIIWVRcVF9pQeBA+p1nTpw+zm4l502ZK0ZnfwcaWYpQjSIukWpUN15zGnMBSGad/z4cXAch4yMDADApk2bYGtrq9OYSN1D61zoL1VdpwQc4O9ipbC/hDFcSaBB37UB5Qf9QcWFChEREYiOjsalS5eqdHzpblGMVTDqjBA95eXlhRUrVmj0nAkJCeA4DlFRUdU+15AhQ3D//v3qB6UB2nitiH6qbn4AXn0BRQO6tWdIsAdOz+iCbeFtcWZGV/w8MljpuIzx265h6LoL6LD4KLZfSqz5QA0U5Qf11bX8QMWFlsgSBwAUFRXpMBJCai8zMzPUr19f12FoFE0CUTeYmJgAoJYLbSs7pe38vo3BIFG6r2zl72RxHpLFeTgbl4pnWdT7wFBRftAdKi60RNZyAVDXKKIdUqkU33zzDXx8fGBqagoPDw98/fXX/P03b95E165dYWZmBnt7e4wZMwbZ2dn8/SNHjkS/fv2wdOlSuLi4wN7eHhEREXwxHBISgocPH2LSpEngOA7cy+4FL168wPvvvw9XV1eYm5sjMDAQ27ZtU4jt22+/VRqbt7c3AKB58+bgOA4hISEqn+O+ffvg5+cHMzMzdOnSBQkJCXL3l232jouLwzvvvAMnJydYWloiODgYhw8fljvGy8sLX331FYYPHw5LS0t4enpiz549eP78Ofr16wc3Nzc0a9YMly9fljvu9OnTeOONN2BmZgZ3d3dMmDCB/9ZZ1WtV0XGyeL788ksMHz4c1tbWGDNmjMrXg9Qe1C1KN95t6YLHog/x1GQmlgzwU7hfyoCRP19Eh8VHMXTdBQxddx45BcU6iLR6KD9oPj+88847sLa2RocOHSg/VICKCy0pXVzQoG7DwhhDTmGOTv5VpgvdzJkzsXjxYsyZMwfR0dHYunUrnJycAJR0tQgNDYWdnR0uXbqEP//8E4cPH8b48ePlznHs2DHExcXh2LFj2Lx5MzZt2oRNmzYBAHbu3Ak3Nzd88cUXSE5ORnJyMoCSYrlly5bYu3cvbt26hTFjxmDYsGG4ePEif94FCxZgyZIlSmOT7Xf48GEkJydj586dSp/fo0ePMGDAAPTp0wdRUVEYPXo0ZsyYUe5rkp2djV69euHIkSO4du0a3nrrLfTp0weJifJdHZYvX44OHTrg2rVr6N27N4YNG4bhw4fjf//7H44fP46GDRti+PDh/O8jLi4Ob731Ft59913cuHED27dvx+nTp/nXU9VrVdFxMkuXLkVQUBCuXbuGOXPmlPscSe1AxYXuSLgXKBDeRDN3a6XdpO6lZPMDwRkDMnKLUFgsfXmb8gNQN/PDBx98gMuXL8Pb2xsjR46k/FAeRsolFosZACYWiyt1XGFhITMyMmIAWGJiopai057CwkK2a9cuVlhYqHyH7GzGSq67JT/rieyCbIb5YJgPlp6dXuH+eXl5LDo6muXl5Sk9R03/yy7IZhKJhKWnpzOJRKIy7szMTGZqasrWrVun9P61a9cyOzs7ll3qd7N3714mEAjY06dPGWOMjRgxgnl6erLi4mJ+n0GDBrEhQ4bwtz09Pdny5csrfB179+7NpkyZwhhjLCMjg5mamrI1a9Yo3Tc+Pp4BYNeuXSv3nDNnzmQBAQFy26ZPn84AsPT0dMYYYxs3bmQ2NjblnqdJkybs+++/5297enqyDz74gL+dnJzMALA5c+bwr/2ZM2cYAJacnMwYYywsLIyNGTNG7rynTp1iAoGAf+8oe63UPa5fv37lPgd1VfTeUfZ+Z6zq17m6rjr5YcuWLQwls6Wqvs7qqQrzA2MGkyN+v/iQNZyxl3lO/5c1nLGXfbL1KvOc/i//r/1X/7GDZ66wlDQxy8orZGk5mTrND4xV/HdO+SGdMab5/MBYyWt/8OBBg8wPsvhVvXdU5QfGKn+to6lotcjY2BjFxcXUckE07s6dOygoKEC3bt1U3h8UFAQLCwt+W4cOHSCVSnHv3j3+W6ImTZpAKBTy+7i4uODmzZvlPrZEIsHChQvxxx9/4PHjxygsLERBQQHMzc3Viq0yz7FNmzZy29q1a1fuMdnZ2Zg/fz727t2L5ORkFBcXIy8vT+GbqaZNm/I/y16LwMBAhW3Pnj2Ds7Mzrl+/jhs3buC3337j92GMQSqVIj4+Hv7+/krjUfe4Vq1alfu8SO1Telxebm4ubGxsdBhN3VV2vQwA+PfGE4UpbJPF+eBypMgr0v8B+JQflNNUfpCN46D8oBoVF1pkYmKCvLw8GnNhYMyNzZE9M7viHbX02EyNpm8zM83MyW5sbCx3m+M4SKXSco/59ttvsXLlSqxYsQKBgYGwsLDAp59+yg8001RsVTF16lQcOnQIS5cuhY+PD8zMzDBw4ECFQXCln7es/6uybbLXIjs7G2PHjsWECRMUHtPDQ/UCXOoeVzrJk7rByMgIAoEAUqmUigsdK7texqIBgfxifGWJjMxxfvhj+DlZwcSoZnuWmxubq7Uf5QflKD/UHCoutMjIqOTlpZYLw8JxHCxMdPfHrE5x4evrCzMzMxw5cgSjR49WuN/f3x+bNm1CTk4Of2E6c+YMBAIBGjdurHYsJiYmkEjkZ1Y5c+YM3nnnHXzwwQcASi6w9+/fR0BAgEJsjRo1UnpOAArnVfYc9uzZI7ft/Pnz5R5z5swZjBw5Ev379wdQcvEuO8ivKlq0aIHo6Gj4+Pio3EfZa6XOcaRu4jgO5ubmyM7OpnEXeqZ0a0ZGdg6Qk8rfx3EczIzNIYApjAVGKCyWwMRIWOOFRnkoPyhH+aHm6M9fQy0k+yOhlguiaSKRCNOnT8dnn32GX375BXFxcTh//jw2bNgAAPjf//4HkUiEESNG4NatWzh27Bg++eQTDBs2jG/mVYeXlxdOnjyJx48fIzW1JMH6+vri0KFDOHv2LO7cuYOxY8ciJSVFLraJEydixowZSmOrX78+zMzMsH//fqSkpEAsFit97I8++ggxMTGYNm0a7t27h61bt/KDCVXx9fXFzp07ERUVhevXr2Po0KEVftOmjunTp+Ps2bMYP348oqKiEBMTg927d8sNvFP2WqlzHKm7ZF1FqLjQP7IpbJs0sIaSMd9ITM/D3aeZeJCag3tPM5GWoz9fIlJ+UI7yQ82h4kKLZM1o1HJBtGHOnDmYMmUK5s6dC39/fwwZMgTPnj0DUPKh5cCBA0hLS0NwcDAGDhyIbt264YcffqjUY3zxxRdISEhAo0aN4OjoCACYPXs2WrRogdDQUISEhMDZ2Rn9+vWTO27atGmYPHmy0tiMjIywatUqrFmzBg0aNMA777yj9LE9PDywY8cO7Nq1C0FBQfjpp5+wcOHCcuNdtmwZ7Ozs0L59e/Tp0wehoaFo0aJFpZ6zMk2bNsWJEydw//59vPHGG2jevDnmzp2LBg0a8Psoe63UOY7UXVRc6L/6ViLYmhuDe1licACMhQK5FmYG4HF6Pj+jlD6g/KCI8kPN4Zg6fTDqsMzMTNjY2EAsFsPa2lrt44qKiuDv74+4uDjs3bsXvXr10mKUmldUVIR9+/ahV69eCv0uAQA5OYClZcnP2dmAnvQJzCnMgeWikrjSp6bD1sK23P3z8/MRHx8Pb29vuQGWuiSVSpGZmQlra2sIBIZX/xty/IYcO1Bx/Kre71W9ztV11ckP+/btw6xZsxAdHY0jR46ga9euWoxUsyrMD0CtyxGu7p4QGJnAxEiAwmIJHqQqDuxu6GAJS5H2e5sb8nXKkGMHanf85X0equy1zvBemRoSGRmJgIAABAcHV/kcsosudYsihJDaQxP5AaCWC0NiYiSApcgIJkYCmBgJlXaVKpJKUVgsRXZ+kV61YhBS06i4UCEiIgLR0dG4dOlSlc9B3aIIIaT20UR+AKi4MFQmRgK42pnxXaVkHqXl6u04DEJqEs0WpUXUckEIIUQVKi4MVz0LU1iaGqOwWApjIw5PxfkQ5xXx98vGYViaGuvVTFKE1AR6x2uRbLYoarkghBBSlmzOfyouDJOsq5SpkRD2FiYK9zMw5BYWUzcpUudQy4UWUcsFIYQQVWRrDOTk6P+qz6R8snEYZWfISUwrKRw5AK52ZqhnYVrToRFS46jlQotozAUhhBBVqFtU7aFqHIaMPk5XS4i2UMuFFtEieoQQQlSh4qJ2KT0Oo1gq5VstZBgYCoulNAaD1HpUXGiRkVHJy0stF4QQQsqiMRe1T8lUtQIUFkuVdpOiwoLUBfQu1yJquSCEEKIKtVzUXqq6SeUXSXQUESE1h4oLLaIxF4S8Mn/+fDRr1qxa50hISADHcYiKitJITIToEg3ort3qWZiisbMVGjpYwM685MvGxLRcZOUX0QxSZVB+qF2ouNAimi2K6LPjx4+D4zhkZGToOhRSCsdx2LVrl67DIDWAWi5qv5Lpao3hamcGS1MjSBlDfGqO3i+0R/lBPxlKfqDiQououCCEaENhYaGuQyAaQGMu6g4Bx6GBrZncNppBimiDPuQHKi60iLpF1T3J4jycjUtFsjhP648llUqxaNEieHt7w8zMDEFBQfjrr78AAIwxdO/eHaGhoWCsZEhhWloa3NzcMHfuXCQkJKBLly4AADs7O3Ach5EjR1Z4XuDVN1pHjhxBq1atYG5ujvbt2+PevXty8S1ZsgROTk6wsrJCWFiY0iJ7/fr18Pf3h0gkwmuvvYYff/xR7v6LFy+iefPmEIlEaNWqFa5du1bh65KcnIzevXvDzMwM3t7e2Lp1K7y8vLBixQp+n4yMDIwePRqOjo6wtrZG165dcf36dbnzrF69Go0aNYKJiQkaN26MLVu2yN3PcRzWrFmDt99+G+bm5vD398e5c+cQGxuLkJAQWFhYoH379oiLi5M7bvfu3WjRogVEIhEaNmyIBQsWoLi4GADg5eUFAOjfvz84juNvy7oMrF+/Ht7e3hCJRPjll19gb2+vcH353//+h+HDh1f4OhHdo5aLuiUpLRc3kjKQmv3qb1Y2g5SmUX5Qrq7nh/79+2Ps2LEVvk7Vxmq5xMRE1rlzZ+bv788CAwPZH3/8UanjxWIxA8DEYnGljissLGQREREMAOvduzd7kpHLzsQ+Z08ycss9Tp39NLVPRfHv2rWLFRYWKt8hO5sxoORfdnaVHkMbsguyGeaDYT5YenZ6hfvn5eWx6OholpeXx2+TSqUsp6Co0v9+ORvPvGf8yzyn/8u8Z/zLfjkbX+lzSKVSJpFIWHp6OpNIJOXG/tVXX7HXXnuN7d+/n8XFxbGNGzcyU1NTdvz4ccYYY0lJSczOzo6tWLGCMcbYoEGDWOvWrVlRURErLi5mO3bsYADYvXv3WHJyMsvIyFDrvMeOHWMAWJs2bdjx48fZ7du32RtvvMHat2/PGGNMIpGwn3/+mZmamrL169ezu3fvss8//5xZWVmxoKAgPv5ff/2Vubi4sB07drAHDx6wHTt2sHr16rFNmzYxxhjLyspijo6ObOjQoezWrVvsn3/+YQ0bNmQA2LVr11S+Lt27d2fNmjVj58+fZ1euXGGdO3dmZmZmbPny5XL79OnTh126dIndv3+fTZkyhdnb27Pnz5+z9PR09tdffzFjY2MWGRnJ7t27x7777jsmFArZ0aNH+XMAYK6urmz79u3s3r17rF+/fszLy4t17dqV7d+/n0VHR7O2bduyt956iz/m5MmTzNramm3atInFxcWxgwcPMi8vLzZ//nzGGGPPnj1jANjGjRtZcnIye/bsGWOMsXnz5jELCwv21ltvsatXr7Lr16+z3NxcZmNjI3dNS05OZkZGRuzw4cNKXxtl73fGqn6dM3S6zA+7du1ie/fuZQDk/i4MQYX5gbFamyM0kR+8ZvzLFu2LZhcepLILD1LZU3EuSxHnsvScggrzA2NMrRxB+UG56uSHFy9eMIlEwrZs2WKQ+SElJYUZGRmxPXv2KH3vqMoPjFX+Wlfri4snT57wb7Tk5GTWoEEDll2JC111ksekSZMYABbYuqPch87fLz5UeszvFx9WuJ+m9mGs/AKkLhcXOQVFzHP6vzr5l1NQpFbiyM/PZ+bm5uzs2bNy28PCwtj777/P3/7jjz+YSCRiM2bMYBYWFuz+/fv8fbIkkJ7+6nVS57yy40p/gJV9SMrLy2MSiYQFBwezjz/+WO4cbdq0kUsejRo1Ylu3bpXb58svv2Tt2rVjjDG2Zs0aZm9vL/e7Wb16dbnJ486dOwwAu3TpEr8tJiaGAeCTx6lTp5i1tTXLz8+XO7ZRo0Zs9erVLD09nbVv356Fh4fL3T9o0CDWq1cv/jYANnv2bP72uXPnGAC2YcMGftu2bduYSCTib3fr1o0tXLhQ7rxbtmxhLi4ucuf9+++/5faZN28eMzY25pOJzMcff8x69uzJ3166dCnz8vJixcXFSl8fKi7k6TI/7Nq1ix0/fpwBYD4+PpU6XtfqcnGh6/zAWMXFBeWHa0pfl+rmhzVr1jCJRMJat27NRo8eLXe/IeSH7777jjVs2JClpaVpvbio9etcuLi4wMXFBQDg7OwMBwcHpKWl8bN0aEuyOB+pRSWzQ9x/kgbnl5NdSxkwc+dNdPJzBADEp+bA28ECj9PzMH3HTf54KQNm7byF15ytkFMogbeDBbLzizFjx01+3mxl+zDGMGPnTTAmv08nP0e42Lzq77n9UiJm7rwJKQMEHLBoQCCGBHsofR5JYjG8HSzkjie6FRsbi9zcXPTo0UNue2FhIZo3b87fHjRoEP7++28sXrwYq1evhq+vr0bOCwBNmzblf5b9jT179gxubm64f/8+xo0bJ7d/u3btcOzYMQAls+PExcUhLCwM4eHh/D7FxcWwsbEBANy5cwdNmzaFSCSSO0d57t27ByMjI7Ro0YLf5uPjAzs7O/729evXkZ2dDXt7e7lj8/Ly8ODBA/6xx4wZI3d/hw4dsHLlSpWvgZOTEwAgMDBQblt+fj4yMzNhbW2N69ev48yZM/j666/5fSQSCfLz85Gbm8t3k1HG09MTjo6OctvCw8MRHByMx48fw9XVFZs3b8bQoUPBccpXCSbydJUfZGjMBdEGyg/KVTc/yLow3b9/Hx999JHc/YaQHzZt2oQRI0bUSH7Q++Li5MmT+Pbbb3HlyhUkJyfj77//Rr9+/eT2iYyMxLfffounT58iKCgI33//PVq3bq1writXrkAikcDd3V2rMcs+uGcnmwIAmKRI7n4pA/637gLiX+TwRYAyEsbQL/IsGAAOgJGAU1iQp+w+DpYmCueUMIaL8WlwtDKFu50ZzsWlKS1kyhYg51I4TPruZIUFSG1jZixE9BehlTrmqTgf3ZedgLTUay/ggMOTO8PZRqT6QCWPzcp7U7yUnZ0NANi7dy9cXV3l7jM1NeV/zs3NxZUrVyAUChETE6Ox8wKvxhQB4C9WUql6fYdlj7Nu3Tq0adNG7j6hUKjWOaoqOzsbLi4uOH78uMJ91tbWlTqXstegvNclOzsbCxYswIABAxTOVTpJKqPsA2/z5s0RFBSEX375BW+++SZu376NrVu3Vuo5GDJDzA+l0ZgLw6PJ/LA7oj3yihSvmV72FrAUKX48MzNW79pI+aHqyssPtra2lTqXPuaHf/75p1LPoar0vrjIyclBUFAQPvzwQ6Uv+Pbt2zF58mT89NNPaNOmDVasWIHQ0FDcu3cP9evX5/dLS0vD8OHDsW7dOq3GmyzO41sOOKOSlgtWXKSw34NU9eY1Z6X+L5Iq/9BZep/n2cpnCZj4e1S5jyNhDOfiXsDZRgRvBwvk5xdi+wOBQitJ2QKkNuI4DuYmlfvTaOhoiUUDAjFr5y1IGIOQ47BwwOto6GhZ6cdXp7gICAiAqakpEhMT0blzZ5X7TZkyBQKBAP/99x969eqF3r17o2vXrgBeLfIokbxa1End81bEz88PFy9e5AcBAsD58+f5n52cnNCgQQM8ePAA//vf/5Sew9/fH1u2bEF+fj5/YS19DmUaN26M4uJiXLt2DS1btgRQ8m1beno6v0+LFi3w9OlTGBkZ8QPiZKRSKTIzM+Hv748zZ85gxIgR/H1nzpxBQECAWs9flRYtWuDevXvw8fFRuY+xsbHc76Qio0ePxooVK/D48WN069YNbm5u1YrRkBhafiiLigvDo8n80NjZBveeZsp9acgBsDU3qdZK3pQflKtufgBKcoSfnx/Onj2LUaNG8dsNIT90794d7u7uyMzMrFac6tD74qJnz57o2bOnyvuXLVuG8PBw/pf8008/Ye/evfj5558xY8YMACWzNfXr1w8zZsxA+/bty328goICudH1sl9CUVERiooUi4SyYp9m8i0HnLDk5WWSIgg48C0A7wS54O+oZIVjR3f0xM9nHkLKSi4wyj5ihrX3xMZz5e/zv9Zu2HYpCSpqEZUm//FqNgSRkQCszMqiEsYQl5IJB3MjoKgIsvq7qKgIUOO1qQmlf0fq/M6KiorAGINUKlX7WxVVBrV0Q0cfezx8kQtPe3O42JhV6Zyy4kIWlzIWFhaYMmUKJk2ahOLiYnTs2BFisRhnz56FlZUVRowYwf8dnDlzBi1atMDUqVMxYsQIREVFwc7ODu7u7uA4Dnv27EGvXr1gZmYGS0vLCs8ri6n0a1b6f8YYPvroI0RERKBly5bo0KEDtm7ditu3b6Nhw4b8vvPmzcOnn34Ka2trhIaGoqCgAJcvX0ZGRgYmTZqE9957D59//jlGjx6NGTNmICEhAUuXLlV47NL8/PzQrVs3jBkzBpGRkTA2Nsa0adP47idSqRRdu3ZFu3bt0K9fPyxevBh+fn548uQJ9u3bh3feeQeNGzfGlClT8N577yEoKAjdu3fHv//+i507d+LgwYNyj6vqNVC1bfbs2ejbty/c3d3x7rvvQiAQ4Pr167h9+za+/PJLACUzghw+fBjt2rWDqakp7Ozs+PeEsuf83nvvYerUqVi3bh02bdpU7ntH9vspKiqS+wZQnWubPjK0/CAj21f2LWZhYSHy8vJgZKT3KRnAq/jLfc6UI+Qoyw8A0MDWDE8y8vhcbmoshJGg/G/5K8oRlB+0kx/69euHli1bYsKECRg1ahSaNWtmcPmhvPeOqvwAVD5HGMaVTIXCwkJcuXIFM2fO5LcJBAJ0794d586dA1DyAo4cORJdu3bFsGHDKjznokWLsGDBAoXtBw8eLLe/m0xGAcBBCAauVMtFIT5tUoxCKQdHEQPwCLte7iPDgcE9Nw7zmgPP8zmYCBiW31LcxyO/4n38ihP4fbKLgE0xis2IXV2kOJbMvTyW8UfL5BdLX24vXWAw/HrwAs5ZAc6CfMhezd//OYgUJoKjiMFWvmW0xuVLXk1nd/ToUYiE5TclGhkZwdnZGdnZ2RqZG9qCAwIcjAEUITOzesk0Kyur3PunTp0KKysrLFq0CAkJCbCxsUFQUBAmTZqEBw8eYPTo0Zg+fTp8fHyQmZmJyZMnY//+/QgPD8fPP/8MKysrzJw5EzNnzkRYWBjee+89/Pjjj+WeNzMzk/+WNSsrCwJBybdrshWGs7OzkZWVhQEDBiA+Ph7Tp09HQUEB+vTpg1GjRuHo0aP8B7LBgweD4zh8//33+Oyzz2Bubo6AgAB8/PHH/D5bt27F5MmT0bJlSzRu3Bhz587F8OHDkZOTo/Lblx9++AGffPIJQkJCUL9+fcydOxe3bt0CY0zuvF999RU+/PBDpKamon79+mjfvj3ftNy1a1csWrQIS5cuxaRJk+Dp6YkffvgBLVq0kHvcvLw8/rasKb90bGVfq3bt2uH333/HN998g2+++QZGRkbw8/PDsGHD+GMWLFiA2bNnY/369XBxccGNGzdQUFAAiUSi9DlzHIc+ffrg4MGD/LeOqt47sg+xJ0+e5Kc3LB1nbaKP+aGs0t+07tq1q0rn0KVDhw6pvE+Yn4+3X/584MABSCro1lFTdJkjlOUHIQAXcyBfAqQVAPlFEqSkZcJMjU9n5eUIyg+azw/m5ubIyspC7969DTY/yN4zyt47qvJD6VjVxTF1+mDoCY7j5PrUPnnyBK6urjh79qzcQJ7PPvsMJ06cwIULF3D69Gl06tRJbmDNli1b5AbVlKbsmyl3d3ekpqaq3R/7zytJmL07Gvkp8Uje+Als6jng+dMnSveRtWZ89U4ABrV00/g+yeJ8hLwcNyEj4IDjUzoBABLTcpGaXYhP/7ih8DxCXKQ4+VSgtAWEAzDg5mEYSYrxR7O3wCD/+Mni/FLf0NRcUskpzIHd0pLBWc8mPoOthW25++fn5+PRo0fw8vKqsE9jTWGMISsrC1ZWVgY5MFff4k9KSoKnpycOHjyIbt26lbuvvsWurh49eiAgIAArVqwoN/78/HwkJCTA3d1d7v2emZkJBwcHiMXiSo870ReGkh+Akm8BDx06hO7du8PKygqMMSQmJsLZ2bmSz1o3ZPH36NFDrg+5nJwcGL8cKFuUng7U0CD5iuhzjkjOzMeL7EKYGgnhW1/162Wo1ylA/2KvTH4A9C9+dcjyw8qVK8uNX1V+ACqfIwy65UIdHTt2rFRTpqmpqcLAJKCk+VrlRbSMoW290dHHAZE/J+IbAJAUKRw7tK03uvg7IyE1F14O5krHMWhiHw8HY6V9PT0crF7eb4VkcR7fbUtGwAFdXKT44n8heCwuhKkRh3dXn5Mb37EjsLvcY0kZ8PmuaCS8yMf60w/kBoJ38nPkZ8bS5pgNY/bqdVbndyaRSMBxHAQCAf8ti67J3q+yuAyNruM/evQosrOzERgYiOTkZHz22Wfw8vJCSEhIhfHoOvbKSk9Px/Hjx3H8+HH8+OOPfLJQFb9AIADHcQp/G+pe22obXeSH0kxMTGBubo6cnBwUFSnmCX1X7vMu+/7Sk+emzznCyUoEcW4RCooleJ5dCAsTIUyMhArjLwztOlWarmOvTn4AdB9/ZZTNDwKBoNz4VeUHoPI5wqCLCwcHBwiFQqSkpMhtT0lJqfY3QJGRkYiMjKzUwJnSXGxE8LEreXmVrTxZso9ZhR+0NbHPkGAPdPJzVFmAuNiYKRQgX77jD4uUG3CxEcHDwQpn41KVju8oiwFYe+oBf1vKgOk7bvKdr0rPOpUszquRgoPULUVFRZg1axYePHgAKysrtG/fHr/99pvBfXBTR/PmzZGeno4lS5agcePG1R43VJvoc34oTVZc1MauaaRyjIQCOFmL8DgjDymZJZ8bOACudmaoZ6HjPse1RF3ODzXJoIsLExMTtGzZEkeOHOGbwqVSKY4cOYLx48dX69wRERGIiIhAZmYmP69yVeIDSprSGWM6bUKrbAHiYG6EfftedZXydrBQbN2QSgBwkKpRvZeedWrGjptISM3BmpO6ad0gtVtoaChCQys3VaShSkhI0HUIekvf84MMzRhFSrMqMwUtA/A4PR+WpsbVmkGKlKD8UDP0vrjIzs5GbGwsfzs+Ph5RUVGoV68ePDw8MHnyZIwYMQKtWrVC69atsWLFCuTk5MhNEaYrpSvhwsJCpc3p+qR0AVJ2ZgDF1g1g4YEfAACzek2EhAFCjsNnbzXGkv13y52pigFYfUK91g1CCFHFkPODjGwSAdmgV1K3FRYrtj4yMBQWS6m4IAZD74uLy5cvo0uXLvztyZMnAwBGjBiBTZs2YciQIXj+/Dnmzp2Lp0+folmzZti/fz+/GmJVaaLZu3RxUVBQoPfFRUXkWjfMAZfFJXOBdNq/DQm54Ltc2Zob80WIACUFg7pdqoDy19SgrlSEEBlDzg8y1HJBSjMxEipMM88BVFgQg6L3xUVISEiFi4qNHz++2s3cZWmi2bt0cZGfn2+ws7CUxrdulPqWzcVaBBeXVzNblO1idfL+c7nxHOq0bkgYQ0JqrlwBIVv5vK6tGE4IUc6Q84MMFRekNBMjAVztzPA4/dXaF3YW1VtUj5CapvfFhSHjOA4mJiYoLCyUm76wLijdxUrZgHJ1WjdKD1EpvfI5ULdWDCeE1F5UXJCy6lmYwtLUGCmZ+UjPLURWfjGkUgaBwDCmPiWEigstE4lEKCwsVDljVF1RdkB5ea0bMp/9dQM/fdACL3IKse1CIsp+QamsdYMQQgwJFRdEGRMjAVxtzZBTUIxCiRTPsgvgbK0fazERUhEqLlTQVJ9a2TiLutZyoQ5VrRs2Zkb46NerSEzLRa9Vp8s9h5sdXWwJITVLk2MuaEA3UUUg4OBiI8LDtFw8zyqAyEgAcxMjUA8pou/oLapCREQEoqOjcenSpWqdR7bKYV1vuVCHi40Z2jWyR0ADG3wzsKnSfcq2Cm+/lKSwj5DZw1QSiJRMKugqsmnTJtja2uo6DEIMiqbyA0AtF6R81mbGMDUSlqzinpaLe08zkZ5bVPGBGkD5gVQVFRdaRi0XVSNVMUhz1XvNsS28LRa80wQA8MOxWOy4moSzcalIFudhx5VkuOZvhHPhIvRcdRnbLyXWZNgGZ8iQIbh//36ljgkJCcGnn36qnYD0CMdx2LVrl67DILUcFRekPEUShoLiVy1kDMCTjDwombFW4yg/qCYUCik/lIO6RWmZrLiglovKUbZon5Dj0NLLjm/hiEnJwq/nEzHlj+tyx3Iva2Ya9F0xMzMzmJnRa1NVhYWF/GKZhFQFFRekPIXFil3vGFAjxQXlh+qpy/mBWi5UiIyMREBAAIKDg6t1Hiouqka2aJ/w5ZRRQo7DwgGvyxUJozs2rPA8skHfNSUpKQnHjh1DUpJidy1NCwkJ4afZtLGxgYODA+bMmSM3NWd6ejqGDx8OOzs7mJubo2fPnoiJieHvL9vsPX/+fDRr1gxbtmyBl5cXbGxs8N577yErKwsAMHLkSJw4cQIrV64Ex3HgOE7lKqDJycno3bs3zMzM4O3tja1bt8LLywsrVqzg98nIyMDo0aPh6OgIa2trdO3aFdevyxeLq1evRqNGjWBiYoLGjRtjy5YtcvdzHIc1a9bg7bffhrm5Ofz9/XHu3DnExsYiJCQEFhYWaN++PeLi4uSO2717N1q0aAGRSISGDRtiwYIFKC4uBgA0bFjy3urfvz84joOXl5fc67N+/Xp4e3tDJBLhl19+gb29vULrZL9+/TBs2DAVvz1iyDSVHwAqLuqSquQH2boXpXFAheMuKD+8fK20kB+aNi3ptk35QTUqLlTQ9JgL6hZVeUOCPXB6RhdsC2+L0zO6KKxp8UScp9Z5Dt5+iqT0XL7rVEUYY8jJyan0vx9//BGenp7o2rUrPD098eOPP1b6HBXN2V/W5s2bYWRkhIsXL2LlypVYtmwZ1q9fz98/cuRIXL58GXv27MG5c+fAGEOvXr0UVmAvLS4uDrt27cK///6Lf//9FydOnMDixYsBACtXrkS7du0QHh6O5ORkJCcnw93dXel5RowYgSdPnuD48ePYsWMH1q5di2fPnsntM2jQIDx79gz//fcfrly5ghYtWqBbt25IS0sDAPz999+YOHEipkyZglu3bmHs2LEYNWoUjh07JneeL7/8EsOHD0dUVBRee+01DB06FGPHjsXMmTNx+fJlMMbk1jo4deoUhg8fjokTJyI6Ohpr1qzBpk2bsHDhQgDAhQsXAAAbN25EcnKy3HUgNjYWO3bswM6dOxEVFYVBgwZBIpFgz549/D7Pnj3D3r178eGHH6r+5RGDpckxFzSg27DUdH4oKsiDrYkUebmv8oOjlalag7opP5TQVH74+uuvAQBHjx4FQPmhXIyUSywWMwBMLBZX6rjCwkK2a9cu1rVrVwaA/fbbb1qKUDtk8RcWFirfITubMaDkX3Z2zQb30pOMXOY941/mOf3VP6/p/zKP6bv5nz3L/POe8S/7/eJD/hx5eXksOjqa5eXl8duys7Nly27U+L/s7GwmkUhYeno6k0gk5T7/zp07M39/fyaVSvlt06dPZ/7+/owxxu7fv88AsDNnzvD3p6amMjMzM/bHH38wxhjbuHEjs7Gx4e+fN28eMzc3Z5mZmfy2adOmsTZt2sg97sSJE1XGJZFI2IULFxgAdunSJX57TEwMA8CWL1/OGGPs1KlTzNramuXn58sd36hRI7ZmzRrGGGPt27dn4eHhcvcPGjSI9erVi78NgM2ePZu/fe7cOQaAbdiwgd+2bds2JhKJ+NvdunVjCxculDvvli1bmIuLC//aA2B///233D7z5s1jxsbG7NmzZ3LbP/74Y9azZ0/+9nfffccaNmwo97upCRW9d5S93xmr+nWurqtufigsLGRr165lAFjfvn21FKXmVZgfGNOLHKFMdkE2w3wwzAdLz06vcP+yfzO6zA/XHzxl1x+ls+SM3ApzBOWHEprOD7JrrCHmB8bKzxGq8gNjlb/WUcuFltGAbu1R1nVqwTuN8Vj0IZ6azMT+ia0wv0+A3DGycRjqtGAYgrZt24Irtdpgu3btEBMTA4lEgjt37sDIyAht2rTh77e3t0fjxo1x584dlef08vKClZUVf9vFxUXhG6WKxMbGwsjICC1atOC3+fj4wM7Ojr99/fp1ZGdnw97eHpaWlvy/+Ph4von6zp076NChg9y5O3TooBC/rJkaAJycnAAAgYGBctvy8/ORmZnJP/YXX3wh97iyb9sq6p7i6ekJR0dHuW3h4eE4ePAgHj9+DKCkO8HIkSPlfjeEKEPdooi67CxK+u9n5hertT/lhxKUH2oeDejWMhpzoV1lF+OzNpNixMEXkAhfwMnaFH7OQoVjJIzh4oM0OFqbws1a8U/A3Nwc2dnZlYrj8ePH8Pf3h1T6apSdUChEdHQ0XF1d1T6Publ5pbtGaZqxsbHcbY7j5J6XpmRnZ8PFxQXHjx9XuK+y0x+Wjll2wVa2TfY8srOzsWDBAgwYMEDuPFKplO/KqIqsG0tpzZs3R1BQEH755Re8+eabuH37Nvbu3Vup50DqJiouDIsu84OJqQgZT7OQXyRBsXHF+2tDXc4PQElX9/J+/5QfSlBxoWXUcqF9pRfjyymU77esbNYpAJi4PQoA4GYtxIpeJRf3wmIpCoslMDESKr1AlMfPzw9r167F2LFjIZFIIBQKsWbNGvj5+VX6+VSmuJCNDZA5f/48fH19IRQK4e/vj+LiYly4cAHt27cHALx48QL37t1DQECAstOpxcTEpMLFw3x8fFBcXIxr166hZcuWAEq+rUpPT+f3adGiBZ4+fQojIyN+QFxZ/v7+OHPmDEaMGMFvO3PmTLXilz32vXv34OPjI7ddKpXy314ZGxtXapG00aNHY8WKFXj8+DG6d++usq8xIaVRcWFYOI7TaX4wNzVCTkEx8tRovKD8UDWq8gPwqgCh/FA+6halgqZmA6FF9HSrbNepso2QjAHpuUVISM3B3aeZeJCag3tPM5GWU/liMCwsDAkJCTh27BgSEhIQFhamgWdQvsTEREyePBn37t3Dtm3b8P3332PixIkAAF9fX7zzzjsIDw/H6dOncf36dXzwwQdwdXXFO++8U+XH9PLywoULF5CQkIDU1FSl31r5+fmhW7duGDNmDC5evIhr165hzJgxMDMz478l6t69O9q1a4d+/frh4MGDSEhIwNmzZ/H555/j8uXLAIBp06Zh06ZNWL16NWJiYrBs2TLs3LkTU6dOrXL8ADB37lz88ssvWLBgAW7fvo07d+7g999/x5w5c+Se55EjR/D06VO5pKfK0KFDkZSUhHXr1tXugXpEo7NF0YDuukFT+cFaVPKNe64an2spP1SNqvwwe/ZsuedJ+UE1Ki5U0NRsINRyoXulZ536fmhzpfuUXaTocXoeCqswkbibmxtCQkLg5uZW1XArZfjw4cjLy0Pr1q0RERGBiRMnYsyYMfz9GzduRMuWLfH222+jXbt2YIxh3759Ck3blTF16lQIhUIEBATA0dERiYnKFyrcvHkznJyc0KlTJ/Tv3x/h4eGwsrLiC26O47Bv3z506tQJo0aNgp+fH9577z08fPiQ7xfbr18/rFy5EkuXLkWTJk2wZs0abNy4ESEhIVWOHwBCQ0Px77//4uDBgwgODkbbtm2xfPlyeHi8mpHsu+++w6FDh+Du7o7mzZW/b0qzsbHBu+++C0tLS/Tr169a8RH9Rit0k6rQRH6wMSvpcFIgASRlm+TLoPxQNaryg6enJ7/Pt99+S/mhHBzTdQdvPZeZmQkbGxuIxWJYW1urfVxRURH27duHo0ePYtWqVZgxYwYWLVqkxUg1SxZ/r169lF9ocnIAS8uSn7OzgUo2E2tLTmEOLBeVxJU+NR22FrZy9yeL89Bh8VG+m5SrlRDzu9RH/QZu4IzkF7txszOHpakR31XKRJ25/zRA1jXH2toaAoHqxwwJCUGzZs3k5gXXB6riT0pKgru7Ow4fPoxu3brpMELV1H3tVenWrRuaNGmCVatWaSG6ilUUf35+PuLj4/k52GWqep2r66qbH3r16oW4uDj4+/vDzs6On2JT31WYHwCDzRFlqfqb0aX7KSXjLtzszFDPwlTpPpQftKM6OULX+QEoP/7y3uuVvdbRmAstk63OSC0X+kHWTWrWzluQMAYBx8HGzEihuxQAPE7Phazy5gC4lnMhJ4qOHj2K3NxcBAYGIjk5GZ999hm8vLzQqVMnXYemcenp6Th+/DiOHz+OH3/8UdfhEANCLReksqxFRsgvkiAzr9hgcxLlh9qNigstozEX+qf0DFOu1kLkpCbDzlqE57kMDAwcABMjAQpKdYsq6SqVD0tT4xprwTB0RUVFmDVrFh48eAArKyu0b98ev/32W7Wa3PVV8+bNkZ6ejiVLlqBx48a6DocYEFlxUVBQwA/2JaQ81iJjPMsqQHZBMbLyimBqXHMt65pC+aF2o+JCy6i40E+yGaby8/MRnwrYmpugnrUJCoulMDESoLBYggep8gMsGRh/vz5QNkWfPgkNDUXPnj11HUaNSEhI0HUIxECVnnkoNzdXbg0BQpQRGQtezoLIEP8iR2nLOuUH/VEX84N+fErSQ5qaDYQGdBsOEyMBLEVGMDESwMRIqNBVStaiUVZhsRTZ+UVVGgBOCDE8mpwtqnTfZuoaRdRRJGFy06vLWtYpBxF9QcWFCpqaDYRaLgxD2XkNTIwEcLUzkyswbM1NFIqLtJwC3KvmFLaE1BSav0MzNDlbFMdxNO7CAOjT306hRLGIkLWsE1JVmnyPU3GhZdRyod9k/TuVJfV6FqZo7GwNe4uSQfmZ+UUoKnVRLyyWICk9jx/0Td8eEX0ne5/Xxn7NhoyKC/1VXo7QFROh4kc3DpzedNklhkmT+YHGXGiZrLiglgv9JBQKYWtri2fPngEoSfKyRXxk6ok4ZOVKUFAowcNnEtiZGwPg8CyrAKxYfiUjBiArRwgL06r/aUmlUhQWFiI/P79K06HqmiHHb8ixA6rjZ4whNzcXz549g62tLQ0a1jNUXOgvdXJETZNKpbAWFEJcWHKbA+BoLYK0uBD5aqzcrUu19RprKJTFr438QMWFllHLhf5zdnYGAD55KFMkkeJ5ZgGeAYgv51wcAGGOCEJB1ZMPYwx5eXlyq5UaEkOO35BjByqO39bWln+/E/1Bq3TrN3VyRE2S/Z3nSY2QVySBpakQ6XkmqHidaN2r7ddYfVde/JrMD1RcaBmNudB/HMfBxcUF9evXR1FRkdJ9nmXl4+Pd5xW2D23tga2XEiHrGzW+qw86Nqre6txFRUU4efIkOnXqZJDdVww5fkOOHSg/fmNjY2qx0FPUcqHf1MkRNUn2d55p54cfjsWhpacdvhloGFOc1uZrrCFQFb+m8wMVF1pGLReGQygUqvzjSnqcjcdZEoXtLRo64Z1WXhj00zkki/ORmsuUruKaLM5DfGoOvB0s4GJjVmEcxcXFEIlEBnnxMuT4DTl2wPDjr6uouDAM5eWImo6juLgYrzWoh8dZ95HzQAxTU1OD+Cbd0K9RFL96DK/DmIGhlovawdvBAmV7Ogk5Dl4O5nCzM8eMnq8BANadeoDMfPlvtn6/mIj2i49i6LoL6LD4KLZfSqypsAkhBoCKC1IVfk6WMBEKkJFbhMQ0eu8Q/UHFhQq0zgUpzcXGDIsGBEL48pshIcdh4YDX+VaIt5s2gG99S2TmF+Pn069GZdx/moUZO29CNsOblAGzdt5Csjivxp8DIUQzNLnOBUDFBakaEyMB/F1KFl28niTWcTSEvELdolSIiIhAREQEMjMzYWNjU+Xz0GxRtceQYA908nNEQmouvBzM5bo3CQUcPu3uh4itV7Hu5AM0cbEBgxSzdt5SOI+EMSSk5lbYPYoQop80lR9kaEA3qaqmbra4niTGzaQM9A1qoOtwCAFAxYXWUbeo2sXFxkxlUdDzdWc4W5viaWYBwrdcVnkOWXcqQggBqOWCVF2gW0lxSy0XRJ9Qtygto25RdUdKVj5SMuV/zxyAyT185cZrzOj5GrVaEEJ4VFyQqgpyswUA3HoshkSqP6uIk7qNigstk7VcFBQUaHRpdaJ/4lNzUPY3zAAEe9njzIyu8HOyBACI83Q/lSEhRH9QcUGqyqe+JcxNhMgtlODB82xdh0MIACoutE7WcgEAhYWFOoyEaFt5M0q52Jhhcg8/AMCvFx4it1DPl1ElhNQYKi5IVQkFHF5vQF2jiH6h4kLLSq95QOMuareKZpTqEeAMT3tzZOQWYceVJF2GSgjRIzSgm1SHbNzFjaQM3QZCyEs0oFvLTExM+J9p3EXtV9GMUh928Ma8Pbex4XQ8hrbxhLBsUwchpM6hlgtSHU354oJaLoh+oJYLLeM4jqajrWNcbMzQrpG90kHbA1u6wVpkhIQXuYg8FkvrXRBCqLgg1SIb1B2dnInCYqlugyEEVFzUCJoxishYmBqhhacdAGDZofu0YjchhIoLUi2e9uawFhmhsFiK+ylZug6HECouagKtdUFkksV5OHn/OX+bVuwmhFBxQaqD4zg0fdl6QV2jiD6g4kKFyMhIBAQEIDg4uNrnKj0dLanb4lNzUHYqctmK3YQQw6DJ/ADQgG5SfU1pUDfRI1RcqBAREYHo6GhcunSp2ueiMRdERtl0tQBw/sELJIvzcDYuFcliep8Qos80mR+AVy0XaWlpSEqimeRI5claLs7GvaCWcKJzVFzUAOoWRWTKTlcrqzNWHolB+0VHMXTdBYR8dxLnUmgWKULqiv/++w8AkJqaCk9PT2zYsEHHERFD8/BFSatXYloujeUjOkfFRQ2gAd2ktCHBHjg9owu2hbfFmRldMLqjNwDwq3tLGbD9gYBaMAipA5KSkjBv3jz+tlQqxdixY6kFg6gtWZyHJfvv8rdpLB/RNSouagC1XJCyZNPVNrA1R1f/+gr3M3BITKNxGITUdjExMZBK5acPlUgkiI2N1VFExNDQWD6ib6i4qAHUckHKo2wcBgcGj3rmugmIEFJjfH19IRDIp2KhUAgfHx8dRUQMjbIcIuQ4eDlQDiG6QcVFDaCWC1KeV+Mw5LdfS8x4OcCbmrYJqa3c3Nywdu1a/rZAIMCaNWvg5uamw6iIIZHlEFkK4QAsHPC60oVcCakJRroOoC6glgtSkSHBHujk54gHz3Ow9kQsTsS8wMQ/bgAABBywaEAghgR76DhKQog2hIWFYfPmzTh16hSWLl2KsLAwXYdEDMyQYA/kFUow/59oNPewpXxBdIpaLmoAtVwQdbjYmKGDjwPm9H4Nr4Z30+A8QuoCPz8/AEB2draOIyGG6jUXawBAem6RjiMhdR0VFzWAWi5IZSSLC/BqktoSNDiPkNrNw6Pkm+bERJpClFRNg5fdoJLFeWCMVbA3IdpDxUUNoJYLUhme9ubgIJ8YOIAG5xFSi1FxQaqrvvXLBXuLpMig1guiQ1Rc1ABZcUEtF0QdLjYiDGkolZv9g6FkFW9CSO1ExQWpLpGxEPYWJgBA6yQRnaIB3TVA1i2KWi6Iuto5MYwb0AmPxYU4FP0UP59JwMydN+FgaQqhgIO3gwXNBEJILVK6uGCMgeO4Co4gRJGzjQgvcgrxNDMPAQ2sdR0OqaOouKgB1C2KVIWLjQgeDlZo7V0Psc9zcPL+cwzbcBEAzSBFSG0jm3o2NzcXaWlpsLe313FExBC52Ihw+0kmtVwQnaoT3aL69+8POzs7DBw4UCePTwO6SXUIBRxm9nxNbhvNIEWIZug6P8iIRCI4OzsDAB4+fKjTWIjhcrYp+TLzKRUXRIfqRHExceJE/PLLLzp7fGq5INWVnluosI1mkCKk+nSdH0qjcRekulz4GaPo8wbRnTpRXISEhMDKykpnj08tF6S6vB0s5AZ4y0gZo1W8CakGXeeH0qi4INXlbE0tF0T39L64OHnyJPr06YMGDRqA4zjs2rVLYZ/IyEh4eXlBJBKhTZs2uHjxYs0HWg5quSDV5WJjhkUDAiEsM8jzf+svYOi6C+iw+Ci2X6IPJKRuqQ35oTQqLkh1udiWfN6gL5yILul9cZGTk4OgoCBERkYqvX/79u2YPHky5s2bh6tXryIoKAihoaF49uxZDUeqGrVcEE0YEuyB0zO6YFt4WywfHCR3H43BIHVRbcgPpVFxQaqrdLcoWkiP6Ipas0XVq1cP9+/fh4ODA+zs7MqdIs/S0hJNmjTBkiVL0LRp02oH2LNnT/Ts2VPl/cuWLUN4eDhGjRoFAPjpp5+wd+9e/Pzzz5gxY0alH6+goECuCMjMzAQAFBUVoahI/UVpZPsWFRXByKjkZc7Ly6vUOXSpdPwqdoBx6X315HmVjreyvzN9Ud5r72BuBAcPa5x/UKxwn4Qx/HHpIfo1bYDE9Dx42pvD5eXgvppU4XtHjxly7EDV46/u89VVjqgN+aG0Bg0aACgZ0K3P70G13meUI7SmvNff3kwIAMgtlCAtKw/WZsYK++hSXb3G6ouayhFqFRfLly/n+6SuWLGi3H0LCgqwb98+jBo1CleuXKlUMJVVWFiIK1euYObMmfw2gUCA7t2749y5c1U656JFi7BgwQKF7QcPHoS5eeVXSD506BBu3rwJAHj27Bn27dtXpbh05dChQ0q3C/Pz8fbLnw8cOACJqOY/xCqTL3nV9ezo0aMQCfUjrqpQ9doDQEYBwEEIBvkPccsPx2H54VgAHDgwDGkohb8tw/N8Do4iBltTLQddSnnx6ztDjh2ofPy5udWbGEAfc4Sh5IfSZC0WMTExBpErynufUY7QPlWvv7mRELnFHP7YewgNKv+2rBF17Rqrb7SdIzimhXazR48eoWXLlhpveuY4Dn///Tf69esHAHjy5AlcXV1x9uxZtGvXjt/vs88+w4kTJ3DhwgUAQPfu3XH9+nXk5OSgXr16+PPPP+X2L03ZN1Pu7u5ITU2FtbX6C9IUFRXh0KFD6NGjB06dOoW33noLTZo0wbVr16rwzGte6fiNjZV885GTA2M7u5J909MBC4sajlC5nMIc2C0tievZxGewtbDVbUBVUOFr/9KfV5Iwe3c0pKxk3YtgT1tcSMhQ2I9DyQrfAg746p0ADGrpprXYAfXj10eGHDtQ9fgzMzPh4OAAsVhcqetcVWkjRxhqfij9e3r+/DlcXV0BAFlZWXyXWn2j1vuMcoTWVPT69/nhLO6mZGPD8Bbo5OuggwhVq6vXWH1RUzlCK4voubu761Wf1sOHD6u9r6mpqdILurGxcZXeSMbGxrC0tARQ8k2aob0ZVT7vUtuMjY3lbuuSMZOPy9Be79Iqin9oW2908XdGQmouvBzMEZ+ag6HrLijsJ/v2QMqA2buj0bqhIyxMhYhPzdHqSt+G/PobcuxA5eOv6eeqTzlC1/mh9HEuLi4QiUTIz8/Hs2fP0LBhw0qfsyaV+7wpR2idqvhdbM1wNyUbz7OL9Pb51dbX3lBoO0dofMxFWlpapQKoDgcHBwiFQqSkpMhtT0lJ4RcjqqrIyEhERkZCIpFU6zzAqwHdNFsU0TQXGzO54kDAlRQRqkgZ8ObyE/w+pVf6Thbnab3gILWTPuYIQ8kPpXEcBw8PD9y/fx+JiYl6X1wQ/eRMa10QHdP4mIuaZGJigpYtW+LIkSN8U7hUKsWRI0cwfvz4ap07IiICERERyMzMhI2NTbXORVPRkpogm6521s5bkDAGAUpaLcrWGqWLDykDpu+4ib8uJ+FyYjoYe1VwdPJzpGKDqEUfc4Sh5IeyShcXhFSFC63STXRMreJixIgRSn8uz+LFi/HRRx/B1ta2SoHJZGdnIzY2lr8dHx+PqKgo1KtXDx4eHpg8eTJGjBiBVq1aoXXr1lixYgVycnL42UH0AU1FS2rKkGAPdPJz5LtKnbz/nC82hByHDzt6Yd2peIXjLj1M53+WFRylx2pQ6wYpj65yRG3ID2V5enoCoOloSdXJiovkTCouiG5oZcwFACxcuBCDBw+udnFx+fJldOnShb89efJkACUJbNOmTRgyZAieP3+OuXPn4unTp2jWrBn2798PJyenaj2uJpu9qeWC1KTSXaXKFhsAsOF0vFzrhayIKKv0WI0ZO27i+qMM/H7pET+AnFo3SHVoIkfUhvxQlmyti4cPH2r83KRu4Ne6yKB1j4huaK240NQkVCEhIRWea/z48dVu5i5Lk83espaLwsJCMMbK7Y9MiKaVHZdRuuuUkOPw2VuNsWT/3XLHajAAWy8+4m/LCg5wkOtKVdK6kY8YMYdkcT48HAx3wBvRLk3kiNqQH8qihfRIdTlTtyiiY1orLsgrolLzexcUFMjdJqSmlW3NcLExg625cYVjNcpieLWTrNi4kSTGtouJkDIhfrxzklo3CKkkKi5IdcmKi6yCYmTlF8FKRF/ykJpFxUUNKD11IRUXRB+Ubc2oaKyGuq0bv1149YGIH7uhtHWDxm4Qokzp4oJauklVWJoawUpkhKz8YqRk5lNxQWocFRcqaLJPrYmJCf9zfn6+xpvRCdGE8sZqVLV1AygpLIBXrRuHolNw5O4zmpmKGCxtjrlwcytZ5DI3NxdpaWmwt7fX+GOQ2s/FRoSs/Gwki/PhU99K1+GQOoaKCxU02aeW4zh+YSSaMYoYCm21bhy+82rxNGUzU33dPxDvt1Zs3aDWDqIvtDnmQiQSwcnJCSkpKUhMTKTiglSJs40Z7qdk01oXRCe0Vly88cYbMDOjDwAypqamyM/PpxmjiEFTp3Vj5s6bkLKSWagANVs3Xv4vZcDMnTfx/ZEYPHmZFDkOeDvQBXtvJsvNVEXdqwwb5QjVPDw8+OKiefPmug6HGCAXaxrUTXSnysXF7du35ZqEhUIhmjRpwt/et29f9SKrZUQiEcRiMbVckFpFWetGO287/LHvGAb36oJz8emVbt0AwBcWQEm3qn9uJPO3pQyYsfMmbiaJsfViYoUFBxUgukE5ouo8PDxw6dIlGtRNqszF9uVaF1RcEB1Qu7g4deoUJk+ejEuXLgEA2rZti9zcXH4aQI7jcODAAXTv3l07kdYwTfeplQ3qppYLUtu52Ijga8PgYiOq0tgNVetulMYY8GuZweMzdtzE4egUHC41nqN/c1f8fe0xtXjUgLqUI7Q55gKgGaNI9fEL6YlprQtS89QuLn788UcMGzZMbtuxY8fg6ekJxhhWrVqF1atX14rEAWi+Ty0tpEfqKm2M3VCGAThUZjzHjquP5W7P3HkTKZn5WHE4RuligG42r2Z2owKkcupSjtDmmAuAVukm1ef88ppF3aKILqhdXFy+fBmff/653DY3Nzf+Ijhs2DD07t1bs9HVIrKWC+oWRUjlZqYSchz6NW+AXdeeVLsAkTJg2aEYudtlB5QP9uaQcyUJs3dH02rklUA5QnOo5YJU16uWCyouSM1Tu7hISkqS+4Zm8+bNcHZ25m/Xq1cPL1680Gx0tQi1XBCiWkWtGy42Zpga2rjcAqSqBQcgP6D89wcC/P4gmr9P1uWK40DjO8pBOUJzZMXFw4cPdRwJMVSyhfTEeUXILSyGuQlNDkpqjtrvNisrK8TFxcHd3R0AMGDAALn74+PjYW1trdnoahFquSCkcsoWHOoUIJpp8VBctIxBcb2OPVFPcDbuBd/i0a+5K3apMb6jthYglCM0R1ZcJCcno6CgQG4hVkLUYWVqBAsTIXIKJXgqzkdDR0tdh0TqELWLizZt2uCXX35BSEiI0vs3bdqENm3aaCoundP0gD1quSBE86rb4qF8MUAGDly5g8oZgDNxr76FlzJgZ5nxHTN23MTZuBfYc/1JnRhgXpdyhLYHdDs4OPBrIz1+/BgNGzbUyuOQ2ovjODjbiBD3PIeKC1Lj1C4uJk+ejO7du8Pe3h7Tpk1D/fr1AQDPnj3DkiVL8Ouvv+LgwYNaC7SmaWtAN7VcEKJdlW3xKD2gvGTMhRRNmwZizu47lV6NvDQGYHfUE/62sgHmM3bexO3HYvx6wfCn1K1LOULbA7o5joOHhwfu37+P3bt3Y9CgQfzK3YSoq4GtGeKe59C4C1Lj1C4uunTpgu+//x6TJk3CsmXLYG1tDY7jIBaLYWRkhBUrVqBr167ajNWgyb7hSk5OrmBPQoi2qRpQ7mpjgmtnjqJXSzd08XfW+IxWZTEG/HJefkrd6TtuYsu5h7j9JBMMJYsI9vCvj8N3nqlRgOQjRswhWZwPDwfj6gVXSZQjNEsoFAIoKdqmTp2KtWvXIiwsTMdREUPibE3T0RLdqNQIn3HjxqFPnz7466+/EBNTMuOKr68vBg4cyPezJYo2bNiA//77DwAwZ84cODs7U5IgRI/Iio2ioiJcK7MNqNkZrQDg1pNM/mfGgIPR8lPsTt9xE79fSkRUopgvQEIDnHAwOgVSJsSPd07yBUhNohyhGUlJSbh79y5/WyqVYuzYsQgNDaUWDKI2mjGK6Eqlpw9wd3fHpEmTtBFLrZSUlIQxY8bwtxljlCQIMUCamNFKkwXItUQx/zNjwP7bKfxtKQNm7byFTn6ONd6FinJE9cXExPCLD8pIJBLExsZS3iBqk611Ef0kE8niPL3sTklqJ5qbTMtiY2MhlUrltlGSIKR2qMqMVtqcUrc0CWNISM2lDxQGyNfXFwKBQC53CIVC+Pj46DAqYmhiUrIAANceZaDD4qM6ac0kdZNA1wHUdj4+PhAI5F9mShKE1B0uNmZo18ie/5Bf9vaQYA+cntEF28Lb4vSMLhjbuREWDQiEkCuZElfIcXi3havc7Zk9X4NAccZcOUKOg5eDufaeGNEaNzc3LFu2jL8tFAqxZs0a+kKKqC1ZnIdN5xL427LWTBp/QWoCtVyooKmpBt3c3LB27VqMHTuWP9f3339PSYIQwtNUl6vS09wuHPA6tVpoibanogWAiRMnYunSpUhKSsLWrVsxePBgrT0WqX3iU3NQpmcdtWaSGkPFhQqanGowLCwMb775Jlq0aIHU1FQEBARoKEpCSG1VlS5XE7s2wh/7jmFwry7wcLDSRdh1granopVp06YNkpKSkJiYWPHOhJTi7WABAQe57pXUmklqCnWLqiHu7u7o1q0bAOD06dM6joYQUhsodrkSwdeG8bPEEMPWsmVLAMCVK1d0HAkxNC42Zlg0IFCu++QX7zShVgtSIypdXNjZ2aFevXoK/+zt7eHq6orOnTtj48aN2ojV4HXs2BEAFReEkNqLcoTmyIqLy5cv6zgSYoiGBHvg5GddYC0qWfPGmb50IDWk0sXF3LlzIRAI0Lt3byxYsAALFixA7969IRAIEBERAT8/P3z88cdYt26dNuI1aB06dAAAnD17Vqt9dQkhRFcoR2iOrLiIjY2FWCyuYG9CFLnZmePdlq4AgN1RT3QcDakrKj3m4vTp0/jqq6/w0UcfyW1fs2YNDh48iB07dqBp06ZYtWoVwsPDNRZobRAYGAgrKytkZmbi1q1bCAoK0nVIhBCiUZQjNMfe3h6enp54+PAhrl69ii5duug6JGKA3mnmio1nEnAoOgU5BcWwMKXhtkS7Kt1yceDAAXTv3l1he7du3XDgwAEAQK9evfDgwYPqR1fLGBkZoV27dgCoaxQhpHaiHKFZrVq1AkDjLkjVBbnZwMveHHlFEhy+k1LxAYRUU6WLi3r16uGff/5R2P7PP/+gXr16AICcnBxYWdFMJcrQuAtCSG1GOUKzaFA3qS6O49C3GXWNIjWn0m1jc+bMwccff4xjx46hdevWAIBLly5h3759+OmnnwAAhw4dQufOnTUbaQ3T1jzmsuLizJkzGj0vIYTog7qQI2pinQsZGtRNNKFvUAOsOhKDk/efIy2nEPUsTHQdEqnFKt1yER4ejhMnTsDCwgI7d+7Ezp07YW5ujhMnTiAsLAwAMGXKFGzfvl3jwdakiIgIREdH49KlSxo9b+vWrWFkZIRHjx7R3OWEkFqnLuQIbeUHZWhQN9EEn/qWeN3VGsVShr03k3UdDqnlqjSqp0OHDvzMR6RyLCws0Lx5c1y6dAmnT5/G0KFDdR0SIYRoFOUIzaFB3URT3glyxa3Hmfj1XAK6+9enNS+I1lSpuJBIJNi1axfu3LkDAGjSpAn69u0LoVCo0eBqq44dO1JxQQiptShHaFbLli3x8OFDXLlyhYoLUmVSlCzXfS8lGx0WH8WiAYEYEuyh46hIbVTpblGxsbHw9/fH8OHD+SbvDz74AE2aNEFcXJw2Yqx1aFA3IaS2ohyheTRjFKmuZHEelvx3l78tZcCsnbeQLM7TYVSktqp0cTFhwgQ0atQIjx49wtWrV3H16lUkJibC29sbEyZM0EaMtY6su8CtW7eQnp6u42gIIURzKEdoHs0YRaorPjUHUia/TcIYElJzdRMQqdUq3S3qxIkTOH/+PD+lIFDSJ3Tx4sXUx1ZNTk5O8PX1RUxMDM6dO4devXrpOiRCCNEIyhGaJysuYmJiIBaLYWNjo+OIiKHxdrCAgINCgWFnbqybgEitVumWC1NTU2RlZSlsz87OhokJTW2mLlmS/e2335CUlKTjaAghRDMoR2iebFA3AFy9elXH0RBD5GJjhkUDAiHkOLntyw/fB2NMxVGEVE2li4u3334bY8aMwYULF8AYA2MM58+fx0cffYS+fftqI8ZaSfbHvHXrVnh6emLDhg06jogQQqqPcoR2UNcoUl1Dgj1wekYXbAtvi02jgmEs5HDgdgrWnnyAs3GpNP6CaEyli4tVq1ahUaNGaNeuHUQiEUQiETp06AAfHx+sXLlSGzHWOklJSdiyZQt/WyqVYuzYsdSCQQgxeJQjtENWXOzbt49yBakyFxsztGtkj5DG9THlzcYAgEX/3cXQdRfQYfFRbL9E62+R6qv0mAtbW1vs3r0bMTExuHu3ZOYBf39/+Pj4aDy42iomJgZSqVRum0QiQWxsLNzc3HQUFSGEVB/lCO1ITU0FABw7dgyenp5Yu3YtvyghIVXRp6kLFiuZQaqTnyOtgUGqpUrrXACAr68vfH19NRmLXomMjERkZCQkEonGz+3r6wuBQCBXYAiFQkq+hJBaozbnCG3mB2WSkpLkWn1krd2hoaH0hRSpsodpijNFyWaQouKCVIdaxcXkyZPVPuGyZcuqHIw+iYiIQEREBDIzMzU+M4ebmxvWrl2LsWPH8slp9uzZlCQIIQapruUIbeYHZai1m2iDqhmkGtiKdBMQqTXUKi6uXbum1sm4MrMQENXCwsIQGhqKDz74ACdOnOCbvAkhxNBQjtAuau0m2iCbQWrWzluQlJox6qu9dzDnbX8kpefB28GCWjFIpalVXBw7dkzbcdRJbm5umD17Nk6cOIEtW7bgm2++gbm5ua7DIoSQSqEcoV2y1u4xY8bwBcZPP/1ErRak2oYEe6CTnyMSUnPxPDsfU/+8gUPRKTgUnQIAEHDAogGBGBLsoeNIiSGp9GxRRLO6du2Khg0bIjMzE3/88YeuwyGEEKKHwsLCcO/ePX6tkLZt2+o4IlJbyGaQ6hvkikX9A+Xukw3ypmlqSWVQcaFjAoEA4eHhAIA1a9boOBpCCCH6ysfHB926dQMA7N+/X8fRkNrIRcl4C9kgb0LURcWFHhg5ciSMjIxw/vx53LhxQ9fhEEII0VNvvfUWACouiHbIBnmXlSLOp4X2iNqouNADzs7OeOeddwAAy5cvx7Fjx2iRJEIIIQpkxcWpU6eQnZ2t42hIbSMb5C0sM/nCp39E0UJ7RG1UXOiJMWPGAAA2bdqErl27wtPTExs2bNBxVIQQQvSJr68vvL29UVhYiOPHj+s6HFILDQn2wOkZXbAtvC1+C2sjd5+UATN33qQWDFIuKi70xGuvvSZ3W7ZIErVgEEIIkeE4DqGhoQCoaxTRHtkgb07Jp0QpAy4npCNZnEddpYhSVV6hm2hWXFycwjZaJIkQQkhZb731Fn766SccOHBA16GQWk7VQntT/4xCYTEDA01XSxRRy4WekC2SVBotkkQIIaSsrl27wsjICLGxsYiNjdV1OKQWKzsGQ8ABzjYiFLwsLACarpYoouJCT8gWSSpdYCxdupRaLQghhMixsrJCx44dAYBaL4jWlR6DcWZGV3w7sKnCPhLGcCUhnbpJEQB1pLj4999/0bhxY/j6+mL9+vW6DkelsLAwPHjwAI0aNQIAPHjwQMcREUJI7WYo+aEsmpKW1CTZGAwXGzP41LdUOl3t+G3XaEYpAqAOFBfFxcWYPHkyjh49imvXruHbb7/FixcvdB2WSp6enli7di0A4Mcff8Tdu3d1HBEhhNROhpYfSpMVF4cPH8aBAwdo8g9SY1RNVytT0k2KZpSqy2p9cXHx4kU0adIErq6usLS0RM+ePXHw4EFdh1Wurl27om/fvpBIJPjkk09o3QtCCNECQ8wPMk2bNoW1tTXy8/Px1ltv0fTlpEaV7ir1w9DmCvdLGHA4OoVmlKqj9L64OHnyJPr06YMGDRqA4zjs2rVLYZ/IyEh4eXlBJBKhTZs2uHjxIn/fkydP4Orqyt92dXXF48ePayL0avnmm28gEAhw+PBhWveCEEKUqKv5AQAeP36MrKws/jZNX05qmqyrVEtPO6XdpObsvo12i45SV6k6SO+Li5ycHAQFBSEyMlLp/du3b8fkyZMxb948XL16FUFBQQgNDcWzZ89qOFLNsrCwAGOv5n6jxEEIIfLqan4AgJiYGLkcAbyavpyQmqRsRqnXG9jI7SNbfO96khgxYg7J4nxdhEpqiN6vc9GzZ0/07NlT5f3Lli1DeHg4Ro0aBQD46aefsHfvXvz888+YMWMGGjRoIPdN1OPHj9G6dWutx11d5SUOmkGKEELqbn4AXk1fLpVK+W00fTnRlSHBHujk54iE1Fx4OZgjPjUHQ9ddkNtHyoBBay6AQYgf75yktTFqMb0vLspTWFiIK1euYObMmfw2gUCA7t2749y5cwCA1q1b49atW3j8+DFsbGzw33//Yc6cOSrPWVBQgIKCAv52ZmYmAKCoqAhFRUVqxybbtzLHlObl5aWQOAQCATw9Pat8zsqoMP6iIhiX3rcGYlJH6Xgr+zvTF9V97+iaIcdvyLEDVY/fUJ9veWpzfgAAJycnrF69Gh999BEYY+A4Dj/++COcnJy0/vtUK37KEVqjr9cpB3MjOHhYAwCKbEyVLr5Xem2MmTtvop23HQDg4YtceNqbw8VGVIMRV56+vvbqqqkcYdDFRWpqKiQSCZycnOS2Ozk58bMsGRkZ4bvvvkOXLl0glUrx2Wefwd7eXuU5Fy1ahAULFihsP3jwIMzNzSsd46FDhyp9jMzHH3+M1atX8wWGmZkZzpw5Aysrqyqfs7JUxS/Mz8fbL38+cOAAJCL9uCDkS141tR49ehQioX7EVRXVee/oA0OO35BjByoff25urpYi0Z3anh+Akufy2WefYcmSJTA1NYWtrS327dtXrXNWRnnxU47QPn2/Tg325rD9gQAMHErKCvmBGVIGfLz+OKIzODBw4MAwpKEU7ZyY0vPpE31/7Sui7Rxh0MWFuvr27Yu+ffuqte/MmTMxefJk/nZmZibc3d3x5ptvwtraWu3HLCoqwqFDh9CjRw8YGxtXfIASvXr1wpQpU3Dr1i1MmDABCQkJ+Ouvv/DHH3+AUzEFnKZUGH9ODv9jaGgoYGGh1XjUlVOYA9ws+blr166wtbDVaTxVoYn3ji4ZcvyGHDtQ9fhl38DXRYaaH2TeeustbN++HQkJCSgsLET//v2rdT51qBU/5QitMZTrVC8A48T5SEzLhchYiMFrLyi0ZNzOeDX0l4HDH/FCjBvQSW9bMAzltVelpnKEQRcXDg4OEAqFSElJkduekpICZ2fnKp3T1NQUpqamCtuNjY2r9Eaq6nEy3t7e8Pb2RoMGDdCuXTvs3r0by5YtQ5s2beDr66v18Rcq4y+1zdjYWO62Lhkz+bgM8Y9fhuLXHUOOHah8/Ib8XFWpC/lB5oMPPsBXX32Fbdu24YMPPqj2+dRVbvyUI7TOEOL3cDCGh0NJb4tFAwIxc+dNSFnJoO/mHna48jBdbn8pA64+EqO+uADeDhZwsTHTRdgVMoTXvjzazhF6P1tUeUxMTNCyZUscOXKE3yaVSnHkyBG0a9euWueOjIxEQEAAgoODqxumRrRs2RKLFy8GUPLtGU1PSwghqtWl/DBs2DAAJd2PyhZThOiLIcEeOD6lE8YHSHB8Sif8MLS50ilsP91+naavNXB6X1xkZ2cjKioKUVFRAID4+HhERUUhMbHkDTd58mSsW7cOmzdvxp07d/Dxxx8jJyeHnx2kqiIiIhAdHY1Lly5V9ylozMCBA+Vu0/S0hJC6jPJDCT8/P7Ru3RoSiQTbt2/XdTiEqORiI4KvDYOLjUitlb5n7LiJKw/TaDE+A6P33aIuX76MLl268Ldl/V1HjBiBTZs2YciQIXj+/Dnmzp2Lp0+folmzZti/f7/CIL7aIC4uTmEbTU9LCKmrKD+88sEHH+DixYvYsmULJkyYoOtwCFFL6SlsX+QUYPzWa3L3MwDvrj7H3xZwoClsDYDeFxchISEK6z2UNX78eIwfP16jjxsZGYnIyEhIJBKNnrc6lM1rDpRMW0sIIXUN5YdX3nvvPUyePBmXL1/G5s2b0a1bN/rSiRgEFxszuNiYIVmcp3T62tKkDJix8ya87S1QzJhej8uoy/S+W5Su6FuzNwC4ublh7dq1EAqFcttnzZqFo0ePUvcoQgipAfqYHxwdHREQEAAAGDlyJI3JIwanbDcpIcdhzBveCvsxBgxee15hXAZ1ndIfet9yQeSFhYUhNDQUsbGxiIuLQ3h4OLZt24Zt27ZBIBBg7dq1CAsL03WYhBBCalBSUhJu3rzJ35aNyQsNDaUWDGIwyq70DQDrT8erbM2QMmD6jps4G/cC/1x/ws9ERV2ndItaLgyQm5sbQkJCSuYOL4UGeBNCSN0UExOj0EVMNiaPEEPiYmOGdo3s+e5SpVszVH1o3R31hC9ApAyYtfMWksV51JqhI9RyoYI+9qktS1UyiYmJoW+qCCFES/QxPygbkycUCuHj46PDqAipvtKtGeYmAvT/8Wy54zIAQMIYPvvrBs7EplJrhg5Qy4UK+tintixZMilr69atePToEY4dO0atGIQQomH6mB9kY/JK54S5c+fSF02kVpC1ZgS52ymMy5jZ8zWl62WcikmVa82YufMmrj9Kp5aMGkDFhQErO8Cbe/nHtn79enh6etJCe4QQUoeEhYXh4cOHeOONNwAAjx490nFEhGjekGAPnJ7RBdvC2+L0jC4Y27mRfNcpDmjjZadwnJQB70SepYHgNYCKCwMXFhaGhIQEHDt2DImJiVi6dCkA8N2laBwGIYTUHW5ubli4cCEAYMuWLXj+/LmOIyJE80qPywDkC44zM7pixfvKV/+WkS3Qt2T/HXRYfJRWBNcwKi5UiIyMREBAAIKDg3UdSoVkA7zd3NzQokULhftpUB8hhGiOvueHDh06IDg4GAUFBfjpp590HQ4hNaK8geDK6gwGYPXxBwoDwanrVPVRcaGCPvapVYeqcRjx8fE0BoMQQjRA3/MDx3GYNGkSgJJCqKCgQMcREVLzSrdm7IpoX25LhoyEMfQrt+tUvpajrh2ouKhlVC209+GHH9IYDEIIqSMGDhwINzc3pKSkYNu2bboOhxCdqMpAcNlEVLLVwL/df5fvOhXy3UmcS1GjSqnjaCraWqj0QnsA0KVLF/4+WliJEEJqP2NjY3zyySeYPn06lixZAg8PD/j5+dF1n9RZZRfoc7Exg625MWbtvAUJY+DwqrCQYQyIPB7H35YyYPsDAQYmiVEoBbwdLPhxH+QVarmopWTjMMqugwGUjME4fPgwkpKSqKsUIYTUUuHh4TAxMcHdu3fRrVs3arkmdV55A8HV7TrFwGHgmgs061Q5qLhQQd8H7KlL1RiM8PBweHh4UFcpQgipJEPJDzk5OSgqKuJv0+yBhCiqStcpGdmsU5/9dZ1mnSqFigsV9H3AnrrKjsEQCoXw9fVFcXExTVdLCCFVYCj5ISYmRqH1mmYPJES1itbQUDXr1B+Xk2jBvlJozEUdUHoMho+PD+7du4fu3bvL7SORSHD27Fk4OjrC19cXTk5OOoqWEEKIJsharqVSKb9NKBTCx8dHh1ERot9kU9nKlB6rYSxgGLTmPJjSMuMV2YJ9AMBxwML+gXi/tQeSxXmIT82p9WM1qLioI9zc3OQG8pVNOAAwZMgQ/r7Vq1dTgUEIIQZM1nI9duxYSCQSAMC4ceNoUDchlSQrOIqKijCkoRR/xAshZSVdpz57qzGW7L/Lt1yUxV62ZGw59xB3kjPBULKK+KIBgejk51griw3qFlUHqZquVkYqlWLcuHFITU2t4cgIIYRoUlhYGBISEvD+++8DAE6cOKHwxRIhRH3tnBiOT+mksuuUqg/W0S8LC6CkZWP6jptov0hxnEZtGBhOLRd1VOmuUs+ePeNbLWQkEglOnTqFrl27IiEhAb6+vvRtFyGEGCA3Nzf88MMP2Lt3L27cuIE//vgD7733nq7DIsRgudiI4OFgxd8u3XXK3ESA/j+elWvJUDbNLaC4pkbcs2ysPx0PKXvVujEk2EObT0UrqOVCBUOZDaQ6ZNPVtm/fXumMUps3b0bDhg1pRilCCCnFEPNDvXr1MG3aNADAzJkzcejQIZrEgxANKm/WqRkVzDoFlHSfWnsqXmFgeLI4z+BaM6jlQoWIiAhEREQgMzMTNjY2ug5Hq8r2yxUKhWjbti3OnDnD7yObUapp06bIzs6Gr6srqB2DEFIXGWp+mDhxIpYsWYKEhAS8+eabEAgEWLt2LcLCwnQdGiG1SkUL9glQ0mqhYpgGT8qAnitOQpxXbFBjNai4IAAUZ5S6c+cO3nzzTbl9JBIJ2rRpA8ZYSVICQCmJEEIMg1gsRk5ODn9b9qVRaGgodXslRMPKm3XKy8EcJ+8/54uN8gaGZ+QV8z/LxmrIlO46pU8zUVFxQXilZ5QqKioCx3EKc6TLrY0BIBQAHj9GzOPHNC6DEEL0WHnrXtC1mxDtK11wVNS6IeQ4DG3jgS3nH6o8n6zY2HntMS7Gp4Ex/WjdoOKCKOXm5oZx48bhp59+gkQiUTp1rQTAMAAn/f0hlUqpiZ0QQvSYsnUvBAIBrXtBiI5U1LoBAL9deKhymluZCw/S+J9lBYdsELl860Y+YsQcksX58HAw1sIzKkEDuolKPXr0QExMDI4dO4bz588rHfR9HOATFa30TQgh+kvZNORNmjSBq6urDqMihJQmGxguKzzKTnNbwbhwntxMVDtuYsK2qwj57iR+iBYi5LuT/NS32kDFBSmXbEap4OBguaQkFArxppL9JRIJzp8/j6SkJBw7dowKDUII0SOydS82b94MIyMj3Lx5E3v27NF1WIQQFYYEe+D0jC7YFt4WZ2Z2xeJ35WeimqnOTFQA9lxPlpuJatbOW1qbfYq6RRG1yQ36btAAaNwYngDKLsc0ZMgQMMZeDfxeuxahoaGIiYmhcRmEEKJjbm5uGD58OO7cuYPFixdj4sSJ6NGjB8zNzXUdGiFEicqM1VB3JioJY0hIzdXKeAxquSCVImvJcHs5Fe1agG/NEAgE8PLyglQqlRv4HR4eDk9PT1ovgxBC9Mjs2bPh4eGBhw8fYsaMGdTaTIiBKN11Cqha64aQ4/hxHZpGLRcqREZGIjIyEhKJRNeh6LUwAKHR0Yh98gQ+Pj6IiYlB165d5faRtWIA8lMfAqDWDEKIwakt+cHCwgIrVqzAgAED8P333+P777+niTkIMVDqtG7M3HmTX/174YDXtTaLFLVcqBAREYHo6GhcunRJ16HoPTdX15LWDDc3fjaS8kgkEkyePJlaMwghBqk25YdWrVrJ3aaJOQipHZS1bhyf0gnjAyQ4PqUThgR7aO2xqbggGlV2NhKBQACOUxxp9OeffyrMMnXp0iVqlieEkBoUGxursE229gUhpHZxsRHB14bBxUak1ceh4oJonGw2kmPHjuHhw4dYt26dXLHRrFkzhWNkq39TSwYhhNQcZa3NtPYFIaQ6qLggWsEP/HZzUyg2/vnnH6Vdp0qPyxgzZgwSEhJoSltCCNEiZWtfeHp6wtnZWYdREUIMGRUXpEaULjaUdZ0qSyqVokmTJvDw8JBrzaBigxBCNEv2BdD27dthZWWF+Ph4LFmyRNdhEUIMFM0WRXSi9JoZFhYWaNu2LT8GQyY3N5f/WSqVYvTo0eA4Tm79DJrRhBBCqs/NzQ2DBw/G/9u787io6v1/4K9hZHFjc0EQRAT3QC2VtNvVEEPtmnZLy1uJBmqCjzCzslu55C3T0p/ll1BaJLPFyqX1aqhg5UoqaqkICBqJaymKqTjz/v1hc+4MMwMzOMzCvJ6Ph4/HnDOfM/Oe45nz5j2f8/mcq1evYuzYsZg1axZatWqFYcOGcTY/IrIKey7IYWq6+/e0adNMblP90qljx46xN4OIyEYeeeQR9OrVCxqNBpMmTeIYOCKyGnsuyCkY3P37r4GEixcvNurN0KfVatGlSxdcvXpV6c1Y8tYS4AKA34HffvsN/p387fMBiIgagN9++w379u1TlvXvTcQeDCKyBHsuyGnUNi7D1JS2V65cMejNSH08Ffh/AN4HYrrGKL+4sXeDiKh2hYWFRj/qcGpaIrIGiwtyWjVNaVvTpVM6Wq0WEyZMwJQpU3jDPiIiC5i7EWpNvchERPpYXJiRnp6Obt26oU+fPo4Oxa2Zm9K2tLQUTz75ZK13AxcRpKen84Z9RGQzDTk/mJqaFrhxV/KCggKeN4moViwuzEhNTcXBgweRl5fn6FBIT10unapOo9Ggb9++iIuLQ7t27ZCeng6Al04RkWUaen7Q/yFn7969CAkJweHDh9GlSxfExcUhKioK2dnZjg6TiJwUiwtyadUvnfq/jP8D/qov1Go1FixYUGPvhohgypQpiIiIMLqnBsCCg4jck+6HnJ49e2LJkiUGz2m1WmRkZPC8SEQmcbYocnm6XgwASByfiNSCVOB3YN+/96F7p+4IDAzEpEmToNFo4OHhYfLa4dLSUuWx7p4an332GbKzs6HVapX7aiQkJKCwsBAdO3bkzClE5BYCAgKM1mm1WhQXFyMiIsIBERGRM2PPBTU8fgAigLZt2wIw7N3YsWOHUU+GuUupNmzYYDBWIzk5mb0bROR2TA3y9vDwQGRkpIMiIiJnxuKC3EJNN+ybP39+rQPDdfSnvZ0wYQLGjRtnNBMViw0iakhMDfJu1qwZfHx8HBgVETkrFhfkdqrPOvX0008bFRy1jdUAbhQa77//vkHvxoQJE0xOe8uCg4hcme68uWbNGoSGhqKiogIjRozAhg0beF4jIgMcc0FuSX+cBmB8h/DQ0FCjsRoiovRcmKPfRldsbN26VSlCOHaDiFyV7rwZERGB2NhY5OXlYciQIcp5LSkpydEhEpETYM8F0V/0p7kFar+Jn6W9G8uXL+fYDSJqMPz8/FBVVaUs6+4fxHMXEQHsuSCqkX4Ph616NwDjsRtfffUVvvrqK6V3Y/LkyYiJiUFpaSl7N4jIqRQVFRmd5zQaDYqKiniuIiL2XBBZo756N7744guD3o309HRERkayd4OInE5UVJTJWfb27dvngGiIyNmw54LoJtmjdyM5ORlvvvkmDhw4ABGBh4cHli1bhiFDhhiN3SgrK+N4DiKqN6GhoUhJScHSpUuh0WigUqkgInjqqafQvHlzREREoGPbtuDZh8g9seeCyMbqo3cDAPbv3290OVVYWBji4uLQrl07zJ49G0uXLuXUuERU7wYPHozCwkJl1r3ExERoNBokJSXdOP907Yp3HR0kETmEWxQX9913HwICAvDAAw84OhRyU/oFR21T4Xp4eBhdcmDuRn86IoI5c+Zg8uTJnBqXyArMD3WnO6/pftzQp9VqMQkAzzBE7sctiou0tDSsWLHC0WEQKWrq3SgqKkJKSkqdb/SnT0SMio17773XooKDBQi5A+YH2ygpKTFapwFQZP9QiMjB3KK4GDhwIJo3b+7oMIhqpF9wVL/kwJLeDQ8PD4sGj+tmpQL+N54jNjbWYHrcxMRE9niQW2B+sI2OHTuaPP8cBpCzZQvPGURuxOHFxffff4/hw4cjJCQEKpUK69atM2qTnp6O9u3bw8fHB7Gxsdi1a5f9AyWyM2vHbmRmZtZagJiza9cug/EcK1asMChAJk6ciOnTp1s0nqOsrAwHDhzgHxN005gfXEdoaKjB+UdnMoC4e+4x+JGCiBo2h88WVVlZiR49euCxxx7DP//5T6PnV61ahWnTpmHp0qWIjY3F4sWLkZCQgIKCArRu3RoA0LNnT1y/ft1o2++++w4hISH1/hmI7KW2makAGKzbsGGDMlOVWq3GvHnzMGPGDKVwAKDM9FITrVaLhQsXGiwnJycr2+pmr1KpVJg4cSK0Wi1mzZrFu5HTTWF+cC365yR1VRX+fvfdynO6G+0lJCTwPEDUwDm8uBg6dCiGDh1q9vlFixZhwoQJGD9+PABg6dKl+Oabb/Dee+9hxowZAID8/HybxXP16lVcvXpVWa6oqAAAVFVVGdyRtDa6ttZs40xqjb+qCp76bZ3kc+rHa+3/mbOw5tgJCgpCUFCQQXv9dWPHjkVcXByKi4sRGRmJ0NBQ+Pn5ISUlRSk4Xn75Zfz73/82KDgsVX32Kn26dSqVSrk5YEZGBsaPH4+ysjKDoqj6sqM0+O9tLds5G+YH51RT/LrzT+769UbPaTQaHD58WDk/OYK75Qhn48qxA+4bv7XtHV5c1OTatWvYvXs3nnvuOWWdh4cH4uPjsX379np5z3nz5mHOnDlG67/77js0adLE6tfLzs62RVgOYy5+9ZUr+Mdfjzds2ACNj4/9gqrBFc0V5fHmzZvho3aOuOrC1sfO/v37sX//fgQFBWHZsmUoLy9HcHAwWrZsicmTJyMjI0MpAgYMGIAtW7Yoy48++ihWrFhh0f059Onf00P3y2VGRgb27dsHEYFKpcLAgQORm5urLKekpGDw4ME4e/asQYz21FC/t+Zcvny5niKpP8wPjldT/H/89hs8AFT/ySIrKwvHjh1z2HebOcI5uHLsgPvFb22OUIm1fy3UI5VKhbVr12LkyJEAgBMnTqBt27bYtm0b+vXrp7R75plnsGXLFuzcudOi142Pj8e+fftQWVmJwMBAfPbZZwavp8/UL1NhYWE4e/YsfH19Lf4sVVVVyM7OxuDBg+Hp6Vn7Bk6m1vgrK+EZEHCj7R9/AE2b2jlC0yqvVSLg9RtxnU47Df+m/o4NqA4cdeyUlZUZ9HBUX16+fLnS42HqZoC68R03e0rx8PBAcnIy3nnnHYMej8GDBxv1bti6x6PBf2/NqKioQMuWLXHhwgWrznP2xPzgPCyKv7ISKwICMAk3Zo3Sv/xS/3JKXW+mvTBHOJYrxw64b/zW5gin7rmwlY0bN1rc1tvbG97e3kbrPT0963Qg1XU7Z2E2fr11np6eBsuO5CmGcTXIfV9PIiIiEBERYXZ54sSJGDZsmNnxHMuWLQOAOt2NXJ9Wq0VmZqbB8qRJk5Q/SFQqFWbOnIlWrVrhiSeeUAqQzMxMJCUlGd2hvC53LHe3Y8eVP+vNYn6ouxrj9/REEoAEAEXffovIW27BokWLsHjxYoPezJSUFAwbNsxul0MyRzgHV44dcL/4rf2sDp8tqiYtW7aEWq3GqVOnDNafOnUKbdq0qdf3Tk9PR7du3dCnT596fR8iV1LTzQCTkpKQlJSEwsJCzJ07F0VFRXW+G7kpuj9IdDcMnDJlitE9PO6//36DKXXHjRvHe3o0UMwPriEUwMC//x1hYWG49957jZ7XaDQoKuLdMIgaEqcuLry8vHDbbbdh06ZNyjqtVotNmzaZ7ba2ldTUVBw8eBB5eXn1+j5Erqz6dLm6ddHR0RbdjVytViMxMdEmBYiIYM2aNQa/ir7//vtGU+omJycbFBz6BUhUVJRyLSoLDufG/OB6zN0Lo2nTpvy+ETUgDr8s6tKlSwa/WpSUlCA/Px+BgYFo164dpk2bhsTERPTu3Rt9+/bF4sWLUVlZaddrNImo7vSnzwVMT6H7n//8x2A5MDCwxsurdH+gWDvDlVarNZhrX1eA6C+/9dZbAGAwwN3clLp1ueSKLMf80LDo7oWh+27rDB48GBcvXjS6xJGIXJPDi4uffvoJd911l7I8bdo0AEBiYiKysrLw4IMP4syZM5g5cyZOnjyJnj17Yv369Q6dyo6Ibk71gqO2AqS28R3m7uFRFyKC9PR0ZVl3Tw8dlUqFJ554AoGBgZgzZw4LkHrE/NDw6H+3vby8kJiYaFBA8n4YRK7P4cXFwIEDax3wOWXKFEyZMsVOEd2Qnp6O9PR0g19XiMh+rL1hYPUeD7VajUceeQQrV660aQEiInjjjTcM1pkqQNLS0hAQEFDnAoRFCfNDQ6X/3V60aJHRWAzdOAx3Pe6JXJ3DiwtnlZqaitTUVFRUVMDPz8/R4RC5veq9G6bWWXvJlakC5JFHHsEHH3xw0wXI4sWLDdaZKkBSUlIQEBCAV155xaAAAaDc6Zy9Is6H+cF2evXqBQ8PD4Pvm+7Gmzk5OTy2iVwQiwsialCsveRKvwAJDw/H/v370b9//xrv6VHXMR/6ql9+pXs9/QJEt676nc7NFSBxcXE4cOAAYmJiDKYRJnJWpsZhiAgGDRoEAByDQeSCnHq2KEfiVINEDVf1Wa6qL48fP16Z5erYsWNGU+pmZmYazHrl4eGh3ERQx8PDo87T7lYnIgazXiUnJyM5OdloKt7IyEi8+OKLiIqKMhi4TrbF/GBb+rPK6U+wAPxvDAZnkSJyHSwuzOBUg0TuzZJ7etRHAaJSqYzaWUK/d4V/kNUv5gfb033fwsLCjJ7TaDQ4ePAgp6slchG8LIqIyAK1jfmwZNC5JbNeWXKnc13xUdNgZw6KJVekuxdG9UsOx40bh1OnTnG6WiIXwOKCiMhGbFGA1KUoMVWAqNVqREVF2eNjE9lM9TEYHh4e8PHxQXl5udKG09USOTcWF0REdmTJrFfV19W1V4R/eJErqn68b9++HaNHjzZow545IufF4sIMzmNORM7EkgIkLi4OH374IR5++GHOFlWPmB/qn/6x3a9fP5OXSh05cgQiwulqiZwMB3SbwQF7RORqQkNDER0dzT+06hnzg33pLpXSTY6gM2nSJMTFxSE8PJyzoxE5ERYXRERE5NT0Z2fLyckxmFFNq9Vi4sSJnEWKyEmwuCAiIiKnp5uutvrkBcCNAuP555/H0aNHOV0tkYOxuCAiIiKXoZuutroVK1YgMjKSl0oRORiLCzN4B1YiIjKF+cGxqo/BUKvV+Ne//mXQRjddbV5eHnsyiOyMxYUZHLBHRESmMD84nv4YjNLSUiQnJxu10Wg0iI2NZU8GkZ1xKloiIiJyOdWnZzY1Xa1ubIauJ+POu+60a4xE7og9F0REROTSql8qZWpMhkajQWFhIXABQAnw22+/2TlKIvfA4oKIiIhcnv6lUjt27DBZYExMmggsBvA+ENM1hpdKEdUDFhdERETUIOimq+3Tp49BT4ZKpYKPjw9O/HYC+GsWW92lUhzsTWRbLC7M4GwgRERkCvODa9DvyTh+/DhWrlxp1Eaj0eD777/njFJENsQB3WakpqYiNTUVFRUV8PPzc3Q4RETkJJgfXIclg74ffvhh5bnMzEwkJSXZNUaihoY9F0RERNTghYaGYslbSwDVXytUhs9rtVpMnDgRZWVlKCsrY28GUR2xuCAiIiK3kDg+EZgKIBF4L+s9o+e1Wi0SEhIQHh7O+2MQ1RGLCyIiInIffgAigL639zU5o9TBgweVS6c46JvIeiwuiIiIyO20bdvWYEYptVqNRx55xKidRqNBTk4OL5MishAHdBMREZFbSkpKQkJCAoqKihAVFQUA+Oijj4wGfY8dOxYAB30TWYI9F0REROS2dPfG0M0sVf3+GPp0g76Li4s56JvIDBYXZnAecyIiMoX5oWHTvz/GJ598YvS8VqtFr1690K5dOw76JjKBxYUZqampOHjwIPLy8hwdChERORHmh4ZP15vRv39/k4O+L168CJEbt/rWH/TN3gwiFhdEREREJlW/TEqtViMtLc2onUajwdixYzmFLRE4oJuIiIjILFODvpcsWWI06DsnJ0d5rOvNiImJwaVLl9CxY0eDO4UTNWTsuSAiIiKqQU2DvtVqNQYPHmy0jUajQWxsLHsyyO2w54KIiIjICqZ6M8LDw416M6qPy0hISAAAFBYWsjeDGiz2XBARERFZyZopbIEbPRlPPfUUx2VQg8eeCyIiIqKbpN+b0bRpU9x+++1GPRmffvqp8pjjMqihYs8FERERkQ3oejP69Olj0JPh4eGBHj16GLXnuAxqiNhzQURERGRjHJdB7oo9F0RERET1oC7jMiZNmsRxGeTS2HNRC90vChUVFVZtV1VVhcuXL6OiogKenp71EVq9qjX+ysr/Pa6oADQa+wVXg8prlcCVG48rKirgoXG9+rnBHztOzJVjB+oev+78pjvfkWWYH2qInznCpFGjRqF///44evQomjRpgkGDBhl977799lvlsVarxcSJE9GhQwdUVlYiMjISrVu3dtnjxy2OfSdmrxyhEmYTk9LT05Geno5r166huLjY0eEQEdW7X3/9lZdgWID5gYjckaU5gsVFLbRaLU6cOIHmzZub7MI0p6KiAmFhYfj111/h6+tbjxHWD1eO35VjBxi/I7ly7EDd4xcRXLx4ESEhIfDwcL3ePkdhfmD8juDK8bty7ID7xm9tjuBlUbXw8PC4qV/yfH19XfIA1HHl+F05doDxO5Irxw7ULX4/P796iqbhYn5g/I7kyvG7cuyAe8ZvTY7gT1RERERERGQTLC6IiIiIiMgmWFzUE29vb8yaNQve3t6ODqVOXDl+V44dYPyO5MqxA64fv7tw9f8nxu9Yrhy/K8cOMH5LcUA3ERERERHZBHsuiIiIiIjIJlhcEBERERGRTbC4ICIiIiIim2BxcRPS09PRvn17+Pj4IDY2Frt27aqx/WeffYYuXbrAx8cH0dHR+Pbbb+0UqaF58+ahT58+aN68OVq3bo2RI0eioKCgxm2ysrKgUqkM/vn4+Ngp4v+ZPXu2URxdunSpcRtn2e8A0L59e6P4VSoVUlNTTbZ39H7//vvvMXz4cISEhEClUmHdunUGz4sIZs6cieDgYDRu3Bjx8fEoLCys9XWt/e7UR/xVVVV49tlnER0djaZNmyIkJARjx47FiRMnanzNuhyDto4dAMaNG2cUx5AhQ2p9XXvte3fH/GD//AAwR9hz3zM/GLNXfqgtfsCxOYLFRR2tWrUK06ZNw6xZs7Bnzx706NEDCQkJOH36tMn227Ztw5gxY5CUlIS9e/di5MiRGDlyJH7++Wc7Rw5s2bIFqamp2LFjB7Kzs1FVVYW7774blZWVNW7n6+uL8vJy5d+xY8fsFLGh7t27G8Tx448/mm3rTPsdAPLy8gxiz87OBgCMGjXK7DaO3O+VlZXo0aMH0tPTTT6/YMECvPnmm1i6dCl27tyJpk2bIiEhAVeuXDH7mtZ+d+or/suXL2PPnj148cUXsWfPHqxZswYFBQW49957a31da47B+ohdZ8iQIQZxfPzxxzW+pj33vTtjfnBcfgCYI+y175kfTLNHfgCcPEcI1Unfvn0lNTVVWdZoNBISEiLz5s0z2X706NFyzz33GKyLjY2VSZMm1Wucljh9+rQAkC1btphts3z5cvHz87NfUGbMmjVLevToYXF7Z97vIiJpaWkSGRkpWq3W5PPOst9FRADI2rVrlWWtVitt2rSR1157TVl3/vx58fb2lo8//tjs61j73bGV6vGbsmvXLgEgx44dM9vG2mPQFkzFnpiYKCNGjLDqdRy1790N84PjMEc4BvPDDY7IDyLOlyPYc1EH165dw+7duxEfH6+s8/DwQHx8PLZv325ym+3btxu0B4CEhASz7e3pwoULAIDAwMAa2126dAnh4eEICwvDiBEj8Msvv9gjPCOFhYUICQlBhw4d8PDDD+P48eNm2zrzfr927RpWrlyJxx57DCqVymw7Z9nv1ZWUlODkyZMG+9fPzw+xsbFm929dvjv2dOHCBahUKvj7+9fYzppjsD7l5uaidevW6Ny5MyZPnoxz586Zbevs+76hYH5w/HmKOcLxOYL5wfH5AXBcjmBxUQdnz56FRqNBUFCQwfqgoCCcPHnS5DYnT560qr29aLVaTJ06FXfccQduueUWs+06d+6M9957D1988QVWrlwJrVaL/v37o6yszI7RArGxscjKysL69euRkZGBkpIS3Hnnnbh48aLJ9s663wFg3bp1OH/+PMaNG2e2jbPsd1N0+9Ca/VuX7469XLlyBc8++yzGjBkDX19fs+2sPQbry5AhQ7BixQps2rQJ8+fPx5YtWzB06FBoNBqT7Z153zckzA+OPU8xRzhHjmB+cGx+ABybIxrVOWpqEFJTU/Hzzz/Xek1gv3790K9fP2W5f//+6Nq1K5YtW4a5c+fWd5iKoUOHKo9jYmIQGxuL8PBwfPrpp0hKSrJbHLbw7rvvYujQoQgJCTHbxln2e0NXVVWF0aNHQ0SQkZFRY1tnOQYfeugh5XF0dDRiYmIQGRmJ3NxcDBo0yG5xUMPlavkBcJ7vpy0wRzgHV8wPgGNzBHsu6qBly5ZQq9U4deqUwfpTp06hTZs2Jrdp06aNVe3tYcqUKfj666+Rk5OD0NBQq7b19PREr169UFRUVE/RWcbf3x+dOnUyG4cz7ncAOHbsGDZu3Ijk5GSrtnOW/Q5A2YfW7N+6fHfqmy5xHDt2DNnZ2TX+KmVKbcegvXTo0AEtW7Y0G4cz7vuGiPnBuc5TzBGOwfxwg7PkB8C+OYLFRR14eXnhtttuw6ZNm5R1Wq0WmzZtMvgFQV+/fv0M2gNAdna22fb1SUQwZcoUrF27Fps3b0ZERITVr6HRaHDgwAEEBwfXQ4SWu3TpEoqLi83G4Uz7Xd/y5cvRunVr3HPPPVZt5yz7HQAiIiLQpk0bg/1bUVGBnTt3mt2/dfnu1Cdd4igsLMTGjRvRokULq1+jtmPQXsrKynDu3DmzcTjbvm+omB+c6zzFHOEYzA83OEt+AOycI6wa/k2KTz75RLy9vSUrK0sOHjwoEydOFH9/fzl58qSIiDz66KMyY8YMpf3WrVulUaNG8vrrr8uhQ4dk1qxZ4unpKQcOHLB77JMnTxY/Pz/Jzc2V8vJy5d/ly5eVNtXjnzNnjmzYsEGKi4tl9+7d8tBDD4mPj4/88ssvdo39qaeektzcXCkpKZGtW7dKfHy8tGzZUk6fPm0ybmfa7zoajUbatWsnzz77rNFzzrbfL168KHv37pW9e/cKAFm0aJHs3btXmS3j1VdfFX9/f/niiy9k//79MmLECImIiJA///xTeY24uDhZsmSJslzbd8de8V+7dk3uvfdeCQ0Nlfz8fIPvwtWrV83GX9sxaI/YL168KNOnT5ft27dLSUmJbNy4UW699Vbp2LGjXLlyxWzs9tz37oz5wTH5QYQ5wp77nvnBcfmhtvgdnSNYXNyEJUuWSLt27cTLy0v69u0rO3bsUJ4bMGCAJCYmGrT/9NNPpVOnTuLl5SXdu3eXb775xs4R3wDA5L/ly5crbarHP3XqVOWzBgUFybBhw2TPnj12j/3BBx+U4OBg8fLykrZt28qDDz4oRUVFZuMWcZ79rrNhwwYBIAUFBUbPOdt+z8nJMXms6GLUarXy4osvSlBQkHh7e8ugQYOMPld4eLjMmjXLYF1N3x17xV9SUmL2u5CTk2M2/tqOQXvEfvnyZbn77rulVatW4unpKeHh4TJhwgSjBODIfe/umB/snx9EmCPsue+ZHxyXH2qL39E5QiUiYl1fBxERERERkTGOuSAiIiIiIptgcUFERERERDbB4oKIiIiIiGyCxQUREREREdkEiwsiIiIiIrIJFhdERERERGQTLC6IiIiIiMgmWFwQEREREZFNsLggohqNGzcOI0eOVJYHDhyIqVOnOiweIiJyDswPZAqLCyInkpubC5VKhfPnz9v0dWfPno2ePXva5LXWrFmDuXPnWtSWiYaIyDaYH8hVNHJ0AETW0Gg0UKlU8PBgXewogYGBjg6BiMgI84PjMT8QwJ4Lukmff/45oqOj0bhxY7Ro0QLx8fGorKwEAGi1Wrz00ksIDQ2Ft7c3evbsifXr1yvbmvoVJj8/HyqVCqWlpQCArKws+Pv748svv0S3bt3g7e2N48eP4+rVq3j22WcRFhYGb29vREVF4d1331Ve5+eff8bQoUPRrFkzBAUF4dFHH8XZs2dv+vOWlZVhzJgxCAwMRNOmTdG7d2/s3LlTeT4jIwORkZHw8vJC586d8cEHHxhsr1Kp8M477+C+++5DkyZN0LFjR3z55ZcAgNLSUtx1110AgICAAKhUKowbNw4AsH79evztb3+Dv78/WrRogX/84x8oLi62KLasrCzMmTMH+/btg0qlgkqlQlZWlsnPp9FoMG3aNOV9nnnmGYiIQZvqvza99dZb6NixI3x8fBAUFIQHHngAwI3u8i1btuCNN95Q3re0tBQajQZJSUmIiIhA48aN0blzZ7zxxhsG76Hran/99dcRHByMFi1aIDU1FVVVVUobRx0DRGQZ5gfmB+YHNyVEdXTixAlp1KiRLFq0SEpKSmT//v2Snp4uFy9eFBGRRYsWia+vr3z88cdy+PBheeaZZ8TT01OOHDkiIiI5OTkCQP744w/lNffu3SsApKSkREREli9fLp6entK/f3/ZunWrHD58WCorK2X06NESFhYma9askeLiYtm4caN88sknIiLyxx9/SKtWreS5556TQ4cOyZ49e2Tw4MFy11133dTnvXjxonTo0EHuvPNO+eGHH6SwsFBWrVol27ZtExGRNWvWiKenp6Snp0tBQYEsXLhQ1Gq1bN68WXkNABIaGiofffSRFBYWyhNPPCHNmjWTc+fOyfXr12X16tUCQAoKCqS8vFzOnz8vIiKff/65rF69WgoLC2Xv3r0yfPhwiY6OFo1GU2tsly9flqeeekq6d+8u5eXlUl5eLpcvXzb5GefPny8BAQGyevVqOXjwoCQlJUnz5s1lxIgRSpsBAwZIWlqaiIjk5eWJWq2Wjz76SEpLS2XPnj3yxhtviIjI+fPnpV+/fjJhwgTlfa9fvy7Xrl2TmTNnSl5enhw9elRWrlwpTZo0kVWrVinvkZiYKL6+vvL444/LoUOH5KuvvpImTZpIZmam0sYRxwARWYb5gfmB+cF9sbigOtu9e7cAkNLSUpPPh4SEyMsvv2ywrk+fPpKSkiIilicPAJKfn6+0KSgoEACSnZ1t8n3nzp0rd999t8G6X3/9VTkp19WyZcukefPmcu7cOZPP9+/fXyZMmGCwbtSoUTJs2DBlGYC88MILyvKlS5cEgPz3v/8VEdP7xJQzZ84IADlw4IBFsc2aNUt69OhR20eU4OBgWbBggbJcVVUloaGhZpPH6tWrxdfXVyoqKky+nn7bmqSmpsr999+vLCcmJkp4eLhcv35dWTdq1Ch58MEHRcRxxwARWYb5wRDzgzHmh4aLl0VRnfXo0QODBg1CdHQ0Ro0ahbfffht//PEHAKCiogInTpzAHXfcYbDNHXfcgUOHDln1Pl5eXoiJiVGW8/PzoVarMWDAAJPt9+3bh5ycHDRr1kz516VLFwAw6irW0W/7+OOPm2yTn5+PXr16mb2m9NChQxZ9Xv3P0rRpU/j6+uL06dMmX1OnsLAQY8aMQYcOHeDr64v27dsDAI4fP25RbJa4cOECysvLERsbq6xr1KgRevfubXabwYMHIzw8HB06dMCjjz6KDz/8EJcvX671vdLT03HbbbehVatWaNasGTIzM5XPotO9e3eo1WplOTg4WNlP9XEMEJHtMD8YYn5gfnAnHNBNdaZWq5GdnY1t27bhu+++w5IlS/D8889j586daNGiRa3b6wbdid41m/rXTOo0btwYKpXKYLkmly5dwvDhwzF//nyj54KDg01uk5+frzz29fU12aa297WUp6enwbJKpYJWq61xm+HDhyM8PBxvv/02QkJCoNVqccstt+DatWs2jc1azZs3x549e5Cbm4vvvvsOM2fOxOzZs5GXlwd/f3+T23zyySeYPn06Fi5ciH79+qF58+Z47bXXDK5NBmreT/VxDBCR7TA/1A3zA/NDQ8CeC7opKpUKd9xxB+bMmYO9e/fCy8sLa9euha+vL0JCQrB161aD9lu3bkW3bt0AAK1atQIAlJeXK8/rn8TNiY6OhlarxZYtW0w+f+utt+KXX35B+/btERUVZfCvadOmJrfRb9O6dWuTbWJiYpCfn4/ff//d5PNdu3at8fNawsvLC8CNgXM6586dQ0FBAV544QUMGjQIXbt2VX4BtDQ2Ly8vg9c0xc/PD8HBwQYn8evXr2P37t01bteoUSPEx8djwYIF2L9/P0pLS7F582az77t161b0798fKSkp6NWrF6Kioqz+tag+jgEisi3mh/9hfmB+cCcsLqjOdu7ciVdeeQU//fQTjh8/jjVr1uDMmTPo2rUrAODpp5/G/PnzsWrVKhQUFGDGjBnIz89HWloagBsn7LCwMMyePRuFhYX45ptvsHDhwlrft3379khMTMRjjz2GdevWoaSkBLm5ufj0008BAKmpqfj9998xZswY5OXlobi4GBs2bMD48eNrPYHWZMyYMWjTpg1GjhyJrVu34ujRo1i9ejW2b9+ufN6srCxkZGSgsLAQixYtwpo1azB9+nSL3yM8PBwqlQpff/01zpw5g0uXLiEgIAAtWrRAZmYmioqKsHnzZkybNs2q2Nq3b4+SkhLk5+fj7NmzuHr1qsn3T0tLw6uvvop169bh8OHDSElJqXFO9a+//hpvvvkm8vPzcezYMaxYsQJarRadO3dW3nfnzp0oLS3F2bNnodVq0bFjR/z000/YsGEDjhw5ghdffBF5eXkW7yPd6zriGCAiyzA/MD8wP7gxRw/6INd18OBBSUhIkFatWom3t7d06tRJlixZojyv0Whk9uzZ0rZtW/H09JQePXooA9N0fvzxR4mOjhYfHx+588475bPPPjMasOfn52f03n/++ac8+eSTEhwcLF5eXhIVFSXvvfee8vyRI0fkvvvuE39/f2ncuLF06dJFpk6dKlqt9qY+c2lpqdx///3i6+srTZo0kd69e8vOnTuV59966y3p0KGDeHp6SqdOnWTFihUG2wOQtWvXGqzz8/OT5cuXK8svvfSStGnTRlQqlSQmJoqISHZ2tnTt2lW8vb0lJiZGcnNzjV6rptiuXLki999/v/j7+wsAg/fTV1VVJWlpaeLr6yv+/v4ybdo0GTt2rNkBez/88IMMGDBAAgICpHHjxhITE2Mwq0dBQYHcfvvt0rhxY+X/9cqVKzJu3Djx8/MTf39/mTx5ssyYMcNgQGFiYqLBe4qIpKWlyYABA5RlRx0DRFQ75gfmB+YH96USqTZJMRERERERUR3wsigiIiIiIrIJFhdERERERGQTLC6IiIiIiMgmWFwQEREREZFNsLggIiIiIiKbYHFBREREREQ2weKCiIiIiIhsgsUFERERERHZBIsLIiIiIiKyCRYXRERERERkEywuiIiIiIjIJlhcEBERERGRTfx/6WaAROpbMOEAAAAASUVORK5CYII=",
"text/plain": [
"<Figure size 800x400 with 2 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"dij = np.r_[:15:1j*ng]\n",
"figure(figsize=(8, 4))\n",
"for i, att_i in enumerate([1, 2]):\n",
" subplot(1, 2, i + 1)\n",
" gij = G(dij[:, None, None], c, s, nug=1e-3, att=att_i)\n",
" axvline(Ds, color='r')\n",
" axvline(Dc, color='g')\n",
" semilogy(dij, gij, '.-')\n",
" semilogy(dij, 1/(1e-3+(dij)**att_i), '.-k')\n",
" ylim([0.02, 1000])\n",
" legend(('source diameter', 'contact diameter', 'extended geometry', 'point geometry'))\n",
" grid(1)\n",
" title(['hyperbolic', 'parabolic'][i] + ' attentuation')\n",
" xlabel('source - contact distance')\n",
" ylabel('log G_ij')\n",
"tight_layout()"
]
},
{
"cell_type": "markdown",
"id": "2ee6aabb-48cf-4b8a-bf92-8e0f4f4a3b92",
"metadata": {},
"source": [
"the extended geomtry produces much less sparse gain matrix and may strongly affect intepretation of spatial specificity."
]
}
],
"metadata": {
"kernelspec": {
"display_name": "jax",
"language": "python",
"name": "jax"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.11.5"
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment