{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# [MPR model from tvbk](https://github.com/Ziaeemehr/vbi_paper/blob/main/docs/examples/mpr_tvbk.ipynb)\n", "\n", "\n", "\"Open" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "- Check if tvbk is available, else install it" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Running cmake --build & --install in /home/ziaee/git/inference/tvbk/build/cp310-cp310-linux_x86_64\n" ] } ], "source": [ "import subprocess\n", "import sys\n", "\n", "try:\n", " import tvbk\n", "except ImportError:\n", " print(\"tvbk not found. Installing...\")\n", " try:\n", " # Clone the repository\n", " subprocess.check_call(\n", " [\"git\", \"clone\", \"--branch\", \"master\", \"https://github.com/maedoc/tvbk.git\"]\n", " )\n", "\n", " # Change directory to the cloned repository\n", " subprocess.check_call(\n", " [\n", " \"bash\",\n", " \"-c\",\n", " \"cd tvbk && \"\n", " \"pip install nanobind 'scikit-build-core[pyproject]' pytest pytest-benchmark numpy cibuildwheel scipy && \"\n", " \"pip install --no-build-isolation -Ceditable.rebuild=true -ve .\",\n", " ]\n", " )\n", " print(\"tvbk installed successfully.\")\n", " except subprocess.CalledProcessError as e:\n", " print(f\"An error occurred during installation: {e}\")\n", " sys.exit(1)" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "import warnings\n", "import numpy as np \n", "import networkx as nx\n", "import matplotlib.pyplot as plt\n", "from vbi.models.tvbk import MPR\n", "from vbi import report_cfg\n", "\n", "warnings.filterwarnings(\"ignore\")" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "seed = 2 \n", "np.random.seed(seed)\n", "\n", "LABESSIZE = 14\n", "plt.rcParams['axes.labelsize'] = LABESSIZE\n", "plt.rcParams['xtick.labelsize'] = LABESSIZE\n", "plt.rcParams['ytick.labelsize'] = LABESSIZE" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "nn = 6\n", "weights = nx.to_numpy_array(nx.complete_graph(nn))\n", "params = {\n", " \"G\": 0.0,\n", " \"weights\": weights,\n", " \"delays\": np.ones_like(weights),\n", " \"num_time\": 1000,\n", " \"num_batch\": 1,\n", " \"noise_amp\": 0.0,\n", " \"dt\": 0.01,\n", " \"tau\": 1.0,\n", " \"I\": 2.1,\n", " \"eta\": -4.6,\n", " \"decimate_rv\": 10,\n", " \"RECORD_BOLD\": False,\n", " \"RECORD_RV\": True,\n", "}" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "obj = MPR(params)\n", "data = obj.run()" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "# obj.p[0,0,:,:]" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "(101, 1, 2, 6, 8)" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" } ], "source": [ "data['rv_d'].shape" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiwAAAGiCAYAAADEJZ3cAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjAsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvlHJYcgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAahRJREFUeJzt/Xl8VPXZP/6/ZjKTSTLZN7KSEBIWWQQEREEWF3BrpUWLeIOidSnWelu9b5F+ehe1tehdW7X6+1YrdwVF6l6tbS1oJSAUkFVBWQKBEEIg+ySZZNZzfn/MnJMEss3MmTlnZl7PxyOPh85MTt5zgOTK9b6u660TRVEEERERkYbp1V4AERER0UAYsBAREZHmMWAhIiIizWPAQkRERJrHgIWIiIg0jwELERERaR4DFiIiItI8BixERESkeQa1F6AUQRBw5swZJCUlQafTqb0cIiIiGgRRFNHW1oa8vDzo9X3nUSImYDlz5gwKCwvVXgYRERH5obq6GgUFBX0+HzEBS1JSEgDPG05OTlZ5NURERDQYra2tKCwslH+O9yViAhZpGyg5OZkBCxERUZgZqJyDRbdERESkeQxYiIiISPMYsBAREZHmMWAhIiIizWPAQkRERJrHgIWIiIg0jwELERERaR4DFiIiItI8BixERESkeQxYiIiISPMYsBAREZHmMWAhIiIizWPAQgRg+/FG/PnLU2ovg4iI+hAxpzUTBeK/3/sKp5s7MaU4DaXZ/R9xTkREoccMCxGAJqsDAHC6uVPllRARUW8YsFDUEwQRHQ43AKC+za7yaoiIqDc+BSw1NTV4/vnnMXfuXAwdOhSxsbHIycnBggULsHPnzkFfp7y8HDqdrs+PNWvW+Po+iPzW6XTL/13fzoCFiEiLfKphefHFF/HMM89g+PDhmDt3LrKyslBRUYEPP/wQH374IdavX4+FCxcO+nqzZs3C7NmzL3h8woQJviyLKCBWh0v+b2ZYiIi0yaeAZerUqSgvL8esWbN6PP7FF1/gqquuwrJlyzB//nyYTKZBXW/27Nl4/PHHfVkCkeI6Hd0yLAxYiIg0yactoe9///sXBCsAcMUVV2DOnDlobm7GgQMHFFscUShY7QxYiIi0TrG2ZqPR6LmgYfCXrKiowPPPP4/Ozk4UFBTgyiuvRH5+vlJLIhqUju5bQqxhISLSJEUCllOnTuGzzz5Dbm4uxo0bN+jPW79+PdavX9+1GIMBP/nJT/Cb3/wGMTExSiyNaEBWbgkREWlewG3NTqcTS5Ysgd1uxzPPPDOoQCMrKwtPP/00Dh48iPb2dpw7dw4ffvghSktL8dxzz+HRRx8d8Bp2ux2tra09Poj80WHvyrC02VywdesaIiIibQgoYBEEAUuXLsWWLVtwzz33YMmSJYP6vDFjxmD58uUYM2YMzGYzsrOzcdNNN2HTpk3IysrC73//e9TV1fV7jVWrViElJUX+KCwsDOStUBTrnmEBmGUhItIivwMWQRBw1113Yf369Vi8eDFefvnlgBeTk5ODm266CS6Xa8C5LitWrIDFYpE/qqurA/76FJ06u9WwAKxjISLSIr9qWARBwJ133onXX38dixYtwpo1a6DXKzM0NzMzEwBgtVr7fZ3JZBp0+zRRf5hhISLSPp+jjO7BysKFC/HGG28oWiArZVaKi4sVuyZRf7rXsAAMWIiItMingEXaBnr99ddxyy23YN26df0GKw0NDTh8+DAaGhp6PL5nz55eX//CCy9g06ZNKCsrw5QpU3xZGpHfmGEhItI+n7aEnnzySaxduxaJiYkYMWIEfvWrX13wmvnz58uj9V966SU88cQTWLlyZY+JtgsWLIDRaMTkyZNRUFAAq9WKHTt2YN++fUhNTR0wECJSkjSHxRijg9Mtoo4BCxGR5vgUsJw8eRIA0N7ejqeeeqrX1xQXFw94FtCyZcuwYcMGbNmyBY2NjdDr9SgqKsJDDz2ERx55BAUFBb4siygg0knNhekJqKy3MsNCRKRBOlEURbUXoYTW1lakpKTAYrEgOTlZ7eVQGLl77W58dugcrhyVjc8P1+HiwlR89OPpai+LiCgqDPbntzKtPURhTNoSKspIAAA0MMNCRKQ5DFgo6klFt8UZZgCeotsISTwSEUUMBiwU9aS25qHeDIvDLaC109XfpxARUYgxYKGoJxXdpifEIjnOU4de325Tc0lERHQeBiwU9aQaFrMpBllJnunJbG0mItIWBiwU9aQalvhYgxywsLWZiEhbGLBQVHO6BThcAgDAHBuDrKQ4AAxYiIi0hgELRbWObmP5E2INyEr0Zlh4YjMRkaYwYKGo1n0sf6xBzy0hIiKNYsBCUU3KsCTEerqDGLAQEWkTAxaKah12KWDxHLaZzYCFiEiTGLBQVLN6t4SkgIUZFiIibWLAQlGtawZLzy2hpg4HnG5BtXUREVFPDFgoqlnP2xJKS4hFjF4HUQSarA41l0ZERN0wYKGo1uktujV7i25j9DpkmGMBcFuIiEhLGLBQVJNrWLxbQgDrWIiItIgBC0U1ua3ZGCM/xoCFiEh7GLBQVLPapQxLt4CF026JiDSHAQtFtY7zalgAZliIiLSIAQtFtQ5HLxkWBixERJrDgIWimrWXDEs2T2wmItIcBiwU1Tq8NSzxsb1kWFjDQkSkGQxYKKr1lmHhlhARkfYwYKGo1l8NS7vdJT9PRETqYsBCUU06rbl7hsUcG4N471yWhjaO5yci0gIGLBTV5MFx3WpYdDqdnGWpa7Opsi4iIuqJAQtFNet5pzVLWMdCRKQtDFgoaomi2GuGBeC0WyIirWHAQlHL7hLgFkQAvQQszLAQEWkKAxaKWlJ2BQASYrklRESkZQxYKGpJLctxRj1i9Loez2V6t4Qa2tklRESkBQxYKGr1dvChxOydy8I5LERE2sCAhaKWtZex/BJpi6j7thEREamHAQtFrf4yLFIRbicDFiIiTWDAQlFLyrB0H8svkbIuVm4JERFpAgMWilqdzn5qWLyPMcNCRKQNDFgoalntvQ+N6/5YNNSwiKIoZ5uIiLSKAQtFrY4+xvIDXVtCnU43BO9wuUj1s78cwMRffopjde1qL4WIqE8MWChqSRmW3rqEum8TSVtHkUgURWz85hwcLgFfVbeovRwioj4xYKGoJWdYeglY4ox66HTS6yI3YDnXakej1TMcr9HKqb5EpF0MWChqSR1A54/lBwCdTod4Y+S3Nn9zxiL/N6f6EpGWMWChqCXPYemlrRnoKryN5Nbmb860yv/dwHOTiEjDGLBQ1OqQu4QuzLB0fzySt4QO1nRlWOrbGbAQkXYxYKGo1bUl1H+GJbK3hLplWLglREQaxoCFopaUOekrwxLp025bOhyoaemU/7+RGRYi0jAGLBS1pGFpfdWwRPq0Wym7InVJNVodET9zhojCFwMWilrSfJWBMiyRWsMidQhdNjwTAOAWRLR0OtVcEhFRnxiwUNSSBscN1CXUEaFbQlKGZeLQVKTEGwEADdwWIiKNYsBCUUsKRBKM0dklJHUIXZSXjMzEWABsbSYi7WLAQlFJEMSuotsBMyyRF7B0OFyobLACAMbkJSMz0QSArc1EpF0MWCgqdT8fyNznHBaprTnytoQO1bZBFIGsJBOyk+KQmeQJWBrZ2kxEGsWAhaKSlDXR6TznBvWmq6058jIsUsHt2LxkAECWN8PCGhYi0iqfA5aamho8//zzmDt3LoYOHYrY2Fjk5ORgwYIF2Llzp0/XEgQBL774IsaNG4f4+HhkZWVh0aJFqKys9HVZRD7pOvjQAJ10yuF5Irmt+ZsaT8HtmLwUAOiqYWHAQkQa5XPA8uKLL+KnP/0pKisrMXfuXDzyyCOYMWMGPvroI1x++eV4++23B32t++67Dw8++CBEUcSDDz6Ia6+9Fh988AGmTJmCiooKX5dGNGhWeSx/7/UrQPe25sjbEvqm1pNhGePNsGTIGRZuCRGRNvW+ed+PqVOnory8HLNmzerx+BdffIGrrroKy5Ytw/z582Eymfq9zqZNm7B69WrMnDkTn376KWJjPb/h3Xbbbbj++uvxwAMPYMOGDb4uj2hQOgYYy9/9uUjbEnK4BBw52wYAGJsvZVi4JURE2uZzhuX73//+BcEKAFxxxRWYM2cOmpubceDAgQGv8+qrrwIAfvnLX8rBCgBcd911mD17NjZu3IhTp075ujyiQbEOMJYfiNwtoYq6NjjdIpLjDChIiwfQtSXEolsi0ipFi26NRs/wKYNh4MRNeXk5zGYzpk+ffsFz8+bNAwBs3rxZyeURyToGGMsPRO6WkDQw7qK8ZLl+p3tbsyhyPD8RaY/PW0J9OXXqFD777DPk5uZi3Lhx/b7WarWitrYWY8eORUzMhT8wysrKAKDfOha73Q67vSt93dra2udric430MGHnuci87Tmb8/0LLgFPO3NgGe7qM3uQnKcUZW1ERH1RZEMi9PpxJIlS2C32/HMM8/0GoR0Z7F4Cv5SUlJ6fT45ObnH63qzatUqpKSkyB+FhYV+rp6ikdwl1E+GJVJrWKQJt2Pzk+XH4owxSDR5gjdOuyUiLQo4YBEEAUuXLsWWLVtwzz33YMmSJUqsa0ArVqyAxWKRP6qrq0PydSkySEFIfB9j+YGu7EskZVgEQcSh2gszLACQIbc2s46FiLQnoC0hQRBw1113Yf369Vi8eDFefvnlQX2elFnpK4Mibe/0lYEBAJPJNGAnElFfBlPDImVYHG4BTrcAY0z4z1k82WiF1eGGyaBHSaa5x3OZiSZUNXawU4iINMnv78CCIODOO+/E2rVrsWjRIqxZswZ6/eAuZzabkZubixMnTsDtvvC3V6l2RaplIVLaYLqE4ru1PEfKeUKV9Z7zg0qzE2E4LwDr6hRiwEJE2uNXwCIFK6+//joWLlyIN954Y8C6lfPNmjULVqsV27Ztu+A5af7KzJkz/Vke0YCkAMTczxyW2Bg9DHpPF02kbAu1dDoBdA2K666rU4hbQkSkPT4HLNI20Ouvv45bbrkF69at6zdYaWhowOHDh9HQ0NDj8XvvvRcA8D//8z9wOLq+QX7yyScoLy/H3LlzUVRU5OvyiAZFHhxn6jvDotPpIq61uaXD828tNf7CLiAOjyMiLfO5huXJJ5/E2rVrkZiYiBEjRuBXv/rVBa+ZP38+JkyYAAB46aWX8MQTT2DlypV4/PHH5dfMmTMHd999N1avXo1JkybhhhtuQG1tLd5++22kp6fjxRdf9PtNEQ1EGs3fX4YF8NSxtNlcEbMl1OrNsKT0FrB4W5vZJUREWuRzwHLy5EkAQHt7O5566qleX1NcXCwHLP155ZVXMG7cOPzxj3/ECy+8gMTERHzve9/DU089heHDh/u6NKJBkzIm8QMELJ5pt/aICVikLaHUhF4CFjMPQCQi7fI5YFmzZg3WrFkz6Nc//vjjPTIr3en1ejz44IN48MEHfV0GUUCscg1L//8EIm9LaOAMS6OVNSxEpD3h36dJ5IdOuYZl4C0hz+sjI8Ni6W9LKJFbQkSkXQxYKCp11bAMlGHxPB8p0267toRiL3hOamu2OtwREaA1ttvxz4O1EASejUQUCRiwUFQazGh+oKsotzNCtoT6K7pNNBlgMni+JURCHcvTnxzGj9btxZs7q9ReChEpgAELRSV5NP+ga1jCP+MAdGtr7qXoVqfT9Ti1Odwd8J6Z9Nevzqi8EiJSAgMWijpOtwCHSwAwuLZmIDK2hARBlGtYepvDAnSfdhvehbeCIOJko2eq7+6qZtS12VReEREFigELRZ3u2ZL+RvMDXTUukbAl1O5wQSrnSO4zYImM4XG1rTbYnJ6gVBSBDd+cU3lFRBQoBiwUdaSCUmOMDrGG/v8JRNKWkMXb0hxn1CPO2HtmKVI6hU54z0yS/PNgrUorISKlMGChqGOVWpoHyK54XhM5bc3SDJbU+As7hCSZSZExPK6yoR0AMConCQCwo7IJTZwvQxTWGLBQ1OkY5Fh+oCuosUbAllB/M1gkXVtC4f3DXTqVeuaILIzJS4ZbEPHpt2dVXhURBYIBC0Ud6yDH8gNdGZZI2BJq6fQEISm9dAhJIqVLqLLBE7CUZJpx3dgcAMAnBxmwEIUzBiwUdbpmsETrllDfAUuG3CUU3gHLCe+W0LBMM64dmwsA2HasQc4yEVH4YcBCUUfKliQMIsMSSZNuB7MllBUBW0J2lxunmzsBACVZiSjNTkRZdiKcbhH/OsRuIaJwxYCFok7HIMfye14TOZNuLf2c1CyRtoQsnU55Vk24qWrsgCgCSSaDPFeG20JE4Y8BC0Wddrt08OHAAUsktTV3Tbntu0soJd4Ig14HAGi0hue2kFRwOyzLDJ3O816uG+fZFtp8tF7+8yei8MKAhaKO1fsDK3GAc4SAri6hSAhYpAxLX0PjAECv18l1LA1t4bktJLU0l2Sa5cdG5SShOCMBDpeATYfr1FoaEQWAAQtFnXap6NaHLaEOhwuiGN6n/g6m6Bbo1tocphkWaWjcsMxE+TGdTicX3/6T20JEYYkBC0UdKcMymC4haUtIEAF7mNZ0SAZTwwIAGWE+7faE1NKcZe7x+LwxQwAAW481hHxNRBQ4BiwUdazeotvEQbU1d70m3FubB9MlBHQdgBiunULSDJZhmT0DltJsT8bF0umUg1YiCh8MWCjqtNkGn2GJ0etg8p43FO7Tbgczmh/o3tocfhmWlg6HPIL//IAlKc4oB6lnW3l6M1G4YcBCUadrS2jgolsgMobH2V1udDo96+9v0i0Q3ic2S9mVnOS4XgPSnJQ4AMBZCwMWonDDgIWijpQpSYobOMMCREankLQdpNN55pP0J5wPQOwquDX3+nyuN2CpZcBCFHYYsFDUkeZwDKZLCOgqvA3nLSFLR1f9it47Z6UvUoalMQxrWPoquJXkyhmWzpCtiYiUwYCFoo4vXUJA92m34ZthaekcXEszAKR5B8s1d4RfwFLZ7Qyh3uSkxAMI7wxLraUTT39yGHVt4fseiPzBgIWiji9dQkBkTLvtnmEZSJrZG7BYnWE3e0aacjs8K7HX53MjoIbl+U8r8PLm41j50TdqL4UopBiwUFQRRVHe2hlshqWrhiV8t4SkDEtKP2P5JWneolyHWwirIE0QRJxs7L+GJScCali2VzYCAP75zVlU1rervBqi0GHAQlGlw+GGlDQYbIYlIQIyLPI5QoPIsMQbY+RWbqlFOBzUttpgcwowxuhQkBbf62vkDEuYtjXXtHTiVFMHAEAUgT9uqVR5RUShw4CFoopUcKvXAXHGwf31j4SApXWQQ+MAzxh7qY5Fmt0SDqRsw9D0BBhiev+zzU32BDJNVgdszvD789zpza6ke7ftPthbg3NhGnwR+YoBC0WV9m4Ft9JJvgORtoQiouh2gBksEul14VR4e6LhwjOEzpccb0C80ROAhmMdyw5vwHLL5AJMKU6Dwy3gT1tPqLwqotBgwEJRReoQGmgWSXcJEdDW3OJD0S3Q9Rt8OAUsXQW3vdevAJ7sUTjPYtlR2QQAmFaSgWWzhwMA3tx5Sp6zQxTJGLBQVGn3saUZiIxJt4M9R0gitzaHUQ1LX2cInS83VapjCa9ZLFL9Soxeh8lFaZgzMhsjhySh3e7Cuh1Vai+PKOgYsFBUkVqafQlY4iNg0m3XltDAXUIAkGb2BDZNYVTDcsI7g6Wkj5ZmSU5yeM5ikepXxuanICnOCJ1Ohx/NLgEAvLbtRFjW5BD5ggELRRVpS2iwHUJA96Lb8N0Ssni3dnzNsLSEyZaQzenG6WZPxmTADEuYzmKR6lemlaTLj904Pg/5qfFoaHfgvT2n1VoaUUgwYKGo0u7jwYdAZHQJWXwuupVqWMIjw1Ld1AFR9NQmZSb2n0UK11ks3etXJMYYPe65YhgAT4tzuA36I/IFAxaKKr6O5QfC//BDQRC7ApZBF916u4TCpIZFmquSlxo/YPdXOGZYzpxXv9LdD6YUwhijw6mmjrALwoh8wYCFokp7AFtC4Vp022Z3QfD+4p08yIAlNczOEzrX6jlZOjvZNOBrwzHDsvNEz/qV7hJiDSjK8GyDVdRx8i1FLgYsFFUC6RIK17ZmaWhcnFGPOOPgtsLSw6xLSBqelpMcN+Brc70HIDa02+FwCUFdl1J2HJe2g9J7fb4s21NofIwBC0UwBiwUVfwrug3vwXHSDJbU+MF1CAHdT2wOjxoWKWAZMoiAJS3BiFjv0QPhMiV2xwmp4Daj1+dL5YClLWRrIgo1BiwUVXw9qRkI/6Lblk7fOoQAINVbw9LpdIdFu2xXwDLwllC4DY8709KJqsbe61ckpRGQYdl0uA5jfvFPvL79pNpLIY1iwEJRJaDBcU43BCH8ujDkoXGD7BACPN02Br2neDUc6li6algGzrAA6BawaH94XH/1KxIpYKmoaw/LTqF2uwsrPjgAq8ON3/+rAnaX9oNkCj0GLBRVuraEfGlr7gpuOsMg23C+ri2hwQcsOp2uq/DWqv1toToftoSArjqWcOgUGqh+BQCGZyVCp/P8WTeGSd1Rdy98dlTu9Gpod+CfB8+qvCLSIgYsFFX8ybDEGfWQOmXDcVvI17H8Erm1WeMZFkEQUdfmybAMZksICK9OoYHqVwAgzhiDwrQEAOG3LXSothV/2nYSADC91PMeX9/OowboQgxYKKr4E7DodDr5hN9wnHbr69A4Sbi0NjdaHXAJInQ6IDNxcAFLuMxiaWy3o6qxAzod+qxfkXTfFgoXgiDi5x8ehFsQcd3YHDz3gwkw6HXYU9WMgzUWtZdHGsOAhaKKP11CQHgPj5PG6w/2HCFJWoKUYdH2lpBUcJuZaIIxZnDf0qT251qNdwmd8B7omJcS32f9ikQKWI6HUcDy3p7T2FPVjITYGPziOxchOzkO143LBQC8wSwLnYcBC0UVfw4/BMK7U0jKsAx2aJwk3Rwes1jq2gbfISTpqmHRdtHtiUGeQA10z7CER2tzs9WBVZ8cAgD89OoR8p/J7ZcVAQA++qombM6yotBgwEJRw+ES4HB7BoUlxvoXsITjLBZ/im6B8NkSkjqEhiQNruAW6KphqWuzw+nW7vC4k42egKU4M2HA14Zba/P/bjiM5g4nRuUkYen0YvnxyUVpGJ2bDJtTwLu7eaAjdWHAQlFD2g4CfDv8EAjvabf+1rDIW0Iaz7BIW0KDbWkGgAxzLIwxOogiUO8t2NWikw0dAIDijMFnWM612tFq0/Y2ns3plk+XfvKmsT228nQ6nZxlWbezKixHCVBwMGChqCEV3MYZ9TAMstZBEs7Tbv3tEgqXabdyhsWHLSG9Xie3QGt5FosvW0LJcUb5Hmg9y3KwxgKnW0RWkglTii8sJr5pQh6S4gyoauzA5op6FVZIWsSAhaKGlB3xteAWAOLDuIbFn9H8QFfAovU6Al9nsEjyvDUTWm1tFkWx25bQwAELED7bQvtOtQAAJham9nq6dkKsAT+YXAiAxbfUhQELRQ2rHy3Nkq6i2/DaErK73PKwO18m3QJAmrfotknjActZH8byd5ej8dbm+jY7Ohxu6HWQZ6wMpCw7CUAYBCzVzQCAiUP7btVePM2zLbTpSJ1m/4wotHz/zk0Uptps3oDFx4JbIHzbmqXtIJ3OM27fF1INS4vGJ912bQn5lmHR+nlC0nZQQVqCfFjjQIaHW4ZlaGqfrxmWacbY/GQcrGnFnqpm3DA+NzSLC9Cftp7A27uqYXO54XAJsLsEuAUR988ejvtmDVd7eWGNGRaKGv4cfCgJ17ZmS0dX/Ypef2HqvT/SllCb3QWHS5udNE63gEarfwGL1jMsvm4HAUBZGLQ2n7XYUGuxQa8Dxhek9PvaCYWpAID93oyM1lk6nHj6n4dx5Fwbqho7UGuxocnqgKXTid//q0KuoyP/+BywrFu3Dvfddx8mT54Mk8kEnU6HNWvW+HSN8vJy6HS6Pj98vR7RYHRtCfnWIQSE75ZQS6d/Lc2AZ26LFONIJz5rTUO7HaIIGPQ6pPs4GE/rByCe8HYIDcsY3HYQ0FXDcrq5U7OnbEvBx6ic5B7ndPVmYqFny0jKyGjdR1/VwOESMGJIIt5fdhk+fmAGNjw0EyWZZlgdbny0v0btJYY1n3/V/PnPf46qqipkZmYiNzcXVVX+F0TNmjULs2fPvuDxCRMm+H1Nor74M5ZfErZbQh3+dQgBQIxeh5R4I5o7nGi2OpHtw5yTUJFPaU4y+ZxBytH4AYgnG3zPsGSYY5GaYERLhxPH69sxJq//DIYapOBjQj/bQRJpy+hAjQUOlzDorTG1vL2rGgCwaOpQXFLUdVjlbZcOxa/+fgjrdpzCbVOH9lpoTAPz+Tv36tWrUVZWhqKiIjz99NNYsWKF31989uzZePzxx/3+fCJfSBmWpDj/t4TCra1ZyrCk+Jh9kKSZYz0Bi0YLb/2ZwSKRMizn2uxwCyJifAx4gs2fLSGdToey7ETsOtmMY3XaDlgmerd7+jMs04yUeCMsnU4cPtuK8QUDf45aDtZY8M2ZVsTG6DF/Qn6P526+pAC/2XAEh2pbsa+6BZP6KTamvvkcrl599dUoKioKxlqIgqrd4X/RbXy4bglJ5wj5kWEBtN/afM7PDiHAc/ZQjF4HtyCioV1bw+MEoauledgghsZ1p+XWZqdbwNc1LQD67xCS6HQ6uY5F69tC7+z2ZFfmjhkid9hJUhNiceP4PADAmztOhXxtSlG7lk3V/FpFRQWef/55rFq1Cm+88QZqari/R8ETSFuzFORYwyzD0urn0DiJ1CnUpNFOoXN+zmABPFteQ5I8gY7WOoXOtdlgcwow6HUoSIv36XOHZ3kLb89pL2A5crYNNqeA5DgDSgaZOZK2hfZXtwRvYQGyOd34cJ/n59fCKYW9vuY/pg0FAPzt6zOa/QWgPx/sPY0bX/xC1ZovVdua169fj/Xr18v/bzAY8JOf/AS/+c1vEBPTf2Gk3W6H3d71W1Fra2vQ1kmRQYkuoXDdEvJ1LL8kTePnCfnb0izJTY3HGYsNtS2d8m/yWiC1NBemJ/g8lblsiHcWS732ApZ9pzwFtxOGpg265kjKxEifq0UbvjmLVpsL+anxmD48s9fXTCxMxUW5yfi2thXv7TmNu68oCfEq/ffZt+fw3+99Dbcg4t3dp/HgVWWqrEOVDEtWVhaefvppHDx4EO3t7Th37hw+/PBDlJaW4rnnnsOjjz464DVWrVqFlJQU+aOwsPeolkgiz2EJaNJtuG0JBZhh0fiJzYFkWICu1uYzGsuwdJ0hNPgOIYm0JXSywaq5gx19qV+RTPDWrZxs7ECTRv8eSsW2N19S0GcgptPp5CzL+p2nIIrhcUbSzspG/Hj9XrgFEd+flI8H5pSqthZVApYxY8Zg+fLlGDNmDMxmM7Kzs3HTTTdh06ZNyMrKwu9//3vU1dX1e40VK1bAYrHIH9XV1SFaPYUrZdqawyvD4u85QhKtnydU58c5Qt3lp3q2W2qatdXa7E/BrSQvJQ7m2Bi4BBFV3utohbSt09/AuPOlJBgxPMvs/XztZVlONXbg38cbodMBt0wu6Pe18yfkI9FkQGWDFduPN4Zohf47WGPB3Wt3w+4ScPXobDyzYLzP3XhK0lSPWE5ODm666Sa4XC7s3Lmz39eaTCYkJyf3+CDqTyBnCYVrW3OLQjUsWt1zP9cWWIZFqg+paelQbE1K8OXQw/PpdDpNTrxttjpQ6X1fvm6/TdDwPJb39nh+WZ5RmomCAY5QMJsM+N5ETwfRup3aPiOpsr4dd/zpS7TZXZg6LB0v3Tapx6naatBUwAIAmZme/T+rVVu/GVD4C2wOS3jWsDR6u18yEv3LQKQmaPc8IZvTLW95DfFzRoycYWnRWIZFmsHiY4eQpFSDhbf7T7cAAEoyzfLfq8HSauGtWxDx7p7TACAf1jgQaVto4zfnUNemra1IiaXDiSX/9yUarQ6MyUvG6jsmI87oe2ZaaZoLWKTMSnFxsboLoYgjbQkFUnTrcAuaqwvoiyh2tetmJvo3hyXdLLU1a29LSNoOMhn0SI73r38g35thOa2hLSFBEFHV5J1y60eGBQBKh3gzLBoqvPVlYNz55IDlVAsEQTu1H/8+3oBaiw2pCUbMHTNkUJ8zKicZEwpT4RJEfPZt/6UPanl79ynUtHRiaHoC1t41Fclx/mVolRbUgKWhoQGHDx9GQ0NDj8f37NnT6+tfeOEFbNq0CWVlZZgyZUowl0ZRKLAuoa7PCZdtoQ6HGzanJ7jK9DPD0tXWrL0Mi7QdlJMS5/fkUCnD0tLhlANatZ2xdHqmusbokZfqW0uzRDq1WUsZFqnLZzDzV843ckgS4o0xaLO7cFxDQVj5kXoAwLVjcmAyDD4DcfXobO/nay9gEUURf/7Ss811/+zhfn/vCAa/Jt1u3boVAHDgwAH5sfLycgDAjBkzcPfddwMAXnrpJTzxxBNYuXJlj4m2CxYsgNFoxOTJk1FQUACr1YodO3Zg3759SE1Nxbp16wZsaybyhSCIcg2LP1tCsQY9DHodXIKITofb75qQUJKyK3FGvZwh8pXUJdRqc2puGqzcIRTAkQFJcUZ5kmpNSydGeFuC1SR1CBWmx/t9v6VOoeP17Zr4cxMEsavg1o/2cUOMHuMLUrDzRBP2nWqRW7fVtu2Y55fxGWW9tzL3ZfbIbDy78Si2HWuA3eX2KdgJtu2VjTjRYEWiyYDvXJyn9nJ68Pk799atW7F27doej23btg3btm2T/18KWPqybNkybNiwAVu2bEFjYyP0ej2Kiorw0EMP4ZFHHkFBQf+V1kS+6nC6IXUR+pNhATytzW02V9i0Nje0e7IimYkmvzMQ0oRcUfR0HKWb/dtaCgb5HCE/O4Qk+anxsHQ6cbq5QxMBy4lG/wtuJYVp8Yg16GF3Cahp7sRQP9qjlVTZ0I42mwtxRj1G5fh3jycMTfUELNXN+EEfw9lCqa7NhsNnPadiX97H7JW+XJSbjKwkE+rb7Nh9shnTS337/GBav9MzifemCXl+/XIXTD5vCa1ZswaiKPb50f2k5ccffxyiKF5wXtDy5cvx+eefo6amBjabDR0dHTh06BCee+45BisUFFK6X6/zZBz8IU+7tYfHllBX/Yr/P9ANMXr57CWtbQvVBTiDRSJ3CmmkjiXQglvA8+cmTZKtqGtTZF2B2OutXxlfkOrzIDyJ1k5ultqSx+Ql+xzI6/U6zBqRBUBb20IN7XZs+OYsAM+BjVqjuaJbomDo3iHkd7ZBavHt1NYP7r40yhmWwLIiXYW32nrfgZwj1J1ceKuRTiF/TmnujbRtUqGB1uavvR1CgUwTlgpvj55rk/89q2lrhXc7yM/syJyRnjqWTd46GC14f89pON0iLi5I0eTBmQxYKCoE0iEkkYaoaS3T0BclMixAt9Zmjb3vswplWKTCW610CimxJQRoq7X5QI3n6JSx+f7/EBySHIf81HgIYlcApBZRFOX6FX+3c2aUZSJGr8OxunZUN6k/B0gQRPz5S892kBazKwADFooSgcxgkaSZpSFq2mvx7U2DPIMlwAxLgjbft9TWnB1A0S0AediXFraEXG5B/uEVeIZFG63NTreAQ7WegGVcAAELAM2c3HyiwYozFhtiY/SYUpzu1zVS4o24xNsxVX5U/SzLjspGnGzsQKLJIJ8srTUMWCgqBNLSLAm3DEtjt6LbQGj1AERpS0g6D8hfBRqaxXKmxQanW4TJoEdugJmjMmna7bk2Vc+tOVbXDodLQJLJgKL0wIp/pW0htQ9ClLIrlxSlyeeM+WP2KG8dy2H161je9GZX5k/UXrGthAELRQUlt4S0VsvRl3qlt4Q09L7b7S5YvfNwspMC7xICPBkpm1PdgmppO6goIyHgM1uKMsyI0etgdbhRq+LhjgdqLACAMfnJAb+nSUWejMTeUy2qBmHbjnkKbn1tZz7f7BGeOpZ/H29U9e9eQ7sdG6Vi26lFqq1jIAxYKCq0B3DwoUSaSdKksa2RvjQqtSUkbYVZtfO+pexKkskQ8G+DqQlGmL2/JZ9RufBWiQ4hSaxBL5/2rOaZQge9AUug20GApyMn1qBHk9WBk43q1H24BRH/Pu7JsFw+PCOga43OTcKQZBM6nW58eaJJieX55T2p2LYwFRflafdcPgYsFBUUqWHR+EGA55PmsGRFYIZFClgCncECeA4L1MqI/kAOPeyNPPFWxYBFyrAEUnArMRliMN57nT1V6mwLHayxoNXmQlKcIeAgTKfTyd1C5Sp1C4miiHd3eybb3jZV/fk2/WHAQlFBkS0hc/jUsDhcAizek5r9PfhQosW25nMKdQhJ5MJblTMs0tj54d4On0DJhbcqzWJxKVhwK7nEuy2kVsCy1Vu/cllJht8zZbqbPVLdeSzfnGnF8XorTAY9rh+Xq8oaBosBC0UFJTIs6QnaPQjwfFJQFaPXydNq/SXNn2nW0PuWptwqFbDIpzarnGGRWpCHZysTsEgj+tXaEjpW3w6bU0CiyaDINhfQrY5FpYDF33H8fZlemgmDXofKBiuqvDVMofTxV2cAAFeNzkaSRg457AsDFooK0TaHRWppTjfHBlzoKHcJaeh9K7klBHQ/tVm9eRhtNqc8W6ZU4YDl6Ll2VYpUD5z2bAddlBd4wa1kkrcV+Ghdm5xFDJVOhxu7T3oCJaXG6SfFGTG52NveHOJtIUEQ8VdvwPLdi/ND+rX9wYCFooLU1mwOoAVRmsPS6XSr3k0yEKWGxgHdtoQ6nRAE9TozulPi4MPu5PH8Km4JHa/3/HadnWRS7HDN4VmJ0Ok850BJNU2hpGTBrSQryYSijASIIuQDFUNld1UTHG4BuSlx8tEHSpDqWP4V4vbm3VXNqLXYkGQyyFtTWsaAhaKCtCWUGEDKM9FkgMH7W6LWZpKcr0GhsfxA15aQWxDRZlN/JDrQtXUjZUYCpYUtoYpznjoTpbIrABBnjEGhtz5HjTOFDgQhYAEgD1wLdR2L1M48vTTT7yM+enPV6CEAgO3HG9BqC13W6KP9NQCAa8fmIM6onROj+8KAhaJC15aQ//8odTpd2BTeNiqYYTEZYuTMlFYCtVPeabDSD+NASYHP2VYbnG5BkWv6SppIq2TAAnQNkDse4joWl1vAt7WBj+TvjVp1LFuPebZs/D0/qC+l2YkozU6E0y1iU4iyLE63gH8cqAUAfHeCNifbno8BC0UFJYpuge6tzdopQO1N15ZQ4BkWoKu1WQsBS7vdJRcAF6Yrk2HJSjTBZNBDEIGzKg1ZkwKKMoUDllJvp1CoW5uP11thcwowx8Youn0CdHUK7TvVDHeItinr2+w46D0T6fLSwOav9ObaMTkAgE8OnFX82r3ZWtGA5g4nMhNNuKxE+fcTDAxYKCooF7CES4bFs75AW5olUv2OFt63dNZOWoJRsa4GnU4nbwtVq1R4KwUUSnUISdQ6BFGecJuXoljBrWTEkCQkmgywOtw4cjY0W11bvOf9jM1PDvj8qt5cO9YTsJQfrUOnI/g1clKx7Y3jcxVpzw6F8FglUYCU6BICwmc8v1Jj+SXSN+i6Nrsi1wuEvB0U4Lk055O2hdSoY7E53XIgpviW0BDP8LhQH4J4UMGBceeL0evkc4X2huhcIemAQmmcvtLG5CWjIC0eNqeAzUE+DLHT4cYG7yj+cNkOAhiwUJSQu4QCDVjkGhatbwlJGRZltoTyUj0Bi9qj64GuDItS9SsSNTuFTjRYIYhAcpwh4MnE55MCoPo2e0gDbbngtiA4o94nDg1dHYvLLcgZlmB10+h0OnlbSAomguVfh8+hw+FGYXo8JnpPwA4HDFgo4jlcAhzeQsrE2MACFulcHS3UcvRHKrpV6odfnga6aCTS+HzFMyyp6o3nl7aDyoYkKdp9AniyinneE61DNUDOLYj49oyyE27PJ0+8DUGG5avTLbB0OpESb8SEIP6Al7aFPjt0Dg5X8Iq/P9ovzV7JU/zvWzAxYKGIJ20HAYEdfgh0G6Km4YBFEEQ0WqW2ZoUGq6WqP6dE0rUlpEzBrUTNLSEpkChVaCT/+YaHeOJtZX07Op1uJMTGYFhmcN7ThMJU6HRAVWMH6oO8VSkNdLuiLDOo9R6ThqYhK8mENptLPmBRaZZOJzZ73084DIvrjgELRTyp4DbOqA/4m004FN1aOp1y54Q09C1QUsByxqJ+wBK8LSH1zhOSOoSUrl+RhPoQxK6C22TEKFxwK0mJN2KE930Fu45FClhmjwxO/YpEr9dh3hjPTJZgbQt9uK8GDreAUTlJGJmTFJSvESwMWCjiWR3KFNwCXd0yWm5rllqaU+KNiDUo809c2hI6a7GFrI20N6Ioyls2Q4O0JXSmpTPk71Ea6ia1ICutLMStzUqe0NyfUMxjqW+zy+9n1ojgT4O9doznAMKN35xT/O+hKIpYu/0kAOC2S4cqeu1QYMBCEa/dpkxLMxAeGRalC24Bz7j4GL0OTrcoB0RqaGh3oNPphk7XFUQpZUhyHAx6HVyCiLq20M1icbkFnGjwjOUP1paQfAjiudC0AAdjJH9vQnFyc/d25qwkZQuie3NpSTpSE4xotDqw62STotfedqwRlfVWJJoM+P6kAkWvHQoMWCjiyTNYAiy4BcKjrVnJc4Qkhhg9crwnI6tRlCqR6ldyk+MUyx5JYvQ65KaG/j2eauqA0y0i3hgjZ3mUJg2jO2OxBf3AQJdbwDdBLriVSAHL1zUW2F3BmV0S7Hbm8xlj9LjaO6r/nweV3RaSsisLJuUrknEONQYsFPGklmZltoQ8AYvV4Q7aN8hANSjcISTpvmWiFuk0ZaU7hCRqnCnUNTDOrPiANUlqQqy8hbYvyPUeX9dY0OFwIzXBiOFByhhJijMSkJkYC4dLkE9RVlIo2pl70729WalTtk83d+Bfh84BAJZcVqzINUONAQtFPKs85Tbww72S4wxyEaFW61gag7AlBGhjFkt1kIbGSdQovA12h5BkcnFoDgzcftxzQOClw9KDFoBJdDodrhoVnGwE0NXOnBxnCGo78/lmlGXCHBuDWosNuxX683pz5ykIIjC9NCNoxd3BxoCFIp5SY/kB7wGICdoZU9+bYGwJAV01I2oGLEofeni+rlksoRvPH+wOIcnkonQAULwu4nw7Kj0BS6jOp7l2XFc2QlC4SFVuZx6RFdLx9XHGGNww3lN8+3L58YCvZ3O68daXpwAAt4dpdgVgwEJRQKmx/BJ5FotmA5ZgZVjUn8VS3SQNjQtOrYc0iyWUNSwVIQpYpngzLPurW4J2InX3rZnLhit7onFfLh+egSSTAXVtduyrblH02lLAMifI7cy9+dGs4dDpgH8drpOH8Pnrb1/XornDifzUeFw1KvTvRSkMWCjitSvY1gx0Hx6nzS2hYGVY5MFqLeqcZgx0HUyodEuzJNTj+QVBxPF6KWAJ7kyM4VmJSIk3wubsKopV2lenW9DpdCPDHIsRQWrRPp/JEIMrR3t+CCs5uyTU7cznK8lKxA3jPFmW/1/5sYCu9Ua3VuZwOeiwN+G7cqJBsiq4JQR0O7lYo51CXQGLshkWtYtunW5B/trBqmGRtppON3UGLQvRXW2rDR0ONwx6HYoygvOeJHq9Tu6q2R2kbaEd3vqVaSUZIR35LhWpfnKwVrEi1c+8BaqhamfuzY/nlAIA/nGg1u8pxfurW/DVaQtiY/S4dUqhkssLOQYsFPGkOSxKZ1haNLolJBXdKp1hyfWeR2PpdMp1QaFU22KDIAKxBr3iHVCS/NR4JJkMcLiFkIyxr/DORSnONMMYgt98g114u71SCljSg3L9vswamYU4ox7VTZ34tlaZ7JFU83HjePVOMx6dm4yrRw+BKAJ/8LOW5dUvKgEAN16ci4wg/bsJFQYsFPHaFTqpWSKf2KzBDIvV7kKn0/N+lQ5YkuKMSI7z3EM1sizSdlBhWnzQuk/0eh1G53lOFw7Wtkl3UlBUFqKuja7C22bFMhESu8stB0KXDQ9Nwa0kIdYgb9tsUKBb6GCNBV+dtsAYo8Mtl6g7YO2BKz1Zlg/318hdcoO18Zuz+PvXtdDpgLumDwvG8kKKAQtFPCXbmgHIXUJabGuWsitxRj0SYpV5v93lq3jeTrBbmiVjvAFLoIWOg9FVvxKagGV8QQpiY/RoaLfLHVdK2XeqBXaXgKwkU9Dnr/RGOun4EwUCljd3nvJeU/2sxITCVFxRlgm3IOLlzYPPsjS027HigwMAgPtmDg/6MQmhwICFIp6SZwkB2h7PX9+t4DYYNQT5Ks5iCXZLs2RMnucb+zdnLEH9OgBQcS60AUucMQZj8z0B2S6FB61tV6l+RXLlqCEw6HWoqGsPaDuv3e7CX/fXAAD+QyPn7TzgrWV5d/dpnLUMXPQuiiIee/8AGq0OjMpJwk+vKQv2EkOCAQtFPCXnsADaHs/f6A1YgvVboZqzWKqDdOjh+S7K9WZYalsV3zbpThRFHAtxhgUAJhd7toX2VClbeLs9xPNXzpcSb8TlpZ5W6kC6hT7cVwOrw43hWWZcOiy0tTh9ubQkA1OK0+BwC3h245EB5828u/s0Pjt0DrExejy3cAJMBuWzrWpgwEIRT/E5LBquYZFmsGQp3CEkyVNhdL2ka0soODNYJGVDEhEbo0ebzSXPfQmGujY7Wjqc0OmAkswQBixyp5ByGRab0439p1oAhL5+pbvrxnYNkfOHKIrydtCiqUNVyRT15aGrRwAA3ttzGv/59n7YnL0fDVLd1IEnPv4GAPDw3BEY7Q3AIwEDFopooijK81JS4o2KXDPdLA2O014Ni9TSnGEOXhcNAJxRYRaLFLAUBHlLyBijx4gcTwDxbW3wtoWkAtWRQ5IQH4R6o75Irc0Vde2KZQn3VDXD4RaQkxyH4iC3Z/fnmouGQKcDvj5t8Wta8f7qFhyqbUWsQY+bVS62Pd/00kz8783jYdDr8PFXZ7B49c4LtqXPWmx4+J39sDrcmFqcjnuuKFFptcHBgIUiWqvNBYfLM09DqVkKUtFtu73r2lohbQllJgU5wxLiLSGr3YVG7zfnoSH4gShtCwWzU0jKcEwpDu22Q0aiCSVZZgDKtTdL9SuXDVenfkWSmWiS76c/Zwut92ZXbhyXi9SE4PwbCsQPJhfi9bumIinOgN1Vzfje/7cNX55owp+2nsDNf/g3pq36F3adbIY5Nga//cHF8rlnkYIBC0W0+jbPD/CkOAPijMr8FpscZ4T0fUBrdSwNQZrBIpEyLGdbbXArfG5Lf6RR+SnxRiTHKZMp609X4W0QAxZvDYk0GyWUpG0hpQpv1a5f6e5677bQn7aekLeDB8PS6cTHX58BAPzHNG0U2/bm8tJMfLDschSkxaOqsQM/eGU7nvzbt/IhiZOL0vDq7ZOD3k2nBgYsFNGkgEXJSZV6vU7+7Utr4/kbglx0m5VkgkGvg1sQUdcWum2hUyGqX5EEu7XZanfJwdDkEGdYgK55LEoU3lrtLnzlPcNHzfoVycIpQ1GQFo8zFht+9+nRQX/eX/aehs0pYOSQJEwaGvog0hdlQ5Lwl/unY+LQVACeIOUXN16E7SuuxHvLLpeLjyONMlWIRBol/VDNVni0dlqCEU1Wh+Zam4M1ll8So9chJyUOp5s7UdPcidyU0AQQUv1KsDuEJKNyk6HTeTJJje12xQPAr6pb4BZE5KXEyVmrUJKyOl+dtsDucgfURbK7qhkuQUR+arwmfquPj43Br+aPxdLXduG1bSfwvYn5A84gsXQ48X/bTgDwnLejpWLbvmQlmfD+jy5Hm92lWH2e1jHDQhGtK8MSp+h1tdraLNV5BGtLCOjaFgplHUvXlNvQ/EBMNBlQnOGp8wjGtpC0FaNGdgUAhmWakWGOhcMl4GBNYIXFG70dOZdrILsimT0yG9+5OA+CCKz44ABc/ZwL5XILeODPe1Hd1Im8lDh8f1J+CFcaGL1eFzXBCsCAhSKcNEhN6bNntNja7HQL8vTdUAQsoewUktqLC0L4G/xFeV3zWJQm1a9MUaF+BQB0uq6DEL884X8dS7vdhQ/3eYasfW+itn7Q/8+No5EUZ8CBGgte317V5+t+9fdD+KKiAfHGGLx6x2QkhaBGivzDgIUiWn2rJ2DJTlZ+SwgAmjW0JSSN5Y/R65AaxN+61BgeF+otIaCrjkXpDIvLLWCvt0DykiL1BpNdUeapc3h/72m/B+T9xTtkrSTLrIn6le6yk+Lw2HWjAAC/3Xik17+vb+6swpp/nwQAPLfwYrnYmrSJAQtFtGBnWLRUdCvVr6SbY4N2OCAQ+tZmURR7HHwYKl2tzcrOYjl8tg1WhxtJJgNG5iQpem1fzJ+YD3NsDI7VtePf3rZkX4iiiDd3eDIXiy8t0mTdx6IpQ3FJURqsDjcefe9rbDvWgFpLJ0RRxL+PN2DlR54Ba/81dwSuHZur8mppICy6pYgWjC4hAEiXuoQ0lGHpGhoX3PkR+WmhzbA0WR3ocLih03V97VCQfts+0WCF1e5S7GiH3Sc920GTitJUnZORFGfEgksK8Pp2T5Zhuo+dJXuqmnH4bBvijHos0NiQNYler8OvvzcON/z+C2w91oCtxxoAAObYGLhFES5BxHcvzsOPvWf1kLYxw0IRra4tWFtCUoZFOwFL1yTY4P5Qlw5ADFWGpbLBCgDITY4L6ZkoWUkmZCeZIIqerIhSpHkZatWvdHf7ZcUAgH8dOif//RmsN7zZle9enKfpws+ROUn44+2X4OrRQ1CSaUaMXgerww2bU8DFhan435vHazI7RBdihoUiltMtyG3HwSu61c6W0IkGzw8cqbslWKRW5jabC602Z9AHuR047dmSuUiF+oKL8pJRd6Qe356xyEWqgRBFEbu8GRY161ckpdmJuKIsE19UNGDdjiqsuH70oD6vsd2OTw54uoMWTysK5hIVceWoIbhy1BAAgMMl4FRTB2paOnFJUZpiAyUp+JhhoYjVvQg1TeEx21LRrZbamk82ejIRxZnBDVjMJgNSve8/FNtCB7xtt+MLQh+wKF14e7q5E+da7TDodZhQmKrINQN1hzfL8tauanQ6ej9Q73zv7D4Nh1vA+IIUjC9IDd7igiDWoEdpdiJmjchS7EBUCg0GLBSxpPqVzETli1DlDIuGalhOerdOhgU5YAG6tzYHP2D5+nQLAGCcKgGL52sq1dosnd0zNj8lpAce9mfOqGwUpsfD0unER/trBny9WxCx/suuYluiUGHAQhGra8qtskPjgK4aljabC85+hlKFisstyOPrg51hAbp3CgV3FkubzSnXsIwbYFppMEidQofPtiny5yxtB01WYHtJKTF6HW6fVgwAWPPvkwO2OG85Wo/qpk4kxxnwnYvzQrBCIg8GLBSxgtUhBHgO4dPJByCqX8dS09IJlyDCZNAjN1n5AO18ocqwfHOmFaLo+XrBHIbXl6HpCUg0GeBwCThe3x7w9XarPOG2L7dMLkCcUY/DZ9vw5Ym+zxcSRRGvbz8JALj5kkLNZIkoOvgcsKxbtw733XcfJk+eDJPJBJ1OhzVr1vj8hQVBwIsvvohx48YhPj4eWVlZWLRoESorK32+FlFv5IAlCD/oYrqNxNZCp9AJbxaiKCMhqDNYJHneTqFgByxSwa0a2RXA0xYrZVkO1gS2LWTpcOJonafbSI0TmvuTmhArT6pd6w1IzicIIp74+FtsOlIPnU7bJxpTZPI5YPn5z3+OP/7xj6iqqkJurv+Ddu677z48+OCDEEURDz74IK699lp88MEHmDJlCioqKvy+LpEkWC3NEi3NYpHqV4LdISTJT/VMnK1pDm7A8rW34FaN+hXJJO/2zaffng3oOntPNUMUgZJMsyrZooHccXkxAOAfB87iobf2odE71wfwdNb859v75amwK2+8CMOzElVYJUUznwOW1atX4+TJk6ivr8ePfvQjv77opk2bsHr1asycORN79+7FM888gzfeeAMffvghmpqa8MADD/h1XaLugrklBHSfdquBgKXRU78SioJbwJPJAYCj59ogCP6NdR+MA96CWzU6hCQ3TfDUaWw6XB9QV9iOE55pskq0RwfDqJxkPHR1GfQ64MP9Z3D17zbjg72n0W534a41u/DxV2dg0Ovwwq0TsHT6MLWXS1HI54Dl6quvRlFRYJXhr776KgDgl7/8JWJju9pNr7vuOsyePRsbN27EqVOnAvoaRMEayy+RzxPSQA2LtCUUioJbwDOMy2TQo9XmwglvO7XSLB1OORBTa0sIAEbnJmNUThIcbgF/P1Dr1zWcbgF/2evpwJk1MkvJ5SnqoatH4C/3T8eonCQ0dzjx8DtfYfrTn2PrsQYkxMbgT0un4KYJ2jrkkKKHKkW35eXlMJvNmD59+gXPzZs3DwCwefPmUC+LIkzQMywJ2mltlmewhGhLyBijl7Me+061BOVrHPSe4TM0PQGpCs/R8ZVU3yGdTOyrz749h7o2OzITTZh7UY6SS1PcxYWp+PgnM/DotSMRa9DD0ulEujkW6++ZhpkjtBtsUeQLecBitVpRW1uLYcOGISbmwgrzsrIyABiwjsVut6O1tbXHB5FEFMWgtjUDXVtCag+Pc7oFnPbWkoRqSwiAPPhsf3VzUK7/9Wn161ckN03Ih04H7DrZ7PMIe6BrjP2tUwoRa9B+c6YxRo/7Z5diw0Mz8dDVZfhg2eWaGXRH0Svk/3IsFs83oZSU3r8JJScn93hdX1atWoWUlBT5o7CwUNmFUlhrt7tgc3rmZmQmBee3864Mi7pbQtVNHXALIuKNMRgSpALj3kwc6qnFCFaG5UCN57rjVdwOkuSkxOHy4RkAfM+yHK/3nIas1wGLLg2vzpphmWY8dPWIkG01EvVH+6F+H1asWAGLxSJ/VFdXq70k0hBpOyjRZEBCbHDGb6ebPTUsDd26KdQgbQcVZSSE9BC3iUNTAXiGqg12pLsvtJRhAYDvTfScSPyXfTUDDlfr7s0dnnq8K0dly/NriMh3IQ9YpMxKXxkUaWunrwyMxGQyITk5uccHkURuaQ5S/QoADE33/NZ5MkhFp4MlHXoYyu0gwHMI4pBkE9yCKJ/3o5Qmq0Pe5hqrgQwLAFw7NgdxRj0qG6xyMDWQTocb7+3x/DIVDocEEmlZyAMWs9mM3NxcnDhxAm73hb+VSbUrUi0LkT/kc4SCGLAMz/IECNVNHbC7lM8wDNbJEHcIdTexUNoWUraORQqASjLNQT8NerASTQa5YPYvg9wW+virM2i1uTA0PQEzy1iwShQIVbaEZs2aBavVim3btl3w3IYNGwAAM2fODPWyKIIEu0NIunaiyQBBBKoafS/EVIqU4RkWog6h7qRtof3VLYpe94CKBx7253uTPN1CH391ZlBnC63b6Sm2ve3SoSGZQEwUyYIasDQ0NODw4cNoaGjo8fi9994LAPif//kfOBxdHRaffPIJysvLMXfu3IBnvVB0C8WWkE6nk7MslQqcM+OvUM9g6U7qHFG68PZrlUfy9+WK0kxkJsai0erA1oqGfl/7VXULvj5tQaxBjx9MZlMAUaD8mnS7dOlSLF26FO++++4Fj61evVp+7UsvvYTRo0fjpZde6nGNOXPm4O6778aWLVswadIkLF++HLfffjvmz5+P9PR0vPjiiwG+LYp2ociwAECJdzz58Xp16ljsLrd8nk9xZkLIv/64ghTE6HU422pDrUW5Mf3SltD4glTFrqkEQ4xePqH4/b2n+33tOm8r8w3jcpFuVneODFEk8Ll9YuvWrVi7dm2Px7Zt29Zje+fuu+8e8DqvvPIKxo0bhz/+8Y944YUXkJiYiO9973t46qmnMHz4cF+XRdRDsKfcSqQMixIn+fqjuqkDggiYY2OC/l57kxBrwMghSfi2thX7T7Ugd1zgXTB1bTbUWmzQ6YAxedorpv/+xAK8tu0k/vZ1LUYOqcADV5Ze0J31z4Nn8dFXZwAAi3lIIJEifM6wrFmzBqIo9vnR/eTmxx9/HKIo4vHHH7/wC+v1ePDBB3Hw4EHYbDY0NDTgrbfeYrBCigh1hqVSpQyL1CFUnGkOaUtzd1Idyz6F6lgOerMrpVmJMJuC05IeiHEFKbhvVgkA4LefHsV/vfs1HC5PPYvLLWDVJ4fwo3V74HAJuHJUNiYN1ebZQUThJmznsBD1pz7IU24lJd1qWHyZzaEUNTuEJNIAuf0K1bFobf5Kb1ZcNxq/nD8WMXod3t97Grf/aSeO1bVjyf99iVc2VwIA7rliGF5ZcolqgSRRpNHery9EAXK5BTR6z/cJdoalOMMMnQ5otbnQ0O4I+tc73wkVO4QkUuHt1zUtcLoFGGMC+z1o18kmANoruD3fkmlFKEyLxwPr92FHZROu/p3n/DNzbAx+c8vFuH5crsorJIoszLBQxGmyOiCKgF6HoBc7xhljUJDmqdtQo1NICxkWz6wUA2xOAUfOtgV0rbpWG7YfbwQAXDVqiBLLC6rZI7Px3rLLkJfiyeSVZifiowdmMFghCgIGLBRxpJbmzEQTYkIw+2K4ip1CUsAyTIUOIYler8PFUntzgHUsf/3qDAQRuKQoDUMz1HtPvhiVk4yPfzIDL9w6AR/9eDpKsxPVXhJRRGLAQhEnVAW3kpJMqfA2tBkWm9ONMxZPrU6xiltCQPeDEAObePvhfs8E2fkT8wNeUyhlJJpw04R8TRYJE0UKBiwUcUIesEiFtw2hzbBI03WT4gyqz/lQYuJtxbk2HKxphUGvw43cUiGi8zBgoYgTqhkskq4todBmWE7I20HqtTRLJngHvFXWW2HpcPp1Del8ntkjs5HGQWtEdB4GLBRx6lq9Lc3JoQpY1DkEUTpDSO3tIABIM8fKp0V/6e3y8YUgiPhov2fQ2vfCbDuIiEKDAQtFnFBnWNQ6BFELHULdzRmZDQBY8+8TPn/urpNNqGnpRJLJgKtGZyu9NCKKAAxYKOJ01bAEd2icRK1DEE9ooEOoux9eMQwGvQ7bjjXiKx9rWaRi2+vH5SLOGBOE1RFRuGPAQhFHPqk5RFtCQOgPQXS5BXmE/agcbZy3k58aj5smeLZz/lB+fNCfZ3O68bevawGEX3cQEYUOAxaKOHKGJYSHAYb6EMSDZ1phdbiRmmDEyCFJIfmag/Ej7xk7G749i2N1g7sXmw7Xoc3mQm5KHC4dlh7M5RFRGGPAQhHFanehw+EpfA3lmPxQH4K4o9IzDXZqcTr0IRiON1hlQ5JwzUVDIIrAK5sHl2WRuoNumpCvqfdCRNrCgIUiirQdlBAbE9IhXqE+BFEKWC4tyQj61/LVstmeE9f/sq8GZ1o6+31traUTm47UAQC+P4nbQUTUNwYsFFGk7aDsEB9CeP4hiMHkcgvYfdIzUXZaifa2UCYNTcO0knS4BBGrv+i7Y6jV5sSdr+2C0y1iQmEqRmhoa4uItIcBC0WUUE+5lYTyEMRvzrSi3e5CcpxBMwW357t/dikA4M9fnkKT9cIAzu5y477X9+Dw2TZkJZnw4qKJoV4iEYUZBiwUUc56h8aFOmABus4UCnan0M4T3vqVYRkhOdzRH1eUZWJMXjI6nW78duMRWO0u+TlBEPHIO19he2UjEk0GrLlzCgrTtdGaTUTaxYCFIsqxujYAXePyQ2l4VmgOQdxR6Zkkq8XtIIlOp5OzLG/uPIWpT32G5e99jT1Vzfj1Pw7hb1/XwhijwytLLsGYvBSVV0tE4YBHi1JEOXzWE7CoUQ8RikMQ3YKIXSekgEV7BbfdXT8uB098dwxe23YCJxs78Pbuary9u1p+/jc3X4zppZkqrpCIwgkDFooYoijiqDdgGZUT+oAlFIcgfnumFW12F5LiDBidq836FYlOp8Mdlxfj9suK8OWJJry9uxr/OFALm1PAiutGcUgcEfmEAQtFjNPNnbA63DDG6FQ5X+f8QxBNBuVHzMv1K8Xpmq1fOZ9Op8OlJRm4tCQDT3x3DM612lCazY4gIvINa1goYhw911W/YowJ/V/t7ocgBmuAXNf8Fe3Wr/QnKc7IYIWI/MKAhSLGEW/AMlKF7SDAk0mYUJgKAPiiol7x67sFETvDpH6FiEhpDFgoYhxRseBWMnfMEADAhm/OKX7tQ7WtaLO5kGgy4CKN168QESmNAQtFDClgUfMwwGsu8gQse081o67Npui1pe2gKcVpMKiw5UVEpCZ+16OI4HQLct2IWltCAJCbEo+LC1IgisBn39Ypem1pO0iL5wcREQUbAxaKCCcbrHC4BZhjY5CfGq/qWuaOyQEAbPz2rGLXFAQRX7J+hYiiGAMWighSwW3ZkCToVW73neetY/n3sUa02ZyKXPPQ2VZYOp0wx8ZgbB7rV4go+jBgoYhwRMWBcecrzU5CSZYZDreA8iPKdAu9t+c0AOCy4RmsXyGiqMTvfBQRtNAh1N3cizzbQhu+CXxbyNLpxDu7PCPtl1xWHPD1iIjCEQMWigjS0DgtZFiArm2h8iP1sLvcAV3r7V2nYHW4MWJIImaW8ewdIopODFgo7HU4XKhq6gAAjNBIwHJxQSqyk0xot7uw/Xij39dxugWs2XYSAHD3jBLodOExjp+ISGkMWCjsHatrhygCGeZYZCaa1F4OAECv1ykyRO4fB2pxxmJDZmIsvjshT6nlERGFHQYsFPbkgXEaya5IpDqWT789B0EQff58URTxf1tPAACWTCtGnFH5wxSJiMIFAxYKe1oruJVMK8lAUpwBDe127Ktu9vnzd51sxtenLTAZ9Fg8bWgQVkhEFD4YsFDYU/vQw77EGvS4clQ2AOAfB3zvFlr9RSUA4PuTCpChka0uIiK1MGChsHdUowELANwwLhcA8Mb2Khw4bRn0551osOLTQ57alx/OGBaUtRERhRMGLBTWWjocONdqBwCUZSeqvJoLXXPREMy9aAgcbgE/Xr8XrYOcfPvK5uMQReDKUdko1eD7IiIKNQYsFNak+pX81HgkxRlVXs2FdDodfnPzxchPjceppg6s+OAARLH/AtwX/1WBt7yD4u65oiQUyyQi0jwGLBTWjmhsYFxvUhKMeOm2iTDodfj717V4c+epPl/7/GdH8dtPjwIA/nveSFw2nAcdEhEBDFgozMkdQhoOWABg4tA0LL92FADgyb99i2/O9KxnEUURv/v0KJ7/rAIA8Nh1o/DjOaUhXycRkVYZ1F4AUSC0NpK/Pz+cMQzbKxvx+eE63L12N2aUZmJIchyyk02orLdizb9PAgD+3/Wjcc9MbgUREXXHgIXClsst4HCtNmew9Eav1+G3t1yMG37/Bc5YbHjXewJzdz+/YTTuZt0KEdEFGLBQ2Np7qgVtdhfSEoxhEbAAQJo5Fn9/8Ap8dugc6trsONdqQ12rHZZOJxZcUoCbLylQe4lERJrEgIXC1ueH6wAAs0ZkIUYfPocCppljccvkQrWXQUQUVlh0S2Gr/IgnYJnjnSZLRESRiwELhaUzLZ04fLYNeh0wsyxL7eUQEVGQMWChsLTJm12ZODQNaeZYlVdDRETBxoCFwtKmw/UAgDkjmV0hIooGDFgo7Nhdbmw71gCA9StERNGCAQuFnS9PNKHT6UZ2kgkX5SarvRwiIgoBvwKWXbt24frrr0dqairMZjOmTZuGd955Z9Cfv2bNGuh0uj4/ysvL/VkWRQmpnXnOyGzodOHTzkxERP7zeQ7Lpk2bMG/ePMTFxeHWW29FUlIS3n//fSxcuBDV1dV45JFHBn2tm266CRMmTLjg8eLiYl+XRVGk/Ii3fmUU61eIiKKFTwGLy+XCPffcA71ejy1btsjBxi9+8QtMnToVP/vZz3DzzTejqKhoUNebP38+li5d6uuaKYqdaLDiRIMVxhgdppdmqr0cIiIKEZ+2hD7//HMcP34ct912W4/MSEpKCn72s5/B4XBg7dq1Sq+RSLbJux00pTgdSXFGlVdDRESh4lOGRaotmTt37gXPzZs3DwCwefPmQV9v3759aGxshMvlQnFxMa6++mpkZGT4siSKMtL8lTkj2R1ERBRNfApYKioqAABlZWUXPJeTk4PExET5NYPx+9//vsf/x8fHY+XKlVi+fPmAn2u322G32+X/b21tHfTXpfDU4XBhZ2UTALYzExFFG5+2hCwWCwDPFlBvkpOT5df0Z9iwYXjxxRdx9OhRdHR04PTp03j99deRnp6Oxx57DC+++OKA11i1ahVSUlLkj8JCHiYX6bYda4TDLaAwPR7Ds8xqL4eIiEJIJ4qiONgXz507F59++ikqKipQWlp6wfP5+flob28fVNDSm2+++QaTJ09GQkICzp07B4Oh7wRQbxmWwsJCWCwWJCdzNkckuvO1L7HpSD3unF6Mld8Zo/ZyiIhIAa2trUhJSRnw57dPGRYps9JXQCJ9UX+NGTMGM2bMQFNTEw4dOtTva00mE5KTk3t8UOQ6VteGTUfqodMBd1xWrPZyiIgoxHwKWKTald7qVM6ePYv29vZe61t8kZnpaVW1Wq0BXYciy/9tPQEAuGb0EBRncjuIiCja+BSwzJo1CwCwcePGC57bsGFDj9f4w+12Y/fu3QAw6FkuFPka2u14f28NAOCemSUqr4aIiNTgU8By1VVXoaSkBOvXr8f+/fvlxy0WC379618jNjYWt99+u/x4bW0tDh8+fMEW0p49ey64ttvtxmOPPYZjx45hzpw5yM3N9fGtUKR6Y3sVHC4BFxekYHJRmtrLISIiFfjU1mwwGLB69WrMmzcPM2fO7DGav6qqCs8++2yPsforVqzA2rVr8dprr/WYaDt58mSMHz8e48ePR35+PpqamrB582YcPXoUBQUFWL16tVLvj8KczenGGzuqAAB3X1HCs4OIiKKUz2cJzZkzB1u3bsXKlSvx9ttvw+l0Yty4cXjmmWewcOHCQV3jkUcewY4dO/Dpp5+iqakJsbGxKC0txc9//nM8/PDDSEvjb9Hk8Zd9NWiyOpCfGo/rxuaovRwiIlKJT23NWjbYtigKH4Ig4prnNuN4vRU/v2E07r6C9StERJEmKG3NRKFUfrQOx+utSDIZsHAKBwMSEUUzBiykSaIo4pXNlQCAW6cW8qBDIqIox4CFNOnjr2ux80QTjDE6LJ0+TO3lEBGRyhiwkOY0WR14/K/fAAAemFOG/NR4lVdERERqY8BCmvPkx9+gyerAyCFJWDZ7uNrLISIiDWDAQpry+eFz+HD/Geh1wDM3j0esgX9FiYiIAQtpSJvNif/3l4MAgB/OGIYJhanqLoiIiDSDAQtpxjP/PIxaiw1D0xPw8DUj1V4OERFpCAMW0oStFQ1Yt+MUAODpBeMQHxuj8oqIiEhLGLCQ6g6fbcWydZ4DMRdNLcTlwzNVXhEREWkNAxZS1ZmWTiz90y602V2YWpyOld8Zo/aSiIhIgxiwkGosHU4sfe1LnG21oSw7Ea/ePhlxRm4FERHRhRiwkCpsTjfueWM3jp5rx5BkE9bcNRUpCRy/T0REvWPAQiFnc7rxn2/tw5cnmpBkMmDNnVM5zZaIiPplUHsBFF2arA7c+/pu7K5qhjFGh1duvwSjc/s+TpyIiAhgwEIhdLy+HXet2YWqxg4kxRnw8uJL2BFERESDwoCFQmL78Ub8aN0eWDqdKEyPx2tLp6A0O0ntZRERUZhgwEJB5RZE/GnrCfzvhsNwukVMGpqKP94+GZmJJrWXRkREYYQBCwXNiQYr/vvdr7C7qhkAcMP4XPz2lovZukxERD5jwEKKEwQRa7efxDP/PAybU0CiyYD/uXE0fjC5EDqdTu3lERFRGGLAQorad6oZT/39kJxVmV6agWcWjEdBWoLKKyMionDGgIUUcbLBit9sOIK/H6gFACTExmDF9aOx+NKhzKoQEVHAGLBQQM612vCH8uNYt6MKLkGETgfcPKkAD88dgdwUDoMjIiJlMGAhvxw524ZXv6jER/tr4HSLAIDZI7Pw2HWjMCqHg+CIiEhZDFho0ARBxL+PN2L11kqUH6mXH59anI7/vLoM00s5BI6IiIKDAQsNqNbSifd2n8Y7e6pR3dQJANDpgGvH5ODemSWYODRN5RUSEVGkY8BCvbLaXfjs0Dl8tP8Myo/UQfDs+iDJZMD3JuXjhzOGoSjDrO4iiYgoajBgIZnN6camw3X4+Osz+PxwHWxOQX5u6rB03DqlENeNzUV8LAe/ERFRaDFgiXIN7XZ8fqgOnx46hy8q6nsEKcUZCbhxfB6+PykfJVmJKq6SiIiiHQOWKCMIIg6esWDzkXqUH63H3lPNEMWu5/NT43Hjxbn4zvg8jMlL5gwVIiLSBAYsUeCsxYZtxxqwpaIeX1Q0oMnq6PH8+IIUXD16CK4ePQSjc5MYpBARkeYwYIlArTYndlY2YduxBmw91oBjde09nk8yGXB5aQZmjcjGlaOykZMSp9JKiYiIBocBSwSwOd3YU9WMbccasO14Iw6cbpG7egBPC/K4/BTMKM3ErBFZmFSUBmOMXr0FExER+YgBSxhyCyIO1liw9VgDth1rwO6qZjhcQo/XFGckYEZZJmaUZmJaSQZSE2JVWi0REVHgGLCEiVpLJ7Ycrcfmo/XYWtGAVpurx/NDkk2YPjwTl5dm4rLhGchP5Tk+REQUORiwaJRbELH3VDM+/fYcyo/U4ei58+pQ4gy4rCQDM8oycfnwTAzPMrNYloiIIhYDFg2xOd3YcrQen357Dp8frkNjt24evQ64uDAVs0ZkYeaILIzPT4GBdShERBQlGLCozC2I2FHZiA/31eCfB8+izd611ZMSb8SVozydPFeUZbIOhYiIohYDFpVUNVqxfucpfLi/Buda7fLjeSlxmDc2B9dcNARTitPZzUNERAQGLCElCCI2H63H69tPovxovTxhNjXBiOvH5WL+hHxMLkqDXs9aFCIiou4YsISAwyXgnd3VePWLSlQ1dsiPzx6ZhdumDsXskdmINTCTQkRE1BcGLEHkdAt4b89pvPT5MdS0dAIAkuMM+MHkQiyeVoTiTLPKKyQiIgoPDFiCQBBEfLCvBi/86yiqmzyBSnaSCffPHo6FU4YiPjZG5RUSERGFFwYsCjtU24r/95cD2HuqBQCQmWjCstnD8R+XDkWckYEKERGRPxiwKMRqd+H5z47iT9tOwi2IMMfG4IEry7D08mJmVIiIiALEgEUB5UfqsOKDA6i12AAA143NwS++cxFyUzgen4iISAkMWALgcAl4duMR/HFLJQCgMD0eT353LOaMylZ5ZURERJGFAYufTjV24Cd/3ouvTlsAALdfVoQV143m9g8REVEQMGDxw8dfncGKDw6g3e5CSrwR/3vzeMwbk6P2soiIiCIWAxYfCIKIZzcewf9XfhwAMKU4Dc/fOhH5qaxVISIiCiYGLIPU6XDjkXf34x8HzgIAls0ejkeuGcETk4mIiELA75+2u3btwvXXX4/U1FSYzWZMmzYN77zzjk/XsNvtePLJJ1FWVoa4uDjk5eXh3nvvRV1dnb/LCoq6VhsW/nE7/nHgLGJj9PjdDy7G8mtHMVghIiIKEb8yLJs2bcK8efMQFxeHW2+9FUlJSXj//fexcOFCVFdX45FHHhnwGoIg4KabbsKGDRswbdo0LFiwABUVFVi9ejX+9a9/YceOHcjKyvJneYr69kwrfrh2F2otNqQlGPHKksmYOixd7WURERFFFZ0oSmcGD47L5cKoUaNw+vRp7NixAxMmTAAAWCwWTJ06FSdPnsTRo0dRVFTU73Vee+013HXXXVi0aBHefPNN6HSeE4pffvllLFu2DPfeey9eeeWVQa+rtbUVKSkpsFgsSE5O9uUt9UkQRFz/+y9w+GwbhmeZ8aelU1CUwfN/iIiIlDLYn98+72l8/vnnOH78OG677TY5WAGAlJQU/OxnP4PD4cDatWsHvM6rr74KAFi1apUcrADAfffdh5KSErz55pvo7Oz0dXmK0ut1+P2iibh2TA4+uH86gxUiIiKV+BywlJeXAwDmzp17wXPz5s0DAGzevLnfa9hsNuzcuRMjR468IBOj0+lwzTXXwGq1Yvfu3b4uT3EjhiTh5SWXICXeqPZSiIiIopbPAUtFRQUAoKys7ILncnJykJiYKL+mL8ePH4cgCL1eo/u1B7oOERERRQefi24tFs9k15SUlF6fT05Oll8TyDW6v643drsddrtd/v/W1tZ+vyYRERGFr7Dty121ahVSUlLkj8LCQrWXREREREHic8AiZUX6yn5I1b6BXqP763qzYsUKWCwW+aO6unrAtRMREVF48jlg6a++5OzZs2hvb++zNkVSUlICvV7fZ41Kf3UyEpPJhOTk5B4fREREFJl8DlhmzZoFANi4ceMFz23YsKHHa/oSHx+PqVOn4siRI6iqqurxnCiK+PTTT2E2mzF58mRfl0dEREQRyOeA5aqrrkJJSQnWr1+P/fv3y49bLBb8+te/RmxsLG6//Xb58draWhw+fPiC7Z97770XgGdrp/vsuldeeQWVlZX4j//4D8TH81BBIiIi8mPSLdD3aP6qqio8++yzPUbzL126FGvXrsVrr72GpUuXyo8LgoDrr79eHs0/a9YsHDt2DB988AGKi4uxc+dOn0bzB2PSLREREQVX0CbdAsCcOXOwdetWTJ8+HW+//Tb+8Ic/YMiQIXjrrbcGdY4QAOj1enz00Ud4/PHHUV9fj+eeew7btm3DD3/4Q2zfvl0T5wgRERGRNviVYdEiZliIiIjCT1AzLEREREShxICFiIiINI8BCxEREWmez2cJaZVUisMzhYiIiMKH9HN7oJLaiAlY2traAIBnChEREYWhtra2fo/kiZguIUEQcObMGSQlJUGn0yl23dbWVhQWFqK6uprdR0HGex0avM+hwfscGrzPoRHM+yyKItra2pCXlwe9vu9KlYjJsOj1ehQUFATt+jyvKHR4r0OD9zk0eJ9Dg/c5NIJ1nwc6NBlg0S0RERGFAQYsREREpHkMWAZgMpmwcuVKmEwmtZcS8XivQ4P3OTR4n0OD9zk0tHCfI6boloiIiCIXMyxERESkeQxYiIiISPMYsBAREZHmMWAhIiIizWPA0o9du3bh+uuvR2pqKsxmM6ZNm4Z33nlH7WWFnZqaGjz//POYO3cuhg4ditjYWOTk5GDBggXYuXNnr5/T2tqKhx9+GEVFRTCZTCguLsZ///d/o729PcSrD2/PPPMMdDoddDodduzYccHzvM+B+ctf/oJrrrkGGRkZiIuLw7Bhw7Bo0SJUV1f3eB3vs39EUcQHH3yAOXPmIDc3FwkJCRg5ciTuu+8+VFZWXvB63uf+rVu3Dvfddx8mT54Mk8kEnU6HNWvW9Pl6X++nIAh48cUXMW7cOMTHxyMrKwuLFi3q9c/KLyL16vPPPxeNRqOYlJQk3nPPPeLDDz8sFhUViQDEZ599Vu3lhZXly5eLAMThw4eLP/zhD8XHHntMXLBggRgTEyPq9Xrxrbfe6vH69vZ2ccKECSIAce7cueLy5cvFuXPnigDEKVOmiJ2dnSq9k/By4MAB0WQyiWazWQQgbt++vcfzvM/+EwRBvPfee+W/1/fff7+4fPlyccmSJeLQoUPFL774Qn4t77P/Hn74YRGAmJubK/7oRz8SH330UXHevHmiTqcTk5KSxAMHDsiv5X0emPQzLDMzU/7v1157rdfX+nM/7777bhGAOGbMGPHRRx8VFy9eLMbGxorp6eni0aNHA14/A5ZeOJ1Ocfjw4aLJZBL37dsnP97S0iKOGDFCjI2NFU+ePKneAsPM+++/L5aXl1/w+JYtW0Sj0SimpaWJNptNfvwXv/iFCEBcvnx5j9dLgc+vf/3roK853DkcDnHSpEnipZdeKi5evLjXgIX32X/PP/+8CEC8//77RZfLdcHzTqdT/m/eZ//U1taKer1eLCoqEltaWno897vf/U4EIN55553yY7zPA/v000/ln12rVq3qN2Dx9X5+/vnnIgBx5syZot1ulx//xz/+IQc9gWLA0osNGzZc8I9BsmbNGhGA+MQTT6iwssgjRey7du0SRdHzm2teXp6YmJgotre393hte3u7mJiYKJaUlKix1LCycuVK0WQyid988414xx13XBCw8D77r6OjQ0xLSxNLSkp6BCa94X323/bt20UA4m233XbBc0ePHhUBiDfeeKMoirzP/ugvYPHnfi5atEgEIG7evPmC682ePVsEIFZVVQW0Ztaw9KK8vBwAMHfu3AuemzdvHgBg8+bNoVxSxDIajQAAg8FzDmdFRQXOnDmD6dOnw2w293it2WzG9OnTUVlZeUGNAHXZu3cvnnrqKaxcuRIXXXRRr6/hffbfxo0b0dzcjPnz58PtduODDz7A008/jZdffhnHjh3r8VreZ/+VlZUhNjYW27ZtQ2tra4/n/va3vwEArrrqKgC8z0rz536Wl5fLz51PqZ+bDFh6UVFRAcDzD+Z8OTk5SExMlF9D/jt16hQ+++wz5ObmYty4cQD6v/fdH+f9753dbsftt9+OCRMm4NFHH+3zdbzP/tuzZw8AICYmBuPHj8eCBQuwYsUKLFu2DCNHjsR//dd/ya/lffZfRkYGnn76aZw6dQqjRo3CsmXLsHz5clx77bVYvnw57r//fjzwwAMAeJ+V5uv9tFqtqK2txbBhwxATEzPg6/1lCOizI5TFYgHQ93HXycnJ8mvIP06nE0uWLIHdbsczzzwj/yUfzL3v/jrq6Re/+AUqKiqwZ8+eXr9xSHif/VdXVwcA+N3vfodJkybhyy+/xOjRo7Fv3z7ce++9+O1vf4vhw4dj2bJlvM8B+ulPf4r8/HzcfffdePnll+XHZ8yYgdtuu03OzPI+K8vX+xmq+88MC4WcIAhYunQptmzZgnvuuQdLlixRe0kRYfv27Xj22Wfx85//HGPHjlV7ORFLEAQAQGxsLD788ENMmTIFiYmJuOKKK/Duu+9Cr9fjt7/9rcqrjAxPPvkkFi9ejJ/97Georq5GW1sbvvjiC9hsNsyePRt//etf1V4ihRADll5IUWJf0WBra2ufkST1TxAE3HXXXVi/fj0WL17c47cmYHD3vvvryMPlcuGOO+7A+PHj8dhjjw34et5n/0n3ZPLkycjLy+vx3NixY1FSUoLjx4+jpaWF9zkAn332GVauXIkHHngAjz32GAoKCpCYmIgZM2bg448/htFoxCOPPAKAf5+V5uv9DNX955ZQL7rvt11yySU9njt79iza29sxdepUNZYW1gRBwJ133onXX38dixYtwpo1a6DX94yZB9rrHGhvNVq1t7fL9yY2NrbX11x22WUAPMPOpGJc3mffjRw5EgCQmpra6/PS452dnfz7HIBPPvkEADBnzpwLnsvJycGoUaOwb98+tLe38z4rzNf7aTabkZubixMnTsDtdl+wHa3U/WfA0otZs2Zh1apV2LhxI2699dYez23YsEF+DQ1e92Bl4cKFeOONN/oszsrLy8O2bdtgtVp7VKhbrVZs27YNw4YNQ2FhYSiXr3kmkwk//OEPe31uy5YtqKiowHe/+11kZWWhuLiY9zkA0g/QQ4cOXfCc0+nEsWPHYDabkZWVhZycHN5nPzkcDgBAfX19r8/X19dDr9fDaDTy77PC/Lmfs2bNwltvvYVt27Zh5syZPa4n/dw8/3GfBdQUHaGcTqdYUlLS7+C4EydOqLa+cON2u+VZILfccsuAsys4AEpZvc1hEUXe50BI84NeffXVHo8/+eSTIgBx8eLF8mO8z/7585//LE9NPX9w3B/+8AcRgDh9+nT5Md5n34Tj4DidKIpiYCFPZNq0aRPmzZuHuLg43HrrrUhKSsL777+PqqoqPPvss/LeKQ3s8ccfxxNPPIHExET853/+p1zZ3938+fMxYcIEAJ4Ifvr06fjqq68wd+5cTJo0CXv37sXGjRsxZcoUbN68GfHx8SF+F+Fr6dKlWLt2LbZv345p06bJj/M+++/48eO4/PLLUVdXhxtuuEHenvj8889RVFSEHTt2ICcnBwDvs7/cbjeuvPJKbNmyBdnZ2fjud7+L1NRU7N27F59//jni4+NRXl4ub8/zPg9s9erV2Lp1KwDgwIED2Lt3L6ZPn47S0lIAnu6ru+++G4B/9/Oee+7B6tWrMWbMGNxwww2ora3F22+/jcTERGzfvh0jRowI7A0EHPJEsJ07d4rXXnutmJycLMbHx4tTp0694NwbGpj0G35/H+dH+S0tLeJDDz0kFhYWikajURw6dKj4yCOPiK2treq8iTDWV4ZFFHmfA3Hq1Clx6dKlYk5Ojmg0GsXCwkLxxz/+sXju3LkLXsv77B+bzSauWrVKnDhxopiQkCAaDAYxPz9fXLx4sfjtt99e8Hre5/4N9L34jjvu6PF6X++n2+0WX3jhBXHMmDGiyWQSMzIyxIULF4rHjh1TZP3MsBAREZHmsa2ZiIiINI8BCxEREWkeAxYiIiLSPAYsREREpHkMWIiIiEjzGLAQERGR5jFgISIiIs1jwEJERESax4CFiIiINI8BCxEREWkeAxYiIiLSPAYsREREpHkMWIiIiEjz/v8huo20ue4E5wAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.plot(data['rv_d'][:, 0,0,0,0])" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "language_info": { "name": "python" } }, "nbformat": 4, "nbformat_minor": 2 }