{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Validation of SIS model (Numerical vs Analytical)\n", "\n", "* This is in reference to Haris et al. (2018)" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "from ler.lens_galaxy_population import OpticalDepth\n", "import matplotlib.pyplot as plt\n", "%matplotlib inline \n", "\n", "from astropy.cosmology import LambdaCDM\n", "cosmo = LambdaCDM(H0=70, Om0=0.3, Ode0=0.7)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "* optical_depth_SIS_hemanta is the numerical solution of the SIS model" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "z_to_Dc interpolator will be generated at ./interpolator_pickle/z_to_Dc/z_to_Dc_3.pickle\n", "Dc_to_z interpolator will be generated at ./interpolator_pickle/Dc_to_z/Dc_to_z_3.pickle\n", "angular_diameter_distance interpolator will be generated at ./interpolator_pickle/angular_diameter_distance/angular_diameter_distance_3.pickle\n", "differential_comoving_volume interpolator will be generated at ./interpolator_pickle/differential_comoving_volume/differential_comoving_volume_3.pickle\n", "velocity_dispersion_gengamma interpolator will be generated at ./interpolator_pickle/velocity_dispersion_gengamma/velocity_dispersion_gengamma_2.pickle\n", "optical_depth_SIS_hemanta interpolator will be generated at ./interpolator_pickle/optical_depth_SIS_hemanta/optical_depth_SIS_hemanta_1.pickle\n" ] } ], "source": [ "od = OpticalDepth(\n", " z_min=0.001, \n", " z_max=40,\n", " optical_depth_function=\"optical_depth_SIS_hemanta\",\n", " sampler_priors=dict(\n", " velocity_dispersion=\"velocity_dispersion_gengamma\",\n", " axis_ratio=\"axis_ratio_rayleigh\",\n", " ),\n", " sampler_priors_params=dict(\n", " velocity_dispersion=dict(vd_min=0., vd_max=600),\n", " axis_ratio=dict(q_min=0.2, q_max=1),\n", " ),\n", " cosmology=cosmo,\n", " directory=\"./interpolator_pickle\",\n", " #create_new_interpolator=True,\n", " )" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZcAAAFzCAYAAADys0SZAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/OQEPoAAAACXBIWXMAAA9hAAAPYQGoP6dpAABiXElEQVR4nO3dd3xTdffA8c9N0qRAaWsFOmS0QNlTkAIOVKpFUKiLIQ8CIvgoqFgXG8FRWTIURUXByVKGIqDYgiiUyp5ltGykZZRuOpJ8f3/wkJ+hKXSn47xfr74g957vzbmnIYe7NaWUQgghhChGOmcnIIQQouKR5iKEEKLYSXMRQghR7KS5CCGEKHbSXIQQQhQ7aS5CCCGKnTQXIYQQxU6aixBCiGJncHYCFZnVauWff/6hevXqaJrm7HSEEKLIlFKkpqbi5+eHTpf39ok0lxL0zz//UKdOHWenIYQQxe706dPUrl07z/lOby5z585l2rRpxMfH07p1az788EM6dOiQZ/yyZcsYP348J06cIDAwkClTptC9e3fbfKUUEydO5PPPPycpKYk777yTTz75hMDAQFtMYmIiL774Ij///DM6nY7HH3+c2bNn4+bmZrecGTNm8Nlnn3Hy5Elq1KjBCy+8wNixY/O9btWrVweu/hLc3d3zNcZisRAXF0eDBg3Q6/X5fq/KQGrjmNQlb1Ibx4pSl5SUFOrUqWP7fsuTcqLFixcro9GovvzyS3XgwAE1dOhQ5enpqRISEhzGb968Wen1ejV16lR18OBBNW7cOOXi4qL27dtni3n//feVh4eHWrlypdqzZ4/q2bOnCggIUFeuXLHFdOvWTbVu3Vpt3bpV/fnnn6phw4aqX79+du/14osvqsaNG6tVq1apY8eOqe3bt6vffvutQOuXnJysAJWcnJzvMWazWcXExCiz2Vyg96oMpDaOSV3yJrVxrCh1ye/3mlObS4cOHdTw4cNtry0Wi/Lz81Ph4eEO43v37q169OhhNy0oKEg999xzSimlrFar8vHxUdOmTbPNT0pKUiaTSS1atEgppdTBgwcVoLZt22aLWbt2rdI0TZ09e9YWYzAY1KFDh4q0ftJcipfUxjGpS96kNo6VRnNx2m6x7OxsduzYwejRo23TdDodwcHBREVFORwTFRVFWFiY3bSQkBBWrlwJwPHjx4mPjyc4ONg238PDg6CgIKKioujbty9RUVF4enrSvn17W0xwcDA6nY7o6GgeffRRfv75Z+rXr8/q1avp1q0bSimCg4OZOnUqXl5eea5TVlYWWVlZttcpKSnA1U1Qi8WSr7pYLBasVmu+4ysTqY1jUpe8SW0cK0pd8jvGac3l4sWLWCwWvL297aZ7e3tz6NAhh2Pi4+MdxsfHx9vmX5t2o5hatWrZzTcYDHh5edlijh07xsmTJ1m2bBlff/01FouFV155hSeeeILIyMg81yk8PJxJkyblmh4XF2d3POdGrFYriYmJxMbG3vBMjMpIauOY1CVvUhvHilKXtLS0fMU5/YB+WWS1WsnKyuLrr7+mUaNGAHzxxRe0a9eOw4cP07hxY4fjRo8ebbdlde3AV4MGDRwe0FdKYTab7f4ncO1/FLVr15YDkNepbLXR6/UYDIabnsZusViIjY2lYcOGlaIuBSG1cawodbm2R+ZmnNZcatSogV6vJyEhwW56QkICPj4+Dsf4+PjcMP7anwkJCfj6+trFtGnTxhZz/vx5u2WYzWYSExNt4319fTEYDLbGAtC0aVMATp06lWdzMZlMmEymXNP1en2uX2B2djbnzp0jIyPDbrpSCqvVypkzZ+TamOtUxtpUrVoVX19fjEbjDeN0Op3Dz5mQ2uSlsHXJb7zTmovRaKRdu3ZEREQQGhoKXN1iiIiIYMSIEQ7HdOrUiYiICEaOHGmbtn79ejp16gRAQEAAPj4+RERE2JpJSkoK0dHRPP/887ZlJCUlsWPHDtq1awdAZGQkVquVoKAgAO68807MZrPtVD2AI0eOAFCvXr0ir7vVauX48ePo9Xr8/PwwGo22L0ulFFlZWZhMpkrzBZpflak2Simys7O5cOECx48fJzAwUHbriPKlwKcKFKPFixcrk8mkFi5cqA4ePKiGDRumPD09VXx8vFJKqQEDBqhRo0bZ4jdv3qwMBoOaPn26iomJURMnTnR4KrKnp6datWqV2rt3r+rVq5fDU5Hbtm2roqOj1V9//aUCAwPtTkW2WCzq9ttvV/fcc4/auXOn2r59uwoKClIPPPBAgdYvr7Mqrly5og4ePKjS09NzjbFarSojI0NZrdYCvVdlUBlrk56erg4ePGj3+b2enBGVN6mNYxX+VGSllPrwww9V3bp1ldFoVB06dFBbt261zevSpYsaOHCgXfzSpUtVo0aNlNFoVM2bN1e//PKL3Xyr1arGjx+vvL29lclkUl27dlWHDx+2i7l06ZLq16+fcnNzU+7u7mrw4MEqNTXVLubs2bPqscceU25ubsrb21sNGjRIXbp0qUDrdrPm4ugLozJ+geZXZazNjT4r18gXaN6kNo6VRnPRlFLKudtOFVdKSgoeHh4kJyfbHdDPzMzk+PHjBAQE4OrqajdGKUVmZiaurq4VftdPQVXG2tzos3KNxWLh6NGjBAYGynGF60htHCtKXfL6Xrue7MQVQohK5nxKJj/sTyLbbC2x95DmIsokf39/Zs2a5ew0hKiQ3lt7mPnbLzF6xf4Sew9pLqJALly4wPPPP0/dunUxmUz4+PgQEhLC5s2bbTHXN4Y9e/bQs2dPatWqhaurK/7+/vTp0yfXKeEVyb333mt3VqMQZcXm2Ius33uCDrpDDO5c9LNf8yIXUYoCefzxx8nOzuarr76ifv36JCQkEBERwaVLlxzGX7hwga5du/Lwww/z66+/4unpyYkTJ/jpp59IT08v5eyFqNyyzVbGr9rPcMMqXjSsxHooGepOLpH3ki2XMkIpRUa2+X8/ln/9veR/8ntOR1JSEn/++SdTpkzhvvvuo169enTo0IHRo0fTs2dPh2M2b95McnIy8+fPp23btgQEBHDfffcxc+ZMAgICbvh+GRkZPPPMM1SvXp26devy2Wef2c0/ffo0vXv3xtPTEy8vL3r16sWJEyds8wcNGkRoaCjvvfce3t7eeHp6MnnyZMxmM6+//jpeXl7Url2bBQsW2C33zTffpFGjRlStWpX69eszfvx4cnJybPPfeust2rRpwzfffIO/vz8eHh707duX1NRU2/v+8ccfzJ49G03T0DSNEydOYLFYGDJkCAEBAVSpUoXGjRsze/bsfNVeiOLw+Z/H4OJR/mtYDYC6rf1NRhSebLmUEVdyLDSb8KtT3vvg5BCqGm/+UXBzc8PNzY2VK1fSsWNHh3cjuJ6Pjw9ms5kVK1bwxBNPFOgsrxkzZvD2228zZswYfvjhB1544QU6duxIq1atyMnJISQkhE6dOvHnn39iMBh455136NatG3v37rVd0R4ZGUnt2rXZtGkTmzdvZsiQIWzZsoV77rmH6OholixZwnPPPccDDzxge/BR9erVWbhwIX5+fuzbt4+hQ4dSvXp13njjDVtucXFxrFy5ktWrV3P58mV69+7N+++/z7vvvsvs2bM5cuQILVq0YPLkq/8rrFmzpu3WNcuWLePWW29ly5YtDBs2DF9fX3r37p3vughRGGcuZ/Bh5BG+MCzABTNpfndSpXH3mw8sJNlyEflmMBhYuHAhX331FZ6entx5552MGTOGvXv35jmmY8eOjBkzhqeeeooaNWrw0EMPMW3atFy38XGke/fuvPDCCzRs2JA333yTGjVq8McffwCwZMkSrFYr8+fPp2XLljRt2pQFCxZw6tQpNm7caFuGl5cXc+bMoXHjxjzzzDM0btyYjIwMxowZQ2BgIKNHj8ZoNPLXX3/ZxowbN47OnTvj7+/PI488wmuvvcbSpUvtcrNarSxcuJAWLVpw9913M2DAACIiIoCrd+I2Go1UrVoVHx8ffHx80Ov1uLi4MGnSJNq3b09AQAD9+/dn8ODBuZYtREmY9PNBHrT8xZ36AyhDFRJufw1K8JR+2XIpI6q46Dk4OeR/13Jk4epaerc4qeKS//PcH3/8cXr06MGff/7J1q1bWbt2LVOnTmX+/PkMGjTI4Zh3332XsLAwIiMjiY6OZt68ebz33nts2rSJli1b5vlerVq1sv1d0zR8fHy4cOECcPUkgdjY2FxPw8vMzCQuLs72unnz5na3TfH29qZFixa213q9nltvvdXu5IIlS5YwZ84c4uLiSEtLw2w25zqf39/f3+69fX1983WCwty5c/nyyy85deoUV65cITs723arIiFKSuShBLYePE6k6VsA1N2vkuPmV6LvKVsuZYSmaVQ1Gv73o//X30v+p6BNzNXVlQceeIDx48ezZcsWBg0axMSJE2845tZbb+XJJ59k+vTpxMTE4Ofnx/Tp0284xsXFJVeNrNar5+WnpaXRrl07du/ebfdz5MgRnnrqqRsu40bLjYqKon///nTv3p3Vq1eza9cuxo4dS3Z2dr5zy8vixYt57bXXGDJkCL/99hu7d+9m8ODBuZYtRHHKzLEw8acDvGpYSk0tGW4NRHUcXuLvK1suosiaNWtme2BbfhiNRho0aFCks8Vuv/12lixZQq1atW54lXBBbdmyhXr16jF27FjbtJMnTxZ4OUajMddDlTZv3kznzp154YUXbNP+vZUlREn4eEMsHpcPMMD0+9UJPWaA4ebHS4tKtlxEvl26dIn777+fb7/9lr1793L8+HGWLVvG1KlT6dWrl8Mxq1ev5j//+Q+rV6/myJEjHD58mOnTp7NmzZo8x+RH//79qVGjBr169eLPP//k+PHjbNy4kZdeeokzZ84UermBgYGcOnWKxYsXExcXx5w5c1ixYkWBl+Pv7090dDQnTpzg4sWLWK1WAgMD2b59O7/++itHjhxh/PjxbNu2rdC5CnEzxy+m89kfsbzj8iV6rNCyN9TvUirvLc1F5JubmxtBQUHMnDmTe+65hxYtWjB+/HiGDh3KRx995HBMs2bNqFq1Kq+++ipt2rShY8eOLF26lPnz5zNgwIBC51K1alU2bdpE3bp1eeyxx2jatClDhgwhMzOzSFsyPXv25JVXXmHEiBG0adOGLVu2MH78+AIv57XXXkOv19OsWTNq1qzJqVOneO6553jsscfo06cPQUFBXLp0yW4rRojipJRiwqr9PMF62uiOoUzu8OA7pfb+cuPKEiQ3rixelbE2cuPKoqnMtVmz7xwTvttApOlV3LUM6D4dOgwF5MaVQgghCiEty8zknw8yxuW7q43Ftw20f6ZUc5DmIoQQFcyciKP4p+3kMf1fKDR4eCboSnfLTc4WE0KICuRwfCrf/HWUnwxXb2uk3TEEbru91POQLRchhKgglFKMX7WfwdpqAnVnoVpNuL/gJ6QUB2kuQghRQazYdZZ/ThzmRcP/Tp9/8F2o4umUXGS3mBBCVADJV3J4b00M4YavqKJlg//d0Mp5N0SVLRchhKgAZvx2mLYZW3hAvxOlc7l6Jb4TT9mXLRchhCjn9p1J5seth/nV+DUAWucXoWZjp+YkWy5CCFGOWa2Kcav2M0K/gtraRfCsC/e87uy0pLmIssnf359Zs2aVyLLl+faiIlm87TQZZ/bzrH7N1QkPTQNjVecmhTQXUUAXLlzg+eefp27duphMJnx8fAgJCWHz5s22mOsbw549e+jZsye1atXC1dUVf39/+vTpk6/nn5SE5cuX8/bbbzvlvYUoTpfSspiyNoZ3XL7ERbNA4x7QuJuz0wLkmIsooMcff5zs7Gy++uor6tevT0JCAhEREVy6dMlh/IULF+jatSsPP/wwv/76K56enpw4cYKffvqpSLfcL4zs7GyMRiNeXl6l+r5ClJQp6w7RNTuSIOMhlEtVtIemODslG9lyKSuUgux05/zk896lSUlJ/Pnnn0yZMoX77ruPevXq0aFDB0aPHk3Pnj0djtm8eTPJycnMnz+ftm3bEhAQwH333cfMmTMJCAi44ftlZGTwzDPPUL16derWrctnn31mN//NN9+kUaNGVK1alfr16zN+/HhycnJs89966y3atGnD/Pnz7W78eP1usY8//pjAwEBcXV3x9vbmiSeeyFc9hHCm7ScS+XX7Ica4fA+A1uVN8Kzj5Kz+n2y5lBU5GfCeHxpQpbTfe8w/YKx20zA3Nzfc3NxYuXIlHTt2xGS6+QOHfHx8MJvNrFixgieeeKJAdzOeMWMGb7/9NmPGjOGHH37ghRdeoGPHjrbHH1evXp2FCxfi5+fHvn37GDp0KNWrV+eNN96wLSM2NpYff/yR5cuXO7z76/bt23nppZf45ptv6Ny5M4mJifz555/5zlEIZzBbrIxbuZ/XDUuooaVAzSbQsWw9vkG2XES+GQwGFi5cyFdffYWnpyd33nknY8aMYe/evXmO6dixI2PGjOGpp56iRo0aPPTQQ0ybNo2EhISbvl/37t154YUXaNiwIW+++SY1atTgjz/+sM0fN24cnTt3xt/fn0ceeYTXXnuNpUuX2i0jOzubr7/+mrZt29qa0r+dOnWKatWq8fDDD1OvXj3atm3LSy+9VICqCFH6vo46iSlhF08ZIq9O6PEBGIzOTeo6suVSVrhUhTH/OOeZJS75P7Pk8ccfp0ePHvz5559s3bqVtWvXMnXqVObPn8+gQYMcjnn33XcJCwsjMjKS6Oho5s2bx3vvvcemTZto2bJlnu/172agaRo+Pj5cuHDBNm3JkiXMmTOHuLg40tLSMJvNuZ4vUa9ePWrWrJnnezzwwAPUq1eP+vXr061bN7p168ajjz5K1arOP9tGCEcSUjKZvT6G71y+RIeC1k+B/53OTisX2XIpKzTt6q4pZ/wUsIm5urrywAMPMH78eLZs2cKgQYOYOHHiDcfceuutPPnkk0yfPp2YmBj8/PyYPn36Dce4uLhcVyINq9UKQFRUFP3796d79+6sXr2aXbt2MXbsWLKzs+3GVKt249191atXZ+fOnSxatAhfX18mTJhA69atSUpKuuE4IZzl3V9ieNS8lha6EyhXT3iwbJ75KM1FFFmzZs0KdOaX0WikQYMGRTpbbMuWLdSrV4+xY8fSvn17AgMDOXnyZKGWZTAYCA4OZurUqezdu5cTJ04QGRlZ6NyEKCmbYy8StecArxqWAaAFT4RqNZyclWOyW0zk26VLl3jyySd55plnaNWqFdWrV2f79u1MnTqVXr16ORyzevVqFi9eTN++fWnUqBFKKX7++WfWrFnDggULCp1LYGAgp06dYvHixdxxxx388ssvrFixosDLWb16NceOHeOee+7hlltuYc2aNVitVho3du6tM4S4XrbZyvhV+xnv8i3VtStwW3u4fZCz08qTNBeRb25ubgQFBTFz5kzi4uLIycmhTp06DB06lDFjxjgc06xZM6pWrcqrr77K6dOnMZlMBAYGMn/+fAYMGFDoXHr27Mkrr7zCiBEjyMrKokePHowfP5633nqrQMvx9PRk+fLlvPXWW2RmZhIYGMiiRYto3rx5oXMToiR8/ucx/C5tpacxCqXp0B7+AHRld+eTplQ+L3IQBZaSkoKHhwfJycl2B5ozMzM5fvy43bUX1zjlgH45URlrc6PPyjUWi4WjR48SGBjo8HTryqyi1ObM5Qwe+uB3VmmvU18XD0H/hSJcMFmUuuT1vXa9stv2hBBCADD554MMsq6ivi4e5eYD9411dko3Jc1FCCHKsMhDCRyK2csIwyoAtG7vgWveWwxlhRxzEUKIMioj28z4Fft517AQk5YD9e+F5o85O618KRNbLnPnzsXf3x9XV1eCgoL4+++/bxi/bNkymjRpgqurKy1btmTNmjV285VSTJgwAV9fX6pUqUJwcDBHjx61i0lMTKR///64u7vj6enJkCFDSEtLs80/ceIEmqbl+tm6dWvxrbgQQtzA7IijtEzdxL36PSi9Ebo79+mSBeH05rJkyRLCwsKYOHEiO3fupHXr1oSEhOR5O/YtW7bQr18/hgwZwq5duwgNDSU0NJT9+/fbYqZOncqcOXOYN28e0dHRVKtWjZCQEDIzM20x/fv358CBA6xfv57Vq1ezadMmhg0bluv9fv/9d86dO2f7adeuXbGtu5xLIW5GPiOVV8y5FBb/eYBJLgsB0O56BWo0dG5SBaGcrEOHDmr48OG21xaLRfn5+anw8HCH8b1791Y9evSwmxYUFKSee+45pZRSVqtV+fj4qGnTptnmJyUlKZPJpBYtWqSUUurgwYMKUNu2bbPFrF27Vmmaps6ePauUUur48eMKULt27Sr0uiUnJytAJScn2003m83q4MGD6uLFi7nGWK1WlZGRoaxWa6Hft6KqjLW5ePGiOnjwoDKbzXnGmM1mFRMTc8OYyqq81sZisarQuX+pBWOfUGqiu1Kz2yqVfaXYll+UuuT1vXY9px5zyc7OZseOHYwePdo2TafTERwcTFRUlMMxUVFRhIWF2U0LCQlh5cqVABw/fpz4+HiCg4Nt8z08PAgKCiIqKoq+ffsSFRWFp6cn7du3t8UEBwej0+mIjo7m0UcftU3v2bMnmZmZNGrUiDfeeCPPW8sDZGVlkZWVZXudkpICXD3tz2Kx2MW6u7tz/vx5lFJUrVrVdmqtUors7GyUUpXmdNv8qky1UUqRkZHBhQsXbKd7Xv8ZusZisWC1WvOcX5mV19p8//cprKe387RxPQCW7jNA5wLFtB5FqUt+xzi1uVy8eBGLxYK3t7fddG9vbw4dOuRwTHx8vMP4+Ph42/xr024UU6tWLbv5BoMBLy8vW4ybmxszZszgzjvvRKfT8eOPPxIaGsrKlSvzbDDh4eFMmjQp1/S4uDjc3NzspimlyMnJ4ezZs7m+KC0WS7k+J78kVabaXGuiKSkppKam5hlntVpJTEwkNjYWXRm+qM4ZymNtEjPMTFtzgsUuX6DTFMn+3Tln9oXrjhsXRVHq8u9j0zciZ4vloUaNGnZbSHfccQf//PMP06ZNy7O5jB492m5MSkoKderUoUGDBnlebGSxWOwecGWxWDh16hR169atNF+i+VXZauPi4pKv9bRYLMTGxtKwYcNKUZeCKI+1GblkD72ta2jmchJV5RbcHptFYDHfP6wodbm2R+ZmnNpcatSogV6vz/Vsj4SEBHx8fByO8fHxuWH8tT8TEhLw9fW1i2nTpo0t5voTBsxmM4mJiXm+L0BQUBDr16/Pc77JZHL4AC29Xp/nL1Cv12M0/v9zGCwWCwaDgWrVqpWbfwylRWqTN51Od8PPWWVWnmqz6cgFdu7dy3rTDwBoD76D3t37JqMKp7B1yW+8U7cTjUYj7dq1IyIiwjbNarUSERFBp06dHI7p1KmTXTzA+vXrbfEBAQH4+PjYxaSkpBAdHW2L6dSpE0lJSezYscMWExkZidVqJSgoKM98d+/ebdewhBCiuGTmWBi3Yh+TXRZSVcuCendCm/7OTqvQnL5bLCwsjIEDB9K+fXs6dOjArFmzSE9PZ/DgwQA8/fTT3HbbbYSHhwPw8ssv06VLF2bMmEGPHj1YvHgx27dvtz1fXdM0Ro4cyTvvvENgYCABAQGMHz8ePz8/QkNDAWjatCndunVj6NChzJs3j5ycHEaMGEHfvn3x8/MD4KuvvsJoNNK2bVsAli9fzpdffsn8+fNLuUJCiMrgo8hYmiX/QVfjLpTOBe3hWeXmmhZHnN5c+vTpw4ULF5gwYQLx8fG0adOGdevW2Q7Inzp1yu6AU+fOnfn+++8ZN24cY8aMITAwkJUrV9KiRQtbzBtvvEF6ejrDhg0jKSmJu+66i3Xr1tnd+O+7775jxIgRdO3aFZ1Ox+OPP86cOXPscnv77bc5efIkBoOBJk2asGTJEp544okSrogQorI5mpDKd5v2se7aNS13h0HNRs5NqojkrsglKL93D/23inIX15IgtXFM6pK38lAbq1XR97OtdD/zAYMMv6G8GqA9vwVcHN8FuzjIXZGFEKKC+2HHGbJO/s3T+qsnC2kPzyzRxlJanL5bTAghKqtLaVlMWbOfb/53TQut+kL9Ls5Oq1jIlosQQjjJu2tieDT7Z5rprl7TQsi7zk6p2MiWixBCOMGWuItE79z9/9e0PPA2FPPFks4kWy5CCFHKsswWxi2/7pqWtv9xdlrFSpqLEEKUsk82xtHo8ka66q9e08LDM8v1NS2OyG4xIYQoRccupPH1hn2s/fdzWmo2dm5SJUC2XIQQopQopRi3cj8vaYvx1pJQXvXh7lednVaJkOYihBClZMWus6Qfi65w17Q4IrvFhBCiFCRlZPP+6v0stLum5V5np1ViZMtFCCFKwftrD9Ez66cKeU2LI7LlIoQQJezv44n8uW1nhb2mxRHZchFCiBKUbbYydvneCn1NiyPSXIQQogR9/ucxGlzaUKGvaXFEdosJIUQJOXkpnS8j9rCmgl/T4ohsuQghRAlQSjF+1QFepOJf0+KINBchhCgBq/eeI/loVKW4psUR2S0mhBDFLPlKDu/8tI8FleSaFkdky0UIIYrZ1HWHeCRzVaW5psUR2XIRQohitO1EIhujd1Sqa1ockS0XIYQoJllmC6N+2FPprmlxRJqLEEIUk7kb4miW+PvVa1r0xkpzTYsjsltMCCGKwZGEVBZt3MVal68B0O5+rdJc0+KIbLkIIUQRWa2KUT/u5U3dt9TQUlA1m8Jdrzg7LaeS5iKEEEX0bfRJqp3ZxBP6TSg0tJ4fgsHo7LScSnaLCSFEEfyTdIU5a3ezwvAFAFrQc1DnDidn5Xyy5SKEEIWklGL8yv08Z11CHd0FlEdtuH+8s9MqE6S5CCFEIf2y7xwXDkfxjH4tANrDs8Dk5tykygjZLSaEEIWQlJHN26v2sNDlc/Sagpa9IfABZ6dVZsiWixBCFMJ7a2J4PHM5TXWnUFW8oFu4s1MqU2TLRQghCmhL7EW27/ibtcYVAGgPTamUt3i5EWkuQghRAJk5FsYs38MUl/mYtBxoGAwtn3R2WmWO7BYTQogCmPX7UTon/0KQ7hDKpVqlvsXLjciWixBC5NOBf5L56c/trHP5HgCt63jwrOvkrMom2XIRQoh8MFusjPphL2/pv8RduwK3tYcOw5ydVpklzUUIIfJh4ZYT1I5fz4P6HSidAXp+CDq9s9Mqs2S3mBBC3MTpxAw+/20nq10WAqDdFQbezZybVBlXJrZc5s6di7+/P66urgQFBfH333/fMH7ZsmU0adIEV1dXWrZsyZo1a+zmK6WYMGECvr6+VKlSheDgYI4ePWoXk5iYSP/+/XF3d8fT05MhQ4aQlpbm8P1iY2OpXr06np6eRVpPIUT5o5RizIp9vKK+oaaWjKrRCO55zdlplXlOby5LliwhLCyMiRMnsnPnTlq3bk1ISAjnz593GL9lyxb69evHkCFD2LVrF6GhoYSGhrJ//35bzNSpU5kzZw7z5s0jOjqaatWqERISQmZmpi2mf//+HDhwgPXr17N69Wo2bdrEsGG595/m5OTQr18/7r777uJfeSFEmbdi11nMcX/Q17AR4H93PDY5N6nyQDlZhw4d1PDhw22vLRaL8vPzU+Hh4Q7je/furXr06GE3LSgoSD333HNKKaWsVqvy8fFR06ZNs81PSkpSJpNJLVq0SCml1MGDBxWgtm3bZotZu3at0jRNnT171m7Zb7zxhvrPf/6jFixYoDw8PAq0bsnJyQpQycnJ+R5jNptVTEyMMpvNBXqvykBq45jUJW9Frc3F1EzV8a2f1PHxgUpNdFdqdVgxZ+gcRalLfr/XnHrMJTs7mx07djB69GjbNJ1OR3BwMFFRUQ7HREVFERYWZjctJCSElStXAnD8+HHi4+MJDg62zffw8CAoKIioqCj69u1LVFQUnp6etG/f3hYTHByMTqcjOjqaRx99FIDIyEiWLVvG7t27Wb58+U3XJysri6ysLNvrlJQUACwWCxaL5abjr8VardZ8x1cmUhvHpC55K2ptJv98gIE5S/A3JKCq+2G9bxxUgDoXpS75HePU5nLx4kUsFgve3t520729vTl06JDDMfHx8Q7j4+PjbfOvTbtRTK1atezmGwwGvLy8bDGXLl1i0KBBfPvtt7i7u+drfcLDw5k0aVKu6XFxcbi55e9OqVarlcTERGJjY9HpnL7XskyR2jgmdclbUWqz/UwGsXu3MMP4CwBn2rxK+qkEIKEEMi1dRalLXsemrydni+Vh6NChPPXUU9xzzz35HjN69Gi7raqUlBTq1KlDgwYN8t2gLBYLsbGxNGzYEL1eTnP8N6mNY1KXvBW2NulZZj5b8Qcfu3yOQbNibfYofvcOLsFMS1dRPjPX9sjcjFObS40aNdDr9SQk2P9PICEhAR8fH4djfHx8bhh/7c+EhAR8fX3tYtq0aWOLuf6EAbPZTGJiom18ZGQkP/30E9OnTweunjFitVoxGAx89tlnPPPMM7lyM5lMmEy5D/Tp9foC/QJ1Ol2Bx1QWUhvHpC55K0xtZkce5qG05bRwOYFy9UTXfSpUsNoW9jOT33inbkMbjUbatWtHRESEbZrVaiUiIoJOnTo5HNOpUye7eID169fb4gMCAvDx8bGLSUlJITo62hbTqVMnkpKS2LFjhy0mMjISq9VKUFAQcPXYzu7du20/kydPpnr16uzevdt2TEYIUfHsPp1E5OYoXjH8AIDWLRzcat1klLie03eLhYWFMXDgQNq3b0+HDh2YNWsW6enpDB58dRP06aef5rbbbiM8/OqzEl5++WW6dOnCjBkz6NGjB4sXL2b79u189tlnAGiaxsiRI3nnnXcIDAwkICCA8ePH4+fnR2hoKABNmzalW7duDB06lHnz5pGTk8OIESPo27cvfn5+tph/2759OzqdjhYtWpRSZYQQpS3HYmXUD3t41zAfVy0H6t8Lrfs5O61yyenNpU+fPly4cIEJEyYQHx9PmzZtWLdune2A/KlTp+wOOHXu3Jnvv/+ecePGMWbMGAIDA1m5cqXdl/4bb7xBeno6w4YNIykpibvuuot169bh6upqi/nuu+8YMWIEXbt2RafT8fjjjzNnzpzSW3EhRJnz2aZjtLr4M51dDqIMVa4+tljueFwomlJKOTuJiiolJQUPDw+Sk5MLdED/6NGjBAYGyv7z60htHJO65K0gtYk9n8rA2T+zxvAqHloGPPgOdH6xlDItXUX5zOT3e03OWxRCVHoWq+KNZXsYr/sSDy0D5dcWgp53dlrlmjQXIUSl99WWE/idXUc3/TaUznD1Fi96px81KNcKVb2IiAgiIiI4f/48VqvVbt6XX35ZLIkJIURpOHUpgy9+/ZufXRYAoN3zOvi0dHJW5V+Bm8ukSZOYPHky7du3x9fXF00OdgkhyimlFKNX7GUM8/HS0lDeLa7eTl8UWYGby7x581i4cCEDBgwoiXyEEKLULNl2Go9jv9DD+PfV3WGhH4PB6Oy0KoQCN5fs7Gw6d+5cErkIIUSpiU/O5ONfoln+7weA+bZ2blIVSIEP6D/77LN8//33JZGLEEKUCqUUY1fs43XrfGpoKahazeCe152dVoWSry2Xf9+M0Wq18tlnn/H777/TqlUrXFxc7GI/+OCD4s1QCCGK2U97/sFwZDWPGLeiNL3sDisB+Wouu3btsnt97QaQ/376oxBClAcX07KYuWory1yuntmq3TUS/No6N6kKKF/NZcOGDSWdhxBClIq3fjrASPN8aupTUDWaoHV509kpVUgFPubyzDPPkJqammt6enq6w9vQCyFEWfHbgXiy9/9EqH4LStOhPfoxGHI/JkMUXYGby1dffcWVK1dyTb9y5Qpff/11sSQlhBDFLflKDlNXRPHutd1hd74Mt7VzclYVV75PRU5JSUEphVKK1NRUuzsMWywW1qxZk+vRwUIIUVa8+8tBXsiaT019MtYajdB1GeXslCq0fDcXT09PNE1D0zQaNWqUa76maQ6fHy+EEM72V+xFLu1cxWPGv1CaDl2vj8HF9eYDRaHlu7ls2LABpRT3338/P/74I15eXrZ5RqORevXq2R60JYQQZcWVHCvhq7fylcsXAGidRkCdO5ycVcWX7+bSpUsXAI4fP06dOnXsHuAlhBBl1YIdl3g2Yz7e+iSsXg3R3TfG2SlVCgW+/Uu9evW4fPkyX3zxBTExMQA0a9aMwYMH223NCCGEs20/cZnUI3/whHETCg1d6MfgUsXZaVUKBd782LRpE/7+/syZM4fLly9z+fJl5syZQ0BAAJs2bSqJHIUQosAycyy8vTyK92y7w4ZD3SAnZ1V5FHjLZfjw4fTp04dPPvnE9nhMi8XCCy+8wPDhw9m3b1+xJymEEAU16/ejDEiej68hEcst9dHfN9bZKVUqBd5yiY2N5dVXX7V77rJerycsLIzY2NhiTU4IIQpj35lkDv+1nD6GjSg06PkRGKs6O61KpcDN5fbbb7cda/m3mJgYWreW21ULIZwr22xl4rItvGv4HIDLjXpD3Y5OzqryKfBusZdeeomXX36Z2NhYOna8+gvbunUrc+fO5f3332fv3r222FatWhVfpkIIkQ/z/ojjiUuf4mdIxOzhz4VWz+Ph7KQqoQI3l379+gHwxhtvOJynaRpKKTRNw2KxFD1DIYTIp8Pxqeza8CMLDFdvtqv1+giVLWeHOUOBm8vx48dLIg8hhCgSi1UxcdlWpuuv7g5THYZBvc5w9KiTM6ucCnWdixBClDWfborjkYRPqG24iNm9Lobgt5ydUqVWqMvsv/nmG+688078/Pw4efIkALNmzWLVqlXFmpwQQuRHzLkUtv3+A/0NEQAYHp0LxmpOzqpyK3Bz+eSTTwgLC6N79+4kJSXZjqt4enoya9as4s5PCCFuKNtsZcLizbyn/xQAdcdQCLjHyVmJAjeXDz/8kM8//5yxY8faXevSvn17uYBSCFHq5kQcpX/iHHy1RMy3NEB7YLKzUxIUorkcP36ctm1zP2/aZDKRnp5eLEkJIUR+7Dp1mRObviVUvwWrpsfw+GdysWQZUeDmEhAQwO7du3NNX7duHU2bNi2OnIQQ4qauZFt4b3EkbxuuPllSd/erULu9k7MS1xT4bLGwsDCGDx9OZmYmSin+/vtvFi1aRHh4OPPnzy+JHIUQIpep62J4IXU2t+jTsHi3Qt8l97V3wnkK3FyeffZZqlSpwrhx48jIyOCpp57Cz8+P2bNn07dv35LIUQgh7GyJu0hW9Jfc57IHi86I/vHPQe/i7LTEvxS4uQD079+f/v37k5GRQVpaGrVq1SruvIQQwqHUzBxmLfmVhYZvAdA/8BbUauLcpEQuhWou11StWpWqVeXgmRCi9Lz7837ezJxJVV0W5rp3YQh63tkpCQfy1Vzatm2Lpmn5WuDOnTuLlJAQQuQl8lACXrs/oZ3LUSwubhge+wTkketlUr6aS2hoqO3vmZmZfPzxxzRr1oxOnToBV++KfODAAV544YUSSVIIIS6nZ/PFDz+xwPADAPruU8GzrpOzEnnJV3OZOHGi7e/PPvssL730Em+//XaumNOnTxdvdkII8T9vrdjJ+OzZGHUWLI16oG/zlLNTEjdQ4O3JZcuW8fTTT+ea/p///Icff/yxUEnMnTsXf39/XF1dCQoK4u+//75pDk2aNMHV1ZWWLVuyZs0au/lKKSZMmICvry9VqlQhODiYo9fdGTUxMZH+/fvj7u6Op6cnQ4YMIS0tzTb/8OHD3HfffXh7e+Pq6kr9+vUZN24cOTk5hVpHIUTh/bznH5oe+pAmutPkuNZA33M25HNXvXCOAjeXKlWqsHnz5lzTN2/ejKura4ETWLJkCWFhYUycOJGdO3fSunVrQkJCOH/+vMP4LVu20K9fP4YMGcKuXbsIDQ0lNDSU/fv322KmTp3KnDlzmDdvHtHR0VSrVo2QkBAyMzNtMf379+fAgQOsX7+e1atXs2nTJoYNG2ab7+LiwtNPP81vv/3G4cOHmTVrFp9//rndVpwQouSdT8lkxcqlDNP/AoBL6Bxwq+nkrMRNqQIKDw9Xrq6u6sUXX1TffPON+uabb9SIESNU1apVVXh4eEEXpzp06KCGDx9ue22xWJSfn1+ey+rdu7fq0aOH3bSgoCD13HPPKaWUslqtysfHR02bNs02PykpSZlMJrVo0SKllFIHDx5UgNq2bZstZu3atUrTNHX27Nk8c33llVfUXXfdle91S05OVoBKTk7O9xiz2axiYmKU2WzO95jKQmrjWEWui9VqVf+dv0GdGt9AqYnuyrL8+QKNr8i1KYqi1CW/32sFPhV51KhR1K9fn9mzZ/Ptt1fPM2/atCkLFiygd+/eBVpWdnY2O3bsYPTo0bZpOp2O4OBgoqKiHI6JiooiLCzMblpISAgrV64Ert77LD4+nuDgYNt8Dw8PgoKCiIqKom/fvkRFReHp6Un79v9/q4jg4GB0Oh3R0dE8+uijud43NjaWdevW8dhjj+W5PllZWWRlZdlep6SkAGCxWPL9VE6LxYLVapWneDogtXGsItdl6fYzdDk+kzqGC2S71UYf8m6B1rMi16YoilKX/I4p1HUuvXv3LnAjceTixYtYLBa8vb3tpnt7e3Po0CGHY+Lj4x3Gx8fH2+Zfm3ajmOsv/DQYDHh5edliruncuTM7d+4kKyuLYcOGMXly3ndcDQ8PZ9KkSbmmx8XF4ebmlue4f7NarSQmJhIbG4tOTrG0I7VxrKLWJT41h02rV/KJYSMKjXMdxnLlVAKQkO9lVNTaFFVR6vLvY9M3UqSLKCuDJUuWkJqayp49e3j99deZPn06b7zh+B5Go0ePttuqSklJoU6dOjRo0AB3d/d8vZ/FYiE2NpaGDRvaPdJASG3yUhHrYrUqpn/+G+G6z66+DnqB2nf2KfByKmJtikNR6nJtj8zNOLW51KhRA71eT0KC/f9EEhIS8PHxcTjGx8fnhvHX/kxISMDX19cupk2bNraY608YMJvNJCYm5nrfOnXqANCsWTMsFgvDhg3j1VdfdfgLMZlMmEymXNP1en2BfoE6na7AYyoLqY1jFa0uX205xuPnPqCmPoVsr8YYgydAIdetotWmuBS2LvmNd+p2otFopF27dkRERNimWa1WIiIibBdoXq9Tp0528QDr16+3xQcEBODj42MXk5KSQnR0tC2mU6dOJCUlsWPHDltMZGQkVquVoKCgPPO1Wq3k5ORgtVoLvrJCiHyJPZ/God8+5yH9NqyaAeOTn4NLwc9EFc7l9N1iYWFhDBw4kPbt29OhQwdmzZpFeno6gwcPBuDpp5/mtttuIzw8HICXX36ZLl26MGPGDHr06MHixYvZvn07n312dfNZ0zRGjhzJO++8Q2BgIAEBAYwfPx4/Pz/bnQaaNm1Kt27dGDp0KPPmzSMnJ4cRI0bQt29f/Pz8APjuu+9wcXGhZcuWmEwmtm/fzujRo+nTpw8uLnL3VSFKgtliJXzxembqFgCg3TsKfFs7NylRKE5vLn369OHChQtMmDCB+Ph42rRpw7p162wH5E+dOmV3wKlz5858//33jBs3jjFjxhAYGMjKlStp0aKFLeaNN94gPT2dYcOGkZSUxF133cW6devsrsP57rvvGDFiBF27dkWn0/H4448zZ84c23yDwcCUKVM4cuQISinq1avHiBEjeOWVV0qhKkJUTnMjj/LMhWm466+Q7dsO413y76280pRS6mZB15/6eyMffPBBkRKqSFJSUvDw8CA5OblAB/SPHj1KYGCg7CO+jtTGsYpSlx0nE1n7+QTGGb7BrK+C4YXNcGuDIi2zotSmuBWlLvn9XsvXlsuuXbvy9ab5vXOyEEL8W0pmDh99v5x5+kUAGLq9U+TGIpwrX81lw4YNJZ2HEKISe3v5NsZemYFJZyanYTdc2g9xdkqiiOSqIiGEU63YdYa2B6fSUPcP2VW9cXn0Y7kpZQVQqAP627dvZ+nSpZw6dYrs7Gy7ecuXLy+WxIQQFd/JS+n8sfILZhk2oNAwPvEZVLvV2WmJYlDgLZfFixfTuXNnYmJiWLFiBTk5ORw4cIDIyEg8PDxKIkchRAWUY7Hy9nfreYtPAVCdX4b69zo3KVFsCtxc3nvvPWbOnMnPP/+M0Whk9uzZHDp0iN69e1O3rjwVTgiRPx/+foihF8Px1NLJ9m6Drus4Z6ckilGBm0tcXBw9evQArl5hn56ejqZpvPLKK7YLGYUQ4kaij11C/fkBQbpDmA1VMfb+EvRycXJFUuDmcsstt5CamgrAbbfdZntIV1JSEhkZGcWbnRCiwknOyOGLRUt4WX/1ybWGhz+Q044roAIf0L/nnntYv349LVu25Mknn+Tll18mMjKS9evX07Vr15LIUQhRQSilmPTDFsZnfYBBZyWn2eO4tO7r7LRECShwc/noo49sjwseO3YsLi4ubNmyhccff5xx42SfqRAib8u2nabL0fepo79AllttTD1nymnHFVSBm4uXl5ft7zqdjlGjRhVrQkKIiunYhTR2rv6E9/VbsKLH1GcBuMoZphVVgY+5rFmzhl9//TXX9N9++421a9cWS1JCiIol22zl/e/WME778uqEe0dBnQ7OTUqUqAI3l1GjRjl8hrLVapWtGCGEQzN/3c/wxHDctEyybuuI7p5XnZ2SKGEFbi5Hjx6lWbNmuaY3adKE2NjYYklKCFFxbI69iHvUVFrrjpHt4oGp9xegkzsUV3QFbi4eHh4cO3Ys1/TY2FiqVatWLEkJISqGy+nZfL/4a57TrwbA+OiH4FHbyVmJ0lDg5tKrVy9GjhxJXFycbVpsbCyvvvoqPXv2LNbkhBDll1KKt5f+yYSc2eg0hbnN09Csl7PTEqWkwM1l6tSpVKtWjSZNmhAQEEBAQABNmzbl1ltvZfr06SWRoxCiHPo++iQPHXsXby2JTM+GGLq/7+yURCkq8KnIHh4ebNmyhfXr17Nnzx6qVKlCq1atuOeee0oiPyFEORR7PpXYNbPpr9+BRXPBtc8CMMpu88qkULfc1zSNBx98kAcffLC48xFClHNXsi1M+2Yls7VvANAemAS+rZyclSht+Wouc+bMYdiwYbi6ujJnzpwbxr700kvFkpgQonx6d+V2XksOx1WXQ5b//Zg6Pu/slIQT5Ku5zJw5k/79++Pq6srMmTPzjNM0TZqLEJXYj9tP03bfZAL1Z8muUgvTE5+CTh54Wxnlq7kcP37c4d+FEOKaowmp7P5pDm/r/8KKDmOfheBWy9lpCScp8H8pJk+e7PDW+leuXGHy5MnFkpQQonzJyDYz4+sfGKctuDrh/vHgf6dzkxJOVeDmMmnSJNLS0nJNz8jIYNKkScWSlBCi/FBK8fYP0byZGo5JyyE7IBjdXSOdnZZwsgI3F6UUmoNbZO/Zs8fujslCiMph6bZT3BkziQBdAlnV/DA++bkcZxH5PxX5lltuQdM0NE2jUaNGdg3GYrGQlpbGf//73xJJUghRNsWcS+Hwzx8wQR+NRTNg6vcNVJX/ZIoCNJdZs2ahlOKZZ55h0qRJeHj8/3MYjEYj/v7+dOrUqUSSFEKUPWlZZj78ejGzdP+7nuXBt6F2eydnJcqKfDeXgQMHAhAQEEDnzp1xcXEpsaSEEGWbUop3lm1mTMYUjJqF7MAeGOV6FvEvBb5Cv0uXLlgsFn744QdiYmIAaNasGb169cJgKNQF/0KIcub76BMEH55Ibf1FMqvXxfXxT+RxxcJOgbvBgQMH6NmzJ/Hx8TRu3BiAKVOmULNmTX7++WdatGhR7EkKIcqO/WeTOfPLNPrrd2HWGXF96lt5XLHIpcCndDz77LM0b96cM2fOsHPnTnbu3Mnp06dp1aoVw4YNK4kchRBlREpmDvO++ZZXdYsA0D00FXxbOzkrURYVeMtl9+7dbN++nVtuucU27ZZbbuHdd9/ljjvuKNbkhBBlh1KKd5ZsYtyVaRg0K9nNnsDYfpCz0xJlVIG3XBo1akRCQkKu6efPn6dhw4bFkpQQouz5ZssxHo6diI92mUyPBhh7zZbjLCJPBW4u4eHhvPTSS/zwww+cOXOGM2fO8MMPPzBy5EimTJlCSkqK7UcIUTHsOZ3E5XXvcY9+Hzk6V1z7fwcmN2enJcqwAu8We/jhhwHo3bu37UJKpRQAjzzyiO21pmlYLJbiylMI4STJGTl8+c1CZup+BMDQcxbUaurcpESZV+DmsmHDhpLIQwhRBlmtismLIhiXNQOdpshu9R+Mbfo5Oy1RDhTqOhchROXwccQhep+cSE1dCle8mlLlkenOTkmUE4W6u1xSUhIzZszg2Wef5dlnn2XmzJkkJycXOom5c+fi7++Pq6srQUFB/P333zeMX7ZsGU2aNMHV1ZWWLVuyZs0au/lKKSZMmICvry9VqlQhODiYo0eP2sUkJibSv39/3N3d8fT0ZMiQIXZ3e964cSO9evXC19eXatWq0aZNG7777rtCr6MQ5U1ETALum8YTpDtEtsGNKk99Cy5VnJ2WKCcK3Fy2b99OgwYNmDlzJomJiSQmJvLBBx/QoEEDdu7cWeAElixZQlhYGBMnTmTnzp20bt2akJAQzp8/7zB+y5Yt9OvXjyFDhrBr1y5CQ0MJDQ1l//79tpipU6cyZ84c5s2bR3R0NNWqVSMkJITMzExbTP/+/Tlw4ADr169n9erVbNq0ye46nS1bttCqVSt+/PFH9u7dy+DBg3n66adZvXp1gddRiPLm2IU0Ni7+gKf167GiYXzyC6ghZ4OKAlAFdNddd6lBgwapnJwc27ScnBw1cOBAdffddxd0capDhw5q+PDhttcWi0X5+fmp8PBwh/G9e/dWPXr0sJsWFBSknnvuOaWUUlarVfn4+Khp06bZ5iclJSmTyaQWLVqklFLq4MGDClDbtm2zxaxdu1ZpmqbOnj2bZ67du3dXgwcPzve6JScnK0AlJyfne4zZbFYxMTHKbDbne0xlIbVxrLjrkpqZo4ZP+VRlTvBSaqK7yol8v1iW6wzymXGsKHXJ7/daobZc3nzzTbv7iBkMBt544w22b99eoGVlZ2ezY8cOgoODbdN0Oh3BwcFERUU5HBMVFWUXDxASEmKLP378OPHx8XYxHh4eBAUF2WKioqLw9PSkffv/v4NrcHAwOp2O6OjoPPNNTk6WZ9aICs1qVUz6PpJx6e9h0sxkNuyOocvrzk5LlEMFPqDv7u7OqVOnaNKkid3006dPU7169QIt6+LFi1gsFry9ve2me3t7c+jQIYdj4uPjHcbHx8fb5l+bdqOYWrXsn+1tMBjw8vKyxVxv6dKlbNu2jU8//TTP9cnKyiIrK8v2+tq1PhaLJd+nZVssFqxWq5zG7YDUxrHirMsnEYfofXwsPrrLXPEIxPjYJ1iUgnJac/nMOFaUuuR3TIGbS58+fRgyZAjTp0+nc+fOAGzevJnXX3+dfv0q5imKGzZsYPDgwXz++ec0b948z7jw8HCHj3qOi4vDzS1/F5xZrVYSExOJjY1FJ0/zsyO1cay46vL3mXQ8N73LHYYjZOqr8c9d75FzKh5w/B+u8kA+M44VpS6OHnPvSIGby/Tp09E0jaeffhqz2QyAi4sLzz//PO+//36BllWjRg30en2u28kkJCTg4+PjcIyPj88N46/9mZCQgK+vr11MmzZtbDHXnzBgNptJTEzM9b5//PEHjzzyCDNnzuTpp5++4fqMHj2asLAw2+uUlBTq1KlDgwYNcHd3v+HYaywWC7GxsTRs2BC9Xp+vMZWF1Max4qjLiUvpHFv8Nm8ZIrCi4fLkl/gHdi3mTEuffGYcK0pd8nv3lQI3F6PRyOzZswkPDycuLg6ABg0aULVq1YIuCqPRSLt27YiIiCA0NBS42lEjIiIYMWKEwzGdOnUiIiKCkSNH2qatX7/e9hTMgIAAfHx8iIiIsDWTlJQUoqOjef75523LSEpKYseOHbRr1w6AyMhIrFYrQUFBtuVu3LiRhx9+mClTpuTrjs8mkwmTyZRrul6vL9AvUKfTFXhMZSG1cawodUnLMvPhV98zVX0BGljvHYuhSbcSyNI55DPjWGHrku/4wp5tUFwWL16sTCaTWrhwoTp48KAaNmyY8vT0VPHx8UoppQYMGKBGjRpli9+8ebMyGAxq+vTpKiYmRk2cOFG5uLioffv22WLef/995enpqVatWqX27t2revXqpQICAtSVK1dsMd26dVNt27ZV0dHR6q+//lKBgYGqX79+tvmRkZGqatWqavTo0ercuXO2n0uXLuV73eRsseIltXGsKHWxWq3qzQXrVMKEukpNdFdXvu2vlNVaAlk6h3xmHCuNs8Wc3lyUUurDDz9UdevWVUajUXXo0EFt3brVNq9Lly5q4MCBdvFLly5VjRo1UkajUTVv3lz98ssvdvOtVqsaP3688vb2ViaTSXXt2lUdPnzYLubSpUuqX79+ys3NTbm7u6vBgwer1NRU2/yBAwcqINdPly5d8r1e0lyKl9TGsaLU5ZOIA2rH+NuVmuiuMmbeoVRm6s0HlSPymXGsNJqLptT/7jopil1KSgoeHh4kJycX6JjL0aNHCQwMlM3460htHCtsXTYeSiD+u+foq99AtqE6xhc2gVf9Esy09MlnxrGi1CW/32ty+oQQldDJS+n8uXgqffUbsKLD2HdhhWsswrmkuQhRyaRnmZn95deMUgsAsN4/ARoG32SUEAUjzUWISkQpxXuLIhidFo6LZiGzcS8Md490dlqiApLmIkQl8tGv+3jy2Chqaslk3NIE18c/kUcVixIhzUWISuKn3WdosPlV2uiOkeXiQdUBi8FYzdlpiQpKmosQlcCuU5c5t3wM3fV/Y9EMmPovAq8AZ6clKjBpLkJUcP8kXeHnhVN5TrcKAK3XR+B/p5OzEhWdNBchKrD0LDMfzp/PaMvVu3ln3/kaujYV8wazomyR5iJEBWW1Kt7/ZhWjU9/DRbOQ0ehRjMHjnJ2WqCSkuQhRQc1ZHcXQU2/irmWQVqsdVZ+cJ2eGiVIjzUWICujH6Fju2v4ydXUXSK9aB7eBS8HF1dlpiUpEmosQFUx03AVMv4ygve4ImfrqVBu8HKrVcHZaopKR5iJEBXLyUjp7v3mTh3VRWNBjfOpbqNnI2WmJSkiaixAVRGpmDos/n8JQfgTA0mMWugb3OjcpUWlJcxGiArBYFXMXfs0rVz4CIP2OFzHecePHcgtRkqS5CFEBLN+8h+Hn38KoWUgO6EG1hyY7OyVRyUlzEaKc+27DLv5zegKeWjpJXq3xeOoL0Mk/beFc8gkUohz7ZftRWm8aSoAugRSTL57P/AAuVZydlhDSXIQor7YcPovnT4NorTtGms6dqoOXg1stZ6clBCDNRYhyaf/pRJK/f5Y7dfvJ1Fw5f/9MtJqNnZ2WEDYGZycghCiY05fSOfjl8/TWtmDGAH2+JZs6zk5LCDuy5SJEOXIpLYvfP32V3modVjSye32CS+D9zk5LiFykuQhRTmRkm1k27y0GZy8CIO3+96jatreTsxLCMWkuQpQDORYrCz+bybDUTwBIvCMM93tecHJWQuRNmosQZZxSii+/XsCzF8LRaYoLTf6DV/cJzk5LiBuS5iJEGfft8hX858RojJqFhDoPUbP3HHkuiyjzpLkIUYat/H0j3fe+RDUti4RbO+I98CvQ6Z2dlhA3Jc1FiDIq4u/d3PHnEG7VUklwa4b3sB/AYHJ2WkLkizQXIcqgjTv24f9LP27TLnLBVIda//0JTNWdnZYQ+SbNRYgyZtOug9T+qQ8NtH+4bKiF13O/oLnVdHZaQhSINBchypDNe2LwXtmbhtpZLhtqUP25dei96jk7LSEKTJqLEGVE1P4j3Lq8N4210yTpb6X6sHUYajZwdlpCFIo0FyHKgL8PxOKx7EmaaKdI1ntRbdhaDLUCnZ2WEIUmzUUIJ9t26DhVlz5JM+0EKTpPqgxdg4u33OFYlG/SXIRwop1HTmBa9AQttGOk6jxwHfoLRp+mzk5LiCKT5iKEk+w6egrdd0/QSoslVeeO8ZnVGH1bODstIYqFNBchnGBv3BnUt0/QRjtKmuaGcfDPmGq3cnZaQhQbpzeXuXPn4u/vj6urK0FBQfz99983jF+2bBlNmjTB1dWVli1bsmbNGrv5SikmTJiAr68vVapUITg4mKNHj9rFJCYm0r9/f9zd3fH09GTIkCGkpaXZ5mdmZjJo0CBatmyJwWAgNDS02NZXiL3HzpLzzRPcrh0mTXPDMPhnTHXaODstIYqVU5vLkiVLCAsLY+LEiezcuZPWrVsTEhLC+fPnHcZv2bKFfv36MWTIEHbt2kVoaCihoaHs37/fFjN16lTmzJnDvHnziI6Oplq1aoSEhJCZmWmL6d+/PwcOHGD9+vWsXr2aTZs2MWzYMNt8i8VClSpVeOmllwgODi65AohKJ2p/HNavetGOGNK1auifXoFr3dudnZYQxU85UYcOHdTw4cNtry0Wi/Lz81Ph4eEO43v37q169OhhNy0oKEg999xzSimlrFar8vHxUdOmTbPNT0pKUiaTSS1atEgppdTBgwcVoLZt22aLWbt2rdI0TZ09ezbXew4cOFD16tWrUOuXnJysAJWcnJzvMWazWcXExCiz2Vyo96zIynttIrbtUTETWig10V2lTrpNXTm+tViWW97rUpKkNo4VpS75/V5z2pZLdnY2O3bssNsy0Ol0BAcHExUV5XBMVFRUri2JkJAQW/zx48eJj4+3i/Hw8CAoKMgWExUVhaenJ+3bt7fFBAcHo9PpiI6OLrb1E+LfVm+Kpv5PT9iuYzE+uxZX/yBnpyVEiTE4640vXryIxWLB29vbbrq3tzeHDh1yOCY+Pt5hfHx8vG3+tWk3iqlVq5bdfIPBgJeXly2msLKyssjKyrK9TklJAa7uZrNYLPlahsViwWq15ju+Mimvtfnhtw3cvXUofrpEEl18cBu6Gv2t9YttPcprXUqD1MaxotQlv2Oc1lwqovDwcCZNmpRrelxcHG5ubvlahtVqJTExkdjYWHQ6p59vUaaUt9oopVgftY2nTo6lhpZCgksdUrrN5XyiBRKP3nwB+VTe6lKapDaOFaUu/z756Uac1lxq1KiBXq8nISHBbnpCQgI+Pj4Ox/j4+Nww/tqfCQkJ+Pr62sW0adPGFnP9CQNms5nExMQ83ze/Ro8eTVhYmO11SkoKderUoUGDBri7u+drGRaLhdjYWBo2bIheLw+F+rfyVBurVfHl4sUMOvkm7loG592aUGPYT9SoVqPY36s81aW0SW0cK0pdru2RuRmnNRej0Ui7du2IiIiwneprtVqJiIhgxIgRDsd06tSJiIgIRo4caZu2fv16OnXqBEBAQAA+Pj5ERETYmklKSgrR0dE8//zztmUkJSWxY8cO2rVrB0BkZCRWq5WgoKLtAzeZTJhMuR/mpNfrC/QL1Ol0BR5TWZSH2uRYrHz+1XwGnx5HFS2b87fcTq3nVoKrR4m9Z3moi7NIbRwrbF3yG+/U3WJhYWEMHDiQ9u3b06FDB2bNmkV6ejqDBw8G4Omnn+a2224jPDwcgJdffpkuXbowY8YMevToweLFi9m+fTufffYZAJqmMXLkSN555x0CAwMJCAhg/Pjx+Pn52RpY06ZN6datG0OHDmXevHnk5OQwYsQI+vbti5+fny23gwcPkp2dTWJiIqmpqezevRvA1rSEcORKtoUvP5vJ0AvvXX3mvffdeA9ZCsaqzk5NiNJVyDPZis2HH36o6tatq4xGo+rQoYPauvX/T8/s0qWLGjhwoF380qVLVaNGjZTRaFTNmzdXv/zyi918q9Wqxo8fr7y9vZXJZFJdu3ZVhw8ftou5dOmS6tevn3Jzc1Pu7u5q8ODBKjU11S6mXr16Csj1UxByKnLxKuu1OZ+SqT6eOkqZJ3goNdFdJXzRT6mcrBJ/37JeF2eS2jhWGqcia0op5cTeVqGlpKTg4eFBcnJygY65HD16lMDAQNmMv05Zrk1sfBLb579IX/NPAFxo1JeafT8GXcnnWZbr4mxSG8eKUpf8fq/J2WJCFFH04TOkLhpMX67euiix4yhqhowCTXNyZkI4jzQXIYpgTdQubls7hCBdHDkYyHp4Ll7t+zo7LSGcTpqLEIWglOL7n3+ly47h1NZdJE3njst/FuNW/05npyZEmSDNRYgCyrFYWfDNQvoeH4u7lkGiqQ6ez65EV7Ohs1MTosyQ5iJEAaRk5rD40/cYnDgbF81Cwi1t8R66HKp6OTs1IcoUaS5C5NOxhCS2zR/JsJwVoEFC3YfxHvAFuLg6OzUhyhxpLkLkw197D2FY/ix92AfA+TYv4t1zMsj9qoRwSJqLEDeglOKHn3+m845XuE27yBXNlezus6l1h5wRJsSNSHMRIg9Xsi38+MX7PBk/E5OWw0VjbTwGLaGKXwtnpyZEmSfNRQgHzly4zJ7Pn+c/2WtBgzM1u1D7ma+hiqezUxOiXJDmIsR1duw7gMuPg+jBEaxo/NNmJLV7TpDjK0IUgDQXIf7HalX8svJ7Ou4ZTU0tmTStGpk9P6V220ecnZoQ5Y40FyGA84nJbPsyjEfSfgANzpoacOszS6nhLRdGClEY0lxEpRe9LQqPX56nB8cBiK3Xlwb9Z6LJM1iEKDRpLqLSysoxs/7baXQ98QFVtGySNXfSu82mYdBjzk5NiHJPmouolE6cPs3Zb4bxcPYW0OBY9fb4Df4KD6/azk5NiApBmouoVKxWReQvi2i+fRx3apfIQc+JVq8SGDpazgYTohhJcxGVxonTpznx3UiCM38HDf7R34axzwICGwU5OzUhKhxpLqLCM1usRCyfz+373+FeLRmr0jhUty+N+09D71rd2ekJUSFJcxEV2uHYo1xc+hIh/zu2ctZQB8Ojc2nWvIuzUxOiQpPmIiqk9MxsNi6dw51xH9BYS8eMjqOBz9Kk92Q0lyrOTk+ICk+ai6hQlFJs2vgbt2waTw91GDQ4ZQrErfc8mjZo7+z0hKg0pLmICiMmNo6zP47h/oxf0WmKDFw50+pFGvV6E/Quzk5PiEpFmoso9xJT0tm8eApdzn5OUy0DNDhU6yH8+06nkVy3IoRTSHMR5daVLDMbf1pA4IFZPMIZ0OC0KZAqvWbQpJkcsBfCmaS5iHLHbLGy4dcV+Gx7n4fUUQCSNXcSO75JwAPPg07v5AyFENJcRLmhlCJm/05Slz1HsGU3AFcwcSrwaQJDx+JR7RbnJiiEsJHmIso8q1Wx9c/f4K8PeDRnKwA5GDha+3EaPPEWjT39nJyhEOJ60lxEmZVjtrD19+VU2zaHzpa9AFiVRkzNEOo98Q7NfAKdnKEQIi/SXESZk3Ylk+1rv8F73zzuVrEA5KDncM1uWFr1p0Xnbuj1clxFiLJMmosoM06cPEHcurk0/+dH7tUuAZCJkSO3PU79nm/StEZdjh496uQshRD5Ic1FOFWO2cLOzb9ijp5P+/Q/8NfMoMFlzYPTAb1p1PN1Wnl6A2CxWJycrRAiv6S5CKeIjT3MqQ1f0ODsTwRx7upEDeJMTchp9yyN7h3ALUZX5yYphCg0aS6i1PwT/w9HNy7CPe5nWmfvpqGmAMjAldgaXal1/wgaNOvs5CyFEMVBmosoUWdOn+L4lh+oFreallm78dP+t2tLgyOurTG36kej+/rTqoq7cxMVQhQraS6iWOWYzRzatZnLe1ZT85+NNLYcpfb/tlDQ4KQhgMSAHvh3GUCj2k2cm6wQosRIcxFFYjabORazi/P7IjCd+YvA9F201NL+P0CD4y4NuFzvIerd/RT16jWnnvPSFUKUEmkuokAuJPzD6QNRZJzYRrXzO6h/5QCNtHQaXQvQIB1Xjrm1wxIYQv1OjxJQqy4BzkxaCFHqdM5OAGDu3Ln4+/vj6upKUFAQf//99w3jly1bRpMmTXB1daVly5asWbPGbr5SigkTJuDr60uVKlUIDg7OdX1EYmIi/fv3x93dHU9PT4YMGUJaWppdzN69e7n77rtxdXWlTp06TJ06tXhWuBxITbnM4V1/Er1qHps/fYld74fwz1sNqflJU27f9Ax3nfqEtpl/46GlcwUjMabWbA94geM9V1Bl7GlavraGNr1exr1WXWevihDCCZy+5bJkyRLCwsKYN28eQUFBzJo1i5CQEA4fPkytWrVyxW/ZsoV+/foRHh7Oww8/zPfff09oaCg7d+6kRYsWAEydOpU5c+bw1VdfERAQwPjx4wkJCeHgwYO4ul49vbV///6cO3eO9evXk5OTw+DBgxk2bBjff/89ACkpKTz44IMEBwczb9489u3bxzPPPIOnpyfDhg0rvQKVgOysLC5fPEvKhbOkXTxL1uWzWJNO45J6BrcrZ6mRc46aXKZxHuPP6Py44NYUi9/t1Gh2D3WbdqSpi7FU10EIUbZpSinlzASCgoK44447+OijjwCwWq3UqVOHF198kVGjRuWK79OnD+np6axevdo2rWPHjrRp04Z58+ahlMLPz49XX32V1157DYDk5GS8vb1ZuHAhffv2JSYmhmbNmrFt2zbat7/66Nt169bRvXt3zpw5g5+fH5988gljx44lPj4eo/HqF+eoUaNYuXIlhw4dyte6paSk4OHhQXJyMu7u+Tsb6tTRfRzZ9ju33HILOp0GXN0SQymUApQFpawoiwVlNYMlB2XJQZmzUOYsMGeh5WSA+Qq6nAwM5nSM5lRM5nSqWlOobk3FXcvIVy6JuJNgrEta9QZo3s1wr9ea2k07UNX91nyNL24Wi4WjR48SGBgot3/5F6lL3qQ2jhWlLvn9XnPqlkt2djY7duxg9OjRtmk6nY7g4GCioqIcjomKiiIsLMxuWkhICCtXrgTg+PHjxMfHExwcbJvv4eFBUFAQUVFR9O3bl6ioKDw9PW2NBSA4OBidTkd0dDSPPvooUVFR3HPPPbbGcu19pkyZwuXLl7nllty3d8/KyiIrK8v2OiUlBbj6i8zv1eVntv1E8JES3P12tV9hVjqSNHeS9LeSbqxJZjVfrO61MXr5U92vId51G+NxSy08HCzCWVfKWywWrFarXKl/HalL3qQ2jhWlLvkd49TmcvHiRSwWC97e3nbTvb2989w6iI+PdxgfHx9vm39t2o1irt/lZjAY8PLysosJCAjItYxr8xw1l/DwcCZNmpRrelxcHG5ubg7X53rpypX9usZo2v+2Wq51g//9adV0gIYVHVZNj1UzYNX0WDQXrDoXLDoTVr0Jq8EVpa+CMlYFY3U0kxv6Kh4Yq92Cq9stuFb3RPe/h2qZ/vdzjRU4dzEZLibnK+fSYrVaSUxMJDY2Fp2uTBwuLBOkLnmT2jhWlLpcf2w6L04/5lKRjB492m6rKiUlhTp16tCgQYN87xarX/9FYmMfomHDhrIZfx2LxUJsbKzU5jpSl7xJbRwrSl2u7ZG5Gac2lxo1aqDX60lISLCbnpCQgI+Pj8MxPj4+N4y/9mdCQgK+vr52MW3atLHFnD9/3m4ZZrOZxMREu+U4ep9/v8f1TCYTJpMp13S9Xl+gX6BOpyvwmMpCauOY1CVvUhvHCluX/MY7dTvRaDTSrl07IiIibNOsVisRERF06tTJ4ZhOnTrZxQOsX7/eFh8QEICPj49dTEpKCtHR0baYTp06kZSUxI4dO2wxkZGRWK1WgoKCbDGbNm0iJyfH7n0aN27scJeYEEKIf1FOtnjxYmUymdTChQvVwYMH1bBhw5Snp6eKj49XSik1YMAANWrUKFv85s2blcFgUNOnT1cxMTFq4sSJysXFRe3bt88W8/777ytPT0+1atUqtXfvXtWrVy8VEBCgrly5Yovp1q2batu2rYqOjlZ//fWXCgwMVP369bPNT0pKUt7e3mrAgAFq//79avHixapq1arq008/zfe6JScnK0AlJyfne4zZbFYxMTHKbDbne0xlIbVxTOqSN6mNY0WpS36/15zeXJRS6sMPP1R169ZVRqNRdejQQW3dutU2r0uXLmrgwIF28UuXLlWNGjVSRqNRNW/eXP3yyy92861Wqxo/frzy9vZWJpNJde3aVR0+fNgu5tKlS6pfv37Kzc1Nubu7q8GDB6vU1FS7mD179qi77rpLmUwmddttt6n333+/QOslzaV4SW0ck7rkTWrjWGk0F6df51KRFeY6FzkvP29SG8ekLnmT2jhWGte5yLl5Qgghip00FyGEEMVOmosQQohiJ81FCCFEsZPmIoQQothJcxFCCFHs5N5iJejaWd75vRcPXD1FMC0tjZSUFDl18jpSG8ekLnmT2jhWlLpc+z672VUs0lxKUGpqKgB16tRxciZCCFG8UlNT8fBw9FCOq+QiyhJktVr5559/qF69uu0W+jdz7U7Kp0+fzveFl5WF1MYxqUvepDaOFaUuSilSU1Px8/O74e36ZculBOl0OmrXrl2ose7u7vKPIQ9SG8ekLnmT2jhW2LrcaIvlGjmgL4QQothJcxFCCFHspLmUMSaTiYkTJzp86FhlJ7VxTOqSN6mNY6VRFzmgL4QQotjJlosQQohiJ81FCCFEsZPmIoQQothJcxFCCFHspLk4wdy5c/H398fV1ZWgoCD+/vvvG8YvW7aMJk2a4OrqSsuWLVmzZk0pZVr6ClKbhQsXomma3Y+rq2spZls6Nm3axCOPPIKfnx+aprFy5cqbjtm4cSO33347JpOJhg0bsnDhwhLP0xkKWpuNGzfm+sxomkZ8fHzpJFwKwsPDueOOO6hevTq1atUiNDSUw4cP33RccX/PSHMpZUuWLCEsLIyJEyeyc+dOWrduTUhICOfPn3cYv2XLFvr168eQIUPYtWsXoaGhhIaGsn///lLOvOQVtDZw9Qrjc+fO2X5OnjxZihmXjvT0dFq3bs3cuXPzFX/8+HF69OjBfffdx+7duxk5ciTPPvssv/76awlnWvoKWptrDh8+bPe5qVWrVgllWPr++OMPhg8fztatW1m/fj05OTk8+OCDpKen5zmmRL5nlChVHTp0UMOHD7e9tlgsys/PT4WHhzuM7927t+rRo4fdtKCgIPXcc8+VaJ7OUNDaLFiwQHl4eJRSdmUDoFasWHHDmDfeeEM1b97cblqfPn1USEhICWbmfPmpzYYNGxSgLl++XCo5lQXnz59XgPrjjz/yjCmJ7xnZcilF2dnZ7Nixg+DgYNs0nU5HcHAwUVFRDsdERUXZxQOEhITkGV9eFaY2AGlpadSrV486derQq1cvDhw4UBrplmmV5TNTFG3atMHX15cHHniAzZs3OzudEpWcnAyAl5dXnjEl8ZmR5lKKLl68iMViwdvb2266t7d3nvt84+PjCxRfXhWmNo0bN+bLL79k1apVfPvtt1itVjp37syZM2dKI+UyK6/PTEpKCleuXHFSVmWDr68v8+bN48cff+THH3+kTp063HvvvezcudPZqZUIq9XKyJEjufPOO2nRokWecSXxPSN3RRblVqdOnejUqZPtdefOnWnatCmffvopb7/9thMzE2VV48aNady4se11586diYuLY+bMmXzzzTdOzKxkDB8+nP379/PXX3+V+nvLlkspqlGjBnq9noSEBLvpCQkJ+Pj4OBzj4+NToPjyqjC1uZ6Liwtt27YlNja2JFIsN/L6zLi7u1OlShUnZVV2dejQoUJ+ZkaMGMHq1avZsGHDTR/9URLfM9JcSpHRaKRdu3ZERETYplmtViIiIuz+B/5vnTp1sosHWL9+fZ7x5VVhanM9i8XCvn378PX1Lak0y4XK8pkpLrt3765QnxmlFCNGjGDFihVERkYSEBBw0zEl8pkp9KkAolAWL16sTCaTWrhwoTp48KAaNmyY8vT0VPHx8UoppQYMGKBGjRpli9+8ebMyGAxq+vTpKiYmRk2cOFG5uLioffv2OWsVSkxBazNp0iT166+/qri4OLVjxw7Vt29f5erqqg4cOOCsVSgRqampateuXWrXrl0KUB988IHatWuXOnnypFJKqVGjRqkBAwbY4o8dO6aqVq2qXn/9dRUTE6Pmzp2r9Hq9WrdunbNWocQUtDYzZ85UK1euVEePHlX79u1TL7/8stLpdOr333931ioUu+eff155eHiojRs3qnPnztl+MjIybDGl8T0jzcUJPvzwQ1W3bl1lNBpVhw4d1NatW23zunTpogYOHGgXv3TpUtWoUSNlNBpV8+bN1S+//FLKGZeegtRm5MiRtlhvb2/VvXt3tXPnTidkXbKunT57/c+1WgwcOFB16dIl15g2bdooo9Go6tevrxYsWFDqeZeGgtZmypQpqkGDBsrV1VV5eXmpe++9V0VGRjon+RLiqB6A3WegNL5n5Jb7Qgghip0ccxFCCFHspLkIIYQodtJchBBCFDtpLkIIIYqdNBchhBDFTpqLEEKIYifNRQghRLGT5iJEOXHvvfcycuTIYou9/smNhw4domPHjri6utKmTZtC5ykESHMRotI6d+4cDz30kO31xIkTqVatGocPHyYiIoKFCxfi6enpvARFuSa33BeilGVnZ2M0Gp2dRq473sbFxdGjRw/q1avnpIxERSJbLkKUsHvvvZcRI0YwcuRIatSoQUhICPv37+ehhx7Czc0Nb29vBgwYwMWLF21j0tPTefrpp3Fzc8PX15cZM2bkWu7HH39MYGAgrq6ueHt788QTT9jNt1qtvPHGG3h5eeHj48Nbb71lN//fu8U0TWPHjh1MnjwZTdO49957GTx4MMnJyWiahqZpucYLcSPSXIQoBV999RVGo5HNmzfz/vvvc//999O2bVu2b9/OunXrSEhIoHfv3rb4119/nT/++INVq1bx22+/sXHjRrunJW7fvp2XXnqJyZMnc/jwYdatW8c999yT6z2rVatGdHQ0U6dOZfLkyaxfv95hfufOnaN58+a8+uqrnDt3jp9++olZs2bh7u7OuXPnOHfuHK+99lrJFEdUSLJbTIhSEBgYyNSpUwF45513aNu2Le+9955t/pdffkmdOnU4cuQIfn5+fPHFF3z77bd07doVuNoo/v3Ap1OnTlGtWjUefvhhqlevTr169Wjbtq3de7Zq1YqJEyfa3v+jjz4iIiKCBx54IFd+Pj4+GAwG3NzcbLvLPDw80DStwj2YTpQOaS5ClIJ27drZ/r5nzx42bNiAm5tbrri4uDiuXLlCdnY2QUFBtuleXl52j+d94IEHqFevHvXr16dbt25069aNRx99lKpVq9piWrVqZbdsX19fzp8/X5yrJUSeZLeYEKWgWrVqtr+npaXxyCOPsHv3brufo0eP5tq1lZfq1auzc+dOFi1ahK+vLxMmTKB169YkJSXZYlxcXOzGaJqG1WotlvUR4makuQhRym6//XYOHDiAv78/DRs2tPupVq0aDRo0wMXFhejoaNuYy5cvc+TIEbvlGAwGgoODmTp1Knv37uXEiRNERkYWW55GoxGLxVJsyxOVizQXIUrZ8OHDSUxMpF+/fmzbto24uDh+/fVXBg8ejMViwc3NjSFDhvD6668TGRnJ/v37GTRoEDrd//9zXb16NXPmzGH37t2cPHmSr7/+GqvVarfrrKj8/f1JS0sjIiKCixcvkpGRUWzLFhWfNBchSpmfnx+bN2/GYrHw4IMP0rJlS0aOHImnp6etgUybNo27776bRx55hODgYO666y674zaenp4sX76c+++/n6ZNmzJv3jwWLVpE8+bNiy3Pzp0789///pc+ffpQs2ZN2wkJQuSHPOZYCCFEsZMtFyGEEMVOmosQQohiJ81FCCFEsZPmIoQQothJcxFCCFHspLkIIYQodtJchBBCFDtpLkIIIYqdNBchhBDFTpqLEEKIYifNRQghRLGT5iKEEKLY/R8l+4Yw5Bs55QAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "zs = np.geomspace(0.001, 2, 100)\n", "tau_sis = od.strong_lensing_optical_depth(zs)\n", "tau_sis_haris = od.optical_depth_SIS_haris(zs)\n", "\n", "# plot \n", "plt.figure(figsize=(4,4))\n", "plt.plot(zs, tau_sis, label='SIS hemanta')\n", "plt.plot(zs, tau_sis_haris, label='SIS haris')\n", "plt.xlabel('redshift')\n", "plt.ylabel('optical depth')\n", "plt.legend()\n", "plt.grid(alpha=0.5)\n", "plt.show()\n" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY4AAAF2CAYAAAB5xX8IAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/OQEPoAAAACXBIWXMAAA9hAAAPYQGoP6dpAABiVElEQVR4nO3deViUVfvA8e/MsO+uDLhBSi654BKE+WYLRWUlZbm8lmamVvZL4y3fNJc0C9Ncy7LeSm1xyUoqM8pwKyUU19yBUNxYFGEEZJs5vz/IyQlQBoFhuT/XxaVznnOeuZ/bcW6e7TwapZRCCCGEqCCtrQMQQghRt0jhEEIIYRUpHEIIIawihUMIIYRVpHAIIYSwihQOIYQQVpHCIYQQwipSOIQQQlhFCocQQgirSOEQQghhlVpROBYvXoyfnx9OTk4EBwezY8eOq/Zfs2YNHTp0wMnJiS5durB+/XqL5Uoppk6dio+PD87OzoSGhpKQkGBevnnzZjQaTZk/O3furJZtFEKI+kJj67mqVq9ezbBhw1iyZAnBwcEsWLCANWvWcPToUZo3b16q//bt27ntttuIjIzkgQceYMWKFbz11lvs3r2bzp07A/DWW28RGRnJ8uXL8ff3Z8qUKfzxxx8cOnQIJycnCgsLyczMtFjvlClTiImJISkpCY1Gc824TSYTZ86cwd3dvUL9hRCitlNKcfHiRXx9fdFqr7JfoWwsKChIjR071vzaaDQqX19fFRkZWWb/gQMHqn79+lm0BQcHqzFjxiillDKZTEqv16s5c+aYl2dlZSlHR0e1cuXKMtdZWFiomjVrpmbMmFHhuE+ePKkA+ZEf+ZGfevdz8uTJq37/2WFDhYWF7Nq1i4kTJ5rbtFotoaGhxMbGljkmNjaWiIgIi7awsDCioqIASE5OJjU1ldDQUPNyT09PgoODiY2NZfDgwaXW+d1333H+/HlGjBhRbqwFBQUUFBSYX6u/dtSOHz+Oh4fHtTf2L0ajkT///JMbbrgBnU5X4XH1neSlfJKbskleylfZ3BgMBvz8/HB3d79qP5sWjnPnzmE0GvH29rZo9/b25siRI2WOSU1NLbN/amqqefnltvL6/NPHH39MWFgYLVu2LDfWyMhIpk+fXqo9IyODS5culTvun0wmE4WFhaSnp199V7CBkbyUT3JTNslL+Sqbm5ycHIBrHn63aeGoDU6dOsVPP/3El19+edV+EydOtNjTMRgMtGrVirZt21q9x5GYmEi7du3kt6QrSF7KJ7kpm+SlfJXNjcFgqFA/mxaOpk2botPpSEtLs2hPS0tDr9eXOUav11+1/+U/09LS8PHxsegTGBhYan1Lly6lSZMmPPTQQ1eN1dHREUdHx1LtOp3O6g+tVqut1Lj6TvJSPslN2SQv5atMbira16b7dw4ODvTs2ZOYmBhzm8lkIiYmhpCQkDLHhISEWPQH2LBhg7m/v78/er3eoo/BYCAuLq7UOpVSLF26lGHDhmFvb19VmyWEEPWazQ9VRUREMHz4cHr16kVQUBALFiwgNzfXfKJ62LBhtGjRgsjISADGjRtH3759mTt3Lv369WPVqlXEx8fz4YcfAiXH5saPH8/MmTMJCAgwX47r6+tLeHi4xXtv3LiR5ORknn766WrZNqPRSFFRUak2k8lEfn6+/JZ0hYaWF51Oh52dnVzKLeokmxeOQYMGkZGRwdSpU0lNTSUwMJDo6Gjzye2UlBSLkzu9e/dmxYoVTJ48mUmTJhEQEEBUVJT5Hg6ACRMmkJuby+jRo8nKyqJPnz5ER0fj5ORk8d4ff/wxvXv3pkOHDlW+XTk5OZw6dcp89dVlSimKi4s5ceKEfGlcoSHmxcXFBR8fHxwcHGwdihBWsfkNgHWVwWDA09OT7OzsUifHjUYjCQkJuLi40KxZM4svQqUUBQUFODo6NpgvyIpoSHlRSlFYWEhGRgZGo5GAgICrXvly+fMUEBDQIPbGKkryUr7K5uZq32tXsvkeR31UVFSEUopmzZrh7OxssexynXZycqr3X5DWaGh5cXZ2xt7enhMnTlBYWFhqb1iI2kwufq5GDeELUFSe3Hsg6ir55AohhLCKFA4hhKhHcguK+flQGvvOVnxGC2vJOQ5R4/z8/Bg/fjzjx4+3dShC1AsnM/PYdDSdXw6n83vSeQqNJoJauvDobdXzfrLHIcwyMjJ49tlnad26NY6Ojuj1esLCwti2bZu5j5+fHwsWLDC/3rdvHw899BDNmzfHyckJPz8/Bg0aRHp6ug22oGbcfvvtUvSETRlNivjjmbwVfYSw+Vv51+xNTP32IFuPZeBpzOR5963cZ7+n2t5f9jiE2YABAygsLGT58uXccMMNpKWlERMTw/nz58vsn5GRwV133cUDDzzATz/9hJeXF8ePH+e7774jNze3hqMXon4z5Bex9VgGGw+ns+loOhfy/r65uKX2PE83PUCYZgf67L1oihS5ppuBcdUSi+xx1AClFHmFxVf8GP/xuvp+KnqbTlZWFr/++itvvfUWd9xxB23atCEoKIiJEyeWO4/Xtm3byM7O5qOPPqJ79+74+/tzxx13MH/+fPz9/a/6fnl5eTz11FO4u7vTunVr853/l508eZKBAwfi5eVF48aN6d+/P8ePHzcvf/LJJwkPD+fNN9/E29sbLy8vZsyYQXFxMS+//DKNGzemZcuWLF261GK9//3vf7nxxhtxcXHhhhtuYMqUKRZ397/22msEBgby2Wef4efnh6enJ4MHD+bixYvm992yZQsLFy40PzXy+PHjGI1GRo4cib+/P87OzrRv356FCxdWKPdClOfPjBw++vVPhnz4Oz1mbOD5FXv4Zs9pLuQV0cUpnXdabmKP/k1+c/g/njR8gE/2HjQolG8Pcn1uqba4ZI+jBlwqMtJp6k82ee9DM8Jwcbj2P7Obmxtubm5ERUVxyy23lDmh4z/p9XqKi4tZu3Ytjz76qFWXH8+dO5fXX3+dSZMm8dVXX/Hcc89xyy230LVrV4qKiggLCyMkJIRff/0VOzs7Zs6cyb333sv+/fvNd1pv3LiRli1bsnXrVrZt28bIkSPNT4iMi4tj9erVjBkzhrvvvts8Zb67uzvLli3D19eXP/74g1GjRuHu7s6ECRPMsSUlJREVFcW6deu4cOECAwcOZNasWbzxxhssXLiQY8eO0blzZ2bMmAFAs2bNMJlMtGzZkjVr1tCkSRO2b9/O6NGj8fHxYeDAgRXOi2jYiowmdh7PZOPhdDYeSefPc1fuuSvubpzBcK8/6JH3Ky5Zx+Dc5WUaaNMbOj4IHR7A5O5LZkICTaopTikcAgA7OzuWLVvGqFGjWLJkCT169KBv374MHjyYrl27ljnmlltuYdKkSfz73//mmWeeISgoiDvvvJNhw4aVeh7KP91///0899xzQMlewPz589myZQtdu3Zl9erVmEwmPvroI3MxWrp0KV5eXmzevJl77rkHgMaNG7No0SK0Wi3t27dn9uzZ5OXlMWnSJKBkKvxZs2bx22+/mR/gNXnyZHMMfn5+vPTSS6xatcqicJhMJpYtW2Z+mM0TTzxBTEwMb7zxBp6enjg4OODi4mIxg7NOp7N4Xou/vz+xsbF8+eWXUjjEVV3ILWTzsXRiDqez5VgGF/OLzcscdIohvucY4LKbjllbsM8+Dnl/LdTagX/fv4pFP3C74lHbRmO1xiyFowY42+s4NCMMKDlslZ9fgJNTzUyt4Wxf8ekGBgwYQL9+/fj111/5/fff+fHHH5k9ezYfffQRTz75ZJlj3njjDSIiIti4cSNxcXEsWbKEN998k61bt9KlS5dy3+vKYqTRaNDr9WRkZAAlJ9wTExNLPYUsPz+fpKQk8+ubbrrJ4iY6b29viznLdDodTZo0sThRv3r1ahYtWkRSUhI5OTkUFxeXmlrhn09A8/HxqdDJ/sWLF/PJJ5+QkpLCpUuXKCwsLHMqf9GwKaVISM/hl8NpbDyczu6UC5iuOKLczEXHiFap3G+3g9bpG9FmnP17oZ0TtAstKRY3hoFzo5rfAKRw1AiNRmM+XKSUQmsqxsmhds6M6uTkxN13383dd9/NlClTePrpp5k2bVq5hQOgSZMmPPbYYzz22GO8+eabdO/enbfffpvly5eXO+af09hrNBpMJhNQMkFkz549+eKLL0qNa9as2VXXcbX1xsbGMnToUKZPn05YWBienp6sWrWKuXPnVji28qxatYqXXnqJuXPnEhISgru7O3PmzCEuLu6q40TDUFBs5Pc/M9l4OI2YI+mcumB5j0UXbyee9EmhrymWJqd+QXPiigtSHNxKikTHh0qKhqNbDUdfmhQOcVWdOnUyP8+9IhwcHGjbtu11XVXVo0cPVq9eTfPmza16uuK1bN++nTZt2vDqq6+a206cOGH1ehwcHDD+41DAtm3b6N27t/nwG2CxdyQansJiE1uOZRC15zSbjqaTV/j3Z8bBTssd/i4MaZJI0KXfcDn+Cxy54ul7zo2gfT/o9FDJ4Sj72jWXmRQOAcD58+d57LHHeOqpp+jatSvu7u7Ex8cze/Zs+vfvX+aYdevWsWrVKgYPHsyNN96IUorvv/+e9evXl7qayRpDhw5lzpw59O/fnxkzZtCyZUtOnDjBN998w4QJE676bPirCQgIICUlhVWrVnHzzTfzww8/sHbtWqvX4+fnR1xcHMePH8fNzY3GjRsTEBDAp59+yk8//YS/vz+fffYZO3fuvObVZaJ+UUqx68QF1u45zQ9/nCXriktmm7s7cv+NLgxwO0jHrM3YJcXAySv2PNz00PGBkj2LNreCrvZ+PdfeyESNcnNzIzg4mPnz55OUlERRURGtWrVi1KhR5pPN/9SpUydcXFz4z3/+w8mTJ3F0dCQgIICPPvqIJ554otKxuLi4sHXrVv773//yyCOPcPHiRVq0aMFdd911XXsgDz30EC+++CLPP/88BQUF9OvXjylTpvDaa69ZtZ6XXnqJ4cOH06lTJy5dukRycjJjxoxhz549DBo0CI1Gw5AhQ3juuef48ccfKx2vqDsS03P4du9povae5mTm38Wgmbsjgzu5MMh9Hy1SY9Ac3gymKx7u5tW6pFB0fAha3gx1ZOJLeR5HJV1t3vr8/HySk5Px9/cvNV12ycnx/AYzfXhFNcS8XO1zciV57kTZbJ2XjIsFfL/vDFF7T7P/VLa53dVBR79OTRjR7Bgd0tehSfgZTH9fKUXT9iWHoDo+CPquUA2fd3kehxBC1BIlEwimsnbPGX5LyDBfDaXTaugb0JQRfucIufgTdoe+gSNZfw/Ud4VO/UuKRbP2Nom9KknhEEKIqyg2mvgt8RxRe07z08E0LhX9fZI7sJUXj3fQcJ/aiuvhNbAl8e+B7r7QdSB0GwzNO9og8uojhUMIIf5BKcX+U9ms3XOadfvPcC6n0LzMr4kLA7t48pjLHpr9+T78+uvfA+1dSs5XdBsM/reBtn4eWpTCIYQQf0k5n0fU3tNE7TltMd1HY1cHHurSnMebJ9P27Bdo4tdB8eWT4Brw/xd0+3fJoahacJ9FdZPCIYRo0C7kFrLuj7NE7TnNrhMXzO2OdlruuUnP4/459MqORvfHV7A39e+BTW8s2bPoMhC8WtkgctuRwiGEaHBMJsWmo+ms3HGSzUfTKf7rLLdGA7e2bcqgTo7cbfwVp4NvQvT+vwc6N4Yuj5YUDN8e1XJFVF0ghUMI0WBczC9iTfwplsce58T5PHP7Tb4ePNq1KQ+77sfr2HLY8Auov06Ca+1LpvwI/De0uxvsHGwUfe0hhUMIUe8ln8tl+fbjrIk/Se5fU394ONkxqFdLnmiZRuuTX0HsWij4+34MWvQq2bPoPABcGtso8tpJCoeocdX5zPHbb7+dwMBAi8fbioZJKcVviedYuu04m46mc/lW57bNXBlzsyfhpl9w2DcB4o//PcizFXQdVFIwmgbYJO66QAqHMMvIyGDq1Kn88MMPpKWl0ahRI7p168bUqVO59dZbgdJf+vv27WPKlCn8/vvvGAwG9Ho9wcHBvPPOOzRv3vwq71Y9vvnmm1Kz24qGJa+wmG92n2bZ9uMkpueY2++4sSkvdDAQmLoCzZa1YPzrElsHN+gUXlIs2txaZ6b9sCUpHMKsLj9zvLCwEAcHBxo3lkMKDdWpC3l8FnuClTtSMPz1MCRXBx1DejRnTJO9NDs0G37e+/cA3x4QNKrkjm4HV9sEXUdJ4RDA388c37x5M3379gUwP3e8PFc+c9zOruSjdPm549dy+Znja9asoVGjRrz66qsMGzbMvPy///0va9eu5dSpU+j1eoYOHcrUqVPNexOvvfYaUVFRPP/887zxxhucOHECk8lU6lDVe++9x/z58zl58iSenp7861//4quvvqpsmkQto5RiR3Imy2NT+PlQqnkKkDZNXBgbaE//4h9x3P8FXPrrMludY8k5i6CnoUVP2wVex0nhqAlKQVHe338vzAetsWYu5bN3qdD71KZnjkPFng2emJjI119/zTfffFPmRG7x8fG88MILfPbZZ/Tu3ZvMzEx+vfIuX1Fn5RcZ+XbvKT7cdIqkzD/N7X3aNuI/N5wmMHU5mm0/A39VEs/WcPNT0H0YuFbXk7gbDikcNaEoD970BUADONfke086U6Hd8Nr0zHGo2LPBCwsL+fTTTy2eCnillJQUXF1deeCBB3B3d6dNmzZ07979mrkQtVeaIZ/Pfz/BirgUzueWnKNwstcypIsHz3rF0fzIZ/Dr34WEtndC0GgIuKfeTv9hC3IWSJgNGDCAM2fO8N1333HvvfeyefNmevTowbJly8od88Ybb5CamsqSJUu46aabWLJkCR06dOCPP/646ntd7ZnjUPJs8FtvvRW9Xo+bmxuTJ08mJSXFYh1t2rQpt2gA3H333bRp04YbbriBJ554gi+++IK8vLxy+4vaa0/KBV5YuYdbZ23knY2JnM8txMfTiUmdLrA38HumHRtA8+3TIfNPcPSE4Gfh+V3wxFpof58UjaqmbOzdd99Vbdq0UY6OjiooKEjFxcVdtf+XX36p2rdvrxwdHVXnzp3VDz/8YLHcZDKpKVOmKL1er5ycnNRdd92ljh07Vmo969atU0FBQcrJyUl5eXmp/v37WxV3dna2AlR2dnapZZcuXVKHDh1Sly5duhyUUgU5ShXkKFP+RZWXlaFM+RfNbdX6YzJZtV3/NHLkSNW6dWvz6zZt2qj58+eX27+goEB16tRJDRs2rNw+Za2jW7duatKkScpkMqnt27crnU6nZs6cqXbu3KmOHTumZsyYoTw9Pc39p02bprp161Zq3X379lXjxo0zvy4qKlIbNmxQL7/8srrhhhtUu3bt1IULF66x1TWj1OekHMXFxerw4cOquLi4hiKrHUwmk/rxj7Oq/7u/qTb/XWf+GfzeFrX7hw+V8X93KzXN4++fxSFK7fyk5HPfwFX2M3O177Ur2fRQ1erVq4mIiGDJkiUEBwezYMECwsLCOHr0aJmXcm7fvp0hQ4YQGRnJAw88wIoVKwgPD2f37t107twZgNmzZ7No0SKWL1+Ov78/U6ZMISwsjEOHDpkflvP1118zatQo3nzzTe68806Ki4s5cOBA9W2oRvP34SKlwKQDB6c6MV2BLZ45XlXPBoeSQ3ChoaGEhoYybdo0vLy82LhxI4888kil4xPVSynFrwnnePvno+YHJDnotAztZM+zbltpfmwl7Egv6avRoTo+hDZ4NLQOqRP/p+oDmxaOefPmMWrUKEaMGAHAkiVL+OGHH/jkk0945ZVXSvVfuHAh9957Ly+//DIAr7/+Ohs2bODdd99lyZIlKKVYsGABkydPNj8n+9NPP8Xb25uoqCgGDx5McXEx48aNY86cOYwcOdK87k6dOtXAFtdetemZ41X1bPB169bx559/ctttt9GoUSPWr1+PyWSiffu6/yCd+mrn8Uzm/HSUHcmZALg4aJnW+QLhxetxTFj/9zQgbnpMPYaT1Og2bugaAvJkxBpls8JRWFjIrl27mDhxorlNq9USGhpKbGxsmWNiY2OJiIiwaAsLCzP/RpycnExqaiqhoaHm5Z6engQHBxMbG8vgwYPZvXs3p0+fRqvV0r17d1JTUwkMDGTOnDnmvZayFBQUUFBQYH5tMBiAkkc0Go1Gi75GoxGllPnnSpdf/7Pd1lxdXQkKCir1zPGnn36aSZMmWcR7ebs6duyIs7NzqWeO/+9//+Pxxx+/6jaWlZvL7Q8++CDjx4+3eDb45MmTmT59eqn8lbcOpRSenp588803vPbaa+Tn5xMQEMCKFSvo1KlTrcj/5TjL+gxdyWg0YjKZrtqnrjt4xsC8DQlsPlZynsvJDl6/MYnwi6uxP/T30QDVujemm5+G9v0woqUoMbFe56WyKvuZqWh/mz1z/MyZM7Ro0YLt27cTEhJibp8wYQJbtmwhLi6u1BgHBweWL1/OkCFDzG3vvfce06dPJy0tje3bt3Prrbdy5swZfHx8zH0GDhyIRqNh9erVrFq1iiFDhtC6dWvmzZuHn58fc+fO5eeff+bYsWPl3kD22muvMX369FLtO3fuxM3Ncv59k8lEcXExrVu3LvOy1uLiYvN9D+JvDS0vBQUFpKSkYGdnh/YqdyubTCYyMzNp3LjxVfvVRSezCvl0bya/Hi85tOmoKeZVn108WvANLrknATDpnMj2u4+sgEcp8GpnHluf83K9KpubnJwcbr75Znnm+D+ZTCYAXn31VQYMGADA0qVLadmyJWvWrGHMmDFljps4caLF3o7BYKBVq1a0bdu2VILz8/M5ceIEjo6O5vMql12u046Ojlbd91DfNdS82NnZ0aZNm1KfkysZjUYSExNp165dmfer1EWnLuSxaGMSa/ecxqTAWVPA66120z/va+wzzwCgnLxQQaNRN4/Gw6Ux//waq495qSqVzc3lIynXYrPC0bRpU3Q6HWlpaRbtaWlp6PX6Msfo9fqr9r/8Z1pamsUeR1paGoGBgQDm9ivPaTg6OnLDDTeUutzzSo6OjmXuPeh0ulL/MDqdDo1GY/4py9WWNWQNKS+Xt7Wsz9A/abXaCvWr7dIv5rN4YyIrdqRQZFS4k8d0n1geuhSFXfpfU9u4eUPI82h6jUDj6H7V9dWXvFSHyuSmon1ttn/n4OBAz549iYmJMbeZTCZiYmIsDl1dKSQkxKI/wIYNG8z9/f390ev1Fn0MBgNxcXHmPj179sTR0ZGjR4+a+xQVFXH8+HHatGlTZdsnhPhbVl4hs348wm2zN7E89gTuxmwWNvueve4v8siFj7HLPw9eraHfPBi3H259Aa5RNITt2PRQVUREBMOHD6dXr14EBQWxYMECcnNzzVdZDRs2jBYtWhAZGQnAuHHj6Nu3L3PnzqVfv36sWrWK+Ph4PvzwQ6DkN7jx48czc+ZMAgICzJfj+vr6Eh4eDoCHhwfPPPMM06ZNo1WrVrRp04Y5c+YA8Nhjj9V8EoSox3IKivnkt2T+t/VPLhYU48N5Xm30C/cV/oTuYn5Jp2YdoE9EyRxSugZ39LxOsum/0qBBg8xTeV++uik6Oto8XUVKSorFiZ3evXuzYsUKJk+ezKRJkwgICCAqKsriaqgJEyaQm5vL6NGjycrKok+fPkRHR1scQ54zZw52dnY88cQTXLp0ieDgYDZu3EijRo1qbuOFqMfyi4x8/vsJ3tucRGZuIX6as8zxiOae4s1oLxWVdPLtDv96CdrfL1OZ1zE2u6qqrjMYDHh6epZ59UF+fj7Jycn4+fnh7Gw5M5VSivz8fJycnBrMsfyKaIh5uXTpEsePH8ff3/+aJ8cTEhIICAio9cfyi4wm1sSfYlFMAqmGfDpqTvCy6w/cUbwdDSUXpuD3L/jXf+CG26/rhr26lJeaVtncXO177UqyX1gNLv9DFRYWliocQlx2ed6s+vDgqWKjie/2nWFhTAInzufRQ3OMuS7fc6tpFxT/1enG++BfEdCq/Kn6Rd0ghaMa2NnZ4eLiQkZGBvb29haH25RS5hsJG8pv1hXRkPKilCIvL4/09HS8vLzq9G/LRpNi3f4zLPwlgT/P5fAv7R/Mc/qOnhwCE6DRwk2PQJ8XQV/+DbaibpHCUQ00Gg0+Pj4kJyeXmmNJKWW+0a2+f0FaoyHmxcvLq9xLz2s7k0nxwx9nWRiTQGJ6Dj00x/jaaTU9OVzSQWsPgf+GW8dBk7a2DVZUOSkc1cTBwYGAgAAKCwst2o1GIydOnKBNmzZ1+jfNqtbQ8mJvb18nt9NkUkQfTGXBL8c4lpbDDZozfOS0hlD+munBzgl6joDe/weeLWwbrKg2UjiqkVarLXXS02g0mtvr4hdHdZG81G5KKX46mMaCX45xJPUizbjAbKe1PMomtBhLDkkFDoU7JoGHr63DFdVMCocQolxKKX45nM6CX45x8IwBN/KY6Liep3TrsTf9dR9G+/vhrqnQvKNtgxU1RgqHEKIUpRSbjqaz4JcE9p/Kxp5iRjtsZLxDFC7FWSUnvlsGwd3ToU1vW4crapgUDiGEmVKKLccymP9LAvtOZqHBxACHHUx2/ppGBadLLq1tEgCh06DDA/LgpAZKCocQAqUUvyWeY/6GY+xOyQLgDvtDvOn+FT55R6CAkskHb38Fug+TqUEaOPnXF6KBi006z7wNR9l5/AIA3exSmNd4LW0NcZAHOLiXXFYb8tzfj0AWDZoUDiEaqNTsfF7/4RA/7D8LgL/dORY0/4GumT+jMaiSezFuHgm3vQyuTW0crahNpHAI0cAUG00s236c+RuOkVtopLHmIu+03EjvzLVoMv+676jzo3Dnq9D4BtsGK2olKRxCNCDxxzOZHHWAI6kX0WLilaaxjCr8DF3GX09+8+9bcqWUb3fbBipqNSkcQjQAmbmFzPrxMF/GnwLgFueTLPb4jCbZB0o6eHcuKRht75IrpcQ1SeEQoh4zmRRfxp9kVvQRsvKKcCePJb4/0vvCWjTZJnD0gDunlJzL0Mod+6JipHAIUU8dPJPN5KgD7EnJAhRjGu/lJbUM+8yMkg6dH4WwN8C9bk60KGxHCocQ9czF/CLmbTjG8u3HMSno5JDG/5quokXmXxMRNmkH978Nbe+wbaCizpLCIUQ9oZTi+/1nmbnuEOkXC3CkkAW+MdybvbrkaimdI9z2Usk9GXaOtg5X1GFSOISoB5Iycpj67QG2JZ4HYKDXUabbfYJz5smSDu1C4f45cnmtqBJSOISowy4VGlm8KZEPtiZRZFS0ssviE/3XBJyLKeng7gv3RkKn/nK1lKgyUjiEqKP2nszixdV7ST6Xiw4jkT6/MSjnc7TnckGjg1ueLZlbytHd1qGKekYKhxB1TLHRxHubk1gYk4DRpAh1O8E812V4XDha0qFlEDwwD/RdbBuoqLekcAhRh6Scz+PFL/ey68QFXMjnA+8o/pX9HWQDzo0gdDp0fwK0WluHKuoxKRxC1AFKKb7afYoZ3x8mt9BIH8cklrh9iFv2Xye/Ax8vufNbJiMUNUAKhxC13IW8Qt7YnMZvJ3Kxp5h5TdbzcN6XaHJN4NESwt+DG/raOkzRgEjhEKIW+y3hHP/5ci9pFwvoqDvFp40+plnOX+cyug6G+94CZy+bxigaHikcQtRC+UVG5vx0lI9/S0aDif+4/sxYVqHNKSw5l/HAArgp3NZhigZKCocQtcyRVAPjV+3lSOpFfDnHZ02W0TZ3d8nCdndD/3dlfilhU1I4hKglTCbF0u3HeSv6CIXFRoY5xzLVbhl2uTmYdE4Q9ibam5+SG/mEzUnhEKIWSDPk89KaffyacI5GGFjZ+HN65v0GRaBa9CI58BX8etwpRUPUCrXiYu/Fixfj5+eHk5MTwcHB7Nix46r916xZQ4cOHXBycqJLly6sX7/eYrlSiqlTp+Lj44OzszOhoaEkJCRY9PHz80Oj0Vj8zJo1q8q3TYhriT6QStiCrfyacI577Pex3ePVkqKhtYM7J2N6cj1F7q1sHaYQZjYvHKtXryYiIoJp06axe/duunXrRlhYGOnp6WX23759O0OGDGHkyJHs2bOH8PBwwsPDOXDggLnP7NmzWbRoEUuWLCEuLg5XV1fCwsLIz8+3WNeMGTM4e/as+ef//u//qnVbhbiSyaSY89MRnvl8F4V5F3nf8zM+1L2Fc+F5aNYBno6B214uKSBC1CI2Lxzz5s1j1KhRjBgxgk6dOrFkyRJcXFz45JNPyuy/cOFC7r33Xl5++WU6duzI66+/To8ePXj33XeBkr2NBQsWMHnyZPr370/Xrl359NNPOXPmDFFRURbrcnd3R6/Xm39cXV2re3OFAMCQX8SoT+NZvCmJ7poEfvOcwn0FP5YsvOU5GL0ZfANtGaIQ5bLprzKFhYXs2rWLiRMnmtu0Wi2hoaHExsaWOSY2NpaIiAiLtrCwMHNRSE5OJjU1ldDQUPNyT09PgoODiY2NZfDgweb2WbNm8frrr9O6dWv+/e9/8+KLL2JnV3ZKCgoKKCgoML82GAwAGI1GjEZjhbfZaDRiMpmsGtMQNKS8/JmRw5jP9/DnuRxG2f/ERLsVaAuKUR4tMD20GPxvK+n4Vy4aUm6sIXkpX2VzU9H+Ni0c586dw2g04u3tbdHu7e3NkSNHyhyTmppaZv/U1FTz8stt5fUBeOGFF+jRoweNGzdm+/btTJw4kbNnzzJv3rwy3zcyMpLp06eXak9KSsLNze0aW/o3k8lEZmYmiYmJaGU+IbOGkpcdp3KZtSUdivL4wOkjwtgOCgytQkm9+RVMxe7wj/NxDSU31pK8lK+yucnJyalQvwZ78PTKvZauXbvi4ODAmDFjiIyMxNGx9NPRJk6caDHGYDDQqlUr2rZti4eHR4Xf12g0kpiYSLt27dDpdNe3EfVIfc+LUoolW/9kbkwqfpxluesiWhtPoLR2qLtfx/Xm0bQt54qp+p6bypK8lK+yubl8JOVabFo4mjZtik6nIy0tzaI9LS0Nvb7sG5z0ev1V+1/+My0tDR8fH4s+gYGB5cYSHBxMcXExx48fp3379qWWOzo6lllQdDqd1R9arVZbqXH1XX3NS15hMRO+2s+6/WcJ0+5kodMHOBnzwM0bzWPL0bQJueY66mturpfkpXyVyU1F+9p0/87BwYGePXsSExNjbjOZTMTExBASUvZ/ppCQEIv+ABs2bDD39/f3R6/XW/QxGAzExcWVu06AvXv3otVqad68+fVskhAWTmbmMeD9WH7cf4qJ9qv4wGE+TqY8aN0bxmyFChQNIWobmx+qioiIYPjw4fTq1YugoCAWLFhAbm4uI0aMAGDYsGG0aNGCyMhIAMaNG0ffvn2ZO3cu/fr1Y9WqVcTHx/Phhx8CoNFoGD9+PDNnziQgIAB/f3+mTJmCr68v4eHhQMkJ9ri4OO644w7c3d2JjY3lxRdf5PHHH6dRo0Y2yYOof2KTzjN2xW40uRmscl7MzeqvS8ZvGVsyBbrO3rYBClFJNi8cgwYNIiMjg6lTp5KamkpgYCDR0dHmk9spKSkWJ3d69+7NihUrmDx5MpMmTSIgIICoqCg6d+5s7jNhwgRyc3MZPXo0WVlZ9OnTh+joaJycnICSw06rVq3itddeo6CgAH9/f1588cVSV2sJURlKKT6NPcGMdYfoqo7xP5dFNDWdB3vXknmmOj9i6xCFuC4apZSydRB1kcFgwNPTk+zsbKtPjickJBAQECDHZa9QX/JSUGxkStQBvow/yeO6X3jN/jPsKIYmATDoc2jewep11pfcVDXJS/kqm5uKfq/ZfI9DiPoi/WI+Yz7bxeGUNObaf8IA3a8lCzo+BP0Xg1PFf8EQojaTwiFEFTiSauCppTuxM5zgW6cFtOcEaHQQ+hr0/j+ZnFDUK1I4hLhOm4+m8/yKPQQW7WGJ4yLcyAXXZvDoUvD/l63DE6LKSeEQ4jp89vsJpn17gAHazUQ6fIwdRmgZBAOXg4evrcMTolpI4RCiEowmxZvrD/Pxb3/yot3XjLP7pmRBl4ElV07Zlb5ZVIj6QgqHEFbKLShm3Kq9bDl8mrftP+JR3daSBf96Ce6cLOczRL0nhUMIK6QZ8nlq2U5SzqSyzGEht2r/KDkJ3m8u9Bph6/CEqBFSOISooENnDIxcvhOVfZpvnOYQQErJTX2PLYMb77F1eELUGCkcQlTAxiNp/N+KPbQqSuYz5zk0U+fBzRv+/aU8cEk0OFI4hLiGZduSmbHuEL01f/Ch00JcVB40bQ+PfwVerW0dnhA1TgqHEOUwmhSvrzvEsu3HeVS3hbfsP0KnjNCmDwz+HJxlQkzRMEnhEKIMRpNi3Ko9rNt/hnG6b3jR/uuSBV0eK5k+RC63FQ2YFA4h/sFkUrz81T6i959krsPHDNBuKVnQJwLunALymFLRwEnhEOIKSikmf3uAdbuP84HDQu7S7gaN9q/LbZ+ydXhC1ApSOIT4i1KK19cd5uu4RD50mM/t2n1g5wQDP4Ubw2wdnhC1hhQOIf7y9s9HWbHtCB/Zz+Vf2gNg7wL/Xg3+t9k6NCFqFSkcQgDvbkxg6aaDLLV/mxDdoZIb+4auAb9bbR2aELWOFA7R4H3065+8//M+ljnMJkh7FBzc4fGvoXWwrUMTolaqVOGIiYkhJiaG9PR0TCaTxbJPPvmkSgIToiZ89vsJFv6wi08d3qKnNgEcPeGJb6BlL1uHJkStZXXhmD59OjNmzKBXr174+PigkZlARR21Jv4kc6Li+MxhFoHaJJSTF5on1kKLHrYOTYhazerCsWTJEpYtW8YTTzxRHfEIUSO+23eGN7/ezhcOb9JFexzl3BjNsCjw6Wbr0ISo9awuHIWFhfTu3bs6YhGiRvx0MJUZq7fyhf2bdNKeQLk0RTPsW9B3tnVoQtQJVt8C+/TTT7NixYrqiEWIarflWAbTVmzmM7uZJUXDtTmaJ3+QoiGEFSq0xxEREWH+u8lk4sMPP+SXX36ha9eu2NvbW/SdN29e1UYoRBXZdSKTyZ/F8KluBjdqT6Pc9GiGfw/NbrR1aELUKRUqHHv27LF4HRgYCMCBAweqPCAhqsORVAMvLN3MR5o3SoqGuy+aJ9dBk7a2Dk2IOqdChWPTpk3VHYcQ1SblfB6jPvqVeaa36KhNQbl6oxnxAzS+wdahCVEnWX2O46mnnuLixYul2nNzc3nqKZkETtQu6RfzefLj7UwtmEuw9gjK0QPNE99I0RDiOlhdOJYvX86lS5dKtV+6dIlPP/20SoISoipkXypi2EdxPGNYxN26XSidE5p/r5YT4UJcpwpfjmswGFBKoZTi4sWLODk5mZcZjUbWr19P8+bNqyVIIax1qdDIyGU7CT//PwbabUFpdGgeWwpt5FJyIa5XhQuHl5cXGo0GjUbDjTeWvgpFo9Ewffr0Kg1OiMooMpp47otddD/1Gc/Yfw+A5qFF0OF+G0cmRP1Q4cKxadMmlFLceeedfP311zRu3Ni8zMHBgTZt2uDr61stQQpRUSaT4qU1+2ic8DWvOvx1v1HodOj+uG0DE6IeqfA5jr59+3L77beTnJxM//796du3r/knJCTkuorG4sWL8fPzw8nJieDgYHbs2HHV/mvWrKFDhw44OTnRpUsX1q9fb7FcKcXUqVPx8fHB2dmZ0NBQEhISylxXQUEBgYGBaDQa9u7dW+ltELanlGL69we5uH8db9l/WNIY8jzcOs62gQlRz1h9crxNmzZkZ2fz9ttvM3LkSEaOHMncuXPJzMysVACrV68mIiKCadOmsXv3brp160ZYWBjp6ell9t++fTtDhgxh5MiR7Nmzh/DwcMLDwy3uKZk9ezaLFi1iyZIlxMXF4erqSlhYGPn5+aXWN2HCBNlTqicWxiRw8PefeM9+IXYaE3QbAne/DjIRpxBVS1lpy5YtysPDQ7Vq1Uo9/PDD6uGHH1atW7dWHh4easuWLdauTgUFBamxY8eaXxuNRuXr66siIyPL7D9w4EDVr18/i7bg4GA1ZswYpZRSJpNJ6fV6NWfOHPPyrKws5ejoqFauXGkxbv369apDhw7q4MGDClB79uypcNzZ2dkKUNnZ2RUeo5RSxcXF6vDhw6q4uNiqcfXd9eZlRdwJdc8r76msqXqlpnko9fljShUXVnGUtiGfmbJJXspX2dxU9HvN6kkOx44dy6BBg3j//ffR6XRAyVVVzz33HGPHjuWPP/6o8LoKCwvZtWsXEydONLdptVpCQ0OJjY0tc0xsbKzFFCgAYWFhREVFAZCcnExqaiqhoaHm5Z6engQHBxMbG8vgwYMBSEtLY9SoUURFReHi4nLNWAsKCigoKDC/NhgMQMm2G43Gim3wX/1NJpNVYxqC68nL9qTzvBe1ma8cZuGpyUO1DMI04GNAC/Ugz/KZKZvkpXyVzU1F+1tdOBITE/nqq6/MRQNAp9MRERFh9X0c586dw2g04u3tbdHu7e3NkSNHyhyTmppaZv/U1FTz8stt5fVRSvHkk0/yzDPP0KtXL44fP37NWCMjI8u8aiwpKQk3N7drjr/MZDKRmZlJYmIiWq3VRwrrrcrm5VR2IRN/SGK53dt4a7LI92xLStAbmI6frsZoa5Z8ZsomeSlfZXOTk5NToX5WF44ePXpw+PBh2rdvb9F++PBhunWrG88yeOedd7h48aLFns61TJw40WJPx2Aw0KpVK9q2bYuHh0eF12M0GklMTKRdu3YWxbehq0xesvIKefa77UxX79FJdwLl0gz74Wtp69mymqOtWfKZKZvkpXyVzc3lIynXYnXheOGFFxg3bhyJiYnccsstAPz+++8sXryYWbNmsX//fnPfrl27XnVdTZs2RafTkZaWZtGelpaGXq8vc4xer79q/8t/pqWl4ePjY9Hn8uSMGzduJDY2FkdHR4v19OrVi6FDh7J8+fJS7+vo6FiqP5TsbVn7odVqtZUaV99Zk5cio4nnV+6jv+Fz7rfbgdLaoxn8ObrGbWog0ponn5mySV7KV5ncVLSv1YVjyJAhQMnVSGUt02g0KKXQaDTXPF7m4OBAz549iYmJITw8HCjZxYqJieH5558vc0xISAgxMTGMHz/e3LZhwwZCQkIA8Pf3R6/XExMTYy4UBoOBuLg4nn32WQAWLVrEzJkzzePPnDlDWFgYq1evJjg4uEJ5ELajlGJK1AEaHf+B8Q7fAKB5cAG0vsW2gQnRQFhdOJKTk6s0gIiICIYPH06vXr0ICgpiwYIF5ObmMmLECACGDRtGixYtiIyMBGDcuHH07duXuXPn0q9fP1atWkV8fDwfflhy3b5Go2H8+PHMnDmTgIAA/P39mTJlCr6+vubi1Lp1a4sYLp+jaNu2LS1b1q/DHPXRx78l80f8r3zlsKSk4ZaxcoOfEDXI6sLRpk3VHgoYNGgQGRkZTJ06ldTUVAIDA4mOjjaf3E5JSbE4udO7d29WrFjB5MmTmTRpEgEBAURFRdG5898T102YMIHc3FxGjx5NVlYWffr0ITo62mJ+LVE3/XIojQ/W/06Uw1ycNYXQ9i64e4atwxKiQdEopZS1gz777DOWLFlCcnIysbGxtGnThgULFuDv70///v2rI85ax2Aw4OnpSXZ2ttUnxxMSEggICJDjsleoSF4OnzUw5P0tfMQMemmPoZq0Q/N0DDh71WywNUw+M2WTvJSvsrmp6Pea1dewvf/++0RERHD//feTlZVlPo/h5eXFggULrF2dEBWSfjGfp5ftZJLpfyVFw9EDzZBV9b5oCFEbWV043nnnHf73v//x6quvWlSyXr16WXXznxAVVVBsZPSnuwjLWfvXFOnakinSmwbYOjQhGiSrC0dycjLdu3cv1e7o6Ehubm6VBCXElWZHH8Xj9BZetf8cAM09M6Fd6DVGCSGqi9WFw9/fv8xZZKOjo+nYsWNVxCSE2eaj6fy87XfetX8HHQoCH4dbnrN1WEI0aFZfVRUREcHYsWPJz89HKcWOHTtYuXIlkZGRfPTRR9URo2igMi4W8MqX8Xxg/w4emjxoGQQPzJPZboWwMasLx9NPP42zszOTJ08mLy+Pf//73/j6+rJw4ULzBIJCXC+TSfGfNfsYVfAp3ez+RDk3KjmvYVf67n0hRM2yunAADB06lKFDh5KXl0dOTo48a1xUuU+2JeOQGM1Ihx8B0IS/D/VsDioh6qpKFY7LXFxcKjQluRDWOHA6m0+jf+M7+7/uDA95HtrfZ9ughBBmFSoc3bt3R1PB48q7d+++roBEw5ZbUMyLK3YyX7cIL00uyrcHmrum2TosIcQVKlQ4Ls/xBJCfn897771Hp06dzBML/v777xw8eJDnnpOrXcT1ef2HwzySvYyedgmYHD3QPrYU7BxsHZYQ4goVKhzTpv39G9/TTz/NCy+8wOuvv16qz8mTJ6s2OtGgbE3OIX3Pz8xx+B4Abf93oZGfbYMSQpRi9X0ca9asYdiwYaXaH3/8cb7++usqCUo0PKcu5LFy+1Hm2b9X0nDzKOjUMOY9E6KusbpwODs7s23btlLt27Ztk9lnRaUYTYqXvtzDLM0immguovRd4Z6Z1x4ohLAJq6+qGj9+PM8++yy7d+8mKCgIgLi4OD755BOmTJlS5QGK+m/ptmT6nFlKsN0RTPauaB9bBvbyS4gQtZXVheOVV17hhhtuYOHChXz+ecncQR07dmTp0qUMHDiwygMU9VvK+Tw2/vwdn+vWljQ8sACatLVpTEKIq6vUfRwDBw6UIiGum1KK6V/v4A3Ne2g1iiy/+3HvPMDWYQkhrsHqcxxCVJU1u07RJ+U9/LVpFLvqSe8RYeuQhBAVcF13jgtRWemGfKLXreETu58A0PR/B5PJ3cZRCSEqQvY4hE28sXYn000ll96aegwveXa4EKJOkMIhatyPf5wlKGE+rbQZFLq3Qhv2hq1DEkJYQQqHqFHZeUWsW/sFQ+1iAHB45D1wlENUQtQlFTrHERFR8ZOW8+bNq3Qwov57+7s4XjUuBg0Ye41C53+brUMSQlipQoVjz549FVpZRWfQFQ3Tbwnn6HrgLXztMsl398Ppnum2DkkIUQkVKhybNm2q7jhEPZdXWMx3az5mtt1WFBqcHvsAHFxtHZYQohLkHIeoER9E7+KlgpKrqIqCnoPWt9g4IiFEZVXqPo74+Hi+/PJLUlJSKCwstFj2zTffVElgov44dMZA852zaa7LItfdH9e7ZU4zIeoyq/c4Vq1aRe/evTl8+DBr166lqKiIgwcPsnHjRjw9PasjRlGHmUyKZV+uYYi25Coq1wHvgr2zjaMSQlwPqwvHm2++yfz58/n+++9xcHBg4cKFHDlyhIEDB9K6devqiFHUYSt+/5MnMxei1SgudRoEfn1sHZIQ4jpZXTiSkpLo168fAA4ODuTm5qLRaHjxxRf58MMPqzxAUXelG/I5/dNCOmlPUGDngXO/N20dkhCiClhdOBo1asTFixcBaNGiBQcOHAAgKyuLvLy8qo1O1GkL125hLKsBsL/3dXBtauOIhBBVweqT47fddhsbNmygS5cuPPbYY4wbN46NGzeyYcMG7rpL5hsSJTYdTefWxLdx0+WT590Tlx6lHzcshKibrN7jePfddxk8eDAAr776KhEREaSlpTFgwAA+/vjjSgWxePFi/Pz8cHJyIjg4mB07dly1/5o1a+jQoQNOTk506dKF9evXWyxXSjF16lR8fHxwdnYmNDSUhIQEiz4PPfQQrVu3xsnJCR8fH5544gnOnDlTqfiFpUuFRn74ejn363ZgQofLwwtBK1d+C1FfWP2/uXHjxvj6+pYM1mp55ZVX+O6775g7dy6NGjWyOoDVq1cTERHBtGnT2L17N926dSMsLIz09PQy+2/fvp0hQ4YwcuRI9uzZQ3h4OOHh4eZDZgCzZ89m0aJFLFmyhLi4OFxdXQkLCyM/P9/c54477uDLL7/k6NGjfP311yQlJfHoo49aHb8obfHPf/BC/gcAGIOeAX0XG0ckhKhSyko//PCDio6OLtX+008/qfXr11u7OhUUFKTGjh1rfm00GpWvr6+KjIwss//AgQNVv379LNqCg4PVmDFjlFJKmUwmpdfr1Zw5c8zLs7KylKOjo1q5cmW5cXz77bdKo9GowsLCCsWdnZ2tAJWdnV2h/pcVFxerw4cPq+LiYqvG1RWHz2arxZOHKTXNQ12adaNS+RcrNK6+5+V6SG7KJnkpX2VzU9HvtUo9c3zWrFml2k0mE6+88gr33XdfhddVWFjIrl27mDhxorlNq9USGhpKbGxsmWNiY2NLTboYFhZGVFQUAMnJyaSmphIaGmpe7unpSXBwMLGxsebDbFfKzMzkiy++oHfv3tjb25f5vgUFBRQUFJhfGwwGAIxGI0ajsWIb/Fd/k8lk1Zi6QinFh1//yCzt9wDYPzAbo50zVGBb63NerpfkpmySl/JVNjcV7W914UhISKBTp06l2jt06EBiYqJV6zp37hxGoxFvb2+Ldm9vb44cOVLmmNTU1DL7p6ammpdfbiuvz2X//e9/effdd8nLy+OWW25h3bp15cYaGRnJ9OmlJ+VLSkrCzc2t3HH/ZDKZyMzMJDExEW09O+4fl5JDeOo7OOiMnG/emwztjfCPc0vlqc95uV6Sm7JJXspX2dzk5ORUqJ/VhcPT05M///wTPz8/i/bExERcXevWpHUvv/wyI0eO5MSJE0yfPp1hw4axbt26Mmf5nThxosWejsFgoFWrVrRt2xYPD48Kv6fRaCQxMZF27dqh0+mqZDtqg2Kjic++XcRw3R8Ua+zxeuwdvBr7V3h8fc1LVZDclE3yUr7K5ubykZRrsbpw9O/fn/Hjx7N27Vratm0LlBSN//znPzz00ENWratp06bodDrS0tIs2tPS0tDr9WWO0ev1V+1/+c+0tDR8fHws+gQGBpZ6/6ZNm3LjjTfSsWNHWrVqxe+//05ISEip93V0dMTR0bFUu06ns/pDq9VqKzWuNluz8zgjcj8CLRhvHo1js3ZWr6M+5qWqSG7KJnkpX2VyU9G+Vu/fzZ49G1dXVzp06IC/vz/+/v507NiRJk2a8Pbbb1u1LgcHB3r27ElMTIy5zWQyERMTU+aXN0BISIhFf4ANGzaY+/v7+6PX6y36GAwG4uLiyl3n5fcFLM5jiIrJLSgm5ef3aKc9wyX7RjjeMcHWIQkhqlGlDlVt376dDRs2sG/fPpydnenatSu33Va5J7lFREQwfPhwevXqRVBQEAsWLCA3N5cRI0YAMGzYMFq0aEFkZCQA48aNo2/fvsydO5d+/fqxatUq4uPjzdOdaDQaxo8fz8yZMwkICMDf358pU6bg6+tLeHg4AHFxcezcuZM+ffrQqFEjkpKSmDJlCm3btr1qcRFlW75xH08bV4EG7O+aBM5etg5JCFGNKjWtukaj4Z577uGee+657gAGDRpERkYGU6dOJTU1lcDAQKKjo80nt1NSUixO7vTu3ZsVK1YwefJkJk2aREBAAFFRUXTu3NncZ8KECeTm5jJ69GiysrLo06cP0dHRODk5AeDi4sI333zDtGnTyM3NxcfHh3vvvZfJkyeXeThKlC/dkI9T7Dwaa3O46N4W95ufsnVIQohqplFKqWt1WrRoEaNHj8bJyYlFixZdte8LL7xQZcHVZgaDAU9PT7Kzs60+OZ6QkEBAQEC9OC47Z+WPjDsyFAeNEfXvNWhurNwvE/UtL1VJclM2yUv5Kpubin6vVWiPY/78+QwdOhQnJyfmz59fbj+NRtNgCoeAY2kX6XJoLg46I9m+t+FZyaIhhKhbKlQ4kpOTy/y7aNi+Wfslr+h2YkKLZ/+3bB2OEKKGWH1V1YwZM8qcPv3SpUvMmDGjSoIStV9sYgb3n3kHgJybhoJ36ZtChRD1k9WFY/r06WXeXZiXl1fmndWi/lFKEffdB3TVJpOvdcXjvmm2DkkIUYOsLhxKqTLvrN63bx+NGzeukqBE7bb58GnCsz4FoDjkBXBrZuOIhBA1qcKX4zZq1AiNRoNGo+HGG2+0KB5Go5GcnByeeeaZaglS1B4mk+KPde9zhzaNXLtGuN32vK1DEkLUsAoXjgULFqCU4qmnnmL69Ol4enqalzk4OODn5yc3zzUA0fuO82juCtCA5rb/gGPFJ3gUQtQPFS4cw4cPB0qm9Lja9OOi/io2mkj+8R3u12Ry0cEb95BRtg5JCGEDVt853rdvX4xGI1999RWHDx8GoFOnTvTv3x87u0rdiC7qiO92JDCoYM1fU4u8AvZOtg5JCGEDVn/THzx4kIceeojU1FTat28PwFtvvUWzZs34/vvvLab+EPVHQbGRjJiFNNUYyHZuhWevJ2wdkhDCRqy+qurpp5/mpptu4tSpU+zevZvdu3dz8uRJunbtyujRo6sjRlELfP3bAYYUrQXA+Z7JoJNDlUI0VFbvcezdu5f4+HgaNWpkbmvUqBFvvPEGN998c5UGJ2qHvMJi8rcswEOTR5ZbO7y6PWbrkIQQNmT1HseNN95Y6kFKAOnp6bRrZ/3De0Tt9+WWPQwy/QCA232vgVYmlBOiIbO6cERGRvLCCy/w1VdfcerUKU6dOsVXX33F+PHjeeuttzAYDOYfUfflFxnRbV+Eq6aATK/O2HV6wNYhCSFszOpDVQ88UPLFMXDgQPNNgJdnZn/wwQfNrzUaDUajsariFDby3fb9PGL6CTTgce8UKGPWACFEw2J14di0aVN1xCFqoSKjibyti3HVFHDevQNN2ofZOiQhRC1Qqfs4RMOwfudRHin+ATTgfs8rsrchhAAq+ejYrKwsPv74Y/MNgDfddBNPPfWUxTQkom4zmRQZGxfjockj08Wfxjf1t3VIQohawuqT4/Hx8bRt25b58+eTmZlJZmYm8+bNo23btuzevbs6YhQ2ELP/Tx4uiALA+c6XQWv1R0UIUU9Zvcfx4osv8tBDD/G///3PPMVIcXExTz/9NOPHj2fr1q1VHqSoWUopTmx4n7s1F8lybIFX90G2DkkIUYtYXTji4+MtigaAnZ0dEyZMoFevXlUanLCNbUdP80DOVyVzUvX9D+hkDjIhxN+sPv7g4eFBSkpKqfaTJ0/i7u5eJUEJ2zr84xL0mgtk2zfHNehxW4cjhKhlrC4cgwYNYuTIkaxevZqTJ09y8uRJVq1axdNPP82QIUOqI0ZRg/YczyDswqqSF7e+AHaOtg1ICFHrWH0M4u2330aj0TBs2DCKi4sBsLe359lnn2XWrFlVHqCoWXuil/OUNoOLOi88e4+0dThCiFrI6sLh4ODAwoULiYyMJCkpCYC2bdvi4uJS5cGJmnXyfC43n/kctJDffSTuDvJvKoQordJnPV1cXOjSpUtVxiJsbONPaxmuTaZQ40CzO56zdThCiFpKLs4XAGRfKqL10U8AONf2EXBtauOIhBC1lRQOAUD05q3codkFgE/Yf2wcjRCiNpPCISgymnDYuQSA0953oGl2o40jEkLUZlI4BL/EH+Q+Y8msx83ukb0NIcTV1YrCsXjxYvz8/HByciI4OJgdO3Zctf+aNWvo0KEDTk5OdOnShfXr11ssV0oxdepUfHx8cHZ2JjQ0lISEBPPy48ePM3LkSPz9/XF2dqZt27ZMmzaNwsLCatm+2kwpxfnN7+OkKSLNrRMON/SxdUhCiFrO5oVj9erVREREMG3aNHbv3k23bt0ICwsjPT29zP7bt29nyJAhjBw5kj179hAeHk54eDgHDhww95k9ezaLFi1iyZIlxMXF4erqSlhYGPn5+QAcOXIEk8nEBx98wMGDB5k/fz5Llixh0qRJNbLNtcnOxLOE5X0PgOvt42XqdCHEtSkbCwoKUmPHjjW/NhqNytfXV0VGRpbZf+DAgapfv34WbcHBwWrMmDFKKaVMJpPS6/Vqzpw55uVZWVnK0dFRrVy5stw4Zs+erfz9/Sscd3Z2tgJUdnZ2hccopVRxcbE6fPiwKi4utmpcdfl08etKTfNQmTMDlCouslkctS0vtYnkpmySl/JVNjcV/V6z6ex1hYWF7Nq1i4kTJ5rbtFotoaGhxMbGljkmNjaWiIgIi7awsDCioqIASE5OJjU1ldDQUPNyT09PgoODiY2NZfDgwWWuNzs7m8aNG5cba0FBAQUFBebXl5+pbjQarXpErtFoxGQy1YrH6p7NukSP1C9BC4Xdn8KIBmwUV23KS20juSmb5KV8lc1NRfvbtHCcO3cOo9GIt7e3Rbu3tzdHjhwpc0xqamqZ/VNTU83LL7eV1+efEhMTeeedd3j77bfLjTUyMpLp06eXak9KSsLNza3ccf9kMpnIzMwkMTERrY2fcbFx2688pz1BAQ5kt7iDC1ecB6pptSkvtY3kpmySl/JVNjc5OTkV6tfg58s+ffo09957L4899hijRo0qt9/EiRMt9nQMBgOtWrWibdu2eHh4VPj9jEYjiYmJtGvXDp1Od12xX4+CYhP+q0uuoEpr8wBtO9t2SvzakpfaSHJTNslL+Sqbm8tHUq7FpoWjadOm6HQ60tLSLNrT0tLQ6/VljtHr9Vftf/nPtLQ0fHx8LPoEBgZajDtz5gx33HEHvXv35sMPP7xqrI6Ojjg6lp4pVqfTWf2h1Wq1lRpXlTbv2M/dpljQgO/d42rFf7zakJfaSnJTNslL+SqTm4r2ten+nYODAz179iQmJsbcZjKZiImJISQkpMwxISEhFv0BNmzYYO7v7++PXq+36GMwGIiLi7NY5+nTp7n99tvp2bMnS5cubXC7ulm//g97jZEzHt2waxlo63CEEHWIzQ9VRUREMHz4cHr16kVQUBALFiwgNzeXESNGADBs2DBatGhBZGQkAOPGjaNv377MnTuXfv36sWrVKuLj4817DBqNhvHjxzNz5kwCAgLw9/dnypQp+Pr6Eh4eDvxdNNq0acPbb79NRkaGOZ7y9nTqk33H0wnN/QE04HabTGYohLCOzQvHoEGDyMjIYOrUqaSmphIYGEh0dLT55HZKSorF3kDv3r1ZsWIFkydPZtKkSQQEBBAVFUXnzp3NfSZMmEBubi6jR48mKyuLPn36EB0djZOTE1Cyh5KYmEhiYiItW7a0iEcpVQNbbVv7NnxON00WBl1jPAIfsXU4Qog6RqMawjdlNTAYDHh6epKdnW31yfGEhAQCAgJsclw2M7eQ5Lf60FN7lLPdX8Sn/2s1HkNZbJ2X2kxyUzbJS/kqm5uKfq81rAP7gs1bNtJTe5RidPjc+YytwxFC1EFSOBoQpRTaPZ8CcNr7TnCv/+dzhBBVTwpHA7LnzzPcWVgyC27T28fYOBohRF0lhaMBORbzGR6aPM7Z++La/i5bhyOEqKOkcDQQhvwiOpz+CoD8ro9DA7tvRQhRdeTbo4HY+utmAjUJFKOjxe0jbR2OEKIOk8LRQKhdywBIaX4HGjkpLoS4DlI4GoCDJ1Lpe2kjAE37yklxIcT1kcLRAByN+RQPTR4Z9r54dAy99gAhhLgKKRz1XH6RkdYp3wCQ22mwnBQXQlw3+Rap536L+51eHMaIllZ3yElxIcT1k8JRz+XGfQZAilcwOq+W1+gthBDXJoWjHjt57iJBhp8AcA8ZYeNohBD1hRSOemznxq/x0WRyUetO057htg5HCFFPSOGop0wmheeR1QBk+D0EdqUfeyuEEJUhhaOe2nEokT7GHQC0uGO0jaMRQtQnUjjqqdO/foajppgzTgE4tgq0dThCiHpECkc9lFtQTEDqOgCMXQfbOBohRH0jhaMe2v77NrpqkjCipeVtw2wdjhCinpHCUQ9dil8BwIlGt6Jxa27jaIQQ9Y0UjnomLTuPXoYNALgHP2HjaIQQ9ZEUjnpmx+bv8NWcJ1fjSrOe/W0djhCiHpLCUc84HvwSgDMt7wN7JxtHI4Soj6Rw1CNHT6bRu2AbAD7/etK2wQgh6i0pHPXIsS2rcNPkk2Hng1tAH1uHI4Sop6Rw1BNKKRr/+S0AWW3DQaOxbUBCiHpLCkc9se/YnwQZ9wLQ+vbhtg1GCFGvSeGoJ079+gX2GiMnnW7E0aejrcMRQtRjUjjqgWKjiVanSqYYKeg4wMbRCCHqOykc9cDufXvpxlFMaGjTV276E0JULykc9cD5378AINmtB/ZeLWwcjRCivrN54Vi8eDF+fn44OTkRHBzMjh07rtp/zZo1dOjQAScnJ7p06cL69estliulmDp1Kj4+Pjg7OxMaGkpCQoJFnzfeeIPevXvj4uKCl5dXVW9SjcovLCYgLRoAbbdBNo5GCNEQ2LRwrF69moiICKZNm8bu3bvp1q0bYWFhpKenl9l/+/btDBkyhJEjR7Jnzx7Cw8MJDw/nwIED5j6zZ89m0aJFLFmyhLi4OFxdXQkLCyM/P9/cp7CwkMcee4xnn3222rexusXv+I12mlMUYkebWwfaOhwhREOgbCgoKEiNHTvW/NpoNCpfX18VGRlZZv+BAweqfv36WbQFBwerMWPGKKWUMplMSq/Xqzlz5piXZ2VlKUdHR7Vy5cpS61u6dKny9PSsVOzZ2dkKUNnZ2VaNKy4uVocPH1bFxcWVet9/il7wjFLTPNSRBQ9WyfpsparzUp9IbsomeSlfZXNT0e81O1sVrMLCQnbt2sXEiRPNbVqtltDQUGJjY8scExsbS0REhEVbWFgYUVFRACQnJ5OamkpoaKh5uaenJ8HBwcTGxjJ4cOUfalRQUEBBQYH5tcFgAMBoNGI0Giu8HqPRiMlksmpMeXLyi+ic+TNowCHwsSpZp61UZV7qG8lN2SQv5atsbira32aF49y5cxiNRry9vS3avb29OXLkSJljUlNTy+yfmppqXn65rbw+lRUZGcn06dNLtSclJeHm5lbh9ZhMJjIzM0lMTESrvb4jhYf2/s4jmnPk4URB486lzuXUJVWZl/pGclM2yUv5KpubnJycCvWzWeGoayZOnGixt2MwGGjVqhVt27bFw8OjwusxGo0kJibSrl07dDrddcV06oc3AUhuejsdOnW9rnXZWlXmpb6R3JRN8lK+yubm8pGUa7FZ4WjatCk6nY60tDSL9rS0NPR6fZlj9Hr9Vftf/jMtLQ0fHx+LPoGBgdcVr6OjI46OjqXadTqd1R9arVZbqXFXysrJI9CwCTTgFTSkXvzHqYq81FeSm7JJXspXmdxUtK/N9u8cHBzo2bMnMTEx5jaTyURMTAwhISFljgkJCbHoD7BhwwZzf39/f/R6vUUfg8FAXFxcueusq/b8uo6mGgMGjTstevazdThCiAbEpoeqIiIiGD58OL169SIoKIgFCxaQm5vLiBEjABg2bBgtWrQgMjISgHHjxtG3b1/mzp1Lv379WLVqFfHx8Xz44YcAaDQaxo8fz8yZMwkICMDf358pU6bg6+tLeHi4+X1TUlLIzMwkJSUFo9HI3r17AWjXrp1V5yts6sA3AJzU381NOnsbByOEaEhsWjgGDRpERkYGU6dOJTU1lcDAQKKjo80nt1NSUixO7PTu3ZsVK1YwefJkJk2aREBAAFFRUXTu3NncZ8KECeTm5jJ69GiysrLo06cP0dHRODn9/TS8qVOnsnz5cvPr7t27A7Bp0yZuv/32at7q65eRdZHuOVtBA01v+betwxFCNDAapZSydRB1kcFgwNPTk+zsbKtPjickJBAQEFDp47Ix333KXbv/j0xtYxpPTgRt3T++WxV5qa8kN2WTvJSvsrmp6PeaXMNWB9kfWgvAmRb31ouiIYSoW6Rw1DGp5y/Q49J2APS95TCVEKLmSeGoYw5tWYObJp90bXOadpDnigshap4UjjrG5WgUAGdb9ZPnigshbEIKRx2Slp5BYH7JtPO+fYbaOBohREMlhaMOObJlNU6aIk7rWtKsXS9bhyOEaKCkcNQhrglRAKS3vl8OUwkhbEYKRx1x+vQpuhXsBqB132E2jkYI0ZBJ4agjEraswF5j5Lh9W5r4dbF1OEKIBkwKRx3R6M/vAMi64SEbRyKEaOikcNQBfyYdo0tRyXPV/e+Qw1RCCNuSwlEHpGz9HK1GkeB4E576G2wdjhCigZPCUcsppfBN+R6ASx0G2DgaIYSQwlHrHf4jnhvVnxQpHe3ueNzW4QghhBSO2i5j22cAHHULwsXL28bRCCGEFI5arbComIC09QBoug2ycTRCCFFCCkcttn/bj/iSQQ4utL9toK3DEUIIQApHrVa4+wsAEpqGYufkauNohBCihBSOWurChQt0zd4EgFfIEzaORggh/iaFo5Y6+MtnuGnyOa31wb/H3bYORwghzKRw1FKeR1cDkHrDAJkJVwhRq0jhqIWOHdxFl+IDGJWGG0JH2TocIYSwIIWjFkrf9AEAh9xCaKT3s20wQgjxD1I4apmLF7PpnPEDAHbBI20cjRBClCaFo5bZv/5/eGlyOKtpTodbw20djhBClCKFoxYxGk34HlkOwJkbh6HR2dk4IiGEKE0KRy2yO2Y1/iqFHJzpcP9ztg5HCCHKJIWjllAmE247FgBw2HcArp5NbBuQEEKUQwpHLbFn4xo6Fh8hX9nTrv9/bR2OEEKUSwpHLVBUVEij7TMB2O/7GI28W9s4IiGEKJ8Ujlpgx6o38DelkIUbHQbNsHU4QghxVbWicCxevBg/Pz+cnJwIDg5mx44dV+2/Zs0aOnTogJOTE126dGH9+vUWy5VSTJ06FR8fH5ydnQkNDSUhIcGiT2ZmJkOHDsXDwwMvLy9GjhxJTk5OlW/btSTuj6VX4mIAkrq+jIdXsxqPQQghrGHzwrF69WoiIiKYNm0au3fvplu3boSFhZGenl5m/+3btzNkyBBGjhzJnj17CA8PJzw8nAMHDpj7zJ49m0WLFrFkyRLi4uJwdXUlLCyM/Px8c5+hQ4dy8OBBNmzYwLp169i6dSujR4+u9u29UkrCPty/GYqjpoj9zsH0CH+hRt9fCCEqRdlYUFCQGjt2rPm10WhUvr6+KjIyssz+AwcOVP369bNoCw4OVmPGjFFKKWUymZRer1dz5swxL8/KylKOjo5q5cqVSimlDh06pAC1c+dOc58ff/xRaTQadfr06QrFnZ2drQCVnZ1dsQ1VSmVfOKfi1y9VGz6brbYtGasuTvVWapqHOj79JpWVmV7h9dRHxcXF6vDhw6q4uNjWodQ6kpuySV7KV9ncVPR7zaZ3mBUWFrJr1y4mTpxobtNqtYSGhhIbG1vmmNjYWCIiIizawsLCiIqKAiA5OZnU1FRCQ0PNyz09PQkODiY2NpbBgwcTGxuLl5cXvXr1MvcJDQ1Fq9USFxfHww8/XOp9CwoKKCgoML82GAwAGI1GjEZjhbY341QiPePG/d2ggcP2N9F0xErcPBpXeD31kdFoxGQyNegclEdyUzbJS/kqm5uK9rdp4Th37hxGoxFvb2+Ldm9vb44cOVLmmNTU1DL7p6ammpdfbrtan+bNm1sst7Ozo3HjxuY+/xQZGcn06dNLtSclJeHm5lbeJlq4cP4C+boOGJWGHEdv8lrehm/3ezlvuMR5Q8K1V1CPmUwmMjMzSUxMRKu1+RHUWkVyUzbJS/kqm5uKnueVOS0qaOLEiRZ7OgaDgVatWtG2bVs8PDwqtpKAAIwhd5KYmEivdu3Q6XTVFG3dYzQaSUxMpJ3kpRTJTdkkL+WrbG4uH0m5FpsWjqZNm6LT6UhLS7NoT0tLQ6/XlzlGr9dftf/lP9PS0vDx8bHoExgYaO7zz5PvxcXFZGZmlvu+jo6OODo6lmrX6XRWf2i1Wm2lxtV3kpfySW7KJnkpX2VyU9G+Nt2/c3BwoGfPnsTExJjbTCYTMTExhISElDkmJCTEoj/Ahg0bzP39/f3R6/UWfQwGA3FxceY+ISEhZGVlsWvXLnOfjRs3YjKZCA4OrrLtE0KI+sjmh6oiIiIYPnw4vXr1IigoiAULFpCbm8uIESMAGDZsGC1atCAyMhKAcePG0bdvX+bOnUu/fv1YtWoV8fHxfPjhhwBoNBrGjx/PzJkzCQgIwN/fnylTpuDr60t4eDgAHTt25N5772XUqFEsWbKEoqIinn/+eQYPHoyvr69N8iCEEHWFzQvHoEGDyMjIYOrUqaSmphIYGEh0dLT55HZKSorFyZ3evXuzYsUKJk+ezKRJkwgICCAqKorOnTub+0yYMIHc3FxGjx5NVlYWffr0ITo6GicnJ3OfL774gueff5677roLrVbLgAEDWLRoUc1tuBBC1FEapZSydRB1kcFgwNPTk+zs7IqfHKfkpFVCQgIBAQFyXPYKkpfySW7KJnkpX2VzU9HvNbmGTQghhFWkcAghhLCKFA4hhBBWkcIhhBDCKlI4hBBCWEUKhxBCCKvY/D6OuuryVcwVndvlMqPRSE5ODgaDQS4hvILkpXySm7JJXspX2dxc/j671l0aUjgq6eLFiwC0atXKxpEIIUTVunjxIp6enuUulxsAK8lkMnHmzBnc3d3RaDQVHnd5Vt2TJ09adeNgfSd5KZ/kpmySl/JVNjdKKS5evIivr+9Vp2OXPY5K0mq1tGzZstLjPTw85MNeBslL+SQ3ZZO8lK8yubnansZlcnJcCCGEVaRwCCGEsIoUjhrm6OjItGnTynwoVEMmeSmf5KZskpfyVXdu5OS4EEIIq8gehxBCCKtI4RBCCGEVKRxCCCGsIoVDCCGEVaRw1KDFixfj5+eHk5MTwcHB7Nixw9Yh1bitW7fy4IMP4uvri0ajISoqymK5UoqpU6fi4+ODs7MzoaGhJCQk2CbYGhQZGcnNN9+Mu7s7zZs3Jzw8nKNHj1r0yc/PZ+zYsTRp0gQ3NzcGDBhAWlqajSKuOe+//z5du3Y138wWEhLCjz/+aF7eUPPyT7NmzUKj0TB+/HhzW3XlRgpHDVm9ejURERFMmzaN3bt3061bN8LCwkhPT7d1aDUqNzeXbt26sXjx4jKXz549m0WLFrFkyRLi4uJwdXUlLCyM/Pz8Go60Zm3ZsoWxY8fy+++/s2HDBoqKirjnnnvIzc0193nxxRf5/vvvWbNmDVu2bOHMmTM88sgjNoy6ZrRs2ZJZs2axa9cu4uPjufPOO+nfvz8HDx4EGm5errRz504++OADunbtatFebblRokYEBQWpsWPHml8bjUbl6+urIiMjbRiVbQFq7dq15tcmk0np9Xo1Z84cc1tWVpZydHRUK1eutEGEtpOenq4AtWXLFqVUSR7s7e3VmjVrzH0OHz6sABUbG2urMG2mUaNG6qOPPpK8KKUuXryoAgIC1IYNG1Tfvn3VuHHjlFLV+5mRPY4aUFhYyK5duwgNDTW3abVaQkNDiY2NtWFktUtycjKpqakWefL09CQ4OLjB5Sk7OxuAxo0bA7Br1y6KioosctOhQwdat27doHJjNBpZtWoVubm5hISESF6AsWPH0q9fP4scQPV+ZmSSwxpw7tw5jEYj3t7eFu3e3t4cOXLERlHVPqmpqQBl5unysobAZDIxfvx4br31Vjp37gyU5MbBwQEvLy+Lvg0lN3/88QchISHk5+fj5ubG2rVr6dSpE3v37m3QeVm1ahW7d+9m586dpZZV52dGCocQtczYsWM5cOAAv/32m61DqTXat2/P3r17yc7O5quvvmL48OFs2bLF1mHZ1MmTJxk3bhwbNmzAycmpRt9bDlXVgKZNm6LT6UpdzZCWloZer7dRVLXP5Vw05Dw9//zzrFu3jk2bNllM26/X6yksLCQrK8uif0PJjYODA+3ataNnz55ERkbSrVs3Fi5c2KDzsmvXLtLT0+nRowd2dnbY2dmxZcsWFi1ahJ2dHd7e3tWWGykcNcDBwYGePXsSExNjbjOZTMTExBASEmLDyGoXf39/9Hq9RZ4MBgNxcXH1Pk9KKZ5//nnWrl3Lxo0b8ff3t1jes2dP7O3tLXJz9OhRUlJS6n1uymIymSgoKGjQebnrrrv4448/2Lt3r/mnV69eDB061Pz3asvNdZ1aFxW2atUq5ejoqJYtW6YOHTqkRo8erby8vFRqaqqtQ6tRFy9eVHv27FF79uxRgJo3b57as2ePOnHihFJKqVmzZikvLy/17bffqv3796v+/fsrf39/denSJRtHXr2effZZ5enpqTZv3qzOnj1r/snLyzP3eeaZZ1Tr1q3Vxo0bVXx8vAoJCVEhISE2jLpmvPLKK2rLli0qOTlZ7d+/X73yyitKo9Gon3/+WSnVcPNSliuvqlKq+nIjhaMGvfPOO6p169bKwcFBBQUFqd9//93WIdW4TZs2KaDUz/Dhw5VSJZfkTpkyRXl7eytHR0d11113qaNHj9o26BpQVk4AtXTpUnOfS5cuqeeee041atRIubi4qIcfflidPXvWdkHXkKeeekq1adNGOTg4qGbNmqm77rrLXDSUarh5Kcs/C0d15UamVRdCCGEVOcchhBDCKlI4hBBCWEUKhxBCCKtI4RBCCGEVKRxCCCGsIoVDCCGEVaRwCCGEsIoUDiFqidtvv93i6W3X2/efT1g8cuQIt9xyC05OTgQGBlY6TiGkcAhRT509e5b77rvP/HratGm4urpy9OhRYmJiWLZsWakpt4WoCJlWXYgqVlhYiIODg63DKDUDalJSEv369aNNmzY2ikjUF7LHIcR1uv3223n++ecZP348TZs2JSwsjAMHDnDffffh5uaGt7c3TzzxBOfOnTOPyc3NZdiwYbi5ueHj48PcuXNLrfe9994jICAAJycnvL29efTRRy2Wm0wmJkyYQOPGjdHr9bz22msWy688VKXRaNi1axczZsxAo9Fw++23M2LECLKzs9FoNGg0mlLjhSiPFA4hqsDy5ctxcHBg27ZtzJo1izvvvJPu3bsTHx9PdHQ0aWlpDBw40Nz/5ZdfZsuWLXz77bf8/PPPbN68md27d5uXx8fH88ILLzBjxgyOHj1KdHQ0t912W6n3dHV1JS4ujtmzZzNjxgw2bNhQZnxnz57lpptu4j//+Q9nz57lu+++Y8GCBXh4eHD27FnOnj3LSy+9VD3JEfWOHKoSogoEBAQwe/ZsAGbOnEn37t158803zcs/+eQTWrVqxbFjx/D19eXjjz/m888/56677gJKisCVD25KSUnB1dWVBx54AHd3d9q0aUP37t0t3rNr165MmzbN/P7vvvsuMTEx3H333aXi0+v12NnZ4ebmZj6E5enpiUajqfcPPBJVTwqHEFWgZ8+e5r/v27ePTZs24ebmVqpfUlISly5dorCwkODgYHN748aNad++vfn13XffTZs2bbjhhhu49957uffee3n44YdxcXEx9+natavFun18fEhPT6/KzRKiTHKoSogq4Orqav57Tk4ODz74oMWT2fbu3UtCQkKpw03lcXd3Z/fu3axcuRIfHx+mTp1Kt27dLB4Dam9vbzFGo9FgMpmqZHuEuBopHEJUsR49enDw4EH8/Pxo166dxY+rqytt27bF3t6euLg485gLFy5w7Ngxi/XY2dkRGhrK7Nmz2b9/P8ePH2fjxo1VFqeDgwNGo7HK1icaDikcQlSxsWPHkpmZyZAhQ9i5cydJSUn89NNPjBgxAqPRiJubGyNHjuTll19m48aNHDhwgCeffBKt9u//juvWrWPRokXs3buXEydO8Omnn2IymSwOZ10vPz8/cnJyiImJ4dy5c+Tl5VXZukX9JoVDiCrm6+vLtm3bMBqN3HPPPXTp0oXx48fj5eVlLg5z5szhX//6Fw8++CChoaH06dPH4jyJl5cX33zzDXfeeScdO3ZkyZIlrFy5kptuuqnK4uzduzfPPPMMgwYNolmzZuaT+0Jcizw6VgghhFVkj0MIIYRVpHAIIYSwihQOIYQQVpHCIYQQwipSOIQQQlhFCocQQgirSOEQQghhFSkcQgghrCKFQwghhFWkcAghhLCKFA4hhBBWkcIhhBDCKv8PgQFkHkFHc/AAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "zs = np.geomspace(0.001, 40, 100)\n", "tau_sis = od.strong_lensing_optical_depth(zs)\n", "tau_sis_haris = od.optical_depth_SIS_haris(zs)\n", "\n", "# plot \n", "plt.figure(figsize=(4,4))\n", "plt.plot(zs, tau_sis, label='SIS hemanta')\n", "plt.plot(zs, tau_sis_haris, label='SIS haris')\n", "plt.xlabel('redshift')\n", "plt.ylabel('optical depth')\n", "plt.legend()\n", "plt.grid(alpha=0.5)\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "* close enough" ] }, { "cell_type": "markdown", "metadata": {}, "source": [] } ], "metadata": { "kernelspec": { "display_name": "ler", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.10.12" } }, "nbformat": 4, "nbformat_minor": 2 }