{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Tutorial" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "/=================\\\n", "| rayXpanda |\n", "|-----------------|\n", "| Version: 0.1 |\n", "\\=================/\n", "\n" ] } ], "source": [ "import rayXpanda\n", "import numpy as np" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We call the functions like so:" ] }, { "cell_type": "code", "execution_count": 14, "metadata": { "collapsed": true }, "outputs": [], "source": [ "u = 0.5\n", "cos_alpha = 0.5" ] }, { "cell_type": "code", "execution_count": 16, "metadata": { "collapsed": false }, "outputs": [], "source": [ "cos_psi = rayXpanda.deflection.deflect(cos_alpha, u)[0]" ] }, { "cell_type": "code", "execution_count": 17, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "(0.5000000000000001, 1.0175735797883876)" ] }, "execution_count": 17, "metadata": {}, "output_type": "execute_result" } ], "source": [ "rayXpanda.inversion.invert(cos_psi, u)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Let's time the vectorised versions:" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "collapsed": false }, "outputs": [], "source": [ "cos_alpha = np.linspace(-0.05, 1.0, 1000)" ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "100 loops, best of 100: 128 µs per loop\n" ] } ], "source": [ "%%timeit -n100 -r100\n", "\n", "rayXpanda.deflection.deflect_vec(cos_alpha, 0.5)" ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "collapsed": true }, "outputs": [], "source": [ "cos_psi = np.linspace(-0.5, 1.0, 1000)" ] }, { "cell_type": "code", "execution_count": 9, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "100 loops, best of 100: 707 µs per loop\n" ] } ], "source": [ "%%timeit -n100 -r100\n", "\n", "rayXpanda.inversion.invert_vec(cos_psi, 0.5)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We can examine the performance of [invert()](rayXpandaAPI.rst#rayXpanda.inversion.invert): " ] }, { "cell_type": "code", "execution_count": 76, "metadata": { "collapsed": true }, "outputs": [], "source": [ "from matplotlib import pyplot as plt\n", "from matplotlib import rcParams\n", "from matplotlib.ticker import MultipleLocator, LogLocator, NullFormatter\n", "from matplotlib import gridspec\n", "\n", "rcParams['text.usetex'] = False\n", "rcParams['font.serif'] = 'Computer Modern Roman'\n", "rcParams['font.family'] = 'serif'\n", "rcParams['font.size'] = 16.0" ] }, { "cell_type": "code", "execution_count": 82, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAv0AAALXCAYAAAAe8f9EAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzs3XlclFXbwPHfEbVcwyU1wQQjkU1BwiXFXBLUjMIyLNcKtbLeyh7fLM14itTSXBLNtLJFn9wqNUVzL7SURHFBnyIDQnBJSSWgVOa8fwzM68g2LMMwcH0/n/ngnPucc1/39Dw3F/ecRWmtEUIIIYQQQlRftWwdgBBCCCGEEMK6JOkXQgghhBCimpOkXwghhBBCiGpOkn4hhBBCCCGqOUn6hRBCCCGEqOYk6RdCCCGEEKKak6RfCCGEEEKIak6SfiGEEHZPKfWaUuoXpZRBKfWgreMRQoiqxm6TfqWUs1JqgVLqR6VUtlJKK6VcLGxbSyn1ilIqWSn1t1LqsFLqIetGLIQQwoq2AQOA720diBBCVEV2m/QDbsAjwJ9ATCnbvglEAFHAQGAfsEYpNagiAxRCCFG40jy4UUq1UUqtVUpdUkpdVkp9pZS6/fo6Wut9WuvfKiN2IYSwR/ac9H+vtW6ptR4ErLG0kVKqBfAvYKbWerbWepfWejywC5hppViFEEKYs+jBjVKqPrAT6ACMBkYCdwK7lFINKiFOIYSoFuw26ddaG8rYNBioCyy/oXw54KOUci1XYEIIISxh6YObsUA74EGt9Tqt9XogBGgLjK+EOIUQolqw26S/HLyAf4BfbyhPyPvpWbnhCCFEzVOKBzchwD6ttemerbVOAvYCD1gjNiGEqI5q2zoAG2gKXNRa6xvKM647bkYpFQG8buW4hBDCqrTWytYxlIEXsL6Q8gRgaGk7k/u5EKI6KMv9vCYm/aWmtY7AOPG3UEqpQv6GEEII2zNcvox76+/5Net+W4dSVk0xjvu/UQbQJP9NXjIfDtwKeCulooBuWutT1zeS+7kQwt4pVbbnNzVxeM+fgKMq+InlP+HPQAghqol1T7+GQzN7fMBfOlrrCK21s9b6Jq1187x/nyq5pRBC1Aw1MelPAG4C7rihPH8s//HKDUcIIaxn7nYYPe72kitWXX9y3RP96xT1DYAQQohC1MSkfwtwFRh+Q/kI4FjeBDEhhLB/KSkcyxjI8JE+to6kPBIwjuu/kSfykEYIISxm10m/UuphpdTDgH9e0cC8snuuq3NNKfVR/nut9TlgDvCKUmqiUqq3Uup9oC/wSmXGL4QQ1mT497+p18iZ2+36QT8bgG5KqXb5BXmbePXIOyaEEMIC9j6R98a1nRfl/fwO6J33b4e81/WmAH8BzwOtgJ+BR7TWG60TphBCVL74n3/GsVV9W4dRpLyHNmD+4OYP4A+t9Xd5ZUuBZ4H1SqmpgMa4q3oq8EFlxiuEEPZMySoF5SerPQghqpzffiO4z3Pccd8nLFp0K0qpKrdkp1KqqBvnd1rr3tfVux2YC/QHFLADeEFrnWyNmOR+LoSoysp6P7f3J/1CCCEK8+mnnOQuIkbeautIimTpLy2t9e/AQ1YORwghqjW7HtMvhBCicFc2b+aPiw/i52frSIQQQlQFkvQLIUR1s38/P9SrR6tWrbn5ZlsHI4QQoiqQ4T1CCFHdLFnCmtvu4NbcxraORAghRBUhT/rLQCkVoZTS+S9bx1Ocu+++u0ad11JbtmzB3d0dNzc3Zs6cWWS9ixcv8vDDD9OhQwc8PDz48ccfSU1NpU+fPnh6euLl5cX8+fPN2ri4uODj44Ovry933XWXtS9FCHNaQ3Y2+39tyTPP1LN1NEIIIaoIWb2nAtTU1R601mitqVXLvv52zM3NpX379mzbtg1nZ2cCAgL44osv8PT0LFB39OjRBAYGEh4ezpUrV8jOziYnJ4fTp0/TuXNnMjMz8ff3Z926dab2Li4uHDhwgObNm1f2pQkBGzdyPjmZbvNuZuvWcNrlrW5fFVfvqYpq6v1cCGE/yno/t69sTZRaw4YNSU5OxsPDg7Fjx+Ll5UVQUBA5OTkATJ48mYULF5rqR0REMHv2bACWL19Oly5d8PX1Zfz48eTm5pKcnIy7uzujRo3C29ub1NRU7rvvPjp16oS3tzerVq0ynTffnDlz8Pb2xtvbm3nz5gEUG1NxunfvTlKScdPktLQ0/P39S2hRUGxsLG5ubrRr1466desybNgw1q9fX6DepUuX+P7773nyyScBqFu3Lo6Ojtx222107twZgEaNGuHh4UFaWlqp4xDCKmbNYunZsxgMA0wJvxBCCCFJfw2RmJjIhAkTSEhIwNHRkS+//BKAsLAwVq9ebaq3evVqwsLCOHHiBKtWrWLv3r3Ex8fj4ODAihUrTH0988wzJCQkcODAAVq3bs3hw4c5duwYAwYMMDtvXFwcy5YtY//+/ezbt4+lS5dy6NChYmMaNGgQ6enpBa7BYDCQkpKCi4sLAEeOHKFjx45mdQIDA/H19S3w2r59u6lOWloabdq0Mb13dnYuNGlPSkri1ltv5fHHH8fPz4/w8HCysrLM6iQnJ3Po0CG6du1qKlNKERQUhL+/P0uWLCniv4gQVpCdDdeucTYzk1atWtk6GrtgT8M1hRCiPGQib2X55BNITq64/lxcYMwYi6u7urri6+sLgL+/P8l5sfj5+XHu3DnS09P5448/aNKkCW3atCEqKoq4uDgCAgIAyMnJoUWLFvTq1Yu2bdvSrVs3AHx8fHjppZd4+eWXGTx4MIGBgWbn3bNnD6GhoTRo0ACAIUOGEBMTQ0hISJExRUdHF3oNJ0+exNXVFaWM32gdOXIEHx8fszoxMTEWfyYluXbtGgcPHmTBggV07dqV559/npkzZ/Lmm28C8Ndff/HQQw8xb948Gjf+/wmTe/bswcnJiXPnztG/f386dOhAr169KiwuIYr06acwahSJG36jdWu5vVtCax0BROS/l8RfCFFdyW+FylKKBN0abrrpJtO/HRwczIbSDB06lLVr13LmzBnCwsIA43j90aNHM2PGDLN+kpOTTQk8QPv27Tl48CDR0dFMnTqVfv36MW3atHLHVJijR4+aJfkHDhxg3LhxZnUCAwPJzMws0Hb27Nnce++9ADg5OZGammo6durUKZycnAq0cXZ2xtnZ2fQU/+GHHzZN+r169SoPPfQQw4cPZ8iQIWbt8vtq0aIFoaGhxMbGStIvKkdMDCemTOHv1bWZNMnWwQghhKhKZHiPICwsjJUrV7J27VqGDh0KQL9+/Vi7di3nzp0DICMjg5SUlAJt09PTqV+/PiNGjGDSpEkcPHjQ7HhgYCDr1q0jOzubrKwsvv766wLfBlgqIyMDR0dHAE6cOMGmTZsKDO+JiYkhPj6+wCs/4QcICAggMTGRpKQkrly5wsqVKwkJCSlwvlatWtGmTRt+/vlnAHbs2IGnpydaa5588kk8PDyYOHGiWZusrCzTHx1ZWVls3boVb2/vMl2vEKWSkgJOTix8/31q1+6NLBwlhBDievKkX+Dl5UVmZiZOTk7cdtttAHh6ehIZGUlQUBAGg4E6deqwcOHCAuOEjx49yqRJk6hVqxZ16tTh/fffNzveuXNnxowZQ5cuXQAIDw/Hz8/PNJSnMIMGDeLDDz+kdevWZuXBwcFERUWRmpqKu7s7zZo1o2XLlqW+3tq1axMVFUVwcDC5ubk88cQTeHl5FXr+BQsWMHz4cK5cuUK7du1YtmwZe/fu5fPPPzctywkwffp0Bg0axNmzZwkNDQWMw4Mee+yxAvMchLCK116D0aM5NG0aDRu247q59EIIIYQs2VkRZIk3IYTNde2K3rePESPCSUv7iN27zQ/Lkp2Wkfu5EKKqkyU7hRCipvrpJ+jWjb1795Kb+wB58++FEEIIE0n6hRDC3i1cCNOn8/nnn3PHHX144glbBySEEKKqkaS/DGRdZyFElXH5MtSuDQ0acOzYMVJSGuHhYeughBBCVDWS9JeB1jpCa63yX7aORwhRg0VEQM+enD9/nnr16pGdbeuAhBBCVEWS9AshhD3bswdGjWLTpk0EBk7kuu0vhBBCCBNJ+oUQwl79+iv4+UGtWuzYsYPz5wfwyCO2DkoIIURVJOv0CyGEvXrnHZg+Ha01Wmv++qsW999v66CEEEJURfKkXwgh7NGVK5CdDc2bs3PnTmrXrs21a8Y5vUIIIcSNJOmvQSIiIpg9e3axdd577z08PDwYPnw4AA3LsK3nxYsXWbRokVnZ3XffXep+KtOWLVtwd3fHzc2NmTNnFlnv4sWLPPzww3To0AEPDw9+/PFHUlNT6dOnD56ennh5eTF//nyzNi4uLqbde++66y5rX4qoKaZPh7z/PX3wwQc8/PALNGli45iEEEJUWfJMSJhZtGgR27dvx9nZucx95Cf9zzzzjKnshx9+qIjwrCI3N5cJEyawbds2nJ2dCQgIICQkBE9PzwJ1n3/+eQYMGMDatWu5cuUK2dnZ5OTk8O6779K5c2cyMzPx9/enf//+Zu137dpF8+bNK/OyRHUXHQ379gFw7do1Dh3qRIcONo5JCCFElSVP+qu5t956i/bt29OzZ09+/vlnU/ny5cvp0qULvr6+jB8/ntzcXJ566il+++03Bg4cyNy5c836Kaw+wGeffUbHjh3p1KkTI0eOBGDy5MmcPHkSX19fJk2aBJh/YzBnzhy8vb3x9vZm3rx5ACQnJ+Ph4cHYsWPx8vIiKCiInJycEq+ve/fuJCUlAZCWloa/v3+pP6PY2Fjc3Nxo164ddevWZdiwYaxfv75AvUuXLvH999/z5JNPAlC3bl0cHR257bbb6Ny5MwCNGjXCw8ODtLS0UschhMUSE6FrV6hVi0uXLtG0aVN274YRI2wdmP2RfVeEEDWFJP3VWFxcHCtXriQ+Pp7o6Gh++uknAE6cOMGqVavYu3cv8fHxODg4sGLFChYvXkzr1q3ZtWsXL774oqmfouonJCQQGRnJzp07OXz4sGlYy8yZM7njjjuIj49n1qxZBWJatmwZ+/fvZ9++fSxdupRDhw4BkJiYyIQJE0hISMDR0ZEvv/wSgEGDBpGenl7g+gwGAykpKbi4uABw5MgROnbsaFYnMDAQX1/fAq/t27eb6qSlpdGmTRvTe2dn50KT9qSkJG699VYef/xx/Pz8CA8PJysry6xOcnIyhw4domvXrqYypRRBQUH4+/uzZMmSIv5rCVEKM2fCv/8NwMqVK3n00Udp0gQaN7ZxXHZI9l0RQtQUMrynsuzebXwBtG8PrVsbf16fBI4bB5mZsHGj8SfAo4/C6dPG93FxxjJ/f+jdGxo1KvaUMTExhIaGUr9+fQBCQkIA2LFjB3FxcQQEBACQk5NDixYtiuynqPqXLl1i6NChpmErTZs2LfFj2LNnD6GhoTRo0ACAIUOGEBMTQ0hICK6urvj6+uZdoj/JyckAREdHF9rXyZMncXV1RSnj7+kjR47g4+NT4DOoKNeuXePgwYMsWLCArl278vzzzzNz5kzefPNNAP766y8eeugh5s2bR+Prsq89e/bg5OTEuXPn6N+/Px06dKBXr14VFpeoYS5ehJwcyPv/2/Lly9m4cUdJtwMhhBA1nCT9ZaCUigBeL1Wj3r2NrxtFRBQsc3cv/H0FrcWntWb06NHMmDGjXPUXLFhQIfHku+m6XYUcHBxKHN5z9OhRsyT/wIEDjBs3zqxOYGAgmfl/QF1n9uzZ3HvvvQA4OTmRmppqOnbq1CmcnJwKtHF2dsbZ2dn0FP/hhx82Tfq9evUqDz30EMOHD2fIkCFm7fL7atGiBaGhocTGxkrSL8ru5Zch739jBoMBg8HA7t11CQ62cVxCCCGqNBneUwb28nVwr169WLduHTk5OWRmZvLNN98A0K9fP9auXcu5c+cAyMjIICUlpch+iqrft29f1qxZw4ULF0zlYBzXXliiDcYkfN26dWRnZ5OVlcXXX39NYGBgma4vIyMDR0dHwDgEadOmTQWG98TExBAfH1/glZ/wAwQEBJCYmEhSUhJXrlxh5cqVpm9FrteqVSvatGljmhuxY8cOPD090Vrz5JNP4uHhwcSJE83aZGVlmT6LrKwstm7dire3d5muVwgADh6Ehx8G4Pvvv6dnz558+CF06WLjuIQQQlRpkvRXY507dyYsLIxOnToxcOBA0/AcT09PIiMjCQoKomPHjvTv35/Tp08X2U9R9b28vJgyZQr33HMPnTp1MiW8zZo1o0ePHnh7e5sm8l4f05gxY+jSpQtdu3YlPDwcPz+/Yq+jqDH9wcHBbNmyheHDh7NmzRqaNWtGy5YtS/sxUbt2baKioggODsbDw4NHHnkELy+vQs+/YMEChg8fTseOHYmPj+fVV19l7969fP755+zcudM0ZyB/SNLZs2fp2bMnnTp1okuXLtx3330MGDCg1DEKARhX63nwQdPblStX8vLLL/Pnn+DqasO4hBBCVHlKa1msoLyUUlo+RyGE1Y0YAR9+CDffDMBTTz3FokWLuf9+2LSp+KZKKaryN5NVhdzPhRBVXVnv5/KkXwgh7MEvvxi3281L+BMSEjhz5gy7diHj+YUQQpRIkn4hhLAHEyfCCy+Y3r777rs888wz7N8PhUxBEUIIIczYbdKvlGqjlFqrlLqklLqslPpKKXW7hW1vV0p9qpT6XSmVo5T6RSkVqZRqYO24hRCi1DIzISMD8pa0BeP+EkFBQWzaBHlbVQghhBBFssslO5VS9YGdwD/AaEADkcAupVRHrXVWMW0bANuBOsBrwO9AAPBv4E4gzLrRCyFEKX3xhdnyvpcvX6Zt27YYDKBklL4QQggL2GXSD4wF2gHuWutfAZRSR4BEYDwwp5i2PTAm98Fa6615ZbuUUk2Bfyml6muts60XuhBClILWsHMnrFxpKtq4cSMPP/ww334LJSx+JYQQQgD2O7wnBNiXn/ADaK2TgL3AAyW0rZv38/IN5Rcxfh7y3EwIUXXMmwfdupkVvf/++/Tu3Zsff4TJk20UlxBCCLtir0m/F3CskPIEwLOEttsxfiPwtlLKUynVUCnVF3geWFzc0CAhhKh0//kPPPus6e21a9fQWlO3bl3+/BMK2TxaCCGEKMBek/6mwJ+FlGcATYprqLX+G+iJ8doTgExgB7AReLawNkqpCKWULupVngsRQogiffsttG9vXKozz5o1a7j33nsxGODQIRvGJoQQwq7Ya9JfZkqpm4FVQAtgJHAPMAnjBN6FhbXRWkdorVVRr0oLvgzuvvvuGnVeS23ZsgV3d3fc3NyYOXNmkfVcXFzw8fHB19eXu+66y6L2lvYtRIk2boQFC8yKdu/ezeTJk/nmGxnPL4QQwnL2OpH3Twp/ol/UNwDXexLoDbhprU/mlX2vlLoELFFKLdZaH66wSG3shx9+sFrfWmu01tSqVfBvR2uet7xyc3OZMGEC27Ztw9nZmYCAAEJCQvD0LHxk2K5du2jevLlF7UvbtxBFSk+HOnWgaVNTkdaa2rVrc/PNN7NiBbzyig3jE0IIYVfs9Ul/AsZx/TfyBI6X0NYH+PO6hD9fbN5Pj3LGVqU0bNiQ5ORkPDw8GDt2LF5eXgQFBZGTkwPA5MmTWbjw/7/giIiIYPbs2QAsX76cLl264Ovry/jx48nNzSU5ORl3d3dGjRqFt7c3qamp3HfffXTq1Alvb29WrVplOm++OXPm4O3tjbe3N/PmzQMoNqbidO/enaSkJMC4Trm/v3+pP5PY2Fjc3Nxo164ddevWZdiwYaxfv75C2pe3byFMRo+G554zK9qwYQMGgwGAli3lSb8QQgjL2WvSvwHoppRql1+glHLBuBznhhLangGaKKXcbijvmvczrYJirFISExOZMGECCQkJODo68uWXXwIQFhbG6tWrTfVWr15NWFgYJ06cYNWqVezdu5f4+HgcHBxYsWKFqa9nnnmGhIQEDhw4QOvWrTl8+DDHjh1jwIABZueNi4tj2bJl7N+/n3379rF06VIO5Q1ELiqmQYMGkZ6eXuAaDAYDKSkpuOTtRHTkyBE6duxoVicwMBBfX98Cr+3bt5vqpKWl0aZNG9N7Z2dn0tIK/8+ulCIoKAh/f3+WLFlSYvvS9C1EkX7/Ha5eBVdXs+LFixczefJkDAbjYVF+N87ZsnU8QghhLfY6vGcpxkm365VSUzFuzvUmkAp8kF9JKdUWOAm8obV+I6/4E2AiEK2Uegvj5lx3YdyoKw7jsp8V7pNPPiE5ObnC+nNxcWHMmDEW13d1dcU3bzdPf39/Uyx+fn6cO3eO9PR0/vjjD5o0aUKbNm2IiooiLi6OgIAAAHJycmjRogW9evWibdu2dMtbQtDHx4eXXnqJl19+mcGDBxMYGGh23j179hAaGkqDBsbNjocMGUJMTAwhISFFxhQdHV3oNZw8eRJXV1dU3m5ER44cwcfHx6xOTEyMxZ+JJfbs2YOTkxPnzp2jf//+dOjQoUL7F6JQL70Es2YVKM7JyaFt27asWweNG9sgrmpIax0BROS/l8RfCFFd2WXSr7XOyltmcy7wOca19XcAL2it/7quqgIcuO4bDa11slKqG8abfCTQHOMfC0uAt7TWBmvEXJoE3Rpuuukm078dHBzMhtIMHTqUtWvXcubMGcLCjBsSa60ZPXo0M2bMMOsnOTnZlMADtG/fnoMHDxIdHc3UqVPp168f06ZNK3dMhTl69KhZkn/gwAHGjRtnVicwMJDMzMwCbWfPns29994LgJOTE6mpqaZjp06dwqmIdQ/zy1u0aEFoaCixsbH06NGjyPal6VuIQmVkgKMj5P3BnS82NpbBgwcDsGIFTJlii+CEEELYK7tM+gG01r8DD5VQJ5lCNtvSWh8HHrFOZPYnLCyMsWPHcv78eb777jsA+vXrxwMPPMCLL75IixYtyMjIKDSZTk9Pp2nTpowYMQJHR0c+/PBDs+OBgYGMGTOGyZMno7Xm66+/5vPPPy9TnBkZGTg6OgJw4sQJNm3aRFRUlFkdS570BwQEkJiYSFJSEk5OTqxcuZL//Oc/BeplZWVhMBho1KgRWVlZbN26lWnTphXb3tK+hSjS4sXw9NMFirdu3crw4cMB4xzfvC/JhBBCCIvYbdIvKo6XlxeZmZk4OTlx2223AeDp6UlkZCRBQUEYDAbq1KnDwoULadWqlVnbo0ePMmnSJGrVqkWdOnV4//33zY537tyZMWPG0KVLFwDCw8Px8/MrdqjToEGD+PDDD2ndurVZeXBwMFFRUaSmpuLu7k6zZs1o2bJlqa+3du3aREVFERwcTG5uLk888QReXv8/Lzz//H///TehoaGAcUOkxx57zDRnoaj2JfUtRLGys2HdOnj11QKHNm3axNSpU7l6Fcowf10IIUQNp7SW4YvlpZTS8jkKIcpt3Djo0wcefdSseN++fcydO5dVq1axaxfEx8OLL5aua6UUVX1fkapA7udCiKqurPdze129Rwghqpdr1yAlpUDCDzB37lxeyVuU//33IW+UjxBCCGExSfrLQJZ4E0JUuDffhCImwTdu3Ni00lVaGrRoUZmBCSGEqA4k6S8DrXWE1lrlv2wdjxDCzuXkwC+/QI8eBQ798ssvtGtn3JLkwgWoV6+ygxNCCFEdSNIvhBC29txz8MADhR569dVXCQ4OBiAmBsaOrczAhBBCVBeS9AshhC1lZ8PhwzBsWKGH09LS6Ny5MwDr18P991dmcEIIIaoLSfqFEMKWPvgAXnut0EP79+83De0BuHgR6tevrMCEEEJUJ5L0CyGErVy5Avv2QUhIoYe/+uor5s6dC8B//wuykqQQQoiykqRfCCFs5emnYeLEIg9fvnyZFnlL9SxeDCNGVFZgQgghqhtJ+muQiIgIZs+eXWyd9957Dw8PD4bnLQTesGHDUp/n4sWLLFq0yKzs7rvvLnU/lWnLli24u7vj5ubGzJkzi6zn4uKCj48Pvr6+3HXXXRa1t7RvUcOcPQsnTkDXroUe3rx5M//884/p/cWL8OCDlRWcEEKI6kaSfmFm0aJFbNu2jRUrVpS5j8KS/h9++KG8oVlNbm4uEyZMYPPmzRw/fpwvvviC48ePF1l/165dxMfHc+DAgRLbl7ZvUYM88QTMmVPk4fnz5zMtb91+reHqVahdu7KCE0IIUd1I0l/NvfXWW7Rv356ePXvy888/m8qXL19Oly5d8PX1Zfz48eTm5vLUU0/x22+/MXDgQNM44uLqA3z22Wd07NiRTp06MXLkSAAmT57MyZMn8fX1ZdKkSYD5NwZz5szB29sbb29v5s2bB0BycjIeHh6MHTsWLy8vgoKCyMnJKfH6unfvTlJSEmBc5cTf37/Un1FsbCxubm60a9eOunXrMmzYMNavX18h7cvbt6imfv0VOnSAbt0KPay15vLly7i4uACwZw94eFRifEIIIaodSfqrsbi4OFauXEl8fDzR0dH89NNPAJw4cYJVq1axd+9e4uPjcXBwYMWKFSxevJjWrVuza9cuXnzxRVM/RdVPSEggMjKSnTt3cvjwYebPnw/AzJkzueOOO4iPj2fWrFkFYlq2bBn79+9n3759LF26lEOHDgGQmJjIhAkTSEhIwNHRkS+//BKAQYMGkZ6eXuD6DAYDKSkppsToyJEjdOzY0axOYGAgvr6+BV7bt2831UlLS6NNmzam987OzqSlpRX6mSqlCAoKwt/fnyVLlpTYvjR9ixpk4UJ49tkiD2/ZsoWnn37a9H7RIujfvzICE0IIUV3Jl8WVZPfu3ezevRuA9u3b07p1a9q3b29KHAHGjRtHZmYmGzduJDMzE4BHH32U06dPk5mZSVxcHAD+/v707t2bRo0aFXvOmJgYQkNDqZ+3xl9I3gohO3bsIC4ujoCAAABycnJMkwULU1T9S5cuMXToUJo3bw5A06ZNS/wc9uzZQ2hoKA0aNABgyJAhxMTEEBISgqurK76+vqZrTE5OBiA6OrrQvk6ePImrqytKGTdFPnLkCD4+PgU+g4q0Z88enJycOHfuHP3796dDhw4V2r+oAXbvhqQkcHUtssrWrVt5/fXXTe9TUooc+i+EEEJYRJL+MlBKRQCvl1Tver1796Z3794FyiMiIgrwE4BNAAAgAElEQVSUubu7F/r+/gralUdrzejRo5kxY0a56i9YsKBC4sl30003mf7t4OBQ4vCeo0ePmiX5Bw4cYNy4cWZ1AgMDTX9AXW/27Nnce++9ADg5OZGammo6durUKZycnAo9Z355ixYtCA0NJTY2lh49ehTZvjR9ixri5ZehiD9kAf755x9iY2NxdHTMe1/s3wdCCCGERWR4TxlorSO01ir/Zet4itKrVy/WrVtHTk4OmZmZfPPNNwD069ePtWvXcu7cOQAyMjJISUkpsp+i6vft25c1a9Zw4cIFUzlAo0aNCk20wZiEr1u3juzsbLKysvj6668JDAws0/VlZGSYEqMTJ06wadOmAsN7YmJiiI+PL/DKT/gBAgICSExMJCkpiStXrrBy5UrTtyLXy8rKMl1XVlYWW7duxdvbu9j2lvYtaohvv4VOnaBZsyKrzJ07lwceeMD0ftUqGDy4MoITQghRncmT/mqsc+fOhIWF0alTJ1q0aGEanuPp6UlkZCRBQUEYDAbq1KnDwoULadu2baH9FFW/W7duTJkyhXvuuQcHBwf8/Pz45JNPaNasGT169MDb25uBAweajevv3LkzY8aMoUuXLgCEh4fj5+dnGspTmEGDBvHhhx/SunVrs/Lg4GCioqJITU3F3d2dZs2a0bJly1J/TrVr1yYqKorg4GByc3N54okn8PLyKnD+v//+m9DQUACuXbvGY489xoABAwCKbF9S36IG0dq4++6aNcVW27hxI7t27TK9/+47eP99awdXc5Xlm1shhLBHSssWj+WmlNLyOQohivXmm9CzJ/TpU2SVs2fP8s477/Duu++ayu6+GypixVulFFX5m8mqQu7nQoiqrqz3cxneI4QQ1nbhAmzYUGzCD7Bu3TrTt0kAe/dC3uJUQgghRLnI8B4hhLC2UaPghr0vCvPxxx+zf/9+0/v164td2VMIIYSwmDzpF0IIa0pNhbZtjUN7ivHTTz/h7OxsVpaWZhzeI4QQQpSXPOkXQghrmjwZoqJKrDZjxgwiIyNN7//6C5o0sWZgQgghahJ50i+EENayYAHccUeJ2bvBYMDR0RFPT09T2QcfwA0P/oUQQogykyf9QghhDVlZ8NFHcPBgiVW3b99Onxsm+W7caBzTL4QQQlQESfqFEMIawsNh+nSoVfIXqm+88QbffvutWVm9etC4sbWCE0IIUdPI8J4yUEpFKKV0/svW8QghqpjYWOPYnEGDSqx69uxZateuTYMGDUxlP/0EefvXCSGEEBVCkv4y0FpHaK1V/svW8RTnbhst/WGr81pqy5YtuLu74+bmxsyZM4utm5ubi5+fH4MHD7aofWn6FtVUVBRMnWpR1cjISCZOnGhWtnkzPP64NQITQghRU0nSX839UBFbeRZBa43BYKj085ZXbm4uEyZMYPPmzRw/fpwvvviC48ePF1l//vz5eHh4WNS+tH2Laui116BdO7jlFouqX7p0iZCQELOyLVuMq3wKIYQQFUWS/mquYcOGJCcn4+HhwdixY/Hy8iIoKIicnBwAJk+ezMKFC031IyIimD17NgDLly+nS5cu+Pr6Mn78eHJzc0lOTsbd3Z1Ro0bh7e1Namoq9913H506dcLb25tVq1aZzptvzpw5eHt74+3tzbx58wCKjak43bt3JykpCYC0tDT8/f1L/ZnExsbi5uZGu3btqFu3LsOGDWN9ETMmT506xaZNmwgPD7eofWn6FtXQ2bPGx/QRERZVj42NLfC/4VOnoFEjK8QmhBCiRpOkv4ZITExkwoQJJCQk4OjoyJdffglAWFgYq1evNtVbvXo1YWFhnDhxglWrVrF3717i4+NxcHBgxYoVpr6eeeYZEhISOHDgAK1bt+bw4cMcO3aMAQMGmJ03Li6OZcuWsX//fvbt28fSpUs5dOhQsTENGjSI9PT0AtdgMBhISUnBxcUFgCNHjtCxY0ezOoGBgfj6+hZ4bd++3VQnLS2NNm3amN47OzuTlpZW6Of2wgsv8M4771DrusmYxbUvTd+iGnrjDVi61OLqkydPJiwszKzso49gzJgKjquclFLd8uYybVFKHVFKJSqlflRKfaKUelwpJTsKCCFEFSer91SSTz6B5OSK68/FpXSJgaurK76+vgD4+/uTnBeMn58f586dIz09nT/++IMmTZrQpk0boqKiiIuLIyAgAICcnBxatGhBr169aNu2Ld26dQPAx8eHl156iZdffpnBgwcTGBhodt49e/YQGhpqmqQ4ZMgQYmJiCAkJKTKm6OjoQq/h5MmTuLq6opRxGsWRI0fw8fExqxMTE2P5h1KCjRs30qJFC/z9/dm9e3eF9Suqqd27jevx+/lZVD0zM5N//vmHVq1amZX/9htMm2aF+MpAKTUa+BfgBWQCh4FEIAdoCnQFRgILlVKrgX9rrZNsFK4QQohi2G3Sr5RqA8wF+gMK2A68oLX+3cL2HsAbQB+gAfA7sEhrPd8a8dr6yd1NN91k+reDg4PZUJqhQ4eydu1azpw5Y3rqqLVm9OjRzJgxw6yf5ORks1VG2rdvz8GDB4mOjmbq1Kn069ePaRZmLMXFVJijR4+aJfkHDhxg3LhxZnUCAwPJzMws0Hb27Nnce++9ADg5OZGammo6durUKZycnAq02bt3Lxs2bCA6Opq///6by5cvM2LECCZMmFBke0v7FtVMRoZx4m4p/uh89913eeqpp8zK/v7bOBVAVYHlAZRSR4Bbgc+AUUC81rrAamVKqVuAwcBw4LhSaozWelWlBiuEEKJEdjm8RylVH9gJdABGY3zSdCewSynVoLi2ee3vAvYDNwHhwCDgXcDBWjFXZWFhYaxcuZK1a9cydOhQAPr168fatWs5d+4cABkZGaSkpBRom56eTv369RkxYgSTJk3i4A0bEQUGBrJu3Tqys7PJysri66+/LvBtgKUyMjJwdHQE4MSJE2zatKnA8J6YmBji4+MLvPITfoCAgAASExNJSkriypUrrFy5ssBESoAZM2Zw6tQpkpOTWblyJX379mX58uXFtre0b1HNPPIIzJpVqmz9l19+YeTIkWZlS5dC8+YVHVyZfQS4aq1f1lofKizhB9BaX9Jar9BaDwK6ARcrNUohhBAWsdcn/WOBdoC71vpXMD2VSgTGA3OKaqiUqoXxydUOrXXodYd2WS/cqs3Ly4vMzEycnJy47bbbAPD09CQyMpKgoCAMBgN16tRh4cKFBYYiHD16lEmTJlGrVi3q1KnD+++/b3a8c+fOjBkzhi55i46Hh4fj5+dnGspTmEGDBvHhhx/SunVrs/Lg4GCioqJITU3F3d2dZs2a0bJly1Jfb+3atYmKiiI4OJjc3FyeeOIJvLy8Sjy/Je1L6ltUQ998A61bQ/fuFjf54Ycf6Nq1a4Hyr76Cr7+uyODKrizfemqtD2McAiSEEKKKUUU8vKnSlFI7gJu11j1uKP8OQGt9TzFt+wI7gF5a6woZAK6UKuohmBCiOsvJgREjYM0ai3bezdenTx9WrFhR4A/L3r2NUwOsQSlFWfcVUUoNwTi857cKDqvKkfu5EKKqK+v93C6H92CcVHaskPIEwLOEtj3zft6slNqnlLqqlDqnlHpPKVWvQqMUQlRvzzwDs2eXKuG/cOEC165dK5Dw79kDoaFFNLK9tUCiUuqiUur7vPvlE0qpzkqpurYOTgghRMnsNelvCvxZSHkGUNLScfm/aVcBWzFOBH4H49j+/xTWIG+pOl3Uq2yXIISwa3Pnws03g6trKZvN5aWXXipQvmMHPPBARQVX4doAD2Cc+3QeuB/4EPgJyFRKyZAeIYSo4uw16S+P/GterrWeprXerbWeDfwbeDBvVR8zWusIrbUq6lWp0QshbC8lBVasgOs2trPUr7/+yoMPPligfMsW41K8VZHWOk1rvVFr/abWeojW2hXj4glLMS7lud+2EQohhCiJvSb9f1L4E/2ivgG43oW8n9tuKN+a99OyRbaFEDVXZCSsXFmqYT0A69evp3///gXKDx40zgW2J1rrk1rrp4BPgXhbxyOEEKJ49pr0J2Ac138jT+C4BW2LYyhTREKImmHZMujSBdzcSt307bffZtiwYQXK58yBf/2rIoKziYXAUyXWEkIIYVP2mvRvALoppdrlFyilXIAeeceKsxn4Bwi+oXxA3s8DFROiEKLa2bMHPv0Uxo4tddPjx4/TpEkTs83t8tWvX6oVPyudUupFpVRvpZRjIYdzgNL/BSSEEKJS2es6/UuBZ4H1SqmpgAbeBFKBD/IrKaXaAieBN7TWbwBorS8opWYArymlLmPc5OsuYBrwaf66/0IIYSYrC55+2pj4l0FUVBRvv/12gfL0dGjRorzBWd1M8n5fKKVSgEMYh/ScBh4Bkm0WWTkppSKA120dhxBCWJtdJv1a66y89fbnAp8DCuPa+y9orf+6rqrCuMvujd9ovIFx8tkzwL8w/uKahfEPByGEKGjSJJg3D265pdRNr1y5wqVLl/D29i5wbMoUeOyxigjQqhoC3hjnPPkBnYH/BRpgvH+OLLpp1aa1jgAi8t/LimxCiOrKXof3oLX+XWv9kNa6sda6kdb6Qa118g11kvNW2Im4oVxrredord201nW11m3zVvK5WpnXUNkiIiKYPXt2sXXee+89PDw8GD58OAANGzYs9XkuXrzIokWLzMruvvvuUvdTmbZs2YK7uztubm7MnDmz2Lq5ubn4+fkxePBgi9qXpm9RRa1aBW3bQr9+ZWo+ffp0nn322UKP/fwzFDK3t0rRWl/VWh/SWn+stX5Oa91Da90IuEVr7aS13mnrGIUQQhTPbpN+YR2LFi1i27ZtrFixosx9FJb0//DDD+UNzWpyc3OZMGECmzdv5vjx43zxxRccP170fPD58+fj4eFhUfvS9i2qoB9/hM8+g5dfLlNzg8HAxo0b6V7IoP20NGjcuLwB2o7WOtPWMQghhLCMJP3V3FtvvUX79u3p2bMnP//8s6l8+fLldOnSBV9fX8aPH09ubi5PPfUUv/32GwMHDmTu3Llm/RRWH+Czzz6jY8eOdOrUiZEjjd/wT548mZMnT+Lr68ukSZMA828M5syZg7e3N97e3sybNw+A5ORkPDw8GDt2LF5eXgQFBZGTk1Pi9XXv3p2kpCQA0tLS8Pf3L/VnFBsbi5ubG+3ataNu3boMGzaM9evXF1r31KlTbNq0ifDwcIval6ZvUQVlZ8P48cY1+cto6dKl9OnTp9Bja9fC9Oll7loIIYSwmCT91VhcXBwrV64kPj6e6OhofvrpJwBOnDjBqlWr2Lt3L/Hx8Tg4OLBixQoWL15M69at2bVrFy+++KKpn6LqJyQkEBkZyc6dOzl8+DDz588HYObMmdxxxx3Ex8cza9asAjEtW7aM/fv3s2/fPpYuXcqhQ4cASExMZMKECSQkJODo6MiXX34JwKBBg0hPTy9wfQaDgZSUFFzydjQ6cuQIHTt2NKsTGBiIr69vgdf27dtNddLS0mjTpo3pvbOzM2lpaYV+pi+88ALvvPMOta5bn7249qXpW1RBo0cbd951LGzRGst8//33vPXWW4Ue++9/wde3zF0LIYQQFrPLiby2VpbVHnbvNr4A2rc3bsTTvj0sWfL/dcaNg8xM2LjR+BPg0Ufh9Gnj+7g4Y5m/P/TuDY0aFX/OmJgYQkNDqV+/PgAhISEA7Nixg7i4OAICAgDIycmhRTHLhxRV/9KlSwwdOpTmzZsD0LRp0xI/hz179hAaGmpatnDIkCHExMQQEhKCq6srvnkZkL+/P8nJyQBER0cX2tfJkydxdXVFKeOmyEeOHMHHx6fAZ1BRNm7cSIsWLfD392d3/n9MUX299RYMHVrmcfxg/EP21ltvpW7dugWOnT8PiYml3t9LCCGEKBNJ+sugLKs99O5tfN0oIqJgmbt74e/vv9/SCIuntWb06NHMmDGjXPUXLFhQMQHluemmm0z/dnBwKHF4z9GjR82S/AMHDjBu3DizOoGBgWRmFhx2PHv2bO69914AnJycSE1NNR07deoUTk5OBdrs3buXDRs2EB0dzd9//83ly5cZMWIEEyZMKLK9pX2LKuazz+DCBePSOuXwyiuv8NFHHxV6bMYMGDWqXN0LIYQQltNay6ucL0wLAlUtcXFx2sfHR2dnZ+vLly9rNzc3PWvWLJ2QkKDd3Nz02bNntdZaX7hwQScnJ2uttW7btq3+448/TH00aNCgyPrHjh3Td955pz5//rypXGutz58/r2+//XazWBo0aGAWU1ZWlv7rr7+0l5eXPnjwoE5KStJeXl6m+rNmzdKvv/56sde3dOlS/fLLL2uttT5+/LiuV6+ePnPmTKk/p6tXr2pXV1f922+/6X/++Ud37NhRHzt2rNg2u3bt0vfdd1+J7cvSt7CxY8e07txZ69zccnVz8uRJ3a9fvyKPd+tW7lOUSt59qiLudzsB54roqyq+qur9XAgh8pX1fi5fLFdjnTt3JiwsjE6dOjFw4EDT8BxPT08iIyMJCgqiY8eO9O/fn9OnTxfZT1H1vby8mDJlCvfccw+dOnVi4sSJADRr1owePXrg7e1tmsh7fUxjxoyhS5cudO3alfDwcPz8/Iq9jqLG9AcHB7NlyxaGDx/OmjVraNasGS1btiztx0Tt2rWJiooiODgYDw8PHnnkEby8vEo8vyXtS+pbVDFXr8KYMbBpU7nH3Tz//PO88847hR47cwbuucduh/b0BurbOgghhBClo4x/MIjyUEpp+RyFsHNaG8fbvPaaccJNOfz1119MnDiRJddP2rnOkiXGUxQ25M9alFJorVUF9GMAOmitf6mAsKocuZ8LIaq6st7PZUy/EEIAjB0Ljz9e7oQfjJvchYWFFXn8k0+gCm9dIYQQohqyzy+XhRCiIk2ZAnXqQN++5e7q77//Zvv27fQrYtWfkyfhllvKfRohhBCiVORJvxCiZtuwAQ4ehM2bK6S7SZMmMaqYZXneeAPypr8IIYQQlUae9Ashaq6ff4bVq40TdytAbm4uP/74I2PGjCmyTr160L9/hZxOCCGEsJg86RdC1EzJyTBpEqxaVWHL6KxYsYKpU6cWefz0abBgDzshhBCiwsmTfiFEzXPhgnG3u/feMz56rwAGg4Ht27ebdr4uzAsvQGhohZxOCCGEKBVJ+oUQNcu1azB4sHHXXReXCuv27bff5s4776RWMd8apKRA3nYZQgghRKWS4T1loJSKAF63dRxCiFLSGkaOhE8/rZClOf+/W82aNWs4cOBAkXW++w7uuKPCTmlL/YHfbR2EEEKI0pEn/WWgtY7QWqv8l63jEUJYwGAwjq35178qNOEH+Oijj7jvvvuKfcq/ZQu8+26FntYmtNY7tNZ/2zoOIYQQpSNP+oUQNcMDD0C/fuDvX6Hdaq3ZunUrq1atKrKOwWCcRtCqVYWeWgghhLCYPOkXQlR/zz0HPXsaZ9JWsM8//5yhQ4eiVNFf+i1dCi1bVviphRBCCIsprbWtY7B7Siktn6MQVdSrr0LHjjBsWIV3bTAY8Pf3Jy4urtihPT17Gof3NGxY4SFYTCmFDEcsmdzPhRBVXVnv5/KkXwhRfT36KLRrZ5WEHyAyMpIhQ4YUm/BfvmzcBsCWCb8opb17Ye1aW0chhBAVSpJ+IUT1NHw4uLpCeLhVujcYDMTGxjJlypRi633xBcyYYZUQqgSlVC+l1E5bx1GhevSAP/6AzZttHYkQQlQYSfqFENXPW2/BnXfC9OlWO0VUVBSPP/54sU/5AQ4cgO7drRZGVXArcI+tg6hwTz9tnH2dm2vrSIQQokLI6j1CiOrDYIDXXoM2beCpp6x2mn/++Yft27ezYcOGYuudOAGnTxuH99gbpdTtFla91aqB2NKIEcYNFnr0gNry61IIYd/kLiaEqB4MBuOSnMOGwfjxVj3V008/TVhYWIn1Xn8dXnnFqqFYUzJgyYxWZWE9+3TPPcbxWSEh4OVl62iEEKLMZPWeCiCrPQhhY1evQu/e8MwzxrH8VpSRkcGAAQOIjY0ttp7WxmE9+/ZZNRyLlXa1B6VUFvA9UNKM1ruAcVprh/LEV1UUej/PzYX//V+YMME4MVwIIWyorKv3SNJfBkqpCOD168vkcxTCRrKz4YknjGOw77H+0PLp06fTu3dv7r777mLrbd4M587B6NFWD8kiZUj6fwAytNaDS6j3ELC6Wif9YPzDcu9eCAiABg0qPzAhhMgjS3ZWIq11hNZa5b9sHY8QNdbJk/DIIxAZWSkJf0ZGBomJiSUm/GBM+kNDrR6SNcUBlm5fbLf3QaVUhFJK57+KrFinjnHDhePH4dKlSoxQCCEqhjzprwAyvEcIG9izxzicZ+NGuN3SOafl07dvXz7++GNcXFyKrXfunPFvkd27KyUsi5ThSb8T4Ka1/s6KYVU5Jd7Ps7Jg0iTjOP9bbqm8wIQQIo886RdC1ByHDhl32v3hh0pL+Hft2kXjxo1LTPgBpkyB55+3fkzWpLVOq2kJv0UaNDB+s/Tqq8ahZUIIYSfkSX8FkCf9QlSitWthyxZYsqRS18Ls2bMnGzZsoGnTpiXW7dHDOPy7Kinrk6GaxuL7+ZkzcP48tG8PdetaPzAhhMgjT/qFENWbwQBDhhiz6Q8/rNSEf9myZYSHh1uU8H/3nXFoj71TSm1QSvmVov7NSqmJSinrbZBQlbRqBU2awM8/w7Vrto5GCCFKZLdJv1KqjVJqrVLqklLqslLqq1JsJnN9P5PzJnDtsUacQogKkJUF3brBoEEwd26lnvrq1ats27aNUaNGWVR/40Z47DErB1U5koF9Sqn9Sqn/UUp1VkqZ7e2ilGqtlHpQKfURcBp4Ejhog1htw8nJ+JR/6lTjH6VCCFGF2WXSr5SqD+wEOgCjgZHAncAupZTFa6kppdoBU4Fz1ohTCFEB0tJg5Ej4+GMID6/00z/99NMEBwdTy4JvFv76yzjN4NZqsEet1vp/AE8gFogAfgL+VkplKKVOK6VygFTgK8ALeAHoqLUufgMDO3Dq1CnLK7u7GzeEe/NN6wUkhBAVwC7H9CulngfmAO5a61/zylyBROB/tdZzLOznW4xPs9yB2lrrnmWMR8b0C2EN778PO3bAJ59Aw4aVfvrff/+dsLAwfvzxR4vqv/QSeHrCk09aObAyKM+YfqVUXaA70BVoDdwMXAD+C3yvtU6psEBtTCmln332WV599VVuu+02yxuePw/JyXDXXVaLTQghoIZtzqWU2gHcrLXucUP5dwBa6xIX7FZKPQbMx5jwf4Uk/UJUHQaD8elpnTrw+eeVOn7/ek899RQTJkzAx8fHovr33AO7dtks3GLJRF7LKKV0dnY2iYmJODk50axZM8sbZ2bCTz9B377WC1AIUePVtIm8XsCxQsoTMH4dXSylVBNgLsZvBTIqODYhRHmkpRkflQcHw4oVNsugDx48SP369S1O+LdtgwceqJoJf1kopUregayaqlevHnfeeSdJSUmcOXPG8oaNGsHly7BsmfWCE0KIMrLXX09NgT8LKc8AmljQfhbwC/CJJSe7ccfGG18WRy2EKN6XXxo3Ppo82aZjZAwGAxMnTuSdd96xuM3GjVVzWE85xOSN3V+ilBqYN8SnxqhXrx6enp5ERkZy9uxZyxs++KBxcu+mTdYLTgghysBek/4yU0oFAqOApy0dk6O1jtBaq6Je1o1YiBrAYIARI4xj9z/91Dg50ob+53/+h4ceeojatWuXXBlIT4cjR6rdBq1OwL+BNsDXwB9KqdVKqUeVUo1tG1rlqF+/Pm+//TYzZ84kNzfX8obDhxs3azhwwHrBCSFEKdlr0v8nhT/RL+obgOt9AHwEnFJKOSqlHIHagEPe+5sqNlQhRLF+/dW4xmWvXvDNN8Zx/DaUnp5OXFwczz33nMVtJk0y7sJbnWitz2itF2utBwK3AuOBXOB9jH8AbFVKPa2Uam3TQK2sQYMGzJkzh5iYGEo1d8vREVq2hN27rRabEEKUhr1O5N0J1L1x4q1SajfGaypyIq8Fw3Fe1FrPK2U8MpFXiLJ49104dgxmzYLmzW0dDWBcovOVV17h9tst2/ZDa+P0g61brRxYOVXURF6lVB2gH/AAEAK0AuKAr7XWM8rbv60Vdz//5JNPuP/++y2f3Ks1TJ8O/foZ95kQQogKUNMm8m4AuuWtsw+AUsoF6JF3rDh9CnkdxjgxuA+wtuLDFUKY+fNP6N0bfv/dOOmxiiT833zzDe7u7hYn/GCcazxihBWDqgKUUhuVUv9WSj0AtNJab9FaP621dsJ4392Jcb+Uam3IkCFMmzaNjAwL139QCl59Fb79Fi5etG5wQghRAnt90t8AY6Keg3FzLQ28CTTCuDnMX3n12gIngTe01m8U099uZMlOISrH1q3GJ/tvvlmlnn5mZ2fTvXt34uLiLB7LD8bJux98AKVoYhPlXKd/B+CLcVilxrhowl5gmdZ6fcVFaXsl3c8vXbrE6tWrGTt2bOk63rsXAgKMk3yFEKIcKv1Jv1KqVVnblpfWOgvoi3EFns+BFUAS0Dc/4c+jAAfs9xsNIaqP8+eNGfKPPxoT/yqU8AMMHz6c6dOnlyrhP3jQeFlVPeEvL611P611M+AO4BFgMcZhPV8ppaJsGlwlu+WWWxg7diy7d+8u3Rj/Hj1g/XpITLRecEIIUYwyP+lXSmUBC4C3tdZ/3nCsLuCgtc4pf4hVnzzpF6IECxbAxx/DkiXGp51VzH//+19mzZrFRx99VKp2/fvD4sVwxx1WCqwClfNJf4O8hy03lg/CuLlhL611bHljrAosvZ9rrfnXv/7FtGnTuMXSZZuuXTPO+n7hBWjbtpyRCiFqKluM6b8H40ZYvymlpuYNucnXF7hcjr6FENVBSgo89xwkJEBcXJVM+A0GA6+88gqLFy8uVbusLGjc2D4S/gpwWSl1Qim1Qin1klKqj1KqsdY6GvgMqF47FC+x25wAACAASURBVFhAKcWUKVN49dVXuXzZwl93tWvD229DdLR1gxNCiEKU50n/ncBbQH/gFuACxvHzDoA3cERr3bWC4qzS5Em/EDcwGIzJvuH/2DvzOJvq/48/D7JvGVuy74SSJYkokaJViaIMkVIqsoTs6ZtspVD8QqFVikjKvpetZF8b24wZY2bMYrb7/v3xGWHMcmfunTn3zn0/H4/7aOac8znnPdP4fF7n/XkvDhgwAGrUsNuiVHnhhRfo3r07DzzwQIbGDR0KXbvC7bdnkWFuxkVPvz/QELgTuB0ohIntP4nJrSoE9AF2i0iIWwy2iYzO5xcuXODIkSM0aNCAggULOv+gtWtNMrulrV4URckYdnj65wNNgQ+BAZgut4WBRsAy4CEX7q0oireyfDk0bAh16sDMmR4t+NetW8e5c+cyLPgTEuDECe8R/K4iInNFpH9SsYOiwG2YJoc/YsR/JeBXIMiyrH/tszT78fPzo27duhw+fJjIyMj0B1zhvvtMc4eL6bWWURRFcQ+uePqjgaeStnevHMsFvAK8D/QSka/dYqWHYVnWaGDUtcfU06/4PJGRxqtfqBCMH2/+68E4HA6aNm3KmjVrKFo0Yw1mR440vZcGDMgi47IAFz39vTG1+PeKSHyyc7Mxu7vPYZw+DUVkmKv22kVmd25DQkIYNWoUEydOpJCzf/vh4aak53vvmVgxRVEUJ8jsfO6K6D8GjBOReSmcexUYIiIVMnVzL0PDexSfJiEBXnkFzp+HqVOhShW7LXKKgQMH8uyzz9KoUaMMj23cGP74A3J5UV0wF0V/YtKXCcB+YBcmtKccJp7/bRGZ7A477caV+TwoKIixY8dmTPhfuGByXpo103KeiqI4hR3hPXOB/1mWlVLc/ilM23ZFUXIy331nknPvvht+/NFrBP/KlSspXLhwpgT/n39C+/beJfjdQHFM88IhwF7gLuAdjOBfDHxkn2nuJz4+/WtSokyZMowaNYqjR49y+fJl5wb5+ZlwuL17ITY2cw9WFEVxAlc8/XmAlUBrYBWwAlMrvwQwEogRkQbuMdOzUU+/4nNs3mxqVd5/P7zwglcp4KioKO6++27++OMP8ufPn+Hxzz5rmgjny5cFxmUhrnj6U7mfhVlDHO66pydgWZaMHy+0bQtNm2buHqGhoRw7dox69epRoEAB5wYFBcGECaZxnXr8FUVJg2z39ItIAtAeGAhUxnh6lmESfEsCb2b23oqieCh//QX33guTJhnR7+/vVYIf4KmnnuKDDz7IlODfutU4Zr1N8GcFYshRgv8Kw4bBmjWmh1xmKFGiBOXKlWPw4MHOe/zLlIHXXzdloRISMvdgRVGUNMi0p/+GG1lWGaA64AD+EpFot9zYC1BPv5LjCQw0yYanTpm4fS9tLLRq1Sq2bdvGyJEjMzW+VStYuBDKl3ezYdmAizH9uYBcSc6eK8cexCTwrhGR3W4y03aunc8PHjShPvXrZ+5ep0+fZtq0aXzwwQdYzpbmPHIERKB6da97oVYUJXuwI5E3rUVgtYjsydSNvRAV/UqO5exZ6N0bqlaFfv2gdm27Lco0gYGBDBw4kAULFjgvwK7hwAHjiM2s99duXBT93wCxIvJ80vd9gRlJp+OBDiLyu3sstZfk8/mFC7BjBzz4YObul5CQwKZNm2jdurXzg44fN/H9tWtrHX9FUW7AjkTer4DPrzGgL/ALpl7/dsuyMlb4WlEUz+HiRROr/8wzMHo0TJ/u1YLf4XDQsWNHJkyYkCnBD/DFF+bX4KM0w+RtXWEQMAfTmPEHYLgdRmUHfn7G8T5xovlvRsmTJw+tWrVi3rx5xMXFOTeoalXT8nns2Mw9VFEUJQVcEf0+uwgoSo7lwAF4+WV4801Ta3/jRlOdx8vp2bMnPXr0oFImw5IuXTIRTrVqudkw76E0cAbAsqzqQBXgYxG5hKnklskAGO+gfXuTyvLll5kbb1kWrVu3ZvDgwcQ7WxqocWNo2xbefz9zD1UURUmGK6LfpxcBRclR7N5tOoS+8Qb07w/z5kGFnNFm4++//yYhIYFXX3010/fo1cu7GnFlARGAX9LXrYEQEfk76ftEIONZ0V5Gs2bQvTusW5e58ZUrV6Z///4sXLjQ+UHNm5vwuj/+yNxDFUVRriGPC2N9fhFQFK9n0yaYOdPU1583z2sTdFMjMjKSUaNG8f3332f6HhEREBCQ+WTOHMIWYKhlWQnAG1y/y1sdOG2LVdmMZRmP/5AhMGIEFCmSsfFVq1alSpUqrFu3zvkYfz8/yJ8ffvkFHnoowzYriqJcwRVP/5VFoCM+tghYljXasiy58rHbHkXJEAkJMGaMCR/46iv4/HMTypPDBL/D4eChhx5i5syZ5M6dO9P3ee01ePttNxrmnQzGOHmWYhw6o6859wyw1QabbCFXLhg0yAj/8+czPt6yLO655x6GDh3qfIx/oUKQO7dPJ5UoiuI6rlTvqYER+tWA48ADInIy6dwa4F8R8XeTnR6NVu9RvIKICPjf/4xSqVPHxO3n4JKAzz33HC1btqRv376ZvkdkpIl4mjPHjYbZhDuac1mW5SciF5Idqw8EikiwSwZ6CM7O55GR8PffUKMGlMpE//njx4/z0UcfMXHiRPI624zrl19Mgu9TT2X8gYqi5BiyvWTnNQ/O8YtAeqjoVzyaX34xyYB58sDHH3t1FR5nWbp0KV9//TWLFi1y6T7TpsEdd0BGqi16Kq6KfsuybsE0Y2yF6bweCqwFpohIoHustJ+MzOcJCSb3PX9+I/4zyokTJwgJCeH22293XvjHxsLOnSbeX1EUn8Q20a+o6Fc8kKgo+PRT2L/fdBd6/30oW9Zuq7KF06dPM2jQIBYtWpTp8pxgNkbatoXt291onI24WKe/JrAJKA5sBgKBskBz4CLQUkSOuMtWO8nofO5wmGi5du3gnnsy/rzz588TGBhIrVq1yOdsq+ewMPjxR+jRI+MPVBTF67GjTj+WZd1iWdYky7L+tCzrWNJ/J1qW5RvqQlE8jd9+g/vvN52Emjc3cSnz5/uM4I+MjMTf3585c+a4JPgBXnkFhmvh4Su8D4QDNUXkPhHpKiL3ATWTjvtsXclcuUwri23bTF58RildujR+fn4MHjyY2NhY5wYVL27i/OfPz/gDFUXxWVyJ6a8JbARuJod7ftJDPf2KrZw+DQsWwLFjEBcH774L5cvbbVW243A4aNmyJe+99x733nuvS/eKjTWFUtascZNxHoCLnv4woK+IfJ3Cua7ADBG52VUbPQFX5vOwMDhyJHOtLQICApg0aRITJ04kf34ni999/bWpJVq5csYfqCiK15Lt4T2WZS0B6gFtryTwJh2vBKwC9onIk5m6uZehol/Jdq6E7Pz8s6kbOG0a3Hab3VbZSp8+fahbty5vvPGGy/d66y0TOVGvnut2eQouiv5o4GkRWZ7CuY7ANyJSyFUbPQFX5/OLF2HRIujXL+Njz549S3BwMLVr13Y+1Oevv0xPjRIlMv5ARVG8EjtEv894ftJDRb+SLTgcpsTmjh0QHQ116xplkceVdhs5g1mzZhEXF0f//v1dvldYGAwcCP/3f24wzINwUfRvwfRmeVhEHNcct4DlQDERyUREu+fhjvl8zRrYssWEh2U0yuzChQscO3aM+vXrU6BAAecfGBqqVX0UxUewI6Y/L3AplXOXks4riuIKDocR+Q88YLIEd+wwAcSffgqvv66CH/j555/5+++/3SL4Afr2Ve2UAmOBB4ADlmWNtSzrZcuyxgD7gLbAGFutc4Gs6Lty//0mPGzPHlPhJyP4+flx6623MnjwYGJiYpx/YHg4LFmScWMVRfEZXPH0+4znJz3U06+4FYfDLN6zZkG1asaj36ULlC5tt2Uex4EDB+jatSvbt293PhwiDc6cgc6dYfNmNxjnYbihZGd7YDzQELAAAXYC74jIr+6x0n7cOZ8HBZm2GNWrg7NO+yucPXuWKVOm8MEHHziflP7jj+Ztww3/FhRF8VzsCO9pD/wMHAO+Ac5hEnmfBmoAHURkVaZu7mWo6FdcJj7exOdv2mTKgNxxBwwbBrfeardlHkt0dDTNmjXjl19+4VY3/Z7efBOee840K85puKM5V9J9CmIKOFwUkWjXLfMs3D2fBwXBuHEmv75YsYyNTUxMZOPGjbTOSKOILVtMMkrRohl7mKIoXoMtdfp9xfOTHir6lUwRFmaScTdsMHX/3nwT2rTJuDLwQRITE3nmmWeYOnUqFSpUcMs9AwNNDHZOi+W/grtEf04nK+bz0FB45x0YNSrjG3Yiwueff06XLl0oVMjJXOmZM80uYbt2GTdWURSPx9bmXDnd85Mcy7JGA6OuPaaiX0kXEVi/HubOhbx54eaboU4d6N5dY/MzgMPhoHPnzowePZp6biyv89hj5n9NTi2C4mIi7xCgvIi8lsK5j4BTIvKBqzZ6AlnlxImONlV1S5XKeNuMoKAgxo0bx/vvv++88J86FW6/3cT7K4qSo7AjvMdnFoH0UE+/kiqhoaaG/tGjsHOnCe7t1w+aNrXbMq/liSee4P777+e1126YejLNmjUwfXrOzoN0UfQfBCaLyOwUzvUE3hKRuq7a6Alk5XweHw/79pmv77gjY2PPnz/P8uXL8ff3d37Qvn2QP7/x+iuKkmOwQ/T7zCKQHir6lf+IjzddcT/5xITvlC4N/fubyjt5taCVq7z11lvcdNNNvPfee2697113wS+/5FwvP7ilTv/DIrIuhXOtgeVap985HA4T43/ffZDRHnIiwrp162jdurVzyb0icPIkREQYr7+iKDkCO0p2VgRS67h7HKjkwr0VxTuIizO18594wsSH9O5tFtiPPzYlYJYsMau7Cn6XmTx5MpUrV3a74F+zxhRHysmC3w1EA6llS5cHYrPRFq8mVy4YORJ27zYbgRnBsixatWrFgAEDuHQptYrZ1w0w3Xp/+cWEFiqK4tO44ukPAV4XkYUpnOsGfCQiPrGMqqffh4iNhW++MfXyQ0NNHb5ixaBXLxOfr2QJ7733HnFxcYwaNSr9izOAwwHdusG8eTn/vcxFT/8SoApwl4jEXnM8H7ANCBCRx9xjqb1k53y+YwfUqJHx3P2QkBBGjRrFe++9R1FnqvSIwMSJxgGhoYWK4vXYEd7jM4tAeqjoz8GcOGE89uvWwaVLpsTLXXcZj36NGnZb5xNMnz6dX375hRUrVrj93gMHGkeoG9MDPBYXRf/twBYgBFgAnMF4/rsBfsA9IvKXu2y1k+yezzdtMl7/jP4NXuncW7duXQoXLpz+ABE4exYuXjQlPRVF8VrsEP22LgKWZVUApmK6QVrA78AbIhKQzrjGQB/gXkyIUgiwERghIicyaYuK/pxAXBzs2gUrVphWmqGhULw49OxpOuJq3etsZ926dUyYMIGVK1eSK5cr0Yg3cvEitG1rvK2+gBuaczUFJgHNMaGhDmATJn8rx/wW7ZjP162DjRthxAgTkeMsERERHD9+nEqVKnHzzTc7N+j4cQgIgIzU/lcUxaOwq06/LYtAUonQvzBxpCMw/QHGAwWBBiISlcbYScDdwEJMC/lbgXeA0sAdInIqE/ao6Pc24uJMNZ0//4Tly009vVy54PnnoWFDuPNOuy30ebZv387MmTOZO3eu8x1JM4C/P7zwgu9oHzc25yrA1RLNMa5b5lnYNZ/v32/+W7Nmxir4hoWFMXz4cMaMGUPJkiWdGzRpkpnjtJynongldtfpz9ZFwLKs14EpQC0ROZp0rAomsXiwiExJY2wpEQlOdqwScAIYLyIjM2GPin5PJjoatm833vuffzaCP3du6NHDxOHXqwfO1r5WsoV58+axdOlSFi9enCWC//BhmDYNZsxw+609FhfDe4oDT2BCOiMw8+WG5HNpTsDO+fzcOeOEb9DApAs5y6VLlxg6dCijR4+mVKlSzg2aMsW89fr5Zc5YRVFsw1bRn91YlrUayC8i9yQ7vh5ARFpl4p5BwM8i0isTY1X0ewIJCab7zebNxm22a5cpoVmsGHTsCK1awa23apiOhzN37lzmzZvH2rVr3R7Sc4Vu3Yzg96U/hUwvEpbVAliCCdu8gAmnvJJ6ugIYeMX5khOwez4/exbefdd8ihd3flxUVBRHjx6lQoUKlHC2FNX27VC/PhQsmDljFUWxhczO597aBvQ24KcUju8Dns7ozSzLqoMJ7zngol1KdhAXZ+JSDxww4TmnT5ta1LlymY43LVrA8OGmKU1G3GWK7axdu5YVK1ZkqeCfP9/8afiS4HeR+cBOoLeInLIs62FgHDAdk8P1p2VZ94jIfjuNzCmUKwfjx5ta/h98YKY1ZyhUqBC1a9dmz549lC1blkqVnKiafdddsGyZmScfeMA1wxVF8Xi81dMfB0wRkaHJjo8HhoqI0y8zlmXlAVYDdTDhQhdTuGY0kGatQG/8PXo0IhAUBIcOmdj78+dN/fvTp00VnZtvhiefNPvgderATTfZbbHiIsuWLWP58uXMnDkzS0J6wLwvNmli/qQyEjedE8iMZ8iyrFuBAKDylXwny7ICgHYicjDp+/8BTUSkjbtttgO7Pf1XEDGl9Vu1ylhyb2xsLEOGDOG1116jmrOdeKdMMR7/tm0zZ6yiKNmKr4X3uFP0zwJ6AR1EZFUm7fGIRcLriI83AawnTsCRIyYkJzTUCPrcuc0xPz+oXh1efdWsgs4mqilexfDhwzl48CDff/99lgl+gMmTzZ/TYz5RTPh6Min6iwOhQBERibIsKy9wGSglIheSrikEnNeOvFnDF18Y38Yddzg/Jj4+ntGjRzN69GhuctYh8t135h9GTm9YoSg5AF8L77mISRxOTomkc06R5KHqA7yQWcGvpEFsLJw5Y7zzp06ZkJyLF02t+4sXoXBhqFYNGjWCo0eNC7ZuXfN9Fgo/xbPo168f0dHRLF68OEufExgI+/aZ2vyKc4hImGVZU4GywDERibMs6w/gJ8uyxgAHMQm+J200M0fTrZuJ7w8PN15/Z7jpppt49913Wb9+PS1btnQuVO7pp00+VP36GvumKDkUb/X0rwHyikiLZMfXYX6mdKdGy7KGY8p8viYiH7toj0d5hrKF6Ggj5q/9XL5sPrlzmxCcoCDjzS9aFB56yHjtK1Uy3vxq1aBUKRX3Ps6AAQO47bbb6NUrw/nzGeaee0x11owkR+Yk3FiyszLwOdAaUy45HOgkImtdvbcn4InzuQgsXWoa6mZUj3/xxRfUr1+fhg0bOjfgk0/M/Ny+fcYNVRQlW/C18J43MP0BaorI8aRjlTElO4eKyOR0xvcHPgSGi8gEN9jjcYtEpoiPN7HzQUHXf6KjzXmHw1TIOX/eiPs8eaBdO6hQAVatMmK/eHFT6/7yZShTxvcCpxWnSExM5JlnnqFXr1489NBDWf68yZPNptOUVIv55nzcJfqvuV95oAxwMK3eKN6GJ8/nO3ca30nlys6PcTgcjBkzho4dO9KkSRPnBn34oWkYkA3/NhVFyTi+JvoLYZpzxXC1Odc4oAimOVdk0nWVgGPAWBEZm3SsC7AI+BUYk+zWEZmpQOGRi0RCgvGoX/lcuHD16/DwG68XMVmOefJAlSqmZESePCYeIm9eE3ZTpoxJmt25EypWNB/dBlYySHh4OG3atGHQoEE888wzWf682Fi4917YutX5Sig5kYwuEpZlNReRLVlpkydiWZbIO++YObFMGRPuUr8+OFsGM4v56SdTrKx/f+c3SkWExYsX88gjj5AvXz7nBh05Yv7x1KuXaVsVRckafEr0A1iWVRGYCrTF1I1eDbwhIievuaYyponMGBEZnXRsHvBCKrddLyKtM2GL+0W/wwGRkUagR0Rc/ST/HsziZAy5+vVNN5mE2BIlzNe5c5sF7JZbYNEiI/AB+vQxoTn//GM8+uXKmTJuRYqY+vaK4kaCg4Pp0qULEyZM4K677sqWZ/btC2+9ZRJ4fZlMiP5E4DywDFOnf7WIxGWVfZ7Cf/O5iNnV/Ocf+Ptv4zAB4+ho1Mh8bJoj160z1YoHDcrYuC1btlC7dm3n6/ifPw8HD5q3ZkVRPAafE/2ehGVZIhcuQFSUEc5RUVc/KX2fmGgGJv/dXyvac+UywrtoUfMpVuzq14ULm+uuiPTAQBO7sGWLuXfjxuaa0qXh22/NfnC1anD33Vfr1ufPn32/IEUB/vrrL8aMGcO8efMomk07REuXms2qt9/Olsd5NJkQ/WWBx4HHgPuAWMwO6RJguYhEZImhNpOuEyc83DT+27nz6q6pn5+Zdxs1yrbeIDExprdW69YZGzdt2jTq1q1Lu3btnBswf75xHD37bIZtVBQla1DRbyOWZYlMmwaFCplPwYJXv07+fcGCN8a5JyaaF4LTp00YTp48RpSXKQOffnp1D7dPHwgLg19/NS8FpUubzC4R8zKQL5/G0CseyZQpU/jxxx9Zs2YNebLpb/TyZbNp9eefWoUQXIvptyyrCNAB8wLwEFAAWI95AfhJRM66zVCbydTO7YULsGOHeRGIiTE7q3fcYbLHS5XKGkMxU//s2abCT0aa6k6fPp0qVarQsWNH5wYsWmTytR5+OHOGKoriVlT028h/i0RcnClFGRlpvPT58sGSJWZBuHTJiPbjx+HYMVObHkwzlMREE4qzf7/xGJUvb5JjCxbU6jaK1+Pv7098fDxffvllltbgT06nTian3Bdr8qeEG6v33AS0wbwAPIop57kTWCIi77l6f7txS7hmfDz89Rds2gQhIUadV69u3PJVqrjFziucO2e69777rllGnOX06dMEBQXRqFEj5wZER5swp2bNMmeooihuQ0W/jViWJXLgAHz/vYmhv/VW4+HJn99M+CVKmJcAFfCKDxEbG0vPnj3p0qULjzzySLY+e+9emDvXt6v1JMdFT39ZEQlM5VwzTBjQoyJS1xUbPYEsydESMb1I1q0zWbiWZZKD77/fLTsBYWEmlK1bt4wlq4eFhbFlyxYedtaDHxYG33wDL72UOUMVRXELKvqzEcuyRgOjrj2mv0dFucqePXvo378/X3zxBZUzUl/QDSQkQNeuppNpNoVXewUuiv4oYDrwvohcTHYuL5BbRGLcYKbtZEs1NofDJAivXQvBwSYss3lzkzDrQr7V+vXQtGnG/u6/++47Ll68SJ8+fZwb8MsvcPgwvP565oxUFMVlVPTbiEeW7FQUm/joo49YsGABK1asoGTJktn+/K5d4amnTHiPchUXRX9jYCTQEpgMTL1Sm9+yrPbAMhG5yW3G2ogt83l8vCnEsH69KZNZoYKpkV+pUoZvNWKE+dt3thcXwPLlywHo0KGDcwNWrzY72rVq6Q62otiAin4bUdGvKKbh1siRI4mMjGTq1KnksqEo/tKlJrFx2bJsf7TH46LorwG8iymRXAy4gOmBkhuoB/wtItlTgzWL8Yj5PCDAeNQDAoznv21b48J34t+Uw2Hi+1u2zFhln8uXL7N7927uvvtu5wYcPWoSlitXVuGvKNmMin4b8YhFQlFs5PDhwwwZMoQ+ffpkS4fdlIiLMw2iV67UirQp4aLo3wKUA+YBYUBe4HmgLvA90FdEQt1kqq143HweHQ2//QZ//GFE/xU1n05JqoAAU2a/cWPnHxUREcGUKVMYOXKkcy/t//xjctlGjVLhryjZiIp+G/G4RUJRspEZM2bw+eef88MPP1CxYkXb7HjpJRg4EGrWtM0Ej8ZF0R8NPCUiK645lgt4BXgf6CUiX7vHUnvx6Pk8IQE2bDAJwQkJJhH4vvuMxz0FIiLg99/hySedf8Tff//NggULmDBhgnPldXfuNFXqxo717ZbXipKNqOi3EY9eJBQli4iNjeXdd98lICCAzz//3JZwnivMmmX6JA0ZYpsJHo+Lov8YME5E5qVw7lVgiIhUcNFEj8CyLHnnnXcoWbIkt912G/Xq1aN06dLZWm7WKRISjPhfu9YkAnfoAE2a3OBx/+UXU2lz8GDnnfFHjhwhKiqKOnXqkC9fvvQHHDhgBH/16qm+gCiK4j5U9NuIin7F11i9ejWzZ8+mR48etG/f3lZbzp6Fjh1NbyR1NKaOi6J/BPAq8JiIbE927jHgGxHJEUFVV+bz4OBg9u3bx759+wgKCgIgT5481K1bl0aNGlG5cmXPeRG4fBlWrDAhQEWLGtd+7dr/nf7jD9P7MSMe/6CgIM6fP0+1atUo6Eznr4AA05Omfn1tEqkoWYyKfhtR0a/4Cg6Hg549e3Lu3DkWL15M4cKF7TaJ5s3hs8+gXj27LfFsXBT9eYCVQGtgFbACOAGUwFT1iRGRBm4y1VbSms/j4+PZt28fu3bt4uTJkwDkzp2bOnXq0KhRI6pWrWr/i0B4uAm3OXgQypaFZ5+F0qVJTITNm6FFC+dfjoOCghg7diwTJkygWLFi6Q84cgQ++QTef980p1QUJUtQ0W8jKvoVX+DQoUMMHz6cNm3a8PLLL9ttDgBjxpiiJjblDnsVrnbkTRL+/YCXgNrXnArHxPuvdtFEjyCj83lCQgIHDhxgx44dHD9+HBGhWLFiNGvWjMaNG1PAzmYRZ8/CokWmF8ADD5DQqg1vDszFyJHO9wQLDQ1lxIgRjB8/nhIlSqQ/ICAAPvwQJk3S5F5FySJU9NuIin4lJ+NwOBg+fDghISFMmTKFIkWK2G0SYKr0bNoE48fbbYl34KroT3avMkB1wAH8JSLR7rivJ+CO+TwsLIzt27ezY8cOYmJiyJ07N/Xr16dly5aUKVPGTZZmAIfDZPSuXs2lIuUYfvJF3hxeiCpVnBseGRnJkSNHqFq1qnMe/4QE2LjRJBkriuJ2VPTbiIp+JaeyZcsW+vXrh7+/P/3797fbnP8ICTGFS3bsSLdyoZKEO0V/TiYr5vOEhAT27t3Lhg0bOH/+PLlz5+aOO+6gVatW+Pn5ufVZ6RIURNz8rzh0UMjXvBE1ezoX7xMXoCRh6gAAIABJREFUF8fOnTspXbo01apVc+5Zc+eaBOPSpV00WlGUa1HRn41YljUaGHXtMf09KjmJhIQE3n//fXbt2sXs2bOd29bPRpo3h48/hjvvtNsS78HFmP5cQC4RSbjm2IOYxlxrRGS3m8y0nexw4iQmJrJnzx7WrVtHaGgoefLkoVmzZrRq1cq5pFk34EhwcOa7LRxf+g+tWjrg+echnRyd+Ph4hg0bRrdu3bj99tvTf0hkJAwdaj7ly7vJckVRVPTbiHr6lZzEt99+y7fffkuPHj3o2LGj3ebcwOzZJlrhpZfstsS7cFH0fwPEisjzSd/3BWYknY4HOojI7+6x1F7smM/j4+PZvn0769atIyYmhlKlSvHggw9Su3btLE8MnjULil0OomvYTOPxf+EFqFQp1esdDgcTJ07kjTfeIL8zXfAuXzZ5BT16aHktRXETKvptREW/khM4deoU3bt355ZbbmH+/Pnk9cC4mc2b4dtvTZ6gkjFcFP3/Ymrxf530/TFgNTAQ+AwoKyI5IoDbE+bz8+fP89tvv3HgwAEsy6J58+bcd999zonsTPD99+DnB/fdGQ5ffGHqez76qNlSS+WlY9OmTTRp0sS5Ov4A69fDXXdpu2xFcQMq+m3EExYJRcksDoeD0aNHExQUxGuvvUY9D619eewYdOpkGoBq/5+M46LojwHaichGy7KqA4eBO0Tkb8uy2gGLRKSkO+21C0+bzxMTE9m6dStr1qwhJiaG2rVr06FDB0qWdO+vOzoa9u411bAsRyIsXWqScVu3No0wUvDSL1u2jIiICJ577jnnHjJ2rCm11aSJW21XFF9DRb+NeNoioSjOsnDhQqZPn86kSZNo0aKF3eakSlyc0QlLl6YZeaCkgYuiPwh4SUR+tCzrRWCCiJROOtcG+ElE7G/a4AY8fT4/dOgQP//8MyEhIVStWpUnn3zSbcnAly6ZMrjvvptUZl8ENmyAZctMAk3nzjc03lq8eDGBgYH069cv/QeImJu3bAmtWrnFZkXxRVT024inLxKKkpxz587x9NNPU6NGDT777DNuuukmu01Kk8GDTVOhRx+12xLvxUXRvwS4BRgP/A/YISI9ks69BLwpIrVTv4P3YFmWbNki1KkDxYvbbU3aHDt2jB9++IHQ0FBq1arF448/TnEXjT51yvTWGj8+2c+/Y4eJratRwyT9XhPWs3PnTooWLUr16tXTz0EQMU28YmNN915FUTKMin4bUdGveAuBgYF8/PHHXLhwgXHjxrk9RCAr+Oor+PdfUwBEyTwuiv4amC681YDjwAMicjLp3BrgXxHxd5etdmJZlqxYIRw4AGFhV49XqgQNG8Jtt3lms9mDBw/y448/EhYWRoMGDXjkkUcy3VMjNNT8mytbFm65JdnJ/ftN3H+lStCz53+/jFOnTnHmzBmaNGlCbmfi74KD4e+/oU2bTNmoKL6Min4bUdGveDrR0dH07duXo0eP8sknn9CwYUO7TXKKhQtNkuGSJXZb4v24o06/ZVl+InIh2bH6QKCIBLtkoIeQ0nwuYkTw7t2wb58JNxMxya8NG8Ltt3vOroCIsHfvXpYtW0ZkZCQtWrTgwQcfJE+ysJz0iImBEydMKH/tlPZwDhwwdfhr1jSe/7x5OXDgAJ999hnvvfeec0nHixZBVBT07p0h2xTF11HRbyMq+hVPxeFw8OGHH7J+/Xq6d+9Op06d7DbJaXbsgL59Ydu2G8KIlUzgqui3LOsWTLWeVkAJIBRYC0wRkUD3WGk/GZnPQ0Jgzx7zCQ+/+iLQrJl5GbC7UI2IsHHjRlauXEnBggV55plnqFGjhtPjExNNjH+bNmmE4O/dazz/t90G3bpx6tw5Zs6cyYQJE5x7yIoV5i3q8cedtktRfB0V/Taiol/xNBITExk2bBj79++nd+/ePOplwfCXLsEzz8CCBeBhfcG8FhfDe2oCm4DiwGYgECgLNAcuAi1F5Ii7bLUTV+fzkBD44w/YtcuErVsW1KplXgSqVk21AmaWExERwbfffsuRI0eoW7cuTz31FIUKFUp3nAjMmGHsb9QojQt374Yvv4SWLZHHHmPd+vXcd5+TVVwTEmDrVpPgqyhKuqjotxEV/YqnICLMnDmT2bNn06NHD15//XW7TcowCQnw7LMweTJUqGC3NTkHNyTy1gPaXonlTzpeCVgF7BORJ91iqM24ez53OODwYaNpjx83x8qWNZ7zunXt6Ve1b98+Fi9eTFxcHE888QSN0lTzhtBQk+SbbiPeNWvgp5/gmWeYsWcPHTp0oJIzJbdiY2HUKLO14IlJE4riQajotxEV/YrdxMfH88EHH7Bnzx569epFu3btsryTZ1bRpg2MHKkV/dyNi6I/DOh7pTlXsnNdgRkicrOrNnoC2TGfnz1rSuDv22deCvz84L77oEGD7H0JiIuL44cffmDnzp00bNiQp556Ks2mfCEhMGWK0eVpFvxyOOCbb4jdvp1hERH0ePNN6jtTqefff2HSJJgwATKZhKwovoCK/mzEsqzRwKhrj+nvUbGD6Oho+vfvz549e/D39+eVV17xWrEPxsPfoIFW6skKXBT90cDTIrI8hXMdgW9EJP1YES/ADidOSAisXWuK2TgcUK0aPPAAVKyYfTbs2rWLxYsXU7hwYV544QXKlSuX4nVHj5qO2O++C0WLpnPT2Fgcs2bx0+bNtHnvPYpWq5a+IcHBJoO4Zk3PyY5WFA9DRb+NqKdfyW4CAwOZMWMGgYGBtG3blqefftpuk1zmiy+M6Jk0yW5LciYuiv4tQATwsIg4rjluAcuBYiJyj3sstRdPmM+PHYPff4eAANN9+kovqzSc8G4jJCSEL774gqCgIB5//HHuvvvuFK6B06ehdGlI5d3gei5eZOewYdzs50fVd95JP3wnOtrU8i9RQmP8FCUFVPTbiCcsEopv8OeffzJs2DDi4uL47LPPqFWrlt0muYWVK+G330wcv5I1uCj62wM/A8eAb4BzmETep4EaQAcRWeUuW+3E0+bz+HgTCrRhgylyU7cuPPxw1ie4JyQksGTJErZu3Urr1q3p2LEjua6JPbp8GU6eNDY1aODcPb+ZNo0LS5fyyoAB0LFj2hfHxsLbb0OvXqYykKIo/6Gi30Y8bZFQchYiwurVq/nyyy+JjY1l/PjxVK9e3W6z3Mb8+fDDDyb3T8k63FCysz2mI29DwAIE2Am8IyK/usdK+/Hk+VzElMdfscIk1pYqBR06mEiYrHumsGHDBpYtW0bdunV57rnnyJfkqXc4TPh906bQrp1z91uzZg1/ffUVbxYuDC+/nLbxiYkmubdLF6hXzw0/jaLkDFT024gnLxKK9xIeHs7QoUPZsWMHffr0oXv37s41vPEifv3VJAVu2mRPFRNfwh3NuZLuUxC4GbgoItGuW+ZZeNN8fv68eQE4fNhEzDz4INx1V9aVBf3nn39YtGgRpUqV4sUXX/yv4++cOaYvgRNFgACTQLxl/Xpa799v6vO+/nrqibsiJpYoMhLq1HHTT6Io3o2KfhvxpkVC8XyOHDnCpEmTOHPmDM8//zydO3e226Qs4eBBGDbMNOXMYe8yHom7RH9Ox1vn88uXTZjc9u1QqBA8+ijUr581LwCnTp1i9uzZFCtWjN69e1O0aNH/Sno2aODcM0WEjz76iB7t21Ns7ly48054+unUBwcEmBinbt3c+8Moihfic6LfsqwKwFSgLWar+XfgDREJcGJsfmAc0A3TbGYPMERENmTSFq9cJBTPISIighEjRrBt2zYefvhh+vbtS9myZe02K8vYuxfGjzfJu1qSO3twMaZ/CFBeRF5L4dxHwCkR+cBVG+0gJ1Zji4yEpUtNYvzNN8OTT0IGGvE6zenTp/nss88oWrQoffr0ITa2KBMnmn/bzvy7Dg0N5Z133mH48OGUO3TIGN2/P1SpkvKAH380zQ7efNO+LmeK4gH4lOhP2l7+C4gFRmBiS8cDBYEGIhKVzviFQAdgEHAc6Ac8BNwtInsyYY+KfiVT7Nmzhw8//JDw8HA6derEc889Z7dJWc6OHSY3b8sW45FUsgcXRf9BYLKIzE7hXE/gLRGp66qNnkBOm89DQ2HJElNqs0wZ6NzZyYo7GeD06dPMnj2bwoUL89BDLzFzZlHGjXMu2TgmJobly5fTsWNH8gN89JEpU9SvX8rNALZsMaU8a9aEPHnc+4Moipfga6L/dWAKUEtEjiYdqwIcAQaLyJQ0xt6O8ez3FJG5ScfyAPuAQyLyaCbsyVGLhJK1/PPPP4wePZpTp07x6KOP0q9fP4r7SD3qv/82u/Pr1mV99RHletxQp/9hEVmXwrnWwHKt0+/5BAbC11+b5mB3320K6KTZZCuDXBH/uXLdzAMP9KJcuSKpOu2Ts3XrVipVqmT6Axw8CDNmQNeuxtDknDkDFy+aHQH1HCg+iK+J/tVA/uR1oS3LWg8gIqn28rQs6x3gHaD4tUlolmWNAYYCRUUkNoP25NhFQnEPwcHBfPfdd2zatImbb76ZHj160KRJE7vNylb27YNx4+Djj6FkSbut8T1cFP0hwOsisjCFc92Aj0QkR7zG+cJ8LgLbtsHy5cap/tRTphSouzhz5gwzZnxG3rxVufvurrRr51yDgenTp1O6dGmeeeYZY+SXX5p6/QMH3tioKyjIVAEYPdo0DFAUH8LXRH8g8JOIvJTs+AxM18hSaYz9GmgoIrWSHe+MqT9dT0T2ZdCeHL9IKBknNjaWadOm8dNPP5E/f34GDx5M+/bt7TbLFnbsME23pk2DHJyq4NG4KPqXAFWAu651iliWlQ/YBgSIyGPusdRefG0+j4qC774zpUBr1DDVMQsXds+9Dx06xGuv/UPlyrcwa1az6+r8p8b333/PxYsX6d27tzkQEmImj7vugieeuP7iS5dMNYAhQ6B8efcYrShegK+J/jhgiogMTXZ8PDBURFIN9LMsaxXGm98s2fEHgN+Ae0VkY7Jzo0mW6JUcb/w9Ku4nPDycSZMmcfbsWRITE3nsscd45JFHyOPDsacbNpjcvHXrbnTWKdmHi6L/dmALEAIsAM4At2KKIfgB94jIX+6y1U58TfRfy6FD8NVX5utnn3Vf/f9p0w6yefOv9OxZi/btH8RKJwk3NjaWLVu2cN999109uGIFrFljBH6pa/x68fEmYaFwYe3eq/gMKvrJOtHvhD0+u0gopvLOF198wapVq4iIiKBLly68+OKLPi30r7B1K7z2mhH87vIeKpnDDc25mgKTgOZALsABbMIk8e5wj5X2o/O5qf6zaJGJrGnZ0jQAy53btXtGRQmzZ2/l6NElvPBC53TDG0WEyZMn8+KLL17NeQoPh/ffN+U9n3rq2otNHsDx48ZYRcnh+JroDwJ+zGR4zzfAHRreo7jC7t27mTZtGkePHqVx48a0bduWDh06pOvB8iW2bTO5eHPmmLhhxV7c2JyrAFebc8W4bplnofP5VURg40YT+1+mDPj7mxKgmSUqCgYNclC9+mKCg3fx2muvmcTdVAgNDWXkyJEMGjSISpUqXT3x66/w228wePD18fzz5hmj/f0zb6SieAG+JvrXAHlFpEWy4+swP1NaibwjMWU+kyfyjgbeRhN5lRRITEzk22+/Zc+ePZw5c4YiRYrQuXPn67eflf+YOdNU1ps/XzvtegouhvcUB57AxPVHACeADSIS7EYTPQKdz1Pm9GmYOxcSEoymrlw5c/e5fBneeQcefTSKDRumUaBAAV555ZVUu43HxsaycuVKWrduTbFixa6eiIiAiRNNB7LOna/W7V+xwmQlV6qktfyVHIuvif43MNvMNUXkeNKxypiSnUNFZHIaYxsCu4AeIjI/6VgeYC9wVEQeyYQ9ukjkQE6cOMHPP//Mzz//TFRUFE2bNmXw4ME5ummWOxgyxMQG//CDCn5PItOLhGW1AJZgYvcvYJohXlFfK4CBV0on5wR0Pk+bS5eMQ/3UKRNh07Rpxu8hYorvnDwJpUodY8aMGbRo0YLHH3881d3SXbt2UahQIWrVqnX9idWrzVbE4MFXqwQcOWIeUqmSdv9TciS+JvoLYZpzxXC1Odc4oAimOVdk0nWVgGPAWBEZe834r4EHMc25TgAvAx2B5iKyKxP26CKRAwgPD+fzzz9n+fLlxMTE0KFDB5o2bcr999/vVNUJxVTnOXYMpk+32xIlOS6I/mMYh0pvETllWdbDmPl2OiaRtxEmkXe/Ww22CZ3PnSMhARYvhj/+gHbtzCejjvWYGPi//zN9uFav/p2ff/6ZXr16Ub9+/RSvX7BgAREREbz88svXvxxERsL//geNGl2t8HPmDLz3ninpqTWClRyGT4l+AMuyKgJTgbYYz9Nq4A0ROXnNNZUxon6MiIy+5ngB4F3gWaA45gViSEqNZ5y0RRcJLyQqKooffviBxYsXU6FCBfLmzUuFChXo1q0bJXWRyBCJiTB0qNny79fPbmuUlMjMImFZ1q1AAFBZRE4lHQsA2onIwaTv/wc0EZE27rbZDnQ+zxgisGqVCbNv1QoeeSRjO3x//AHffgtjx0LevAn83//9H2fOnOGNN96gRAod/NauXcvevXvp37//jTdbvtwkIbz9NhQrZrYlRowwuwC33urCT6konoXPiX5PQhcJ7+DChQv89NNP7N69G4C9e/fStGlT/P39qVOnjs3WeS9hYdC2Lbz5pinzp3gmmRT9xYFQoIiIRFmWlRe4DJQSkQtJ1xQCzmtHXt9GBNavh2XLTEn9Tp2cr/hz9iz8849pvlukiJmrp06dSs2aNenevfsNIT8Oh4P169ennFMVHGw8/I8/Dvfea7Ykjh2DggW1pKeSY1DRbyO6SHgeIsKJEydYsWIFR48eZceOHRQoUIB77rmHIUOGUKBAAbtNzBGcPGk8e599ZhZsxXNxIbxnMjBDRI4lfb8NSADGAAcxCb4vicht7rTXLnQ+d52tW03oT4MGxhHgTPVihwN27zba/IoPZuvWrXz11Ve89NJL3HbbjX9eEydOpFOnTlSrVu36E1e6+QYEwKBBJq7/8GFTUuz5593wEyqKvajotxFdJOwnODiYzZs3s3HjRnbt2kVCQgLVq1enS5cu3HfffeTVmpFuZ9s2mDoVZs1yrYyfkj24sWRnZeBzoDUmnyoc6CQia129tyeg87n72L3b1Pu/Iv7T8/yLmDK/hQpBjx7mWEJCArNmzSIsLIw333yTQoWubijFx8czevRo2rZtS+vWrW+84YkTMGUK9O5tjFi+HP7+21Qb0DwtxYtR0W8jukhkL+Hh4fz111/89ttv/PXXX4SGhlKnTh2aNGnCo48+SpEiRa5bGBT3M3iwCZf95BNdO70Fd4n+a+5XHigDHBSRKHfd1250Pnc/O3YY8d+8OTz5ZPpzxqpVJg/32nL7Z86cYerUqbRs2ZLHHnvsv+MiwtatW/Hz87uxsg+YhKOPPzZZxv36wd69cNNNULUq6I6v4qWo6M9Gkmr6j7r2mP4e3Y+I8O+//7Jnzx42bNjA8ePHCQ0NpWHDhpQvX54WLVpQsmRJqlevrk2xsomEBHj0UbP9PjnVwriKJ+Ju0Z9TUdGfdWzcaEr5tmsH7dunXe1HxHTybt36+utWrlzJqlWrGDhwILdek5x79uxZ9u7dS5s2bVLuhv7XX6ZT4NChJt4oMNCU+CxTxm0/n6JkFyr6bUQXCdcQEYKDgzlw4AD79+9n7969nD17ltDQUB555BGKFy9O1apVKV++PDVq1NDymTZx/DiMHGlEf+fOdlujZBQV/c6h83nWcqXaz8qVprrmvfemff3o0SYp+NoqnlFRUUyePJly5crRs2fP/9aEQ4cO8cknnzB27FiKFy9+482io2HCBGjSBNq0geHDoU8fSCFfQFE8GRX9NqKLhHOICKdOnfpP3AcHB3Pq1CmCgoLInTs33bt3p2LFigQGBlKvXj1q1qypAt9DmD/f9MAZN870u1G8DxX9zqHzefYgYrz+W7bAiy9eTd5NTmKiabxbu/bVEvxX+OOPP1iwYAGvvvoqNWvWBCA0NJQ5c+YwaNCg1HeAf/zRxBwNGWKai3TqZLr4KoqXoKLfRnSRuIqIEBgYyNGjRzl69CgBAQHExsYSExPDiRMnSExMpE6dOrRp04Y777yT/fv3U79+/RTrMSv243BAz54mfv+77zR+35tR0e8cOp9nL3FxpkHX2bPw6qupR9scO2aqcTZrlnx8HNOnTydXrly8+uqr3HTTTQCsW7eOZs2akT9//pRvePq0aejVp48J8wkONsJfQ0UVL0BFv4342iJx6dIlAgICCAgI4N9//+XcuXOICA6Hg8jISBITE7l48SIlSpSgXbt2+Pn5UbVqVS5dukTVqlXVe+8lnDtnHGGNGsHrr9ttjeIqKvqdw9fmc08hLMwUBsiXD155xZTuTE5MjCkF+uijULTo9ecOHDjAjBkz8Pf358477wTgq6++IioqihdffDHlhyYmwkcfmYTeJ5+EDz4wXcI0wVfxcFT020hOWiQSExMJDAz8T9AHBAQQFXW1MMfly5eJjY2lQYMGFCxYkKNHjxIbG4tlWTz33HPkyZOHcuXKkT9/fk2u9WLmzYPffjM9bipWtNsaxR2o6HeOnDSfeyP//gszZ5pwn27dbizzGRICY8bAyy/fGJHjcDj47LPPuHDhAgMGDKBAgQKsWbOG1atXM378+NTXpO3bTXmhnj1N05Fhw7SDr+LRqOi3EW9YJC5fvkxgYCCBgYGcO3eOc+fOERwcjMPhAP77AyJXrlzcfPPNxMXFER8fT8GCBbn33nspXbo0n376KUWLFqVatWo8+OCDWJalTa5yGJcvm/BWPz8j/HVTJuegot85vGE+9wV27IAFC+CxxyB5493ERHPu8cehWLEbxwYEBDBlyhSeeuopWrRoQXh4OHv27KFly5ap7zRHRMC770LLlqacp4gm+Coei4p+G7FjkbgSQhMSEkJISAgXLlz47+vo6Ogbrs+XLx9ly5albNmyFC5cmPLlyxMfH09AQACbNm0CoGnTphQsWJAqVapw+PBhqlWrRsWKFf+LkVRyNseOmfDW/v3NQqvkLFT0O4eKfs9BxITzbN9u4v2TFxHYtcuU3K9T58auvyLCF198wcmTJxk8eDD58+dn+fLllChRgubNm6f+0EWL4OhR4/XfuhWeftr9P5iiuIiKfhvJ7CIhIkRGRhIeHk5YWBjh4eHXfR0WFsbly5dTHJsrVy5KlChByZIl8fPzo2TJkpQsWZJixYqRL18+duzY8V/cfffu3blw4QJr164lNDSUm2++mSeeeILcuXPj5+enoTg+jsMB778PR46Ybpip5b0p3o2KfudQ0e95REdf7a/16qvXh9wfPGjOjRwJpUvfOPbUqVNMmjSJrl27ctdddzFjxgwKFChAz549U3/gsWOm3Xj9+ibB9+23028nrCjZiIp+G7EsS1auXElUVNR1n0uXLiEi/4XOJBsDQOHChSlWrBjFihWjePHi1/23WLFi14XPhIWFcerUKc6ePUvZsmUJCwujSJEiLF26FIBy5crRuHFjKlasyO7du6lYsSIVKlSgYEoZUYqC8aC99JLx8L/yit3WKFmJin7nUNHvuZw8aQR+8+amfOcVX9WlSyYyZ8wYkwicHBFh7ty5nDt3joEDB7J9+3bKly9PhQoVyJs3b8oPi483tUJz54Y77oAHHrhxO0FRbEJFv41YliWbN2+mUKFC130KFy6cbqWaK4mz58+f58KFC9SrV4/4+HgWLlxITEwMIoK/vz8nT54kPj6eiIgIypUrR40aNfDz89NKOEqmcDhMstwPP8A330DJknZbpGQ1KvqdQ0W/57N6NSxbBr17Xx92v22bCfUpWjTlypsnT55kypQpdO/enVtuuYXg4GDKlClDuXLlUn/Y77+bTmLt2kG9epDWtYqSTajoz0YsyxoNjLr22JXfo4j8Fyt/6tQpzp8/T7169QgODsayLNasWQNArVq1KFeuHLVq1WLZsmWULl2aUqVK0aRJE0SEfCm5KxTFDWzebAT/ww/Ds8/abY2SXajodw4V/d5BfDzMng3nz8Obb15N6N21y/QUGTky5cqbDoeDOXPmcOHCBfr27cuECRN44okn0o7zDw42nQnj4434T94lTFGyGRX9NmJZlowaNQoR4aWXXiIoKIjw8HAASpcuTcWKFSlUqJDGzSu2Eh0Nzz1n/vvVV6D90HwLFf3OoaLfuwgKgilT4M47oXNn4+E/c8aUGx42LHXH/LFjx5g2bRr+/v5s27aNDh06ULFixdTXaRGYNQv+/BNatDCJvopiEyr6bUQXCcXT+eYb+Pprsx3+8MN2W6PYgYp+59D53DtZuxZ+/NEk+taoYZzyZ86Yjr81a6Y8xuFwMGvWLKKioujUqROXL1+mYsWKFC5cOPUH7dkDc+aYOqJPPqkdfBVbUNFvI7pIKJ7KoUMwcCC0bWtKcer65Luo6HcOnc+9l7g4k+gbH2+6iOfPbyJzPvoIhg9PvTLZ4cOH+eijj+jSpQvffvstL7/8MnXq1En9QVFRMH48nD1rHpTUAVhRsgsV/Taii4TiacTEgL8/nDtnmmxVqWK3RYrdqOh3Dp3PvZ8TJ4zQf+ghE4IfEGDKEg8YANWqpTwmMTGR6dOnk5iYSGxsLC1btqRly5ZpP+jrr80E27mzhvso2YqKfhvRRULxFERg4UL49VezDj3yiN0WKZ6Cin7n0Pk8ZyACS5bAli1G7JcsacT/hQtw112pj/vnn3+YNWsWffv2JTw8nObNm6edj3fypMkafugh6NrV7T+HoqSEin4b0UVC8QR++snUqR40CLp00VAe5XpU9DuHzuc5i8hIk+hbsqTpR5KQAJ9+Ci+/DKmX6I9n2rRpFCxYkCNHjjBgwAAqVqyY+kMSEuB//zNJBMOGQYUKWfPDKEoSKvptRBcJxU7OnDF6alCKAAAgAElEQVSNtQoWhP/7P/NfRUmOin7n0Pk8Z/L336bEZ8+eULw4TJoEb72Vdujjzp07mTNnDrlz5+axxx6jbdu2aT9k5UrjeXnrLejUyb0/gKJcg4p+G9FFQrGDY8fMwpWYaBpHFi9ut0WKJ6Oi3zl0Ps+5JCaawjvnz0O/fqYjeY0aUL166mMuX77MpEmTKFCgAI0aNaJ169ZpPyQkxLxZ9OwJHTtqF18lS1DRbyO6SCjZSWAgvPiiKSDx8cfXd6RUlNRQ0e8cOp/nfE6fNiE/Dz4IdevC7t1w//2QVqXOLVu28PXXX1OqVCm6d+9O5cqVU79YBCZMMC8AAwZouI/idlT024guEkp2EBUFI0bA5cvwwgvQrJndFinehIp+59D53HdYsgS2boXnnzehP927Q+PGqV8fFRXFhAkTOHjwIM8++yyd0gvh2b7dxF6+9ZYm+SpuRUW/jegioWQlZ86YrejISJgxI/VGM4qSFir6nUPnc98iPNyESVatCg6HqXhWsiTkypX6mDVr1jB79mx69epFmzZt0q7uEx1tBH/HjibkJ3du9/8Qis+hoj8bsSxrNDDq2mP6e1TczcmTJv50wwZTGKJ5c7stUrwZFf3OoaLfN9m82XQu79ABypQxn1tuSf36iIgIxowZQ+HChenatSu1a9dO+wFz5sCOHTB0KKQVGqQoTqCi30Z0kVDcyY4d8PbbcPvtJnY/vbVEUZxBRb9z6Hzuu8TFmaZe8fEQFgb33mteAtJi2bJlTJo0iS5dutC3b9+0vf4HD8Jzz5l6oS++6F7jFZ9CRb+N6CKhuIM//oCBA6FECZg61Ww3K4q7UNHvHDqfK4cOwSefQNmy8NhjZi4uUCD160NDQ3n55Zd54IEHePrppymeVim1hASTRFC2rIkrSiuOSFFSQUW/jegioWSW+HgYO9bE7VerZhprpdYwRlFcQUW/c+h8roApwDN/vnHOd+xo+p/ceWfaY+bPn8/WrVt55JFH6JDeFsGGDab82vvvp90sQFFSQEW/jegioWSUsDB44w3Yv9/s9r72mjp8lKxFRb9z6HyuXEtQkHHIR0VBxYpmN/amm1K/PjAwEH9/f2655RZmzpxJvnz5Ur84NNRsJbRqBePGaRt1xWlU9NuILhKKM4jAwoWwcaPZ4X3zTahXz26rFF9BRb9z6HyupMQvv8B330HDhtCrV9qdz0WEadOmsWfPHkaOHEm1atXSuhiGDIG9e+Gnn3SrV3EKnxP9lmXlAoYALwFlgUPAWBFZnM64osAbQHugFpAb2A9MFJEfM2mLLhJKqkRGmsTc7dtNbf3x46FoUbutUnwNFf3OofO5khqRkaapV0wMPP44NGmS9g7tqVOnGDZsGA0aNGDgwIHkSuvio0dNPf9Ro8ybhaKkgS+K/neBt4DhwE6gC9Ab6CgiK9IYVw/4HZgLbAAcQFfgBeBVEfkkE7boIqFch8MBn38O27YZr/5rr0GjRnZbpfgyKvqdQ+dzJT127YIPPjBfv/de2hU4RYR33nmHNWvWsHDhQqqkFb+fmAi9e0NwMHz7bdrZw4pP41Oi37Ks0sAp4H8iMuqa46uBUiLSII2xhQARkehkx1cDNUSkYibs0UVCAeDIEROaefiwCdMcMQKKFLHbKkVR0e8sOp8rzpCYaJol/vyzcfCUK5d2SP4///zD66+/zttvv51+Q69ly2D4cFiwABqkKmcUH8bXRH934Augpogcuea4P/A5UFVETmTwnu8DA0UkTybs0UXChzl/HkaPhogI09Bl6FAoVcpuqxTlelT0O4fO50pGOHPGOHcaN4aHH067EI/D4WDUqFFcvHiRIUOGUKFChdQvjokxi0nFivD665Anw9JEycFkdj731nohtwGxwNFkx/cl/bduJu55L3DQFaMU3yEyEqZPh06doEsXaN0avvwSJk9Wwa8oiuIr3HorzJ1rPP3PPw8TJphdgJTIlSsX48aNw9/fn06dOjFhwgRSfcEsUAA+/NDUcm7c2DRyURQX8VZP/2fAoyJSNtnx6sAR4HkR+TID9+sDfAp0E5GFKZwfDYxKfvxavPH3qGSM8HAj6n//HW67zSRyPfywVllTvAP19DuHevqVzBIVZXqt5MtnSnuWL5/6tQ6Hg4EDB5KYmMiAAQOonFZiQEgIPPUUtGihpT0VwMvDeyzLegD4zYlL14tIa3eKfsuyWgMrgW9E5IWMWf7fPXSRyKEEBcHSpbB2rfn+nnugT5+06zQriieiot85dD5XXGXPHuMgql/flGZOa704ePAgEydOpGLFiowePTrtG69cCf/3f6ZxQKVKbrVZ8S68XfQXBJxJoI0WkYCk+PvXgQLXzs6WZTUFtmMq+Cx34rlNgNXAJuAxEYnPpP26SOQgjh41Zdn27IFChaB/f+jQQZtnKd6Nin7n0PlccQeJifDuu7Bihan007Jl6teKCCNGjOD333/n888/57bbbkv94shIs81cpQp8+qkuTD6KV4v+jGJZ1vPAfEy1naPXHO+BKcWZbiKvZVn1gfXAP8CDIhLjgj26SHgxiYmmAsOcOVC6NJQoAf7+UKeO7qIqOQcV/c6h87niTs6cgX794P774eWX0/b6BwQEMGjQIJo2bcqAAQPSrvAzYwbMmwfff2+SfRWfwtdEf2ngNPCuiIy55vjvQBkRqZ/O+BrAxqR73C8iES7ao4uEl3H4MKxaZXKjwsNNEtYbb0CtWnZbpihZg4p+59D5XMkKFi82n7v/n737Dq+qyP84/p70QkhCCYQuTapYEFEsiODaEUHsCgqLPxd7B12qihULlhUriopiwbIWQGQRsIAi0hSkJiQQCAkhPbnz++OEXYwJ3IQk55bP63nuk+Tk3JvvxTjzyZw5MyfC6NEHH1CaPXs2r776KnfccQenn3565Sfm5jqbwPTo4bxoaGjNFy4+KahCP4AxZgrOzrpjgJ+AS3B2573AWvvpAefNB1pba9uXfZ0E/ADUB64Cdpd76Z+ttYVVrEWdhI/buxfmz3cGRnbtgoYNYcwYZ8Mszc+XYKDQ7x2151Jbioqc1TdXrnSuLHfuXPm52dnZXH755YSHhzNz5kxiY2MrP/ntt2HKFGeu/4ABNV+4+JxgDP2hwL04u/A2BX4DJlprZ5c77xugjbW2TdnXfYEFB3npI6y1m6tYizoJH5OV5UzZ+egjSE93Fj3o3NmZm9+okdvVidQ9hX7vqD2X2rZhg3Ojb1KSs8fLwUb9P//8c958802GDBnCoEGDKj9x3z647DLncvXkyRAVVeN1i+8IutDvS9RJuG/nTueGqc8/dzbI2r97+SWXwNFH614nEYV+76g9l7ryzTdw770wYgRcd13l55WWlvLCCy8wc+ZMXn75ZTof7BLBqlUwcaKznvSwYTVdsvgIhX4XqZOoW6WlsHq1s17+hg3O9Jx166BbNyfk9+ypkC9SnkK/d9SeS10qLoYbboC1a52r0wkJlZ+7efNmrr76arp27cqzzz5LSGUdncfj3D38yy8wYwa0b187xYtrFPpdpE6idqWkwHffwbffQkmJMx+ycWNnzfz/+z9nIxSFfJGDU+j3jtpzccP27c6of48ezl4w9epVfu7ChQt58sknufvuu+ndu3flJ27c6IT/c891/rJQRxkwFPpdpE6iZlgLmzY56+P//DPs3g1r1kBMDJx0kjOK3769ltEUqQ6Ffu+oPRc3ffWVE/7PPRcmTKi8vyspKeHmm29mxYoVTJ8+nS5dulT+ojNnwpw5cNVVcP75tVO41CmF/jpkjBkPjDvwmP4dq6agwJmis2KFMxixc6fzsV49OPtsZ6pOt27OvUi6H0nk8Cn0e0ehX9xmrbOx17JlcMcdzkIUldm0aRPDhg2jf//+3H777cTExFR8Yl4eXH017NnjTPlp3rx2ipc6odDvInUSlSssdNbEX7MG1q93wn5KirOiTny8cwNTVJSzo3iDBge/pCki1afQ7x215+Irioqce3LXrnVW+el+kB2IlixZwnPPPUfnzp0ZM2ZM5Rt7rVgBzz4LTZs6lxI05ccvKfS7KNg7CWudqTgbNjgBf9Mmp7EC+O035/Jko0Zw553OUpphYU7Ij493t26RYKLQ751gb8/F92zZ4qzuExUFr7/u7DNTEWstEyZMYM6cOUyYMIELLrig8hddvBhuvx1GjXK2oBe/otDvomDoJDweSEtzgv0ff8C2bc4xcDYFTEx0jsfEOFcNL7kE4uKc3W7bttVggojbFPq9EwztufinxYudmTmtWzt5PTKy4vNyc3O54447KCgoYMqUKTRp0qTiE0tK4PrrnUsJH3zgrHctfkGh30WB0EkUF0NqqhPm9z/y8pzvWessk5mU5AT5evWcFXRCQ6FTJ+eSY7t2zjmaniPimxT6vRMI7bkEtu+/dxbj6dcPHn648kG19PR07r//ftLS0nj55ZcrD/8ZGTB2LDRrBjfeWPmlBPEZCv0u8vVOorTU+X86JQW2bnUC/e7df14VIDzc+SM/KsqZmjNggDNNZ+tWZ3Q/MtK56b9JE2e5TK2gI+JfFPq94+vtuch+Tz8Ns2fDbbfBwIGV98tLlizhjjvu4KSTTuLBBx8kIiKi4hN/+MEJ/T16wLRpUNl54jqFfhe50UlY60ydSU//8yMr66/nhoY6c+pbtIDkZGfE/tdfnek42dnOKl7WwjvvOH/ot28Pxx7r/AGgcC8SGBT6vaPQL/7EWif8f/89nHiik9krs3LlSiZNmkSvXr247bbbCA0NrfjE1193/poYMgSuuaZ2CpfDotDvosPtJDweJ3zv3g27dv354969lT8vIcG5Af/AR3i4M5Jfrx7s2OFM0Zk71wnvHTrAEUdAly7ODbft2zsr5ijYiwQ+hX7vKPSLPyopcdb3nz8fRo+Ga6+t/NznnnuO6dOnM2jQIO67776Kd/a1Fl59Fd5/H668Ei67rPaKlypT6HeRMcZu22bZu9cJ7+U/7tv3v5teK36+E+AbNnRG5A/8WL/+/0K5x+MseRkS4mzgsWWL88fB3//ujNpv2uTs6te6tbOZVYMGzvNFRBT6vaPQL/6spMTJ6l99BcOGOZt8VcTj8fDQQw/x+++/07dvX4YNG1bxMp+FhfB//+dsrPPYY3DKKbVav3hHod9Fxhj74ouW+vWdZSjLf6xX79Cr1xQUOAF+9WonuDdq5Nw026qVc7e+tRAbC1dc4Uyz27YNWrZ0ztNIvYgcikK/dxT6JRB4PM7a/p9+6gzS33FHxVnB4/EwZ84cJk+ezJAhQ7j77rsrHvnPzobp051VPP7xDzjhhFp/D1I5hX4XVdZJFBU5YX3lSmfH2Z07nd1mU1Jg+XIn3ANceKET+uPjnTn5zZo5U3W0E62I1BSFfu8o9EsgKS11wv+WLc79ubfeWvEgpMfjYeLEiXz11Vfce++9nHvuuRWH/9xcZ55/Who8+qgzrUDqnEJ/HTLGjAfGHXhswQJLYiJ8+KET9JOSoG9fZ7Wb775zvk5KckbnK7t3RkSktij0e0ehXwKRtfDll/Dgg85S2088UfES29Zapk2bxrJly2jatCmTJ08mPDz8ryempTk3DzRt6nzs3Ln234T8l0K/i9RJiIivU+j3jtpzCXSzZsGTTzozD665xrkPsLz94X/GjBmccMIJPP7440RWtBtYSgq89JKzYsj998NZZ9X+GxCFfjepkxARX6fQ7x215xIsVq2C115zphuPHQv9+1d83g8//MCTTz7JEUccwahRo2jVqtVfT8rMdEb8d+925hOdeGJtlh70FPpdpE5CRHydQr931J5LsElPdzb42rTJWfbz3HMrnob8/fffM2bMGEpLS5k0aRKnVLSST3ExTJrkLCnYu7dz0++hVjKRKlPod5E6CRHxdQr93lF7LsFq/9L8S5fCnj3w+OMVT/1JTU3lpZdeYs2aNQwZMoTBgwf/9abf4mIYN85ZO/Rvf3M+1w6/NUah30XqJETE1yn0e0ftuQh8+60zS8cYZ6r+ySf/dcDe4/HwxBNPMGvWLI488kgeffRRkpOT//piS5bAs886a4yPHu3sFCqHRaHfReokRMTXKfR7R+25yP8UFjo3/W7YADk58MAD0K7dX8+bO3cun376KZmZmQwfPpx+/fr99aSff4Z77nGW/RwzBs45p/bfQIBS6HeROgkR8XUK/d5Rey5SsUWLnNBvLVx3HQwaBOVX8ywsLGTs2LEsWrSIDh068Mgjj9CsWbM/n7RnD7z1lvOCXbs68/4bNKi7NxIAFPpdpE5CRHydQr931J6LHJzH48z9//FH+O03uP56GDr0rzv+zps3jy+//JK0tDTOOussLrvsMkLL3yH8wQfOpYSQEJgwAU47re7eiB9T6HeROgkR8XUK/d5Rey7ivYwMmDgR1q+Hjh2dQfsjj/zzOSUlJTz66KN8/PHHxMbGMm7cuL+u/JOR4Wwg8O23zi6mN93kfJQKKfS7SJ2EiPg6hX7vqD0XqZ7Nm2HmTPj0U2jYEG6+GQYM+PM5qampfPLJJyxatIiIiAhGjRpF7969/3zS5587Wwbn5DiXES67DCraGCyIKfTXIWPMeGDcgcf07ygivkyh3zsK/SKH7/ff4cUXIS3NGbA/6ii49NI/rwC0cuVKHnzwQbZs2cIZZ5zBBRdcQK9evf53QkkJfP01vPkmZGfDhRc6Wwhr3X+FfjepkxARX6fQ7x215yI1Kz0dHn4Yvv8e4uJg+HA4+2yIj//fOWlpaUyfPp3vvvsOay3XX389559//v/W/9+zB6ZMcW7+7dkT+vSBiy8O2j8AFPpdpE5CRHydQr931J6L1J7iYnj/ffjpJ1i8GJo2hWuvdVbv3H8j8KZNm3jkkUcoLS0lPz+f0047jcGDB5OYmOickJcH774L773nTAG66CIYNQqio917Y3VMod9F6iRExNcp9HtH7blI3fnlF3jpJSgqcjbsLS6GESOcwXxwbgJ++eWXmTVrFoWFhfTp04eLLrqIE044AWOME/qnTYPUVOdqQLt2zlJC3bq5+8ZqmUK/i9RJiIivU+j3jtpzEXd4PPDFFzBjBqSkwLHHQvPmziqevXo5M3m2bdvGjBkz+OWXX0hNTaVTp06MGDGCE0880dlAYO5ceP552LkTevSA00+HgQOdvygCiEK/i9RJiIivC/TQb4xpB7wOJAG5wEhr7bJqvI7acxEf4PHAmjXOMv6//eYcO/lkZz+v886D+vU9LFiwgMWLF7N+/XpCQ0MpLi5m8ODBnH/++YQXFsLs2TBvnrO0UHy8cyPwRRdBWJir7+1wKfS7SJ2EiPi6IAj9c4F3rbXTjTEDgGlAp6o2zmrPRXxXVpZzNeCjj2D7dudYr17Oo39/y/r1P/DWW2+RkpJCTEwMjRo1IiIigiFDhnBcw4aEzJ/v3AwcEuJcTjj1VLjiCmdakB8JutBvjAkB7gZGAU2B34CJ1tr3q/g6bYFVQDTQwVq7oRq1qJMQEZ/ma6HfGNMCpw3vCfTAaYOPsNZuruDclsBUYABggHnALdbarWXfbwxsBBpYa4vLjv0OXF7V0X615yL+JT3dmdWzeLFzZQCgWTPnvoCePfNZvfoTFi/+iIyMDNq2bUtMTAwFBQX069WLs7Ozqbd6NRQUOFcCduyAvn3h/POhVStX39fBBGPofwC4AxgLLAcuBUYC51lr/12F1/kCp8NpikK/iAQoHwz9fYFZOO13KHAmFYR+Y0wM8AtQCNwHWGAyEAMcZa3NNcYcB7xtre14wPO+Al6w1n5QxbrUnov4ueJiWLIE1q2DDz+E3FyIiYE2bSA5uZTff99Ebu5CMjM/o2fP1uzbt4/09HSSIyM5vrSUswoKaNm6tbMp2KpVcPTRcMop0L8/1Kvn9tsLrtBvjEkCtgFTrLXjDjg+H2hsrT3Ky9e5HGf06KGyjwr9IhKQfDD0h1hrPWWfjwCmU3Hovxl4Ajhyf/tsjDkCWA/cZa19QqFfRLy1cyf8+9+wYoVzr0CnTs4iQCkppWzfvpfCwp20bv0H+fkbiYpKZW/maiLyNtA9PIyuycl0DwujYU4OR3o8hDVu7Kw32ro1dO8OTZrUyd4BwRb6rwJmAB2ttesPOD4ceAVoa63ddIjXSATWAnfijDK9ikK/iAQoXwv9BzpE6J8PRFlr+5Q7vhDAWnuapveISE2w1pnqn5YG69fDxo3OpmI5OZZ69fIpKtqDtXv45RdDcXExoaaItk3WU690N3+klxBevIP6Ybs4u0UmEZGGBnElNMpJpVmzJI448kiSBgwgunVr57JDs2YQFVWtOqvbnvvr7ctdcS71lg/oq8s+dgEOGvqBR4B11to3jDHDarY8ERGpIV2BORUcXw1cDGCtzTDG/AAMA/bfyGtwpg6JiHjFGGjZ0nn06vWn7+DMKIwBmv/3aGEh7N7di717YflySE+3bN6ci6m/j982FrJpUz7bdxuKU4po+GsGhe8UEeJZT/qeJLDbiQrJJyE2k3omjz9yOjp/dVhLu/gN5JXGsqswgaKSGAyQHLmNopBImjXIqvb789fQ3wDIqmA4JvOA71fKGHMKcDVwjDc/zBgzHhh3iHO8eSkREamaBsCeCo5nAokHfH098Lox5k4gD7iioiF7teci4ob0Co7tLYWdFWT4nypo8bLznI8b91W/Bp8I/caY/sBcL05daK3te5g/KwL4FzDVWrvGm+dYa8cD4w/ymtaHL5vXWm2H+9pVfX5Vzj/UuYfz/Yq+p9+Bunu+t8/x5ryq/nc+1Pd8/ffA7RpqU9lUz5O8OG88as9r/LXVnteNQPodqMpz1J7/WXXbc58I/cASoLMX55X9ncMeIMH8dfLl/hH+TCp3C87o0NPGmISyYzFlH+OMMXHW2hwv6xYRkdq1hz+P6O9X2RUAERGpgE+EfmttHrCuCk9ZDUQC7fjzvP4uZR8PNoLfBWd5ztQKvvcTztJwR1ehFhERqT2rceb1l9eFg7f1IiJygNpfV6h2fAEUA1eUO34lsOoQK/dMAU4v93j4gOePqNlSXTfBh1+7qs+vyvmHOvdwvl+b/6a1IZB+B6ryHG/Oq+5/Z3/7HQD/rBngY6B32UaKABhj2gB9yr4XTALp/2W159UTSL8DVXmO2vM/q1bNfrlkJ4AxZgrOVJ0xOCP0l+DsznuBtfbTA86bD7S21rY/yGsN4/CX7PTJeV9SN/Q7IKDfg6oyxgwp+/QMnBtxbwAygAxr7cKyc2JxrsDm87/NuSYBcTibcx3GbW0V1qT/hkFOvwMCgfl74BPTe6ppLLAPuBlnus5vwNADA3+ZUPz7fYqIBKr3yn39XNnHhUBfAOvsuNsPZwPFN3DWzpsP3FLTgV9EJJD57Ui/LwnEvwalavQ7IKDfg0Cg/4ai3wGBwPw98Nc5/b7GH+eDSc3S74CAfg8Cgf4bin4HBALw90Aj/SIiIiIiAU4j/SIiIiIiAU6hX0REREQkwCn0i4iIiIgEOIV+EREREZEAp9AvIiIiIhLgFPpFRERERAKcQr+IiIiISIBT6BcREb9njGlnjPnWGPO7MeZnY0xPt2sSEfElCv0iIhIIXgBet9Z2BO4CZhpjjMs1iYj4DIV+ERGpc8aYFsaYZ4wxS40xecYYa4xpU8m5LY0xs40x2caYvcaYD4wxrQ74fmOgN/AagLV2LmCA42r9jYiI+AmFfhERcUN7YCiwB1hU2UnGmBjga6ATcA1wFdABWGCMiS07rRWQZq0tPuCpm8uOi4gIEOZ2ASIiEpT+Y61tAmCMGQGcWcl5I4G2wJHW2g1l568E1gOjgCfqoFYREb+nkX4REalz1lqPl6deAHy3P/CXPXcTsBgYWHZoK5BsjAk/4Hltyo6LiAga6a8WY8x4YJzbdYiIVIW11h9vbO0KzKng+GrgYgBrbYYx5gdgGDDdGDMAZ07/8vJPUvstIoGgOu25Qn81WGvHA+P3f22MsdZa1+oRETkUP17IpgHOvP/yMoHEA76+HnjdGHMnkAdcUVHDXL79Lk/tuYj4uuq25wr9IiLi96y164GT3K5DRMRXaU6/iIj4sj38eUR/v8quAIiISAUU+kVExJetxpnXX14XYM3hvrgxZnzZHgHWGKN5PSISsBT6q0GdhIhInfkY6G2Mabv/QNkmXn3KvndYrLXjrbVm/+NwX09ExFcZ3bB0+HTjl4j4OmOMz63eY4wZUvbpGTg34t4AZAAZ1tqFZefEAr8A+cB9gAUmAXHAUdbafTVck9pzEfFp1W3PFfprgDoJEfF1Phr6K2s4F1pr+x5wXitgKrB/Kc75wC3W2s21UZPacxHxZQr9LlInISK+zhdDvy9Sey4ivq667bnm9IuISNDSPVoiEiw00l8NFe3oqH9HEfFlGun3jkb6RcTXaXqPi9RJiIivU+j3jtpzEfF1mt4jIiIiIiIVUugXEREREQlwCv0iIiIiIgFOoV9EREREJMAp9IuISNDSkp0iEiy0ek8N0GoPIuLrtHqPd9Sei4iv0+o9dUgjQyIiIiLiTzTSXwM0MiQivk4j/d4xxtjce+4hplkzOPJI6NgRWrWCEI2RiYhvqG57HlYbxYiIiPirJ2ISKNq8GZYvh927aVpaSvfGjemelERC8+bQvTscdRQ0bux2qSIiXtNIfw3QSL+I+DqN9HvHGGMnTbIUFztfW2uJj99LePg69u37joLMbbBzJ2G7dnFkdDRHJSXRoVEjwtq3h2OOgS5dIEzjaSJSe6rbniv01wCFfhHxdQr93infnns8kJICv/4Kq1ZBbi4YA6GhpdSvn4q1v5CV9ROeXTux6em0LCqiZ5MmdG/ShIiOHaFnT+jUCUJDXXxXIhJIFPpdpNAvIr5Ood873rbnhYWwdi2sXAkbNzp/HHg8loSEPYSFrSBrz2JKd6VDWhpHlJRwfNOmdElKIvSYY+CkkyA5uQ7ejYgEIoV+Fyn0i4ivU+j3zuG059bCluRs+2MAACAASURBVC2wbBmsXg0lJQCWxMRdwDKyMpdid6QRlZ7O8XFx9G7enPpt2jh/BBx1lK4GiIhXFPpdpNAvIr5Ood87Nd2eezzwxx/w44+wbp3zdXh4EfHx68jJmU/h7i2YlBS6Wkvfli1p2rIlnHYaHH207g0QkQop9LtIoV9EfJ1Cv3fqoj3ft8/5I+D77517BMCSkLCVwsK55O1ei9m2je4eD31btiSpdWs480zo3Nm5mUBEgp5Cv4sU+kXE1yn0e8eN9tzjcW4U/uYbyMgAYywJCVsoKPiS3J1riEpJ4bToaE5q2ZLwE06Afv2gfv06rVFEfIdCv4sU+kXE1yn0V8wYMx4Yd+Axt9vz0lLnBuGFC50/AkJCiomJWcG+vf+G9M0csXcvZzdvTvNWrWDAAOd+AF0FEAkaCv0uUugXEV+n0O8dX2zPCwth8WL49lvn89DQXYSEfEHx3p+ITEmhb1QUJzZvTnjv3s5UoOhot0sWkVqk0O8iX+wkREQOpNDvHX9oz9PT4bPPnBuEw8OLadBgObt3fU7xlk30Ki7mrObNie7dG84+G2Jj3S5XRGqYQr+L/KGTEJHgptDvHX9rz3NzYe5c58ZgsDRq9DuZu9+jJPUPehYUcHbz5sQcdxyce67uAxAJEAr9dcgX54CKiByMQr93/C30H6i0FJYsga+/hoICS3T0JgoL3sWzYz3H5OVxbrNmxB5/PAwcqClAIn5Mod9F/txJiEhwUOj3TqC059bCmjUwZw7s3Wtp3Hgje7PfoXDbBk4vLqZ/y5aEnnMOnHyybgIW8TMK/S4KlE5CRAKXQr93ArE9t9bZE+Dzz6GkxNKgwc9kpL9LbFoKQ6Ki6NyuHQwdCu3auV2qiHhBod9FgdhJiEhgUej3TqC35yUlMH8+/Oc/YG0h0dGfs2/XAo7YtYtLGjYksVcvGDIEoqLcLlVEKqHQ76JA7yRExP8p9HsnmNrzvDyYPRtWrYKGDTPIz3+dwm1rudBaerVpg7niCo3+i/gghX4XBVMnISL+SaHfO8Hanq9dC+++C4WFpcTHL2D39k/olJnJpQ0bEtOvn7P6T2io22WKCAr9rjLGWFtSogZRRHyWQr93gjX071dYCB99BD/9BLGx2ynIf4nIHZu5IjSU9t27w7BhWvpTxGUK/VVkjEkE3gA6AvnADuAGa+2GaryWXXDNNZzQtCnRiYlw0klw/PGaEykiPkOh3zvBHvoPtG4dvPUWGFNIdPS7ZG76lgtLSjipdWu49lpo0cLtEkWCkkJ/FRljEoCe1tp5ZV/fBFxkre1bjdeyEyb8TFHRfEzRDiLS0+kFnNK0KTFNmkC/ftCjB4SE1OybEBHxkkJ/xbTvyqFlZsLrr0NamoeGDb9h59b3OSUvjwuSkgi57DI4+mi3SxQJKn4f+o0xLYC7gZ5ADyAaOMJau7mCc1sCU4EBgAHmAbdYa7cexs/vCcy21rapxnPtd99ZFi2CrCzweEqIifmdgoJ/Y/O2k5iRwRmRkfRo3JiQDh1gwABo2bK6pYqIVJlCv3c00l+54mL48EP47jtLUtIadu98hU5Ze7gqPp6Iyy93rnCLSK0LhNDfF5gFLAdCgTOpIPQbY2KAX4BC4D7AApOBGOAoa21uNX/+m8Bua+3N1XjunzqJ0lL4+WdnV8TsbDAmj+jo78jL/QaydtFh3z4GJCTQPCEBTjwRTj9dU4FEpFYp9HtHof/QrIUFC+CzzyA5eRs706fRfs8uromNJfKyy5x+TURqTSCE/hBrrafs8xHAdCoO/TcDTwBH7p9/b4w5AlgP3GWtfaLs2DygsmuOA621iw94zXHA34D+1tq8atR+0E4iJwe++QZ+/NHZGCUmZifWfkFB3joiduzgDGM4sXFjQtu3d1ZISE6uagkiIgel0F8xTe85PN99B++/DwkJqezZ/QwdcnYzLCqKyMsvV/gXqSV+H/oPdIjQPx+Istb2KXd8IYC19rQq/qz7gPOBM6212dWst0ojQ5s2wVdfwdatAMXExy8jO/sLbOZOji0s5G+JicQ1bgxnngnHHKMt0kXksCn0e0cj/dWzciW8+SbEx6eStecZjtybybDYWMJGjICuXd0uTySgBFPoTwfmWGtHlTv+HHCxtbZxFX7OOOAcDhH4KxoJKq+6/44FBc5l0iVLnKsAiYkbycv7kKKcNJpnZ3NhdDTNGzaEv/0NTjhBNwOLSLUo9HtHof/w/PyzE/6TkraQlvIEpxYXMighAfN//wetW7tdnkhACKbQXwQ8Ya29p9zxycA91towL39GV2AV8Aewr+xwibW2ZzXqrZFOwlpYvhy+/BJycyEubg/WfkjOnrUkZ2VxUWQkLRo2hP79nWVBtS+AiHhJod87Cv01Y+lSZ7OvRo1WsiPlBQYbOC05GW66CeLj3S5PxK8p9Fcx9Nek2uokNm92NknZsQPq188CPmBv5lqSc3K4KDycFomJcN55zooJmgIkIgeh0O8dhf6aNXcufPqppXnz+aRvnM0NQPuePWH4cA1ciVRTMIX+HcBHNTG9p6bURSexebOzVNqOHZCQkIW175O1aw3t8/IYEhVFYpMmMHgwdOhQq3WIiH9S6K+YbuStfaWl8MYbsHp1CTExrxC+6xdu9niIGzrUWb1ORKokmEL/10CEtfbkcse/wXk/VbqRt5r1jcfFTmLTJucPgJ07oXHjdHJzZ5K7eyu9i4s5JyaGyLZt4dJLoUGDOqtJRHybQr93NNJfe/btg2nTYN++vRQWPEL3giyuionB3HEHJCW5XZ6I3wim0H8L8BjQ0Vq7sexYG5wlO++x1j5ep8XibiexcqUzBSg/35KcvI6dO9/BZGdyIXBsTAzmjDPgjDN0GVUkyCn0e0ehv/Zt2gRPPw2tWq1l07onuSEsjE7HHutM+dFiFSKHFBCh3xgzpOzTM4DrgRuADCDDWruw7JxYnM258vnf5lyTgDiczbn2lX/dOqjb9U6itNTZDOzrryEiopjExK/YnvoNrQsLuTw0lMSkJLj8cq2eIBKkFPq94wvtebCYMwcWLiylfv1Xidn9CzeXlhJ5ww3QrZvbpYn4tEAJ/ZUVs9Ba2/eA81oBU4EBgAHmA7eUvypQV3ytk9i7F2bNgvXroVmzDHJyXicvM5WzreWU8HBMv37OEqAaUREJGgr93vG19jzQ5ebCk09CUVEGmbsmMyQMTktKgttug+hot8sT8UkBEfr9lS93Er/84uyW6PF4SE5eyOZNn9GutJQrrCWuVSu45hpo2NDtMkWklin0e8eX2/NAtmYNPP+8JSnpS/J3fMqdRUUkjhzprE4nIn+i0O8if+gk8vKc0f/Vq6FNmzR27pwOOXu4AjiyXj0n/Ldv73aZIlJLFPq94w/teaDyeOCll2Dbtlxy9k7gvJBi+jdp4oz6R0S4XZ6Iz1Dor0Nur95zOKyF77935lLWq1dERMTb7Ej5iQvDwuhjDObCC6FPH637LxJgFPor5s/teaD64w+YOtXSvPmXFKR+yt0lJcSMHq25/iJlFPpd5K8jQzt3wiuvQE6Oh2bN5rJxwxf0CQ9nYHExoSeeCBddBGF1vteZiNQChX7v+Gt7Hmg8HvjXvyA1NYs9meMYFhXK8W3awOjRuh9Ngp5Cv4v8vZPIzYUZM5xl1Nq3/4X169+kfUgIVxUVEdOunbOMWmys22WKyGFQ6PeOv7fngWb9enjqKUvjxu8Ss2cpt5WUEPrPf2pdfwlqCv0uCpROoqQEPvjAmf7ToUMKKSkvEldcyIiiIhq2bg0jRyr8i/gphX7vBEp7HkhKSuDxx6G0NIXUTRO5LzyM5KFDoW9ft0sTcYVCv4sCrZOwFr75Bj77DFq0yGbPnmlE52czqqSExP3hPybG7TJFpAoU+r0TaO15IFmwAD7+uAjsJM6LLOKMRo2cm3y1+aQEGYV+FwVyJ7FyJbz5JjRrlk1m5lPUK8xlVEkJ8W3awHXXKfyL+AmFfu8EcnseCDIy4IEHIDHxU6L3LuD2ggJCx43TdB8JKgr9LgqGTmLFCpg5E1q2zCQj42kSSov4e1ERcW3awIgREBXldokichAK/d4Jhvbc33k88NxzsGvXdnZtn8C4EEPjkSPhuOPcLk2kTij0uyiYOokff4R33oG2bXexfftTNMbDqNxconv3hqFDtaqCiI9S6PdOMLXn/m7hQpgzp5CiwjGMiIKjjzsOLr/c7bJEap1Cfx3Sus6weDHMng0dO6azceMTHB0TzWV79hAyeDCcdprb5YlIOQr93lHo9y8pKfDQQ5b69V/gOLONIZGRMHaslpuWgKbQ76Jg7iS+/trZ6Ouoo1ay8peXGNK4MadkZTk3+3bq5HZ5IlJGod87wdye+6uCApg4ESIj/0PE3jncXVhIyAMPQHy826WJ1AqFfhcFeydRWurc7Lt6taV160/YuP5rro+Pp4O1cMst0KCB2yWKBD2Ffu8Ee3vur6yFl16C7dtT2JM2mQesJXbMGGjd2u3SRGqcQr+L1Ek4cnNh2jQoKiomMnI6+bu3cGtpKfW7dYOrr9Z8fxEXKfR7R+25f/vqK5g3L4e9e+5kQngITa69Fnr2dLsskRql0O8idRJ/lprqhP+mTbPZvn0Kx8TFcsmOHZjrroOjj3a7PJGgpNBfMd2jFXhWrYIXXijCUzqGm2Kh06mnwgUXuF2WSI1R6HeRQn/FfvjBmfbTo8eP/LJiBv9o3JgjS0udzVQ011KkTin0e0fteWBIT4cJEzxERj7KRTHZnNqxIwwb5nZZIjVCod9F6iQqV1oKL74I6eklREf/C/alcXN+PtEnnwyDBrldnkjQUOj3jtrzwJGXB/fcY4mLe41ent8ZmJwMN93kdlkih02h30XqJA5t+3aYOhXatctg3boHuaBJEv127YK77oImTdwuTyTgKfR7R+15YCkuhvvug6ioj+hY8iNXREc7S3oa/a8g/kuh30XqJLw3dy589pmlffuP2bF1MXcBcT16OBuqqBEWqTUK/d5Rex54PB6YPBk8nm9ILvqKUQCTJkFoqNuliVSLQr+L1ElUTV4ePPYY1K+fzdatEzmrWTJnpqTAHXdAixZulycSkBT6vaP2PDBZC08/Dbt3Lyc+/x1uMwbz0EMK/uKXFPpdpE6iepYsgVmzoFOnT9n6xzfcExpKfM+ecPHFbpcmEnACOfQbYxKBN4COQD6wA7jBWruhGq+l9jyATZ8OaWmriMp+lbvCwuDBBxX8xe8o9LtInUT15efDo49C/fo5bN78T4a0aM7JO3c6cy61wo9IjQnw0J8A9LTWziv7+ibgImtt32q8ltrzAPfSS5CaupqYnFe5MyQENOIvfqa67bl2SxJXRUfDP/8JvXrF4fE8wY+hzZgSGUnx/ffDokVulyci1WCMaWGMecYYs9QYk2eMscaYNpWc29IYM9sYk22M2WuM+cAY06oqP89am7U/8JdZAlT480RGjICWLbuSEzOMx6yFe+5xlpoTCXAK/eITTjoJHnjAsGPH5TRsdT23lpayYfFiePhhNcYi/qc9MBTYA1T617sxJgb4GugEXANcBXQAFhhjYg/j598CzDmM50uAu/ZaaNeuG1mR1/B4SAjce69zx69IAFPorwZjzPiykStrjNF14BoSFwdTpkBCQnMSEp9mZqFhJsAttzg7rYiIv/iPtbaJtfYc4L2DnDcSaAtcaK39yFo7B7gAaA3OIisAxph5xphdlTz6HPiCxphxZa95b42/Kwko11wDHTt2IzPsKqYCjBvn3PErEqA0p78GaA5ozdu6FR55BLp3/5YtG+Zwf2kp0eeeC2ec4XZpIn7JrTn9xpgRwHTgCGvt5nLfmw9EWWvLB/eFANba06r4s+4DzgfOtNZmV7NetedB5tVXYfPm5TQvfI+/JybC3Xe7XZLIQWlOvwSUVq3gqacgO/tkklvdwx0FBaxftsy561eXYEUCRVdgVQXHVwNdqvJCZSP8hwz85a/Uln9UqXoJCMOHQ8OGx7Et/Exm7d0Lzz3ndkkitUKhX3xWaKizYW+HDg2Jq/8Mb+SW8n5oKNx6K2RluV2eiBy+Bjjz/svLBBK9fRFjTFdgPNAQWGiMWWGMWVbRudba8dZaU9mjGu9BAsBNN0FUVD9+Lj2WL9PS4M033S5JpMYp9IvPO+ss+Mc/Qtm9eww7o7vxQL16lN53H6xZ43ZpIuIDrLWry0J7e2vt0WWPnm7XJf5lzBgICxvMl7mtWbpuHXzyidslidQohX7xCy1bwtSpkJZ2Jo1a/Z1bjSHrrbfgo4/cLk1Eqm8PFY/oV3YFQKTWGAMTJoAx1zEzK4E1CxbA8uVulyVSYxT6xW9ERMDEiRAV1ZpmzR/m/uy9rN240flrQDfeifij1Tjz+svrAtTJpTytxiYHCg2Fhx4yWHsbz+RFsOPll52VJUQCgEK/+J1rroHTToshIuJJ3s4q5AtwrssWF7tdmohUzcdAb2NM2/0Hyjbx6lP2vVpXfo5/XfxM8W0RETB5cggeO5FxxYb8Bx6AvXvdLkvksGnJzhqgJd7ckZoKDz4Ibdu+TWLJVq5NSYFJkyAhwe3SRHxOXS/ZaYwZUvbpGcD1wA1ABpBhrV1Ydk4s8AuQD9wHWGASEAccZa3dV1f1HlC32nMBYNMmePTRHKLs3TwWFkrIE09AeLjbZYlUuz1X6K8B6iTck58PY8dChw6L2LP9K+7JzCTkzjuhTRu3SxPxKS6E/soaxYXW2r4HnNcKmAoMAAwwH7il/Jr+dUXtuRzo++/hnXdSSQqdyr3R0c7AkojLFPpdpE7CXR4PTJ4MTZpsYP26Z5kcEkLUJZdAr15ulybiM9zanMvfqD2X8j78EJYs+Ykeoe9y5ZFHOgv7i7hIm3NJ0AoJgX/+E6xtT8s293F7fj67//1v+Pxzt0sTER+nG3nlUAYNguTkY/mp6Ch+XLMGli51uySRatFIfw3QyJDvmDMHfv21gPT0O7mzWTKtO3SAiy92uywR12mk3ztqz6Ui1sJtt1ny88YxMWw3SWPHQrNmbpclQUoj/XVII0O+a+BAGDAgiujoqTyxYze/bd4Mr7zidlkiIuLHjIFJkwweex/jiw3FEydCYaHbZYlUiUb6a4BGhnzP+vUwbZqHsLD7uKplEkcbAzff7HZZIq7RSL931J7LwWzcCI8/nklSxCTGxcY6N5SJ1DGN9IscoEMHuP32EAoLH2DWjn0sCQmBBx7QJl4iIlJtbdvCoEENyDbDebuoCGbNcrskEa8p9EvAatUK7rvPkJ09lrmZJcyLj4fx4xX8RUSk2vr3hxYtjmJpQSfWLV0Kv/3mdkkiXlHorwbN6fcfTZvC5MmGHTtu5cecSD5LSIAJExT8RQRQey7Vc8st4PEMY2puKPlTpzqbxoj4OM3prwGaA+r7cnLg3nshOfkVjquXw1lZWc46n0ZTnCU4aE6/d9Sei7cyMuD++3NoEDaGB+PrO1NIReqA5vSLHERcHDzyCKSlXcuynBi+SkyEiRM14i8iItXSuDFcdVUc+0JG8U5JCbz3ntsliRyUQr8EjZgYePhhSEsbyffZEcxv0MDZUl3BX0REqqFPH2jXrhuL9nVkw8KFkJLidkkilVLol6ASGwtTpkB6+vV8m2n4pkEDePRRt8sSERE/deONAMN4NN9Q8tBD4PG4XZJIhRT6JejExcFDD8HOnaNZkFHM0vBweOEFt8sSERE/FBIC994bSmjkP3kqIgKef97tkkQqpNBfDVrtwf/Vrw8PPggZGbfwUXoeq/PyYOZMt8sSERE/1KIF9OvXmE3FZ7Dsjz/g11/dLknkLxT6q8FaO95aa/Y/3K5Hqic+3tlWPTd3DNO37mRzaip88onbZYmIiB8aMgRiY8/lX3sjyZs2DQoK3C5J5E8U+iWoNWwId99t8Hge4NEtKexYsQIWLnS7LBER8UNjxxpMyFimREbB44+7XY7Inyj0S9Br2RJuuCGUmJjHGJ+6nezPPoO1a90uS0RE/Ez9+nD11fXYVXIFX+3cCcuXu12SyH8p9IsAnTrBxRdH0KjxY4zNzaXo2WchPd3tskRExM+cfDIkJ/didl4D9r34IhQVuV2SCKDQL/JfvXrBqafG0jR5POOiorATJkBurttliYiIn7n9dsDczqPhEfD0026XIwIo9Iv8yYAB0Lp1Yxo0GcHUpk1hzBgoLXW7LBGpJVqNTWpDTAxceWU90osvZklKCqxe7XZJIgr9IuVddRWUlnaiJPpkZrVu7ezaKyIBSauxSW059VRISjqFV3NiKJg2TQNI4jqF/mrQyFDgu+su2LHjDNblJPBtXBy89prbJYmIiJ+56y4DIXfzZGQk/OtfbpcjQU6hvxo0MhT4QkLggQdgx47hfJC6ly3btsG337pdloiI+JG4OBgyJJ4/cs9i5erVsH272yVJEFPoF6lEVBRMmGAoKrqfKanp5H34IWzd6nZZIiLiR/72N0hIOJPnCiLwPPGE2+VIEFPoFzmIxo1h1KgwEhs8xD89HuyUKZCX53ZZIiLiR26/PYTSkNt5o7AQ5s1zuxwJUgr9IofQvTv07Fmfhkk38FTTpjB+PFjdyiEiIt5p2hR69mzBoryOZM6eDYWFbpckQUihX8QLF10ExnQgL+w4Pm/UCF5+2e2SRETEj4wYASGhI3jchMK0aW6XI0FIoV/ES86KPucyb0cpW1NSYNkyt0sSERE/ERoKf/97NNsLLuXHzZthyxa3S5Igo9Av4qWQEGfJ/vz8O5myI4OiGTMgM9PtskRExE/07AktWpzMv3KiKH3mGbfLkSCj0C9SBfXrOzf2JiRO5MGYGJg4ETwet8sSERE/cdttBht6K2/k5MDixW6XI0FEoV+kinr0gA4dGhJW73zeb94cnn7a7ZJERMRPJCZCr17N+E9RJ3LefFMDR1JnFPpFqmHYMMjJ6cN3O0PYkJcHixa5XZKIVIN2WBc3XHcdhISO5BkMzJzpdjkSJBT6RarBGLj/fsjNvYXHt6aS/957kJXldlkiUkXaYV3cEBYGl1xSj9/y+rNl8WLYt8/tkiQIKPSLVFO9ejB6dCgJiZN5OCoKJk/W+v0iIuKVAQMgPv58nvKEw7PPul2OBAGF/mrQ5WDZr0sX6NIlEU/UOXzZqBG89ZbbJYmIiJ+48cZwMvKvYcnGjbBjh9vlSIBT6K8GXQ6WA115JeTl9eWz7QVkrFoFGze6XZKIeEmDOOKmDh2gbdvjmFFcH6sNu6SWKfSLHCZjYOxYsPZeHtiXi33sMSgudrssEfGCBnHEbTfeaMgvHcVn6enwxx9ulyMBTKFfpAYkJsKgQZHUTxzNK82awZNPul2SiIj4gUaNoFu3dnxYnEzpCy+4XY4EMIV+kRrSrx+Eh3fk95xG/JaTAytWuF2SiIj4gRtuMJRyPW9nZanvkFqj0C9Sg+66C3JyRvL0rixKX3wRiorcLklERHxcbCyceGIz5hW2p+jVV90uRwKUQr9IDYqMhH/8I5SY2Ht5rmFDmDrV7ZJE5CB0I6/4imuvBRM6gpfz8+Gbb9wuRwKQQr9IDevaFZo3T2ZLflt+y82Fn392uyQRqYRu5BVfER4OZ53VkKX5x7Lvvfe074vUuKAP/caY4WUjPBe6XYsEjtGjIS/vap7K2EPp9Oma5iMiIod08cUQEXkFL+TmabRfalxQh35jTBtgJPCdu5VIoAkLgxtuCKVe3D0837ixVvMREZFDCgmBgQPjWFFyCjmzZ7tdjgQYnwj9xpgWxphnjDFLjTF5ZSPvbSo5t6UxZrYxJtsYs9cY84ExplU1fmYI8BJwI1B4WG9ApALdukGTJs3ZnNuS9VlZsG6d2yWJiIiPO+88iIoawosa7Zca5hOhH2gPDAX2AIsqO8kYEwN8DXQCrgGuAjoAC4wxsVX8mbcBi621y6tVsYgXbroJcnOH81RWjrPbosfjdkkiIuLDjIGBA+vzc/HJztx+kRriK6H/P9baJtbac4CD/YaPBNoCF1prP7LWzgEuAFoDo/afZIyZZ4zZVcmjjzGmGzAYmFyL70mE8HC4/vpQomJu5Y0GDWDGDLdLEhERH3feeRAZNYTpeVrJR2qOT4R+a623w58XAN9Zazcc8NxNwGJg4AHH+ltrG1XyWAycArQB1htjNgO9gReNMaNr6C2J/FePHpCY2JZlWXHs+vVX2LnT7ZJEpIyW7BRfZAxceGE8ywpPIldz+6WG+ETor4KuwKoKjq8Gunj7Itba5621ydbaNtbaNjg38v7dWjutovPLdwrlH9V5IxJcbr0VSktv5HFr4bHH3C5HRMpoyU7xVc5o/8W8uC8XFi92uxwJAP4W+hvgzPsvLxNIrK0fWr5TKP+orZ8rgSMmBs49N4rikIEsjI2F+fPdLklERHyYMTBoUDw/Fp5E3rvvul2OBAB/C/21wlrb11r7kdt1SGA75xzweE5lVnoOhR98AMXFbpckIiI+7LzzICp6MK/v3g2rKproIOI9fwv9e6h4RL+yKwAiPuX22w0RkffybGwsvPii2+WIiIgPMwb692/AEs+xlLzxhtvliJ/zt9C/Gmdef3ldgDV1VYRu/JLqat4cjjyyIRvz2pP6xx+6qVdERA5q6FAg9DLe274dtm51uxzxY/4W+j8Gehtj2u4/ULaJV5+y79UJ3fglh2PkSCgtHcZTJaXaqVdERA4qLAz69Elmru2Affllt8sRP+Yzod8YM8QYMwQ4ruzQ2WXHTjvgtOnAZmCOMWagMeYCYA6wDfhXnRYsUk1hYTB0aAQFDGYpwLJlbpckIiI+7JprIL/0Cr7asgV273a7HPFTPhP6cTbleg+4vuzr58q+nrD/BGttLtAP+B14A5gJbAL6WWv31Wm1Iofh9NMhPPwU3thTgOe117RTr4iIVCo6Go49ti3vlzYFjfZLNflM6D/Ikph9y5231Vo72Fpb31obZ6290Fq72Z2qRapv9GhDq58R8wAAIABJREFUWPitvB4fDzNnul2OiIj4sJEjDXsKLueH1auhsNDtcsQP+UzoFwk2RxwBzZq1ZFlWHNnLl0NentsliYiIj0pIgK5du/JWSSxo3X6pBoX+atDqPVJTbrwR4EaeCgmB5593uxwREfFhw4eHsjX/YlIWLQKr+CFVo9BfDVq9R2pKbCycfHIsO4tOYFNKCmRkuF2SSFDRII74k9atoUWLE3klJw8WL3a7HPEzCv0iLrv0UvB4LuL5Eg9Mm+Z2OSJBRYM44m8uuSSKX4vPJu/9990uRfyMQr+Iy4yBK64IZ2/xIH7KzYUNG9wuSUREfNRJJ0G9uLN4Kz0dNm92uxzxIwr9Ij6gTx+Ijj6V1/JKsC+84HY5IiLio4yB/v0bstBzNPb1190uR/yIsYe4EcQYc/Xh/ABr7YzDeb4/MMbYQ/07ihzK2rXwzDO/MrDeG/xt8GA44QS3S5IAYoxB01cOTe25+IOiIrjyyj/4e+R4+j//PNSr53ZJUoeq256HeXHOa1Uv578sEPChX6QmdO4MCQnd+DATBrz9NiG9ejlDOiLiFWPMcOAVYJC19iO36xGpLRERcPzxR/DBT43o//bbMHKk2yWJHzjk9B5rbchhPELr4k2IBIrRow2hoaN5KywM5s1zuxwRv2GMaQOMBL5ztxKRujF8eAhpBZew4fvvtXyneEVz+qtBS7xJbWnWDJo3b8XifXEUffKJGnLxO8aYFsaYZ4wxS40xeWXtZJtKzm1pjJltjMk2xuw1xnxgjGlVjZ8ZArwE3Ahoq1IJCo0aQbt2xzAjKw9++MHtcsQPVCn0G2M0co+WeJPa9Y9/gOUGXg8NhX//2+1yRKqqPTAU2AMsquwkY0wM8DXQCbgGuAroACwwxsRW8WfeBiy21i6vVsUifmro0EhWlZxHgZbvFC9UdaT/K2PM9bVSiYgAEB8PnTs3ZlluEgVffKHRfvE3/7HWNrHWngO8d5DzRgJtgQuttR9Za+cAFwCtgVH7TzLGzDPG7Krk0ccY0w0YDEyuxfck4pOOPx7i4vvz/ubNkJnpdjni46oa+rsBqRV9wxjzsjFmzOGXJCIjRwJmJK+Gh8NHuh9R/Ie11uPlqRcA31lr/7sxhbV2E7AYGHjAsf7W2kaVPBYDpwBtgPXGmM1Ab+BFY8zoGnpLIj7LGDj11CS+Np1h5ky3yxEfV9XQnwBkGmNCy0ZZOh3wvQXA5TVXmkjwiomB7t0b8fO+ZPIXLACPtzlKxG90BVZVcHw10MXbF7HWPm+tTbbWtrHWtsG5kffv1toKt7cuf09W+Ud13oiImy69FHYXDmLVTz/pyrAcVFVDfwZO8I8AGpR9vt82nNEWEakB110HxozgpbAweO9gsyRE/FIDnHn/5WUCibX1Q8vfk1X+UVs/V6S2xMZC585dmbnPAwsXul2O+LCqhv65wHU4l2UBjjrgew2BkpooSkQgOhqOPjqRVbktyV28WKP9Il6w1vbVGv0SbK68Mpy1RQPJ/fhjt0sRH1bV0D8RZ77kW8B6oLMx5lxjTHvgzrJjIlJDhg8HGM50a2HOHLfLEalJe6h4RL+yKwAiUomuXaFho9N4Z8sW3dArlapS6C+7yaozcB7QB3gSeBb4DTgBmFrTBYoEs6goOO64BFYXtHPm9mu+pgSO1Tjz+svrAqypqyK074oEiv79G/IfjoV33nG7FPFRVd6cy1qbba393Fq7y1q7BeiOM92nt7X2rRqvUCTIDRsG1l7Nax4PfPWV2+WI1JSPgd7GmLb7D5Rt4tWn7Ht1QvuuSKAYNAiyS89j9XJtVyEVO+wdea21Odbaz6y12g5OpBZERMCxxzbgp4JWFGmzLvEDxpghxpghwHFlh84uO3baAadNBzYDc4wxA40xFwBzcBaF+FedFiwSAKKioFOnzryTmQcrVrhdjvggYw8xXcAYc/Xh/ABr7YzDeb4/MMbYQ/07ihyO/Hy46aYdnBQ5ieFDh8Kpp7pdkvgZYwx1NZJ9kGkyC621fQ84rxXOtNABgAHmA7dYazfXdo2VUXsu/uznn2HCP2fyXo+1hE/WfnWBqrrteZgX57xW9XL+ywIBH/pFalt0NHTu3ISla5O4+sMPCVXoFx/mbWdkrd2Ks5uuiNSAY46BevF9+ffvcxhYWAiRkW6XJD7kkNN7rLUhh/EIrYs3Udd045e44brrwGOvZXZJCfz0k9vliAQEtecSaHr3bsbnJc3hk0/cLkV8zGHP6Q9GuvFL3BAfD+3atWBBYX3srFlulyMSENSeS6C57DLDtsJBpC9Y4HYp4mMU+kX8yMiRUFQynC+ysmDrVrfLERERH9OwIRxxxDG8lZoGaWlulyM+RKFfxI80agQtW7bj37Y+vPaa2+WIiIgPOv/8OH6wp+mqsPyJQr+In7nuOsOevItZsX07ZGe7XY6IiPiY/v2hwJzJT1q6Uw6g0C/iZ1q1gqZNj+MtG6bRfhER+YvQUOjevT0f7MqHdevcLkd8hEK/iB+65JJQtueeRepvv0FxsdvliIiIj7n44nBW2wsonT3b7VLERyj0i/ih44+HhIQzeMVj4b333C5HxG9pyU4JVEcdBVFxffj6119BG84JCv0ifmvAgGg25Pdk39KlatBFqklLdkog69mzFZ/mxcHSpW6XIj5AoV/ET513HkREDOb1oiJYtMjtckRExMdcckkI60suJE8bdQkK/SJ+KzQUjjsugZ+L21CiBl1ERMpp2RKSmvTkE93/JSj0i/i1q66CEs8VfLp7N6Smul2OiIj4mFNOacK84rbw5ZdulyIuU+gX8WOxsdC+fUvmmkbwxhv/3969B1ld5ncef3+5igwiiCDIpbkrKMhVkZGLAqIzoqnBSSWbyUwqmcwktVvJJlup2d3MxtlsJVuVzSZbVk2STSq7WxlnZ1bHEWVURFpBEFBRFHEQkOamqKCA3OT67B/nMPa03dDdnO7n9DnvV9Up+vzOc37nc+zH53z793vO88sdR5JUZhYvDvacvocDy5fnjqLMLPqlDu4b3wg+OvYrbKmrg1OncseRJJWR3r2hpmYiD29/B06cyB1HGVn0t4JLvKmcDB4MAwdO4aGzwE9+kjuOJKnMLFx4JWvTdKf4VDmL/lZwiTeVm0WLurHjxGyOrF+fO4okqcwsXAgfnlvAvhUrckdRRhb9UgWYMwd69LiLh44cgVdfzR1HklRGLrsMRo26np+8U+cUnypm0S9VgAiYMqUPL58bQfKS65KkBhYs6M1LOMWnmln0t4Jz+lWOfuM34NNT9/P8hx/CoUO540iSysiCBcUpPrW1uaMoE4v+VnBOv8pRr14wfPgonkiXwQ9/mDuO1CF4EEfVojDFZ5xTfKqYRb9UQX791zvx7tF72LdpEyTrF+liPIijanLnnVew/tzNTvGpUhb9UgUZNw6uuuqL/ODoMXj55dxxJEllZMEC2H9uvlN8qpRFv1RhZs/uyZtM4txjj+WOIkkqI927w8hR1/OT7Tvg5MnccdTOLPqlCnPffXDqzCKee+89OHIkdxxJUhm5++4rWHfmZvBof9Wx6G8Fv/ilcta9O4wePZyfRS/40Y9yx5EklZH582E/8/lw+fLcUdTOLPpbwS9+qdz92q91Ys/RuziwcWPuKJKkMtKtGwwfMZbH394K587ljqN2ZNEvVaDrr4errrqtcIXe117LHUeSVEbuuKMPL568Hl56KXcUtSOLfqlC3XJLL15L40k//WnuKJKkMrJwIexJCzm6dGnuKGpHFv1Shbr/fvj09L2s3rkTTp3KHUeSVCZ69YJrh4zjmZ//PHcUtSOLfqlC9ewJI0aMZunZrvDUU7njSJLKyMyZ17D8YF/YsiV3FLUTi36pgi1e3Jntny7k+MqVuaNIksrIokXBjvRlTi9ZkjuK2olFv1TBJk2CK66YxWP79sFHH+WOI5Udl2BWtRowAPpeM4EXXOWtalj0t4IfEuooImDixAGs7jIUfvzj3HGksuMSzKpmU6YM5Wf7zsK+fbmjqB1Y9LeCHxLqSH71V2Hf0S/x7ptv5o4iSSoj993XmZ9zD+nxx3NHUTuw6Jcq3MCBMGjQZH708UF4663ccSRJZWLUKOjeZyqvrV+fO4ragUW/VAVuu+0LbOo8hfTII7mjSGXF6ZqqdhMmjGDJOx/AyZO5o6iNWfRLVeDee+HoybvYsGMHnD2bO45UNpyuqWr3pS91ZzN3gKu8VTyLfqkK9OgBI0eO5bHTQG1t7jiSpDIxdSoc6XYrHyxbljuK2phFv1Ql7rmnC1tP3c6pFStyR5EklYlOnWDEqNE8+fbW3FHUxiz6pSoxcyZ06z6XZTt3OndTKnJOvwRz5/Zl9eEa2LYtdxS1IYt+qUpEwPjx17Li3JXw9NO540hlwTn9Etx5Z7CLOzm9dGnuKGpDFv1SFVm8uBM7Pr2bY6tW5Y4iSSoTvXtDv2vHsWbDhtxR1IYs+qUqMno09O07naV79sCxY7njSJLKxKRJQ3hy13E4ejR3FLURi36pykyYMIAXYiB4BUZJUtE993RlCwtd4a2CWfRLVeb++4Pdx+/mkFdglCQVXX89nO45lZ3Ll+eOojZi0d8KrvagjmzIEBhwzVR+uncvfPxx7jiSpDIQAWPGjuSJt9+BZGlTiSz6W8HVHtTRTZ58FS91HQ2PPpo7ipSVB3Gkzyxc2JuXjk5w6c4KZdEvVaGvfAXePbaQD15/PXcUKSsP4kifmTMH3o3ZnHzyydxR1AYs+qUq1L8/DB58E4/s3QsHDuSOI0kqAz16wDXDruPFl1/OHUVtwKJfqlLTpvVmY9cbYcmS3FEkSWXippsGs2znMTh9OncUlZhFv1SlfuVXYN/xO9m3cWPuKJKkMnHXXV3Zcm4uvPRS7igqMYt+qUpdeSUMGXIDj+3ZA4cO5Y4jSSoDN9wAx3pM44Dz+iuORb9UxWbM6M2rnW+EJ57IHUWSVAYioGbkSJ7d/FbuKCqxqi36I6JbRPxtRGyLiE0R4eVJVXXuuQfe+3Q++195JXcUSVKZmD27P7X7rvQscIWp2qIf+AugGzA2pXQj8LuZ80jtrk8fGDLkRpbs2gVHjuSOI0kqA/PnB3XMJ9XW5o6iEiqLoj8iBkfEgxGxNiKOFy+SUtNE2yER8UhEHI6ITyLi0YgY2sLXu5xCkf+dlNI5gJTS+5f6PqSOaOrUPrzSZSw4f1OSBAwYAJcPuIG3ly/PHUUlVBZFPzAK+CpwEHihqUbFYr0WuA74OvA1YDTwXET0bOHrHQS+ExEvR8SaiPhya8NLHdm998KeE3dycN263FEkSWVizNjhPLltd+4YKqFyKfpXpZQGpJTuBh6+QLtvAiOA+1JKj6WUlgCLgGHAt843iohnI+JAE7eZQBdgKLA9pTQN+G3gnyNiRFu9QalcXX01DBo0gSd27YITJ3LHkSSVgYULe7HhyI1QV5c7ikqkLIr+81NsmmERsC6ltL3ec+uANcC99bbNSyn1a+K2BtgNJOAHxfZbgNeBySV6S1KHMmlSP9YxFFasyB1FalcR8UBxSmmKiJQ7j1QuZs6EfZ1mceqZZ3JHUYmURdHfAuOBNxvZvhkY19ydpJQOAMuAhQARMRC4EdjUWPuGHwoNby1+F1KZue8+2H3yLo6sXJk7itSuUkoPpJTi/C13HqlcXHYZDKgZy9rVq3NHUYl0tKK/L4W5+A19DPRp4b5+D/iDiNhE4Q+AP04pvd1Yw4YfCg1vLXxdqewMGgRX95/IM3V1cK65J94kSZVs4sTBLK/7BJLHNytBRyv6SyaltDOldEdK6caU0oSU0kO5M0k5TZgwgFUn+4Br9kuSgAULurPl05th+/aLN1bZ62hF/0EaP6Lf1BkASc20aFHwzpk7OfP007mjSJLKwMSJ8FGXqZxwXn9F6GhF/2YK8/obGge02/Wi/eKXKtHIkfCF3pN4cdu23FEkSWWgc2e4dtQYXnzxxdxRVAIdreh/HLil/tKaxYt4zSw+1i784pcq1dixw1j20QmXaFOHEhHdIuJvI2JbRGyKiHb7PJAq3aRJg1hed9x5/RWgbIr+iFgcEYuBKcVNdxW3za7X7B+BncCSiLg3IhYBS4A9wD+0a2CpAi1c2IWtnReQnngidxSpJf4C6AaMTSndSOGK65JKYMGCbmw9dSts3Zo7ii5R2RT9FC7K9TDw7eL97xfvf+98g5TSMeB2YCvwL8BDQB1we0rpaLumlSrQ9OlwqtMM3t68OXcUdVARMTgiHoyItRFxvDgNsqaJtkMi4pGIOBwRn0TEoxExtIWvdzmFIv8756/5klJ6/1Lfh6SC8ePhULfJHH3qqdxRdInKpui/wJKYcxq0251S+kpK6YqUUq+U0n0ppZ15UkuVpXNnGDlyJE/sfQ8OHcodRx3TKOCrFBZXeKGpRsVivRa4Dvg68DVgNPBcRPRs4esdBL4TES9HxJqI+HJrw0v6ZZ06weBRo1mzfn3uKLpEZVP0SyoPs2dfzuYuM8BVfNQ6q1JKA1JKd1M4W9uUbwIjgPtSSo+llJZQuOr6MOBb5xtFxLMRcaCJ20ygCzAU2J5Smgb8NvDP9b/7JenSTJs2kGfqTjuvv4Oz6G8FV+9RJZs3Dz48NZv969bljqIO6PwUm2ZYBKxLKf1iAfCUUh2wBri33rZ5KaV+TdzWALuBBPyg2H4L8DowuURvSap68+d3ZfvpL8KWLbmj6BJY9LeCq/eokvXsCUOHXc/SHTvg7NnccVS5xgNvNrJ9M4VlmJslpXSAwlXVFwJExEDgRmBTY+0bHrRpeGvxu5CqwNixcLTHTRx+8sncUXQJLPolfc6UKX15+cxQeOml3FFUuZq6qOLHNH4Rxgv5PeAPImIThT8A/jil9HZjDRsetGl4a+HrSlUhAoaMGsULfiZ0aBb9kj7ny1+GujMLOeVVGNUBpJR2ppTuSCndmFKakFJ6KHcmqdLMmHENy+twXn8HZtEv6XMGDoQrrxrPGtdlVts5SONH9Js6AyApo9tv70LdqZth587cUdRKFv2SGnX99UN49sPj8MEHuaOoMm2mMK+/oXHAW+0VwoUZpOYZNQqOXnYDx1asyB1FrWTRL6lRCxZ0YVvn+bBsWe4oqkyPA7fUX1qzeBGvmcXH2oULM0jNEwHXjhrN+uefzx1FrWTRL6lR06bB0U7T2e0Xt9RCEbE4IhYDU4qb7ipum12v2T8CO4ElEXFvRCwClgB7gH9o18CSmmXiTYN4dueJ3DHUShb9reDpYFWDzp1h+PDRPLmjzqU71VIPF2/fLt7/fvH+9843SCkdA24HtgL/AjwE1AG3p5SOtmtaSc1yxx3d2XpkIhw4kDuKWsGivxU8HaxqMWNGb149N8alO9UiF1gSc06DdrtTSl9JKV2RUuqVUrovpbQzT2pJFzNhAuzvPIkzq1bljqJWsOiX1KQ774TdZ+dxavny3FGkNuGZW6n5OneGASNGsfHpp3NHUStY9Etq0tVXw1X9x7H67UavcyR1eJ65lVrm+vFDeXb7h7ljqBUs+iVd0NixxaU7P3SQl6RqN29eT974aAwcP547ilrIol/SBS1Y0IVtneZBbW3uKJKkzKZNg/eYRlq/PncUtZBFv6QLmjYNjnW+md1r1uSOIknK7LLLoPeQkWx94oncUdRCFv2SLqhzZ6ipGcnyHXWQ/J6jJFW7kWOGU/vWO7ljqIUs+iVd1IwZfXj5+BDYvj13FElSZvPmXcn6fdfCuXO5o6gFLPolXdS8ebAz3cG5ZctyR5FKyiU7pZabOTPYe2Y6bN2aO4pawKJf0kUNHAg9rxrHGxs35o4ilZRLdkot17s3dO47nA88ENShWPRLapYRI2pYvnMvnD2bO4okKbPBY0azdrULPHQkFv2SmmXOnMvZnKbAK6/kjiJJyuzWW/vz3K4uuWOoBSz6W8E5oKpGs2bBu+dm8eny5bmjSJIyu+22Luw8OgGOHMkdRc1k0d8KzgFVNerVC/oPGsuLW7bkjiJJymz0aDjUdRxnvEhXh2HRL6nZrrvuWmr3HoZjx3JHkSRlFAH9hw9n0+OP546iZrLol9Rsd9zRlW2d5sALL+SOIknKbPS4YTz38725Y6iZLPolNdv06XCw0y18VFubO4okKbN5867g1feHeLX2DsKiX1KzdekCQ2tG8tw7Xn5dkqrdtGnw3pnJsGNH7ihqBot+SS0ycWJ/1uzvAgcP5o4iScqoVy/o0nco+595JncUNYNFv6QWmT+/EzuY57x+VQSXYJYuzeAxo1n7/MrcMdQMFv2SWmTsWDjZYxLvPv987ijSJXMJZunSzPziNdTWdc4dQ81g0S+pRSKgZngNz7/jHE5Jqna33daFuk/Gw8mTuaPoIiz6JbXY1KlXsfZALzh0KHcUSVJGo0fDwc7XcWbDhtxRdBEW/ZJabO7cYGfMdV6/JFW5CLi6Zhhv/exnuaPoIiz6JbXYiBFwpucE9j73XO4okqTMRowbzupXt+SOoYuw6JfUYhEwrKaG57e7Xr8kVbu5c/vw0nv9c8fQRVj0t4JLvEkwfbrz+iVJMH16sO/4eDh2LHcUXYBFfyu4xJtUmNe/q9PtzuuXpCrXrx+cunwwn65blzuKLsCiX1KrDB8OZ78wgT21tbmjSJIyGzBiJG8sXZo7hi7Aol9Sq/xiXr/r9UtS1Rs7YRgrN+3JHUMXYNHfCs7plwqmTbuKdQd6weHDuaNIkjKaO/cKXv9gcO4YugCL/lZwTr9UMHdusCvNgbVrc0eRWsWDOFJpTJkC+z4d7+IOZcyiX1KrDR8OZ3pNZO+KFbmjSK3iQRypNHr1AnoN4oiLO5Qti35JrRYBQ4YNY7Xr9UtS1btm1Che9cq8ZcuiX9IlmT69H6v3doXTp3NHkSRlNHHqEGrf+jh3DDXBol/SJZk1qxN7On0RXnstdxSpxZzTL5XO7NmXs+XA8Nwx1ASL/lbwQ0L6zJgxcLTbBA4/+2zuKFKLOadfKp2JE+GDk9fBRx/ljqJGWPS3gh8S0mci4Nrhw1m3aVPuKJKkjC67DLr1HcBHK1fmjqJGWPRLumTjxw9k5e5jkDzxJUnVbNCo0bz81NO5Y6gRFv2SLtns2V155/Q02L49dxRJUkY3zxpK7eZjuWOoERb9ki7Z5MnwYaepnPKUrjoYv6MlldaMGd2pO1iTO4YaYdEv6ZJ16wb9h43ktTVrckeRWsTvaEmlNW4c7D85Eo4fzx1FDVj0SyqJkSMH8/yuD3LHkCRl1K0bdO3dj0MeBCo7Fv2SSmLWrMt568hY2L8/dxRJUkYDRoxg45NP5o6hBiz6W8E5oNLnzZgBe7mV5NEdSapqE2+uYeWr+3LHUAMW/a3gHFDp83r3hi9cM4a3vUiXJFW1L972BbYcGJI7hhqw6JdUMkNrhrGqblfuGJKkjCZMgA+Oj4LTp3NHUT0W/ZJK5tZbe/PK+30d6NVhOF1TKr2ePYGe/TixcWPuKKrHol9Sydx2W7A3zYQ338wdRWoWp2tKbaP/sBo2LV2aO4bqseiXVDKDB8Opy8dwqLY2dxRJUkZjJo9g9botuWOoHot+SSV17cgRrN+wIXcMSVJGc+ZeyWvvDcgdQ/VY9EsqqXHjruGF3cdyx5AkZTRpUvD+sTGQ/KpMubDol1RSs2Z1451jN8DBg7mjSJIy6dMHTnW9mtPbt+eOoiKLfkklNWkSvM9kzq1blzuKJCmjftdey5ann84dQ0UW/ZJK6rLLoPfgEWx75pncUSRJGQ2/cSTratfnjqEii35JJTd42DBe3L4jdwxJUkZz5g9gbV2v3DFUZNEvqeRuvbUP69+7wi9wSVIVmzSpEx8cG5k7hoos+lvBKzhKFzZjRvDuuWngF7gkqWoNGgSfnO5POnEidxRh0d8qXsFRurCaGjjW/XqOrVyZO4okKZMI6NW3Lx+88ELuKMKiX1IbiIBBI0fyyqpVuaNIF+SZW6lt9R8+nNefWZ47hrDol9RGRo8ZyKrdrtWv8uaZW6ltTZ1Tw8p1+3LHEBb9ktrIrFk9ePvQaPj009xRJEmZ3DKjJ+8cHJo7hrDol9RGpk6FfeemkN54I3cUSVIm48fD/mM1uWMIi35JbaRXL+jer4b3VqzIHUWSlEmPHnCuay9O792bO0rVq9qiPyLujohXI2JjRGyOiG/nziRVmkEjRvDSq6/ljqEK53gulbd+gwbx9lNP5Y5R9aqy6I+ITsAPgW+klG4C5gN/HRGD8iaTKsvkyVezZo8LoqjtOJ5L5W/ouOG8vGJt7hhVryyK/ogYHBEPRsTaiDheXDqtpom2QyLikYg4HBGfRMSjEdHSb4icX6HhyuK/vYAjwNFWvQFJjZoxowu7jo2D48dzR1E7cTyX1NCsO6/lxbe75o5R9cqi6AdGAV8FDgJNXsEhIi4HaoHrgK8DXwNGA89FRM/mvlhK6SxwP/CTiNgFbAC+nVL6pNXvQNLnjB8PH3IDaePG3FHUfhzPJf2SyZM78/6xEbljVL1yKfpXpZQGpJTuBh6+QLtvAiOA+1JKj6WUlgCLgGHAt843iohnI+JAE7eZEdEF+I/AV1NKw4BbgL9rxREmSRfQrRtcMWgodcuW5Y6i9uN4LumXDBkCh09eA2fP5o5S1cqi6E8pnWtm00XAupTS9nrPrQPWAPfW2zYvpdSvidsa4CZgUErpuWL7N4E3gZtL9Z4kFQyqGc7Lm97MHUPtxPFcUkMR0OOKXny0YUPuKFWtLIr+FhhPYTBvaDMwrgX72QMMiogboDAHFZhQ3I+kEpo+/SrW7O2eO4bKj+O5VEX6Dx3KG8ueyR2jqnXJHaCF+lKYJ9rQx0Cf5u4kpfRBRPwO8MOISEBn4D+klN5qrH1EPAD82YX2GeHV26ULeTD+b+4IKi+O51KV+cGTwH/6bu4YVaujFf0lk1L6EfCjZrZ9AHigqccjIqWUyvJToi2zXeq+W/r8lrS/WNuoCvwUAAAJIElEQVRLebyxx+wD7ff85j6nOe1a+nu+2GPl3g9yZ2grjuf59+143j4qqQ+05DmO57+steN5R5vec5DGjwA1dcRIklSeHM8lqR11tKJ/M4V5oA2NAxo9lStJKkuO55LUjjpa0f84cEtE/GKx1+JFX2YWH9Pnfa+M993S57ek/cXaXsrjbfnftC1UUh9oyXOa0661v+eO1geg/DI7nrdcJf2/7HjeOpXUB1ryHMfzX9aqzJFSeUzzjIjFxR/vAL4N/D6wH9ifUlpZbNMTeB04AfwpkIA/p3AFxgkppSxXYCzneV9qH/YBgf3gPMdzdWT2AUFl9oNyKvqbCrIypTSnXruhwN8A8ylcfn0F8IcppZ1tnbEpldgx1DL2AYH94DzHc3Vk9gFBZfaDsin6O7JK7BhqGfuAwH5QCfwdyj4gqMx+0NHm9JerjjgfTKVlHxDYDyqBv0PZBwQV2A880i9JkiRVOI/0S5IkSRXOol+SJEmqcBb9kiRJUoWz6G9nEdEnIpZGxNaIeD0inomIUblzqX1FxHeLfeBcRNyXO4/aXkSMjIjVxd/7axExNXcmXRrHc4HjeTXqqOO5RX/7S8DfppTGpJQmAkuBf8qcSe1vObAQWJU7iNrN3wP/J6U0BvgT4KGIqKjl4KqQ47nA8bwadcjxvOqL/ogYHBEPRsTaiDgeEal4KfjG2g6JiEci4nBEfBIRjxYvLtNsKaVDKaVn6216EWj09dQ+2rsPAKSU1qWUdlxqdrWdUvaLiLgauAX43wAppeUULkY1pc3fSBVxPJfjuRrjeF5Q9UU/MAr4KnAQeKGpRhFxOVALXAd8HfgaMBp4LgqXk2+tPwSWXMLzdely9wGVp1L2i6HAvpTS6XpP3VncrtLJ/f+y43l+ufuAypPjOdAld4AysCqlNAAgIn4HWNBEu28CI4CxKaXtxfZvANuAbwH/vbjtWeCmJvZxb0ppzfk7EfFnxX3+bgneh1ovWx9QWStpv1C7cDyX47ka43iORT8ppXPNbLoIWHe+ExSfWxcRa4B7KXaElNK85uwsIv4UuBtYkFI63rLUKqVcfUDlrcT9YjcwMCK61js6VFPcrhJxPJfjuRrjeF7g9J7mGw+82cj2zcC4luyoeEToHgofEIdLkE3to2R9QBXlov0ipbQfeAn4BkBEzKcwB3RD+0RUA47ncjxXYyp6PLfob76+FOaCNfQx0Ke5O4mI8cADwFXAyojYGBGvlCSh2lpJ+gBARDwQEXuBGcA/RcTeiBhcgoxqf83tF98GfisitgJ/BfyrlFJqh3z6PMdzOZ6rMRU9nlf99J72llLaTOEvQlWxlNIDFIoFVYmU0jbg1tw5VDqO5wLH82rUUcdzj/Q330Ea/+u/qb8KVXnsA2qM/aLj8Xcm+4AaU9H9wqK/+TZTmOvV0DjgrXbOojzsA2qM/aLj8Xcm+4AaU9H9wqK/+R4HbomIEec3FC/sMLP4mCqffUCNsV90PP7OZB9QYyq6X0QH+N5Bm4uIxcUf76Dw5YzfB/YD+1NKK4ttegKvAyeAP6Vw+fU/B3oBE1JKR9s7t0rHPqDG2C86Hn9nsg+oMfYLi34AIqKp/wgrU0pz6rUbCvwNcH55phXAH6aUdrZ1RrUt+4AaY7/oePydyT6gxtgvLPolSZKkiuecfkmSJKnCWfRLkiRJFc6iX5IkSapwFv2SJElShbPolyRJkiqcRb8kSZJU4Sz6JUmSpApn0S9JkiRVOIt+SZIkqcJZ9EsNRMTgiHgwItZGxPGISBFRkzuXJKllHM+lz1j0S583CvgqcBB4IXMWSVLrOZ5LRRb90uetSikNSCndDTzc2p1ExM6IeKB0sSRJLeR4LhVZ9EsNpJTO5XjdiJgYET+NiI8i4kREvB0R/75Bm4XF09QnIuJwRDwWEWMbtBlT3M+HEfFpROyOiIcjokvx8Ssi4lxEfKt4/7eK93sV76+IiCfb631LUltxPHc812e65A4gCSJiOvA8sB34t8BeYDQwoV6bhcDPgFrgV4EvAP8ZWB0RN6WU3i02/RmFU9m/BxwArgXu5rM/8icBAbxavD8Z2JZSOlLv8e+X/E1KUhVwPFe5suiXSiAiAujcyEOdzh+RKUoppbONtPtvwEfALSml48VttQ3a/BdgB3BXSulM8XXXAluBPwb+KCL6UZjDem9K6fF6z/1hvZ8nA2eAN+rdf7W4v+FAHz77AJGkquJ4rkrl9B6pNGYDpxvchgHfbbBtRcMnRsTlwEzgoXofEA3b9KQwmP/4/AcEQEqpDlhTfH0ofNDsAP5rRHwzIkY3srvJwFsppZMR0QmYyC8fJQI/JCRVL8dzVSSP9EulsQGY1mDb48BS4H/W23aEz+tD4Q/wvRfYfx8Kp3D3NfLY+xQ+kEgppYiYDzwA/CVwVUTUAX+VUvq7YvvJwEvFn8cCPfnsQ2ES8HFKaecFskhSJXM8V0Wy6JdKoDh/8pX62yLiFPBeSumVxp/1CweBcxTmal6oTQKuaeSxa4CP62XZAfxm8RT1ROBfA9+PiJ3AMmAM8P+KzRseCZoPvHiRvJJUsRzPVamc3iNlVjwFvBr4jYjo0USbYxSOPt0fEb+YaxoRw4BbKXxprOFzUkppI/BHxU03UDi6BDCo+O9kYGdK6WBE3A1MB/7XJb8pSapCjucqZx7plxoREYuLP04p/ntXROwH9qeUVrbBS/47YCWwNiL+msKp4RHATSmlf1Ns810KKzksjYjvU1jt4XvAYeCvi7knAP8D+DGFlSM6A9+g8EWv2pTS2Yj4CfDbxQ+bW4ATEfGPwG8Cf59SerQN3p8kZeF4LhVY9EuNa3gRl/NLnq0E5pT6xVJKL0fETApLtj0IdAd2Ue8oTUrp6Yj4EvBnFE7nnqJwROhPUkrvFZu9D+ymcDRoMPApsAn4ckppQ7HN1ykcZVoMjC+2B5iTUlpb6vcmSZk5nktApJRyZ5CUQUSMpHD06K6U0tO580iSWsfxXM3hnH6permcmyRVBsdzXZRFv1S9JgPvppQ+zB1EknRJHM91UU7vkSRJkiqcR/olSZKkCmfRL0mSJFU4i35JkiSpwln0S5IkSRXOol+SJEmqcBb9kiRJUoWz6JckSZIq3P8H8QwYM8GcETsAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig = plt.figure(figsize=(12,12))\n", "\n", "gs = gridspec.GridSpec(2, 2, height_ratios=[1,1], width_ratios=[1,1],\n", " hspace=0.075, wspace=0.25)\n", "ax_rays = plt.subplot(gs[0,0])\n", "ax_error = plt.subplot(gs[1,0])\n", "\n", "ax_deriv = plt.subplot(gs[0,1])\n", "ax_deriv_e = plt.subplot(gs[1,1])\n", "\n", "ax_rays.set_ylabel(r'$\\cos\\alpha$')\n", "ax_error.set_xlabel(r'$1+\\cos\\psi$')\n", "ax_error.set_ylabel(r'$|\\varepsilon|$')\n", "\n", "ax_deriv.set_ylabel(r'$\\partial\\cos\\alpha/\\partial\\cos\\psi/(1-u)$')\n", "ax_deriv_e.set_xlabel(r'$1+\\cos\\psi$')\n", "\n", "ax_rays.yaxis.set_major_locator(MultipleLocator(0.2))\n", "ax_rays.yaxis.set_minor_locator(MultipleLocator(0.05))\n", "ax_error.yaxis.set_major_locator(MultipleLocator(0.2))\n", "ax_error.yaxis.set_minor_locator(MultipleLocator(0.05))\n", "\n", "###\n", "\n", "def add_u(u, color, lw=0.5):\n", " \n", " cos_psis = np.linspace(-0.999, 0.999, 1000)\n", " cos_alphas, deriv = rayXpanda.inversion.invert_vec(np.copy(cos_psis), u)\n", " cos_psis_defl, deriv_defl = rayXpanda.deflection.deflect_vec(np.copy(cos_alphas), u)\n", "\n", " ax_rays.plot(1.0+cos_psis, cos_alphas, color=color, ls='-', lw=lw, label=r'inversion: $u=%.3f$'%u)\n", " ax_deriv.plot(1.0+cos_psis, deriv, color=color, ls='-', lw=lw)\n", " \n", " ax_rays.plot(1.0+cos_psis_defl, cos_alphas, color=color, ls='-.', lw=lw, label=r'deflection: $u=%.3f$'%u)\n", " ax_deriv.plot(1.0+cos_psis_defl, deriv_defl, color=color, ls='-.', lw=lw)\n", " \n", " ax_error.plot(1.0+cos_psis, np.abs((np.arccos(cos_psis) - np.arccos(cos_psis_defl)))/np.arccos(cos_psis), color=color, ls='-', lw=lw)\n", "\n", " ax_deriv_e.plot(1.0+cos_psis, np.abs((deriv - deriv_defl))/deriv, color=color, ls='-', lw=lw)\n", " \n", "###\n", "\n", "add_u(5./8., 'r', 0.5)\n", "add_u(0.5, 'k', 0.5)\n", "add_u(0.4, 'b', 0.5)\n", " \n", "###\n", "\n", "ax_rays.set_xscale('log')\n", "ax_rays.set_xticklabels([])\n", "ax_rays.set_xlim([1.0e-2,2.0])\n", "ax_rays.set_ylim([-0.4,1.0])\n", "\n", "ax_error.set_yscale('log')\n", "ax_error.set_xscale('log')\n", "ax_error.set_ylim([1.0e-8,0.1])\n", "ax_error.set_xlim([1.0e-2,2.0])\n", "\n", "locmaj = LogLocator(base=100.0, numticks=100) \n", "ax_error.yaxis.set_major_locator(locmaj) \n", "\n", "locmin = LogLocator(base=100.0, subs=[0.02, 0.04, 0.06, 0.08, 0.1, 0.2, 0.4, 0.6, 0.8], \n", " numticks=100) \n", "ax_error.yaxis.set_minor_locator(locmin) \n", "ax_error.yaxis.set_minor_formatter(NullFormatter()) \n", "\n", "ax_deriv.set_xscale('log')\n", "ax_deriv.set_yscale('log')\n", "ax_deriv.set_ylim([1.0,10.0])\n", "ax_deriv.set_xlim([1.0e-2,2.0])\n", "ax_deriv.set_xticklabels([])\n", "\n", "ax_deriv_e.set_yscale('log')\n", "ax_deriv_e.set_xscale('log')\n", "ax_deriv_e.set_ylim([1.0e-8,1.0])\n", "ax_deriv_e.set_xlim([1.0e-2,2.0])\n", "\n", "locmaj = LogLocator(base=100.0, numticks=100) \n", "ax_deriv_e.yaxis.set_major_locator(locmaj) \n", "\n", "locmin = LogLocator(base=100.0, subs=[0.02, 0.04, 0.06, 0.08, 0.1, 0.2, 0.4, 0.6, 0.8], \n", " numticks=100) \n", "ax_deriv_e.yaxis.set_minor_locator(locmin) \n", "ax_deriv_e.yaxis.set_minor_formatter(NullFormatter()) \n", "ax_deriv.yaxis.set_minor_formatter(NullFormatter())\n", "\n", "for x in [ax_rays, ax_error, ax_deriv, ax_deriv_e]:\n", " x.tick_params(which='major', colors='black', length=8)\n", " x.tick_params(which='minor', colors='black', length=4)\n", " x.xaxis.set_tick_params(which='both', width=1)\n", " x.yaxis.set_tick_params(which='both', width=1)\n", " plt.setp(x.spines.values(), linewidth=1, color='black')\n", " [i.set_color(\"black\") for i in x.get_xticklabels()]\n", " [i.set_color(\"black\") for i in x.get_yticklabels()]\n", " \n", "_ = ax_rays.legend(loc=2, frameon=False, fontsize=10.0)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The fractional error is defined for the left panels as $|\\varepsilon|\\mathrel{:=}|\\psi-\\tilde{\\psi}|/\\psi$. This choice is made because the map $\\cos\\psi\\mapsto\\cos\\alpha$ encoded by [invert()](rayXpandaAPI.rst#rayXpanda.inversion.invert) is defined for primary images, returning for $\\cos\\psi\\in[-1,1]$ a *positive* value of the derivative $$\\mathcal{D}=\\frac{1}{1-u}\\frac{\\partial\\cos\\alpha}{\\partial\\cos\\psi},$$ with range $\\cos\\alpha\\in[a,1]$ where $a\\leq-1$.\n", "\n", "However, the map encoded by [deflect()](rayXpandaAPI.rst#rayXpanda.deflection.deflect) has a domain $\\cos\\alpha\\in[b,1]$ for primary images where $b\\sim a$ is not known without additional calculation. Thus if $\\cos\\alpha\\notin[b,1]$ (i.e., $\\cos\\alpha