{ "cells": [ { "cell_type": "markdown", "id": "a1668795-158d-4a33-b444-f561e49fce02", "metadata": {}, "source": [ "# Gaussian process regression" ] }, { "cell_type": "code", "execution_count": 46, "id": "43ab5992-63e0-45ee-bed0-ef0511be5687", "metadata": {}, "outputs": [], "source": [ "import matplotlib.pyplot as plt\n", "import numpy as np\n", "\n", "from pyinla.model import *\n", "from pyinla.spde import *\n", "from pyinla.utils import set_pardiso_license" ] }, { "cell_type": "code", "execution_count": 51, "id": "4ec4ecc1-f8b1-49bc-b484-fccc39034df8", "metadata": {}, "outputs": [ { "data": { "text/html": [ "\n", " StrVector with 1 elements.\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
\n", " 'SUCCESS: PARDISO IS INSTALLED AND WORKING'\n", "
\n", " " ], "text/plain": [ " [RTYPES.STRSXP]\n", "R classes: ('character',)\n", "['SUCCESS: PARDISO IS INSTALLED AND WORKING']" ] }, "execution_count": 51, "metadata": {}, "output_type": "execute_result" } ], "source": [ "set_pardiso_license(\"/Users/js5013/.pardiso.lic\")" ] }, { "cell_type": "code", "execution_count": 119, "id": "9e33c82d-9d08-45b9-89b6-af1758932b8b", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 119, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiIAAAGdCAYAAAAvwBgXAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAA9hAAAPYQGoP6dpAABirElEQVR4nO3deXhU5fk38O9MMpksJBNCIBMgQNgJgQRQBEE2gyC4gNYFtXUrVgutArWAvyJStYj6FqtSUVxoawE3EAHFsqgIhp0AISAkhC1kgCRk32fO+0c4cTKZ5ZzZzizfz3VxXSQ5M/PkZJJzn+e57/tRCYIggIiIiEgBaqUHQERERMGLgQgREREphoEIERERKYaBCBERESmGgQgREREphoEIERERKYaBCBERESmGgQgREREpJlTpAdhjMplw8eJFREdHQ6VSKT0cIiIikkAQBFRUVKBjx45Qq+3Pefh0IHLx4kUkJSUpPQwiIiJywvnz59G5c2e7x/h0IBIdHQ2g6RuJiYlReDREREQkRXl5OZKSkpqv4/b4dCAiLsfExMQwECEiIvIzUtIqmKxKREREimEgQkRERIphIEJERESKYSBCREREimEgQkRERIphIEJERESKYSBCREREimEgQkRERIphIEJERESKYSBCREREimEgQkRERIphIEJERESKYSBCREQUhC6X1+K3/9qPjUcuQhAExcbh07vvEhERkWdsOFKIrccv4Wp1PW4b2FGxcXBGhIiIKAitzyoAANyZrlwQAjAQISIiCjqnr1TiyIUyhKhVmDQgUdGxMBAhIiIKMl8dvggAGNkzHvFttIqOhYEIERFREBEEAV9lNQUiUwYpuywDMBAhIiIKKtkF5ThdVIVwjRrjU/RKD4eBCBERUTD58lqSaka/BLTRKl88y0CEiIgoSBhNAjZcyw+5M72TwqNpwkCEiIgoSOw5XYzLFXXQRWgwund7pYcDgIEIERFR0Fh/LUl10gA9wkJ9IwTwjVEQERGRR9U1GvF1diEA4I4031iWARiIEBERBYXvf76CitpG6GPCcUNynNLDacZAhIiIKAiITcxuT0uEWq1SeDS/YCBCREQU4KrrG7H9+GUAwO1pyjcxM8dAhIiIKMBtPX4ZNQ1GdG0XiQGddEoPpwUGIkRERAFu47VlmdsGJkKl8p1lGYCBCBERUUArr23A9z9fAeB7yzIAAxEiIqKA9r9jl1BvNKFXhzbokxCt9HBaYSBCREQUwDYeEZdlOvrcsgzAQISIiChglVTVY+epIgDAbWmJCo/GOgYiREREAWpztgGNJgEpiTHo0b6N0sOxioEIERFRgBKXZXwxSVXEQISIiCgAXa6oxe7TxQCaynZ9FQMRIiKiAPTNUQNMApCeFIukuEilh2MTAxEiIqIAtOGw7y/LAAxEiIiIAs7F0hrsP3sVKhUweYDvLssADESIiIgCzqYjhQCA67vFQa8LV3g09jEQISIiCjAb/KBaRsRAhIiIKICcKarCkQtlCFGrcGuqXunhOMRAhIiIKIBsOtq0LHNjj3aIb6NVeDSOMRAhIiIKIM3VMgN9f1kGYCBCREQUME5dqsAJQwU0ISpM6O/7yzIAAxEiIqKAseFatcyoXu2hi9QoPBppGIgQEREFAEEQsNFPmpiZYyBCREQUAHIKy3G6qAraUDUyUhKUHo5kDESIiIgCwIbDTcsy4/p2QBttqMKjkY6BCBERkZ8TBAEb/aiJmTkGIkRERH4u63wpLlytQWRYCMb26aD0cGRhIEJEROTnxGWZ8SkJiAgLUXg08jAQISIi8mMmk4BNR/2riZk5BiJERER+bN+ZElwqr0N0eChu6h2v9HBkYyBCRETkx8Sddif210Mb6l/LMgADESIiIr/VaDThm6MGAP5XLSPyaCCyePFiXH/99YiOjkaHDh0wZcoU/Pzzz558SSIioqCReboYxVX1iIsKw4092ik9HKd4NBD54YcfMGPGDOzevRtbtmxBQ0MDbrnlFlRVVXnyZYmIiIKCuNPural6hIb45yKHR1uvbd68ucXHK1euRIcOHXDgwAGMGjXKky9NREQU0OobTdic7d/LMoCHAxFLZWVlAIC4uDirX6+rq0NdXV3zx+Xl5V4ZFxERkb/58dQVlNc2okO0Ftd3s35d9Qdem8cxmUx45plnMGLECKSmplo9ZvHixdDpdM3/kpKSvDU8IiIiv7LxSFMTs8kDExGiVik8Gud5LRCZMWMGsrOzsWbNGpvHzJ8/H2VlZc3/zp8/763hERER+Y3aBiP+d6xpWeY2P2xiZs4rSzMzZ87Exo0bsWPHDnTu3NnmcVqtFlqt1htDIiIi8lvfnbiMqnojOsVGYHCXWKWH4xKPBiKCIOAPf/gD1q1bh++//x7JycmefDkiIqKgIC7L3JaWCJXKf5dlAA8HIjNmzMCqVauwfv16REdHw2BomkbS6XSIiIjw5EsTUQAzmgTszS/B5YpadIgOx9DkOL9eIyeSo6quEdtOXALgn3vLWPJoIPLOO+8AAMaMGdPi8x999BEeeeQRT740EQWozdmFWLQhB4Vltc2fS9SFY+HtKZiYmqjgyIi8Y+vxS6htMCE5Pgr9O8YoPRyXeXxphojIXTZnF+Kpjw/C8i+LoawWT318EO88NJjBCAW8DYeblmVuH+j/yzIA95ohIj9hNAlYtCGnVRACoPlzizbkwGjiDRAFrrKaBvxw8jIA4DY/bmJmjoEIEfmFvfklLZZjLAkACstqsTe/xHuDIvKy/x0zoMEooE9CNHonRCs9HLdgIEJEfuFyhe0gxJnjiPzRBrFaZmDgLEF6tcU7EZE99qphOkSHS3oOqccR+Zviyjrsyi0CEDjLMgADESLyEY6qYYYmxyFRFw5DWa3VPBEVAL2uKXghCkSbjxlgNAlI7RSD5PgopYfjNlyaISLFidUwljkgYjXM5uxChKhVWHh7CoCmoMOc+PHC21PYT4QC1obDFwEERu8QcwxEiEhRcqphJqYm4p2HBkOva7n8oteFs3SXAtrl8lrsuZaIPTmA8kMALs0QkcLkVMMM79EOE1MTMT5Fz86qFFQ2HS2EIACDu8Sic9tIpYfjVgxEiEhRzlTDhKhVGN6jnaeGRORzmpdlAihJVcRAhIgU5c5qGO5BQ4HowtVqHDxXCpUKmDQgsJZlAAYiRKQwd1XDcA8aClSbrvUOuSE5DgkxgVeezmRVIlKUO6phpFTdEPmrjc1NzAJvWQZgIEJEPsCVahjuQUOBLL+oCkcLyhCiVuHWVL3Sw/EILs0QkU9wthpGbtUNkT/ZeC1JdUTPeLRro1V4NJ7BQIQoyPhyQqcz1TByqm58+XsnsmZjAO4tY4mBCFEQCcSETqlVNz+evILFX5+AoTxwvncKbCcvVeDnSxXQhKgwoX9gLssAzBEhChqeSug0mgRk5hVjfVYBMvOKvZqLYTQJMJkExEZoHB77+cGCFkEIwGRW8m3isszo3u2hk/Ae91ecESEKAo4SOlVoSugcn6KXtVSh5AyLtdeWy5XvnciTBEHAhmvLMoHYxMwcZ0SIgoCchE6ppM6wSJ0xkTOzYuu1neHM907kacculiO/qAraUDVu7peg9HA8ijMiREHAmTbq9kidYTGZgBc3OZ4xkTOzYu+1XSH1eyfyhg1HmpZlbu7XAW20gX2p5owIURBwZxt1QPoMy+9XOZ4xkZu74ui1nSX1eyfyNEEQsPHwtWWZAG1iZo6BCFEQENuo28qAUKFpBsJRG3WRK7MH5k3G6htNspuRuXvmQu73TuRph86XoqC0BlFhIRjbt4PSw/E4BiJEQcAdbdTNuTp7IM6Y/CfzjOzcFXfOXDjzvRN5mrjT7viUBIRrQhQejecxECEKEq60UbfkaIZFqrMl1ZKO25V7pTmJdUjXtpLKdaVw5nsn8iSjSWje5C5Q95axFNgZMETUgrNt1C2JMyxPfXwQKqDF0orlx/Z0jYuUdNzb3+U1/z9RF46berXDhiMGh4+LiwpDSVV988f6GC2mDe2CbvFR7KxKPmnfmRJcrqhDTHgobuodr/RwvIKBCFGQcaaNujXiDItltYteF44Fk/vhxU3HYSirtRqUqK4d9+vh3fD+znybx1ljKKvFhiMGRIWFoKreaPUY8fl/eHYsDpy9ypbu5De+urYsM6G/HtrQwF+WARiIEJEL7M2wqNUqmzMmALBgcj8cOHsVt6bq8eGuM5JnUsTyYE2oGqp6Y6vHmOd9hIWqudEd+Y36RhO+Ptq0LHNneieFR+M9DESIyCW2ZljszZjckZaIFzcdb/F5lQoQJE6LCABKqxswK6M31uw71+r5uX8M+aMfT11BaXUD4ttogyqAZiBCRFa5Y6daazMmV6vqMWPVwVYzGWKF7uMjuiEiLBRvf5fr8Pm7xUdi59xx3FGXAsL6rKZlmdvTEoPqPcxAhIhaceceMuYzJkaTgJFLtttcglEB+DrbgNfvSZMUiHSIDndbzguRkqrqGrEl5xKA4FqWAVi+S0QWPLlL78pd+ZL6hkBoCnwcuWpWEUPkz7Yev4SaBiO6totEWmed0sPxKgYiREHKcpO5+kYTdp0qwrwvjsrqdCrF5uxCjFyyHS9uOi7p+KKqOiyY3M/hcS9ukj8WIl8kLsvcmd4JKlXwLMsAXJohCkrWll7Uql/yNGwx73QqdTlEnGGREy5I7Z4qdyxEvqikqh47Tl4BANyRFhxNzMxxRoQoyNhaepEzseCOXXqtMd/3xd07BhP5qq+PFqLRJCC1Uwx6dmij9HC8jjMiREFEbmBgi7t26TVnue+Lu3YMdkf1D5EnfSUuy6QFV5KqiIEIURCRExhYI3YsHdK1LTLzih1e3OXMVsRGarD4rgHNVTnifjaOurPa2zXXndU/RJ5QUFqDvWdKoFIBt6UF53uSgQhREHFlGUMMM+5IS8To176TdHGXs1Pu1eqGFh872s8GsL9rrq3cFLH6h5vdkS8Qd9q9ITkOiboIhUejDOaIEAUROYGBJb0uHE+MSsZ7O1qX4For7TWaBJhMguSdclVoXZHj7I7B9pagXKn+IXI382qZYMUZEaIg4mi5w5rYCA2WPTgY13eLw+jXvrN5cRcDifEpemzJMbRaEnHEVkWO2J11d14xMk8XAWhqYDasu+1KGUdLUM5U/xC528lLFTheWA5NiAq3puqVHo5iGIgQBRF7yx2WxAWPv01NhVqlwpvbTkq6uM/9/DA+P1jg9BitLR9ZBjZvf5drN9eDFTfkD8Qk1dG9OyA2Mkzh0SiHgQhRkBGXO+atPYpSi7wMc7Y2p3PElSAEAL7JNmDb8cswmgQIEFBR04gfc4taHVdYVosnPz6Ifz4wCJMGtuy94K6KGyJPEQQB6w83/a7cmR58vUPMMRAhClJldoKQKG0IbhuYiPd25Ltc6ivX5myDrONnrj6EN01Au2htcxXPkK5tXa64IfKkQ+dLcb6kBpFhIcjol6D0cBTFQIQoyEjpJVJVZ8SKH/O9NibR9d3aYnDXtmjfRgtNiBr5RVVY+dMZu48xCcDMNYdafC7x2mzOezvynaq4IfI0cVlmQn89IsJCFB6NshiIEAUZV3uJeEJCjBaL7ujfIt/DaBKwdMtJp57PUFaL93bk44lRyfjqcGGL71fvZB8RNkYjd2k0mrDxSFMgckeQL8sADESIgo7SCZqhahVempKKRF0ESqrqUFJVj7g2WugiwmA0CQhRq6w2IpNDrOL56nAhfnh2LA6cvepSAMHGaOROO3OLUFRZj7ioMIzsGa/0cBTHQIQoyLg7QfNXgzvJSlB9bERX3D+0CzZnF+LVb39ucXGPjdBgZK94bDpS6HJuiljFsy+/BCN6Of/H/usjhfj9qoOtPs/GaOSsdYeafl9uH5gITQjbefEMEAUZsZeIO6jQdHenjwmH1DmGFT+eweKvc6xuvFda04CNbghCzE3/9358c/SiU4/9+shFzFzdOggB2BiNnFNZ14hvjzUlZN81uLPCo/ENDESIgozYS8QdBACG8jpMG9pF1mNW/Oi9apzqBiOe+u8hrD1wHpl5xVifVYDMvGKHwcPm7EL8ftUhu7sSmzdGI5Lim6OFqG0woXv7KAzsrFN6OD6BSzNEAc5akuXE1EQsl9BLJCxEhXqj45ChW3wkXpzSHy98lYNGCbMDSkwgzP7sSIuP46I0mJreCRkp+lZ5I2JlkVRS826Y8Eprry1j3j24M1Qq/uwBBiJEAc1RkuX4FD3e3p6Ld3fkobre2OrxUoIQoKknyZvbT6HRJCAyLMTqc/makqoGfLDrDD7YdaZV4qncyiIpeTdMeKWC0hrszi8GwCZm5rg0QxSgxN1n7W1QF6JWYea4nogJl7YxnSUVgLhIDV7++jiKKuuR1DYCYSGeu8tTAdDHaPHf396Af9yfjtXTh+GfDwyWvLGeLebnxGgSsPPUFcmPTZTQGE3Kz4IC35eHCiAIwLDucejcNlLp4fgMzogQBSBHu8+ab1C3N78EhnL5ZbJio7Dy2kY0mgSkdoxB9sVySY9Vq5xbnhEAvHBHf4ywKHnURWrw4Pt75D+h2fOqAMxbexTzvjiK0hrby1WWHDVGk/Oz4DJN4BIEobla5q5BTFI1xxkRogAkZ/dZqfkNlrMOukgN1Cqg0STg1lQ9iirrHD6H6tq/6TclS3rNVmOI1GB8SutdSod1b4dEnfTKHWsEAKXVDbKCEFvjMSfnZ0GB62hBGXIvV0IbqsatA4J3p11rGIgQBSA5u89K7Suy7IHBWD19GP5xfzr+PKEPKmsbYRKAqYM64cEbusJQ7jgQiYsKwzsPDcb8SSlY/tBgxEbKW1IprW6wesEWK4G8nQNrazzmuBMwAb8kqd7SX49oJ5dCAxUDEaIAJGf3WbGviK3ZBBWa8iCG9WiH4T3aIVEXgX9sa0pMvTO9I16/Jw3FVY6DEAD4y+R+zYmZE1MTceAv4zEro7esHA9bF+yJqYmYldFL8vO4i6MAgjsBU4PRhA2Hm3rZ3DW4k8Kj8T0MRIgCkNTgQiwfFfuKWB5vuUFczsVyPL5yH+oaTRjXtwNevycNIWqVrIuteS8PAHg6oxcOLBiPBZP7SX4OW7rFR0l6Dndy9L3L+VlQYNpx8gqKq+oR30aLm9jSvRUGIkQBSE5wATTNJrzz0GDoLTqu6nXhzS3ML5fX4vF/7UNFXSOGJsfhnw8Obm5PLeViGxupwZzPDmPait14ek0Wpq3YjZFLtjdX7zwyItnlC7Y3ZxWkBhByfxYUeMRlmTvTOyKULd1bUQmC4LO9icvLy6HT6VBWVoaYmBilh0Pkd+T2rrDVcKu2wYj73tuNw+dL0aN9FNb+fgR0EZoWx58pqsYbW5t2yzX/o6Ky+BgWXwPQHOyIZa7WnsP8OFuMJgEjl2yHoazW5mtGhoUgXBOCkqp6m8/jiNTxmGMfkeBUVtOA61/eivpGEzb+YSRSOwVHN1U5128GIkQBztVunoIg4A+rD2HjkUJEhoXg5SkDcEd6R2zJMbS6sIrJp+bdWvUxWtQ2mmx2cFWhaeZl59xxNnfelXPBlhLMjE/RY+WufLy46bi0k2DB2QCCnVWDz+q95zB/7VH0SYjG5mduCppuqnKu3+wjQhTgQtQqDO/RzunHP72mKQgBgOp6I2Z9moVFG49ZDSzKqhsgAJiV0Qvd4qPQITocJkGw2+PDvHx1eI92zR1fnb1gi8tMlsGM3iJ4iI/WSj4HsRGheHREcvP35GwA4erPgvzPFwcuAACmDu4UNEGIXB4NRHbs2IHXXnsNBw4cQGFhIdatW4cpU6Z48iWJyI2WbvkZXx1u3fXT1uyG2KBrzb7zzTMc67MKJL2WefWJqxdsKcGM1HyS/5vUD4+NTJYUeHDGg8zlXanE/rNXEaJW4a5BrJaxxaOBSFVVFdLS0vDYY4/hrrvu8uRLEZGbXSytwVvbc2U/znKGQ6nyVUfBjJhgay+fRHXtOCnBBHNAyNJn+5tmQ8b0bo8OMSzPtsWj6bu33norXnrpJUydOtWTL0NEbtZoNOGxlftc2iVXnOHw1fJVe9UsIgHA7/5zwGGvEO4lQ5YajSZ8cbApELnnuiSFR+PbfKqOqK6uDuXl5S3+EZH3vf6/kzhhqHDpOcQZDl8uX7VVtpyoC8fSe9PQs0MbGMpr8eR/DqCu0fqOwo72kgGa9pIxuhLVkd/ZceoKrlTUoV1UGMb17aD0cHyaTyWrLl68GIsWLVJ6GERBLTOvGO/uyHP68WIVjPkMh9QEUiXYyicBmgKJ/1uXjYPnSrHgy2wsuXtgq4RDOXvJMFE1eHy6r2k2ZMqgTggL9al7fp/jU4HI/PnzMXv27OaPy8vLkZTEKS0ibymvbcCcT7MgCMC91yXhx1NX7OZQ2GJthsPVahhPsswnsZbv8en+C1ABWPKrtBaP5V4yZKm4sg5bj18C0PR7RPb5VCCi1Wqh1UovqSMKZEpUYCxcfwwXy2rRJS4SC29PwY+nrjT35JAiNlKDV+4aYHOGwx/KV8V8D2vB1yf7L6BT2wj88ebezZ/jXjJkad2hAjSaBKR11qGPPlrp4fg8nwpEiKiJEhUYG49cxLpDBVCrgKX3pSNKG9q8pPLcuqMoqbJesmtu2bTBGNHLf/fSsJfvIXpj6ynce12X5rwSR9U31paqKHAJgoDPr/UO+RVnQyTx6MJVZWUlsrKykJWVBQDIz89HVlYWzp0758mXJfJrSlRgXKmow/+tywYAzBjbE0O6tm3+2sTUROyen4G4qDCbjzffodefOcr3AACTADy2ch8ajSYAvp2MS953tKAMJwwV0IaqcUdaR6WH4xc8Gojs378fgwYNwqBBgwAAs2fPxqBBg/D888978mWJ/JZSFRgvfHUMZTUN6N8xBn+8uVerr4eFqvG3qalQIbAvtlLzOHIKy7H02r46gLRNAyk4fLr/PABgYqoeugiNwqPxDx5dmhkzZgx8eCsbIp+jRAXG/44ZsOlo0w64S+4e2LyjriVfrnxxFzl5HMu+y8N1XeMw9lpppi8n45J31DYYsT7rIgDgniFclpGKOSJEPsTbFRjltQ1YsL5pSeaJUd0d7gwa6BdbqfkeN/ftgI/3nMOsT7Ow6Y83oVNsBAD/SMYlz/n2mAEVtY3oFBuBG/k+kIzFzUQ+xNsVGIu/PoFL5XVIjo/C02ZLMkaTgMy8YqzPKkBmXnGLpSDxYntneicM79EuYIIQQHq+x3OT+6F7fBRKqxvw+Mp9qG80tXoue+eQAtOavU3LMncP6Qx1AP1eeBpnRIh8iDcrMHafLsbqvU2J46/cNQDhmhAA3DPF0RIUANz8/35o/toJQwUGv7gFr98zsPn8BPs5DEanr1Qi83Qx1Crgvuu5LCMHAxEiHyLekT/18UGogBbBiNSkUCn9R+obTfjLl01LMg/c0AU3dG+aRrbVQ0Os2AmkxEvL8zSka1scOHu1+eMfnh3b4uOhyXHYkmOwen4q6xrx5McHsfyhwQAQNOcwWEj5nRKD+jF9OjQv1ZE0DESIfIwrSaFS78Q/2pWP3MuVaBcVhrkT+wJwXLGjQlPFzvgUfas/wko0X3OFtfOkVqHFJn/iebszvWn7dik9Rp5fnw21Su3UOSTfJOV3qq7R2Nw75IGhXRQZpz9jIELkg5xJCnU0m7HsgUFoG6XFyUvl+PuWptLTebf2bS4xdLZix9+WIWydJ8sUDssZDCk9Ri5X1Nv9Oved8S9SZwg3ZxtwtboBibpwjOnTXpGx+jMGIkQ+Sk4FhpT+IzNXH2pxsdWEqBClDWn+2JmKHX9bypEyqyGynMFw514x3HfG98mZIVy1p2lZ5t7rkhBqo/ydbOMZIwoAUjuCmmswCpjx30PNnVrlVuwo1XzNFVLOkznzGQx37hXDfWeU56iqSeoM4bpDBdiTXwK1Crh/KJNUncEZEaIA4ModtnhXJ7diR4nma65y9jxdrqjFbQM72j0/QNMsU4PRduDFfWd8g73lRHFJ9BuJWylsPNLUwGxc3w5I1DFJ1RmcESEKAM7eYZsHC3L3TPF28zV3cPY8dYgOd3h+VGg6P9pQ239WBQALJvdjoqqC7O3l9OTHBzHkpS2YtmI3/p15VtLz7T9zFUBT9Rk5h4EIkR8Tp5cNZTWIiwprdYGUSgwW5OyZ4u3ma+4wNDkOcVHS9/8QN/MTZzBsnZ+EGC2eyeiF6HANHh7eze5zvrjpuEc2LiTHpCwnllY73mUaaHpvxEZoUFnXiI66cIzu3cFdwww6XJoh8lPWppedZR4sSK3Y8WbzNXcJUaswNb0TPth1RvJjLPu2WJ6fM0XVWL33HJZuPdV8jDZUjTor3VYB303kDQZyc4RsEd8N8dFalNY04L7ru3CWywWcESHyQ7aml62x9+fR8o5fJKWNu9ylHF+RkaKXdFy7qDCbwYJ4frSharyx9SQM5S1/DraCEMB3E3mDgbuWCfW6cCy4rR9yL1eyk6obMBAh8jNSSlDjojRYel86Vk8fhkFdYq0e445gQc5Sjq8QZ3LsfcdxURpkzr/Z7vjllAJbMs/NIe9xdZnwN8O7YvX0Ydg5dxzyrlQBAG5J0bd6/5M8XJoh8jNSppdLqhqgj2lKsDx4rhQqANERoSivaWw+RkqnVin8bUdeKW30/zZ1AMLsJJ0C7pnm96VE3mDgaDnRkXZRWgzv0Q5lNQ1Ye7AAAPDwjd3cOsZgxECEyM9IvXj971ghPr7WaEkAUF7TiLgoDaamd0LGtXJddwULcpqv+QJX2uiL3BFEeDuR199a8bubvSBUijX7zmHmuJ74/MAF1DQY0SchGsO6+04OlL9iIELkZ6RevD76qXX54dWqBny46wyuD7ILkDWuzuS4EkQokcjrb634PcVWEBqlDUFVndHuYwvLarH7dDH+k3kGAPCbG7tCpQru3yN3YCBC5GekTC/butvjxmstuTKTI3Wa39bPwpuJvP7Wit/TrAWhhvJazPoky+Fjfzh5BWeKqxETHoqpgzp5frBBgMmqRH5GSrWKvQsjEyXdw/znYMttAxORENNy5iQyLMSrF35/bMXvDZaVYfoYaTNc4u/NvdclITKM9/LuwECEyA/Zq1Z5eLi0Do9MlHQPXaTtBmkbjxQCEDAroxeeGt0dAFBdb0TbyDAvjU5eK/5gJqWaqn2bMGSdL4VKBfx6eFevjS3QMRAh8lMTUxOxc+44rJ4+DP+4P725rNAkeD7HgX5Z7nDUifNSeR3e2HoKaUmxuP9av4k/f3EENfX28xHcxR9b8XubmMR7a6re7mxiwrXAf2yfDujaLso7gwsCnFci8mOWOQ6VdY3YdMRx+/DYCA1MggCjSQj6PBFnyOkhYp6X883TN+GHk1dwtrgar337M553sLTjDv7Yit+b5HQozi4oB8CSXXfjjAiRj3O0Xbm5D3fmo6S6Hh2itQBsd1UtrWnAg+/vwcgl27nviRPk9hARlz+OF1Zg8V0DAAAf/ZTvleUQR0sOtrrrBgM5HYpFIWoVbuzuP6Xq/oCBCJEP25xdiJFLtmPait14ek0Wpq3YbTN4uFpVjxU7TgMAFtyWgn8+MBhto+znIohVEwxG5HF2GeNyRS3G9OmAe6/rDEEA/vz54VZLNHICTyn8tRW/pznbGddoErD/7FWPjClYMRAh8lH2tiu3Fjws/yEPFXWNSEmMQYgaeHFTDkqq6u2+RjBXTbjC2WUM8XH/NzkF+phwnCmuxuv/+7n563ICTzlsJTe3jdLgsRHdoIsIC7qfvyudcYM5n8YTmCNC5IMclVxa9gK5UlGHf11rsjS2b3vM+O8hyXd65lUT/tQdVUlyW4VbNjDTRWjw8tRUPP6v/fhgZz5KquqgVqnwxbW24eYKy2rx5McH8c8HBmHSwI5Oj9m8d8aWHAO+zLqIkqp6fLDrDD7YdSbompu5EkwEaz6Np3BGhMgHyS25fG9HHmobTBjYWYdVe845tY8G7/KkESss+iVGSw5CgJbLH5uzCzHns8PNx6w7dNFqEGJu5upD+FpCIrI9IWoVymrq8dGuM61my4Jpmc5oElBUUefUY4M1n8aTOCNC5IPklFwWVdbhP7ub2rmfvlKJSgdtqm3hXZ5jciosRJb712zOLsSTHx+U/domAfj9qoNYrna+GZrcmbZA5MzPUKRCcObTeBoDESIfJKfkcsWO06htMAGAU0GIEvue+CNbbdKtiYvSYMFt/aGPabl/jdEk4IWvjrk0jufWHUVNg6nVc0shZ6YtEJfp5PwMLcW3CcNLU1KDZunKmxiIEPkgRzkIYvDQo30UHlvZenM7qYK5akIOuRUWJVUN0MeEt7qY780vgaHcuSUB8+cW90SRm9cRbM3NzHcbjo/S4oWvpP8MHx/ZFedLavG/nEtIT4rFF0/dyN8RD2EgQuSD7G1Xbh48fLjrDGoanO/QKWfb+2DmTIWFtYu5uy/wcjetC6bmZq4swQDAhzvPIjSk6bft2Ql9GIR4EAMRIh9la7tyMXgYmtwOsz89bOcZrJuV0Qvd4qNkb3sfzJwJIKxdzN19gZeb1yF1ps3fl+lcWYIRCQAajAK6tYvEjQG4TOVLGIgQ+TBr25WLwcOSzSdQXW9EcnwU8ouqHD5Xu6gwvDyVa9zOkBNAmF/MzZcGOkSHY0jXttDHaF1enjEnJ69D6kybPwenzjYqs6WspgEmAQjx31Pi8xiIEPk4y/1kjCYB/ztmwIc78wEAc8b3xvNfHbPbvCwuSoPM+TcjLJQV+86Q2jfE/GK+JcfQajYrUReOO9M74t0d+W4fo9RZG0czbf4eqLrSqMyaq9UNAZu86ysYiBD5EWvr3n9Ybbt5mXhh/NvUAQxCXGBvJsGceDEHYHVpwFBWi/d25ON3o5Lxyf4LDnfulUPOrI29mTZfYTmbJHV8nki0DZTkXV/FQITIT9ha97Z3hx4od7m+wNZMQlyUBlPTOyEjRd+cWzFyyXa7vTq+OlyIvc9lYF9+CTJPF8EkAG0jwxAfrYU+Jhxr9p3D+qyLSNSF45mbe2HJtz/jalW9W/M6LGfafIm1gFtqhZAnEm0DIXnXlzEQIfIDzqx7x0Vp8MOzYzkT4kZSZhIy84ol9eo4cPYqRvSKx4he8a2OSUmMwe7TTc9z4lIFHh7eDUu3nmx1XKDkdZizFXBLrRASl9GkLs/Ym+FSAUiI0cIkCFifVeCTM0eBgH+hiPyAM+veJVUNOMBdQt1OnEm4M70Thvdo1+qi5I5eHbpIDV65eyAA4KNdZ6wGIUDTTIjU0l1/4KjzK+B4g0bz3YbtUaFpluXNaYNsHiMAqG004cH397h1E0JqiYEIkR9wZdt58i539eqoc9AfZlZGL+ycOy5gghBAeufXlbvysT6rAJl5xVaDkompiVj+0GDERmqsPo/5TNLtaR2x/KHB0MdoWxwjPtYyjyeY9uTxFgYiRH7A1W3nyXvEpQFbk/finbi9nA5xZsAWFYA1+863OD4zr9juxdkfSA2cX9x03OoMhfl50EWEYe9zGZiV0RuxES0DEsuZpImpifjtTd0BNAUgKx+5HuGhIVZfW+rMDEnHHBEiP3C1qs7uWralQGlM5Y/c0atDzp4wZTX1Tid2mnO2SsWdnAmcxRmKJ0Yl46vDhVbPw4EF4+1+b9X1jVj+Qx4AYO7EvtBqQmAoD949ebyNgQiRDzK/KJwpqraZI2BNICYw+htXe3VInRnYkmPAR7vOOJ3YKXKlSsWdpPZrMSceZ603i9Tz8K+fzqKosh5d4iLxqyGd8fVRacsuXPp0DwYiRD7G1T0yWLLrGywrbOKjtIAKKKqsQ2Zesd0ZB6kzA19mXbSb2PnCV8cctn53tUrFnaT2a5FKSgv8itoGvLujaTbk6Zt7QROilnz+46O0jg8ihxiIEPkQZ/fIWDC5H+KjtSwv9DFihc3m7EL86fPDkmccpMwMxEVp7HbTBQBDeR3e3p6LpzN6Wf26oyoVOfvYuIut2SRnOVpG+XDnGZRWN6BH+yhMGdQJgPSZmTmfHcYLdzDodxWTVYl8hCt7ZMRHa22Wk5KyxODS8qJqr/rCvATV1k8zJTFa0usv3XrSZoWHnFwUb5qYmoidc8dh9fRh+Mf96VgwuZ/Lz2ltGaWkqh7v/3gaAPBMRu/m3x0p5x8ALpWzgsYdGIgQeZmtCgexgZUzWB3jm1zpiyHODOh1LX+2Idf+au/MlR4c2HoNd/Q88RTzfi2PjEi2W4kkhbXfkTe3nUJFXSNSEmMweUDLWQ3x/CfE2F5+YQWNe3BphsiLbCUF3pGWiE/2XXDqOeOiNEhPikVmXrHP7hsSrOTMOFhubLg3vwR1jSbcd10S/r37bPMyjNEkfxy2libc1fPE01zNHWkXFdaqgiy/qAof7z4LAPi/yf2gtvL7MjE1EdHhGjz4/h6bz80KGtcxECHyElv5H4VltS7txlpS1YD+CzfD/IZMiYoHas2ZGQdXk5Vt2ZV7pVWAerWqzuHjHPU88RZXckfuTO/YKjBf8s0JNJoEjO3THiN6tm6zLyqqdHyOAFbQuIKBCJEXuJL/IYXlrLASFQ/UmtwZB2eTlaV4+7s8fHGwoDlANZoEvLjpuMPHLZjsO2XglpVIRRV1kr6H8Sl6AL/MNGWeLsLmYwaoVcD8SfbzT/xl1sifMRAh8gJn9oqxpo02FJV1jQ6PU6rigVpyVH1h3njO08Eq0DJA1UWESXpPto0K89h4nGmiZr5rsNEk4P2d+XarW8QZHWszTeGaEJy+UoneCbYTf+X8DMk5TFYl8jCjScCu3CKXniM2UoNZGb0kBSEipSoe6Bf2qi8sG8+5K1i1xzy50l7nUHOeWnLYnF2IkUu2Y9qK3U5vKOfo/KrQdH635BisVi5V1xsdVr3I+RmScxiIEHmQ+Mf27e9yXXqeZdMGo1t8lFOP5dq1smxVv1jud+Ktn5MYoJZIzH3wxJKDMyXNtjg6v+NT9A5nmhxVvUj9GZJzuDRD5AQpU8ruWu9P1IVjWI92Ts9scO1aeZa5DdbeM878nNSq1vlBUsVFhTm95ODKvjSeaKJm7/xm5tkvi5da9SLlZ0jOYSBCJJOUfTncud4vTvvK3YeDa9e+xTy3wRo5P1/x0vf2tEFoG6VFdkEp/vb1CVnvN70uwqnN+ay9/2MjNHh0RDfMHNfL4YXZ2ZJmR2ydX3f2SnH0MyTncGmGSAapU8ruWO+PDg/FcrNpX6ndHs2/zrVr/yHn5ysuCUwa2BHDe7TD9FE98MSo7pJfS0zglLvkYOv9X1rTgKVbT2HIS1scLqt4u4kaq158H2dEiCSSM6Xs6h/RsBA1djw7tlXFgq1eCpZT9Nz4zj/Z3LU3RotpQ7ugW3yUzSWB+ZP6IfdyBbaduOLwdcxLcqUuOUiZ5SutbnBYNu7pwMBy2WhI17bQRYSirMZ6ojdnDpXHQIRIIjlTyvFtXNuVc3iPdjbLJq1dOIZ0bYsDZ69y7ToAuJKL8MiIZEmBiOV7S8qSg9RZPgH2czw8WQ5rbdmoQ7QWNfXW29Fy5tA3MBAhkkjqLMeWHAO+PuraJlh3D+7U4mNryYGWFw6uXQcOZ3MRHO3GK3Jmxk7OYwrLarFyVz4eGZHc6gJvr127K4GBreTwyxVN1UHto8MQolLBUP5LtRBnDn0DAxEiiaROFX+464xLrxOuUWPywI7NH0tJjiUCPLvsIfcxL246jvd35lt9n9pcgnLyfS1l2UgQgB1zx3Hm0AcxECGSSEpVg5RyythIDUqrG2xu3jU1/ZfZEFt3eWzhTta4a9nD2gyc+NxykrDtvU/l5KY4OkbKslFRZT0OnL3KmUMf5JVAZNmyZXjttddgMBiQlpaGt956C0OHDvXGSxO5jaMpZQHSejosmzYYFXUNNjfvWr3vPL4/eQULJvfDi5uOu7XfAgU28/eoLY6WPezNwC28PQVP2nluS47ep46WoKTOBnq7Eofcy+Plu5988glmz56NhQsX4uDBg0hLS8OECRNw+fJlT780kdvZK3d8fEQ3Sc9RVFWHiamJ2Dl3HGZl9LZ6jKGsFr9fdUhyciyRSHyPJlq8R0PUKvz93jS7M2iOytMBYPlDgxEbqZE8Hmffp3K6r7JE1795fEbk73//O6ZPn45HH30UALB8+XJs2rQJH374IebNm+fplydyO1tTynvzS/CBhPwQ8z+Ga/ads3qMnMZUvMsjS+bv0dzLFVi69RRKquqx9lABJg/s2Ly3jfn7F4Ck8vSdc8dhfIoeb2/PxUe78lFa0yBpTJcram0us1grubU3FgB4bt1RjOubgLBQNTem83MeDUTq6+tx4MABzJ8/v/lzarUaGRkZyMzMbHV8XV0d6up+yWguLy/35PCInGZtSlnuH0N3bXLGuzyyRnyPDu/RDoO6tMW972bix1NF+PUHe3C2uKpF9UiiLhz3X58kq+Pp0xm9MHNcT6zclY8XNx13OJ4zRdUYuWR7q2WWO9IS8dXhwhafj4vSoKTKfoBTUtWAYYu34W9TUzExNdHmshFLdH2fR5dmioqKYDQakZCQ0OLzCQkJMBgMrY5fvHgxdDpd87+kpCRPDo/IreTu0unqTIYKv3TIJLIntZMOyx4cDLUK2JNf0iIIAZqWO5ZuPSXpuczftyFqFR4ZkYxEXbjNbrAqNCVov7H1ZKtAp7CsFu/uyG/1eUdByC/H1Tcv00zor8cNVn4XuDGd7/Opqpn58+dj9uzZzR+Xl5czGCG/Iqcssb2Mpmfu7LdAwWlUr/aIDrfeYVTOUqDlDJyUJG65ryHXog05KCitxZ78EoSoVXju1r6Ij9ayRNdPeDQQiY+PR0hICC5dutTi85cuXYJer291vFarhVbrWkdKIqVJLUt09MdRXMpZMDkFL25yT78FCl5780tstjmXwl6ehb0A/P7rkyTPtjhDXDJ6eVMOAGD+rX3x+E3S990h5Xk0EAkLC8OQIUOwbds2TJkyBQBgMpmwbds2zJw505MvTaQoKZ0x39+Z3/x/ezMeE1MTMSGV249TEyl9NayRsxTozAycrQB845GLkl/XFSYBuDO9Ix4fmeyV1yP38fjSzOzZs/Hwww/juuuuw9ChQ/HGG2+gqqqquYqGKBjlF1Vh6/GmmcIXbk/BuztO253x4PbjBLjWZVdqUvMzN/fCv3efbdEuXuoMnLX3qavJ1NHhoaiolTaTM7Z3e6hU1gMlZwM48jyPByL33Xcfrly5gueffx4GgwHp6enYvHlzqwRWomDywc7TEATg5r4d8MiIZPx6eDf+kQwizlwUXe2yK6UzcGykBmv2nW8RhERpQ3DfdZ0xPqX1croUUl7XGnEpaPucMRixZLukfXRmf3YYYaEhmDSw5XngNgm+TSUIgidziFxSXl4OnU6HsrIyxMTEKD0cIrcoqarH8MXbUNdowurpwzjTEWScuSgaTUKr0ldz4kV759xxDrumio3J5P7hj43U4JW7Bjh14Zb7uuJ3IAZXm7MLZXV0XW4WlNkK4Cxfg9xLzvXb451Viailj3efRV2jCamdYjCsO0tvg4mcbqHmHPWcERM2V+7Kh9HOPgO2OgNHahxfCkqrG+yO0Z6JqYlY9sBgtI0Ka/H5RF04fjcquVUXWMuS24mpiZI7FwNNVTRGk2B3MzzB4lhSjk+V7xIFutoGI/6deQYAMP2m7jbXsynwOLoo2tuTRWqiqb0db0XWkkqNJhMe+mCvw+cX7IzRns3ZhXhxU06L5ZW4KA0WTO6HSQM74s8T+zlcqspI0UvqXAy0bCkvp0kbKYMzIkRetD6rAEWV9eioC8ekAZwODiZSZzWs7ckiJ+HT0ewK8EtS6Z3pnTC8RzsUS8i/EMndN8bWLNDVqgbMWHUIm7MLW43HWpDTo30UNCHSg5/LFbXcDM9PMBAh8hKTScCKH5tKdh8dkQxNCH/9gokrF0Ux4VPKZdiZJQe5lS1Svxd3LY2cvlKJe9/NRINR+hJKh+hwbobnJ/iXkMhLfjh5BbmXK9FGG4r7hrJjcLBx5aJob/sAa+TueCsGOlJJ/V52ny6WNAu0O6+4+XNGk4DMvGKszypAZl4xNh6+iDve3oUzxdWIiwrDr4d1hb1VIfOtDxwFcNwmwTcwR4TIS1b8eBoAMG1oEmLCpW+jToHB1R1ibXUvtcfejrfmxEBHSmWK1Av35uxCzPviqKRxzlh1EK/cPQAA7H5/JVX1+M/us4iNCEWplS6x1hqv2Ws/b3ksKYPlu0RekF1Qhtve2okQtQo7/jwWnWIjlB4SKcBWGaucUlKjSZC84+2sjN5Ys++c5FLhzdmFmLf2KEqrbW86N65vB/zj/nRE2wmmbZXMuosYVMRGalqM1db3xj4i3ifn+s1AhMgLnllzCF9mXcQdaR3x5rRBSg+HFGTrorhgcj+0jdJKanIm9hWxN7uii9SgrLpBdv8Mo0nA29tz8dGufJTW/HKRD9eoUdtgAtAUAEy/qTsevrEb2mhDWz3eXs8Td1EBSIjR4v/dm46iyjpJ54xNA72HgQiRFzn6A3extAajXv0OjSYBG2aOxIDOOgVHS77A8j1ztaq+1caGju7Y7c2uWJstgMUxjhqgWXtf/3jqCv66IQeni6oAAG20oZg2tAseuKELkuOjAACZecWYtmK3rPPhCjYF9E1yrt/MESFygZQp35U/nUGjScCw7nEMQghAyz1ZNmcXYsYq+a3bXdnxVkr/DGv7xtQ2GFFd/0tuRmVdI1b8eBorfjyNgZ11yOiX4PXmYCy99X8MRIicJGXvjxE947F6zzkATQ3MiMy50uQMcH3HWzkXcUd5H0culOHIhTLJz+cutip4uBTjPxiIEDlB6gXksRHJqKhrRI/2URjbp4OXR0m+Tk6TM/PZCUcXWXf2zzCaBOzOK8a8L47aTT6NCAtBTb1R0uu6g70qIyan+hcGIkROkHoBeXdHHgDgtzd1h1rC3Rjv4oKLM03OpFxkXS0Vtvdatng7CAGsl966uksxeV/QBiJGkwAVIOniQMHJXlAg9QJSVFmPdlFhmDqok8NjeRcXfOTOXEi9yIp9QVzpn+GJEtyosBAAQJWdoEWtAszTTGIjm8qEzRNv9TZ+L1xd6iJlBGUgUlhWg6fXZGF8vwRMH8V1e2rNUVAgpyX0b4Z3Q7gmxOHr8S4u+MiZuZB7kbWXzOoouLX3Wq6oqjfid6OS8d6OfJtlxW9PG9SqjBmApJlCZ5e6SFlBGYj88PMV7M0vwcGzVzE0OQ5pSbFKD4l8iJSgYHyK3u4FRKQNVeOhYV3svh7v4oKXnJmLzDxp7dLNL7K2klkdvY8cXdBd8dXhQix7YHCrcmVHAZKUwIGb3PmnoAxE7rs+CTtOXcHXRw34w+pD2PTHkXa7BFLwkBMU2LuAiB//akhntGujtfuavIsLblJnLpy9yForw5X7HPZI6VtirrCsFm2jwrBz7jibAZKzuVLc5M4/BWUgolKpsPiugTh8vgznSqrx3LpsvHl/OlQq3m0GOzlBga0LSLsoDYqqGqACMKRLWxhNgt0/oryLIykzF968yMp5DjFgAiBprxqg6b1sK0ByJVfKXUm65F1Bu/uuLkKDN6cNQohahQ2HL+Kz/ReUHhL5ALlBwcTUROycOw7/ffwGzBzbA7em6lFR15SIJwCY/dlhjFyyHZuzCwG03lnUaBJ4F0cAfpm5uDO9E4b3aNcqePXmTrKOXgsAtKEq/Oexodg5dxwmpiZiYmoiZmX0kvT8tt7L4rKo5c2AuCwq/h7ZYm+XYm5y57uCNhABgCFd22LOLb0BAM9/lY3cyxUKj4iU5kxQsCXHgD99fhhvf5eHb7INqGs0tThW/CO6+OscjFyyHdNW7MbTa7IwbcVujFyyHVer6rhVOTnkzYus+Fr28p/qGgX8Yc0hbMkxNH9u5rhe0MfY/h2y9152tCwKNC2LOurcKs5U6nUtx6HXhTPp20cFdSACAE+O6oGRPeNR22DCzFWHUNfovVp48j1y7zpt3cGZE679e3dHvtU7vRmrDuGOtMTm57d8PYB3cdTEmxfZ8Sn65tJZW0qrG1rMVISoVXjhjhSoIP+9LGdZ1BFxpnL19GH4x/3pWD19WPPMDfmeoA9E1GoV/n5fGtpFheGEoQJLt9jen4ECn5y7TneUOIqPFSsJeBdHjnjrIrs3v0RS8qmAljMVzgZL7s6VcrTURb4jKJNVLXWIDsff7hqA3/3nAN7dkYeb+3XA9d04DR6spFYxuKvEUbzTc1RJQCRyphJGLjnJ0e4oG2auVPBiIHLNhP56/GpIZ3x+4AJmf5qFb54ehTZanp5gYVkuOD5F3/yH1FBWg5KqesS10UIXEdZcBePuKhZ7lQRE3ib3gu9q2bCzFS/cFsH/8UprZuHtKcjMK8b5khq8tDEHr9w9UOkhkRfYKxcEgFe//dnq19x9Z8Y7PfIlYmAgddbP1fevM23puS1CYAj6HBFz0eEa/L9706BSAWv2ncfWnEtKD4k8zF654JMfH8STdkoJxWoXV7EqhnyReb6UPe58/8rJL3G11Jd8h0oQBHdvJ+A25eXl0Ol0KCsrQ0xMjNde9+VNOVjxYz7i24Th22dGOeyMSf7JaBIwcsl2p/I8xGniBZNT8PtV9ps4JerCcUdaIt7bkQ/A+p0eE1LJV23OLsS8tUetJq566v3raLnF0e+u+Pu5c+44LtMoRM71mzMiVsy5pQ/6JESjqLIeC786pvRwyENcSTY1TzBN66yzeVzbyFDcd11npHTU4ZmM3kiIYVUM+ZeJqYk48JfxmJXRG7ERLct5PfX+dVTx4s5SX1Iec0SsCNeE4PV70jDln7uw8UghbhtowMRUvdLDIjdzR7Lp0YJSHCkos/n1q9WNeGNbbvPH+hgtZmX0Qrf4KCbWkd8IUavwdEYvzBzX02OJoXKSTrktQmBhIGLDgM46/G5Ud/zz+zz85ctsDOseh9jIMKWHRW7kjuTQ3adLIAhAqFqFRgcdHwHgUnkd3th6Cu88NJjVMeR3PFXVJTfplKW+gYVLM3b88eZe6NmhDYoq6/DXjTlKD4fcTMp+GraoAMREhGL7icsAICkIAeS1qiYKBs4knXpz3x3yPAYidoRrQvDqrwZCpQLWHizAd9cuOhQYpHRRtfU1AUB5TaNTr8v1a/IV1jZh9PbrO7O/DDe3CywMRBwY3KUtHh+RDACYv/Yoymsdtzwm/2GvXHD5Q4Ox3MrXEmK00EW4vqrJ9WtS0ubsQqubMHqz7NWVpFNubhc4mCMiwZxb+mDL8Us4W1yNxV8fx+K72OgskDhqR235NZMg4MH397j8uly/Jm8yTwY9U1SNN7aebDUTIS6HOHMhd6bDqatJp860kiffw0BEgoiwECy5eyDuf283Vu89j8kDOmJkr3ilh0VuZC8Jz/JrXxw479Jr2WpVTeQp1pJBrRHQ9P5ctCEH41P0ki/oznY4dUfSKbdF8H9cmpFoWPd2+M3wrgCAeWuPoLreufwA8n8/X6p0+rFcvyZvs5UMaovcHCZXOpwy6ZQABiKy/HliX3TUhePC1Rq8sfWU0sMhBRhNgqTW/7ERGjxzcy/o2cCMFGQvGdQRKcsmUpNN6xtNVpNimXRKAJdmZGmjDcWLU1Lx+L/24/0fT+OOtI5I7WS7qyYFni05BpwuqkJkWAiq6402N+d65e4BmJiaiD/c3Ivr16QYV7oHS1k2kZpsOmzxNpRU1Td/3nzZRkw6tVza0XPzuqDBQESmm/slYPLARGw6Uoj5a49i3e9vRGgIJ5aCQaPRhNe+/RkAcEuKHrekJODFTfb/eDpav+YW5uRJzlRmyclhkvr85kEI0DoplkmnwY2BiBMW3p6CHSev4GhBGVb+dAa/vam70kMiD9ucXYjn1h1FSVVT+faXWQXYk1+MBZP7oW2U1qk/ntzCnDxNbmWW3OUQZyu/rCXFMuk0ePFW3gkdosPx3KR+AID/97+TOF9SrfCIyJPEZDwxCBEZymoxY9UhlNXU29ycy9Fzcgtz8iS53YPl5jC50p2Yjf1IxEDESfddl4Sh3eJQ02DEgvXZEAS26w4EYqfJdQcv4IMfT+OLAxfw3Lps2Z0fHb2GM90kieSSkgw6K6MX/nF/OlZPH4adc8fJmo2z9/xSsbEfcWnGSWq1Cn+7awAm/eNHfP/zFWw4Uog70joqPSyyw1E+htReC+bM7+qkTivL6SbJqWpylaeTQW09f1yUptUsojVs7EcMRFzQs0MbzBjbE0u3nsRfNxzDqF7x3KHXRznKxxCXSpydg5Ba6rg3vwTfSFx24Z0iuYunk0GtPf+Qrm0x+rXvYCirtfp7xcZ+JGIg4qInx3THhiMXkXu5Eou/PoElv2L7d19jK8gQ8zGWPTAIL2467nQQAji+q3NmtoV3iuROnk4Gtfb8C29PwVMfH7RZ5s4eIQQwR8Rl2tAQLL5rAADgk/3nkZlXrPCIyJyUfIy/rM92uteClM6PcjtbspskBQpuTEdScEbEDa7vFocHbuiCVXvO4bl1R/HN0zchXBOi9LAI0vIxpKxjWyPlrk5uZ0veKVKgYY8QcoQzIm4yd2JfdIjWIr+oCm9vz1V6OHSNJ/MspNzVye1syTtFCkTiso3cMncKDpwRcRNdhAaL7uiPp/57EMt/yMPtaR3RRx+t9LCCntQ8i7iosFbdH60fp8GC2/pDH+Pebc5/M7wrbk1N5J0iEQUdzoi40cRUPTL6JaDRJGD+2iMwsQ+E4qTu7vnSnakO+yCoAPxt6gBMHST9rk5qIHRraiLvFIkoKDEQcSOVSoW/3tkfUWEhOHiuFP/dc1bpIQU9qbt7ThrYlFSXqLMeOCSaLZmITc8sdxK1htucExHZpxJ8uCVoeXk5dDodysrKEBMTo/RwJPvXT2ew8KtjaKMNxZbZo5Coi1B6SEFP6r4uRpOAb48Z8MfVh9BoEvD7MT1wU6/2zUsmzuwPI1bNANZLGJkTQkSBRs71m4GIBxhNAu5+5ydknS/FLSkJeO831yk9JIL0nW5f3JiDD3bmY0jXtvj8yeFQqZqOsdWPREpAwQ3uiCiYMBDxAScM5bjtzZ1oNAlYzjtev3Glog43vbodtQ0m/OuxoRjduz2ApiBm5JLtNitgxC6RO+eOs1vKyxJGIgoGcq7fzBHxkL76GPxudHcAwPPrj6G81rleFeRd73yfh9oGE9KSYjGqV3zz5+XsD2MLSxiJiFpjIOJBfxjXC8nxUbhcUYdXN59QejgBR07SqBSFZTX4+FqC8ZzxvZuXZADpZbjcH4aCjbt/Dyn4sI+IB4VrQvDy1FQ8sGIPPt59DlPSO+G6bqyOcAdP5Fy8tT0X9Y0mDO0Wh5vMZkMA6WW43B+Ggglzn8gdOCPiYTf2iMe913UGAMxbexR1jUaFR+T/bO3dIm5it1ni7rbAL3dzH+w8jU/2ngcA/GlCnxazIQDLcIksufP3kIIbAxEveG5SP8S3CUPu5Uos//600sPxKXKndaVsYrdoQ46k6eHN2YUYuWQ7pq3YjRc3HodREKANVaOkqq7VsVL7kTDvg4KBO38PiRiIeEFsZBiev70/AGDZd7nIvVyp8Ih8g3kg8PSaLExbsRsjl2y3eyfljqRR8bWt3c3VNZps3s1xJ1GiJu76PSQCPJgj8vLLL2PTpk3IyspCWFgYSktLPfVSfuH2gYlYe/ACvv/5Cp5bexRrnhgGdRDfPdvqySFO69q6sLsjaVTKjriLNuRgfIq+1QwHdxIlYvI2uZfHZkTq6+txzz334KmnnvLUS/gVlUqFl6akIkITgr1nSvDJ/vNKD0kxrkzrSk0GjY/S2vyaq3dzLMOlYMfkbXInjwUiixYtwqxZszBgwABPvYTf6dw2EnNu6Q0A+NvXx3G5PDjvFqQGArvzilvljzhKGhXN+eyw1eUVo0nArtwiSePk3RyRdUzeJndijoiXPToiGQM761BR24hFG3KUHo4ipF7gZ6w62Cp/ZEuOwWbSqLlL5a0z98WclLe/y5X0+rybI7KOydvkTj4ViNTV1aG8vLzFv0ATolZh8V0DEKJWYdPRQmzNuaT0kLxO6gW+tKZlN1oxfwRo2tclIcb28otw7d+8L45iV24Rvj5iPTnVGt7NETnG5G1yF1nJqvPmzcOSJUvsHnP8+HH07dvXqcEsXrwYixYtcuqx/qR/Rx1+e1My3v3hNBasz8YN3eMQHa5RelheMzQ5DrGRGpRWy2t7L6ApSFi0IQc7545DdLgGD76/x+5jSmsa8OD7e6BWwW5yquXr8G6OyDEmb5M7yApE5syZg0ceecTuMd27d3d6MPPnz8fs2bObPy4vL0dSUpLTz+fLnrm5N745asC5kmos/uYE/jY1eHJptuQYZAchIvNE0qLK1v0+bJHTziA2UoPxKXr5gyMKQmLyNpGzZAUi7du3R/v27T01Fmi1Wmi1tqfbA0lEWAiW3D0Q01bsxqo95zApNREjLdqKByKjScALX7meGyPefXlCaXUD9uaX8I8rEZEXeCxH5Ny5c8jKysK5c+dgNBqRlZWFrKwsVFaymZdoeI92+M3wrgCAuV8cQWVdo8Ij8ry3t5+CwQ3VQmeKqiVX0DiDFTNERN7hsUDk+eefx6BBg7Bw4UJUVlZi0KBBGDRoEPbv3++pl/RLcyf2RVJcBApKa7D46+NKD8ejNmcXYunWU255rje2npRcQeMMVswQEXmHxwKRlStXQhCEVv/GjBnjqZf0S1HaUCy5eyAA4L97zknuceEKJbbtFpuYuZPY/dRa5r6zWDFDRORdHmvxTtLd2CMevx7WFf/ZfRZ//vwIvp01Cm20nvnRKLVtt6MmZuYSdeG477okvLHN9uyJedKqmLm/O68YM1YdbFX2KxX7HxAReZ9P9REJZvNu7YvObT27RKPktt1yci4W3p6C5PZRkp/XaBKaqmiq6vDoiGRnh4jYSA37HxAReRlnRHxElDYUr/5qIB5YsQf/3XMOkwYkYkRP56toxIuzWF0ypGtbu/u7iP05rG305g5Scy5mZfTGxNREZOYVSzr+TFE1Ri7Z3iK4itKGoKrOKHuM2lA1y3aJiLyMgYgPMV+i+dNnh7H56VHQRcpvdGZt+SUuSoOSKttLFuZLHZ4oWxUrXAxltTYbi+ljtJg5rqek41UAdJEavLH1ZKuvOxOEAIChvI5lu0REXsalGR8z79a+6NYuEoVltViwPlv2420tv9gLQszZW0JxJcnV0d4UKgAv3NG/eTbG0fHiK7s7zZZlu0RE3sVAxMdEaUOx9L50hKhV+OrwRazPKpD8WLEyxZWLs60lFHHDOMtN6OTklcjdm8Le8bMyejndndUeW9+/EpVGRETBgEszPmhQl7b447heWLr1JP7yZTaGdG2Lzm0jHT5OTmWKJRWaLvDWylbFWRbLS6+Y5ConwVPu3hS2jt945KLM79A+R9+/EpVGRETBgDMiPmrG2B4Y1CUWFbWNmPPpYUl34M4uK5iXrQJocedf32iym+QKNCW5yl2mGd6jHe5M74ThPdo5TI61drychmO/G9Udq6cPw2MjugGQt225kpVGRETBgDMiPio0RI037kvHrf/4EXvyS7Dix9N4cnQPu49xthuo/trdPYBWFShKJ7naIiX5FQC6xEXimYzeiAgLwfAe7TA0Oa7V7IbexuyGvaUub1QaEREFAwYiPqxruyi8cHt//PmLI/h///sZw7u3Q1pSrM3jxYuzlOWZuCgNFtzWH/qYpuWILTkGq8sv7khy9QQxmfWpjw+2SF619PLUVESEhTR/LGdpyNFSl1JBGBFRIOHSjI+757rOuDVVjwajgBmrDqLMTtdQ80oTR0qqGqCPCW++gHoqydWTbCWzhl4LKqYO6oSberXeLVrq0pDU4IqVNkREzuOMiI9TqVR45e6ByL5YhvMlNfjz54ex/KEhUKl+uXhaNi979MZu+OinMw6fW7yAeiLJ1XJM9hJSXWE5w7HjZBG+OHgBsZEa/GVyP5eeW2pwxQ3yiIicx0DED+giNFj2wGDc/c5P+PbYJXy06wweG9nUytxW8zIpxAuoO5JczYMMb1eZiDMcRy+UYc6nhwEAf70zFe3aaF16XilN1WxV2hARkTRcmvETAzvH4v8mNd3hL/7mOLLOl7rUvCw2UtN8AXUlydWydFepKpPaBiPmfJaFRpOAyQMScftA1wMeR03VAG6QR0TkKs6I+JGHb+yGPfkl+CbbgJmrDqKh0eR0XkdpdQO25BgwMTXRpSRX84uwklUmS7eexMlLlYhvE4YXp6S2WLpyhZiHIrXShoiI5GEg4kdUKhWW/Gogjl0sx7mSateeCy2DgoW3p+DJjw86fJxlkqs5b1aZmOegFFfW470dpwEAi+8aiLioMJee25LcJmxERCQdAxE/ExOuwfKHhmDKsl2oN5qcfh7LoGBiaqLsJFepn5dynJzkVms5KAAwvHscxqckSBqDXGIeChERuRcDET+U0jEGT47ugTe3n3L5ucSgYHN2IdYflravja2cEqm5JkUVdVifVdAccGzJMUhObrXVbh4AMk+XYHN2odPLJd6q9CEiol8wEPFTT2f0wge7Tju95b2oQ3S43Yu7JX2MFkOT46xetKV0O1WrgBc3HW/+ODZSY3XzOmv72Dja1E8F4IWvjiE6XIOiyjpZwQT3kyEiUoZKEASf3Ua0vLwcOp0OZWVliImJUXo4PufrIxfx+1WHWn1e7DQaG6lBWXWD3dLTH54di9GvfSe5j0hspAb3XdcZXx0utHrRBoCnruWauPrGEse4c+44hKhVyMwrxrQVu2U9h5RgwlYgJoYvcjb1IyIieddvlu/6sUkDO+Lv9w5s7iQq0uvCsfyhwXjlrgEA7JeeHjh7VVYzs9LqBry7I99meS4Aq91OnVnhMM9jAZzrd+KobNhRpQ8gf1M/IiKSjoGIn7trcBK2zRkNXURTE7MBnXT47k9jMDE10WYLdPP+H+5qTy5c+/fcuqMY1zcBO+eOw+rpw/CP+9OxYHI/uHIdF8foTL8TR8GEnEofIiJyP+aIBICu7aLw8eM34L73MnG0oAx/+TIbr/1qIFQqlcPSU3e3Jy+pasCQl7bg3iGdkZGix20DO2LjkYsuPac4xuu7tUWEJgQ1DfLyYuyVDXM/GSIiZTEQCRADOuuw7IHBePxf+/D5gQto1yYM8yb2hUqlslt6KiXBVK6K2kZ8sOsMPth1Bom6cNx/fZJTz2PZQn3lT2dkByHmrAUT3E+GiEhZXJoJIGP7dsDfpjblhbz7w2m8uS3X4WPstTF3B0NZLZZuPYXYSI2s57dsob4+qwAvXau2GdUr3qmmZWeKWjeBEwMxW2NToSnhlfvJEBF5BgORAHP/0C5YcFtTYLF060m8+0Oew8fYyiVxB7G1u8hW4mxsZMuN+szzWN7YehJPr8lq/tqOU0Uoqapv/rhtZCh0EY4n99bsO9cqT4T7yRARKYvluwFq2Xe5eO3bnwEAf5ncD7+9qbvDx5j3BjlTVI03tp4E4HoZrmhWRm+s2XfOatmvrTyWd3/IxeJvfrb7vGK5shSrpw+zukzFPiJERO4j5/rNHJEANWNsT9Q2GPHW9ly8tOk4ymsbMSujl93N4CxzSfro27S6OMdFaSTt7mtNt/hI7Jw7zmbirGWAkHW+FK84CEIAeYGSraRT7idDRKQMBiIBbPb43tCGqvH6/07izW2nUF7TgOdvS4Fa4sV1fIoe0VoNMk8XAWgKUq7vFierAZq5DtHhkvdsOXqhDA+u2O222RjzMdjC/WSIiLyPgUiAsLVPysxxvRATocHz649h5U9nUFJVj1d/NRDhmhC7z2dtqWLV3rOYmt4J91/fBW9sPSk5SLCsfnEkM68YT/xnP6rqXWtf78oYiIjIOxiIBABrQUNshAaPjuiGmeN64TfDuyE6PBTPfnYEXx2+iHMl1XjvN0Nszg7YanleUtWAD3adaXr+a8ml1vaJMSc34XPD4YuY8+lh1BtN6KuPxglDhcPHWHtNweJjOWMgIiLvYdWMnxODBsulktKaBizdegpDXtqCzdmFmDqoM/792FDoIjTIOl+KKW/vQtb50lbP52hjOVFZdQNKqxswK6MXlt6bhgWT++HREd1aldWaV7/YIwgC/vl9Lv6w+hDqjSbcmqrH2qdutFtaa04ss/3nA/Y7yRIRkW9h1YwfM5oEjFyy3WG+hgpN+7+MT9FjfVYBXvnmBC5X1CFUrcLciX3x+Mjk5rwRORvLWW5KJ45JbsJnWU0D/vTZYWzJuQQAeOTGblhwW9PshRhoAbaTUi03p3NmDERE5D5yrt8MRPyYnKAhNlKD8NAQGMpbBy2je7fH4rsGoGNsBNZnFbTo2SGFrZJYKQ6du4pnPsnC2eJqhIWosfCOFDwwtEuL6h5rS0/mWGZLRORbWL4bJOTsf9KUy2E9n+OHk1eQ8fcf8Kdb+qBPQrRHxyGqbTBi6ZaTWPHjaZgEoFNsBN55aDAGdo5tdaxlaW18lBZQAUWVdZzxICLycwxE/Ji79j/RhKhQXW/EXzfmoK8+Gm0jNbjqIAnV2XEIgoANRwqx5JsTKCitAQBMHdQJC29PQWyk7bbtIWoVhibHtVhyuW1gRwYgRER+joGIHxP3SXGmp4e5BqOAx0cm49P952VVqcgpiTWaBPzvmAHv/JCHIxfKAAD6mHC8OCUV41MSWhxnLb+DnU+JiAITAxE/Ju6T8uS1ZE5XDOysw8yxPbF8Rx7+9dMZ1DaY7B4vtST2XHE11h68gNX7zuFSeR0AIEITgqfG9MD0m7ojIuyXfia2go070hLx3o78VsmqhrJaPPXxQVbEEBH5MQYifm5iaiKWPzQY89YeddjTw54O0eFoGxWG+bf2w+MjkvH+znx8su88ymqsP2dcVBiendAbo3q3h8kkoK7RhIraBhSU1uDkpQpsPX4Zh8+X4nJFXavHRoeHondCm1ZBiLXeJYVltXh3R77VMYgb6i3akIPxKXou0xAR+SFWzQQIo0nA29tz8dGufJSaBQ/6GC1qG00oq26wWv5qrQRXVFNvxMYjF7HpaCF2nipCo8k9bxVr5bZSypDtcaVyh4iI3ItVM0EoRK3C0xm9MHNcz1Y5FltyDHjq44OyO45GhIXgnuuScM91SSirbsC+MyXYd7YE2QVlOFtcjYulNTCPTVQqQEpYKx7ywlfHmqthXM1zcaZyh4iIlMdAJMBYVpfszS/B+BQ93nlocKv8C72MZE9dpAYZKQnIMEssFZdkahqMCFGrMGHpDzCUt16KscVQXoe3t+eivKZe3jdphbsqiIiIyLsYiAQYe9UlO+eOc2vHUbVahYiwEESEhSAzr1hWECJauvUkorT2N+Czh5vZERH5NwYiAcRWwqc3qktcWRqpqpO2yy43syMiCjzc9C5A2NusTvzcog05MLop4dSSp5dGxvVtz83siIgCEGdEAoSjhE8BTaWwe/NLPFJdIjZXM5TVOty51xnTb+rRqrMqW7sTEfk/BiIBQurSiKeqS8Tmataqc1yVqPsl6GCJLhFRYOHSTICQujTiySWUiamJeOehwa2WUOxNWqjQtDOwCr/kfJh/TQXmgBARBTLOiAQIR0sj3qousdwpt0N0OK5W1WPGqqY29NaSTV+5awAAuFReTERE/omBSICwtzRiXl0CAJl5xR7Ns7C2hPKO2nEfE8sAhjkgRESBjy3eA4y9PiJA61kHb+5ga2tnXSIiCixyrt8MRAKQtQu+2Obd8odtue8LERGRq7jXTJCzXBpx1GOEO9gSEZFSWDUTBOT0GHE3o0lAZl4x1mcVIDOv2GMN1YiIyD9xRiQIeKPHiK3lICVzUoiIyPcxEAkCUnuHFFXUYX1WgexEUmsJsrGRGpRWN7Q61hv73hARkf9gIBIEpLRfV6uAFzcdb/5Y6syFrY32rAUhgPScFFbYEBEFBwYiQUBK+3XL1A0pMxf2kmDtcbTvjb0SZM6iEBEFFiarBgm57del7NjrKAnWEWs5KeIMi+XzioHR5uxCp1+PiIh8D2dEgohl+/WiiroWyzGWHM1cuLqBnmXuCsuMiYiCD2dEgozYY+TO9E6Ij9ZKeoytgMPZDfRU+GVHXXNKlhkTEZEyPBaInDlzBo8//jiSk5MRERGBHj16YOHChaivr/fUS5JM8W2kBSK2Ag4xCVbO3IT5vjeWsxreKDMmIiLf4rFA5MSJEzCZTHj33Xdx7NgxLF26FMuXL8dzzz3nqZckGTZnF2LOp1l2j7E1cyESk2DFYy0fCzSV8ZrT68JtJsBKnWFxdiaGiIh8j1f3mnnttdfwzjvv4PTp05KO514znmGr5NacnD1o7FW5yNlR12gSMHLJdptlxio0BTI7545jjggRkQ/z2b1mysrKEBdn/e4aAOrq6lBXV9f8cXl5uTeGFVSkltzq7ZTLWvb4GJ+itxtwWEt0tcZembG9JR0iIvJfXgtEcnNz8dZbb+H111+3eczixYuxaNEibw0pKEktub33uiSrQYine3yIZcaWr2EvMCIiIv8le2lm3rx5WLJkid1jjh8/jr59+zZ/XFBQgNGjR2PMmDF4//33bT7O2oxIUlISl2bcaH1WAZ5ek+XwuNhIDQ78ZXyL2QdbSzpylnGkYmdVIiL/5dGlmTlz5uCRRx6xe0z37t2b/3/x4kWMHTsWN954I9577z27j9NqtdBqpVVykHOkJnqWVje06B/i7R4fYpkxEREFNtmBSPv27dG+fXtJxxYUFGDs2LEYMmQIPvroI6jVbFuitKHJcYiN0KC0xvpeMObMy2Tl9PhgAEFERFJ5LDIoKCjAmDFj0KVLF7z++uu4cuUKDAYDDAaDp16SJAhRq/DoiG6SjjWfPWGPDyIi8gSPJatu2bIFubm5yM3NRefOnVt8zYsVw2TFzHG98NFPZ2zukCuWyZr3D2GPDyIi8gSPzYg88sgjEATB6j9SVohahVfuGmC1I6qtMllHXVQdNT8jIiKyhkkbQUosk0202I3XVudTKV1U2eODiIjk8mpnVbnYWdXz5JbJerqPCBER+T85128GIiQbe3wQEZE9PtvinQIDe3wQEZG7MEeEiIiIFMNAhIiIiBTDQISIiIgUw0CEiIiIFMNAhIiIiBTDQISIiIgUw0CEiIiIFMNAhIiIiBTDQISIiIgU49OdVcXu8+Xl5QqPhIiIiKQSr9tSdpHx6UCkoqICAJCUlKTwSIiIiEiuiooK6HQ6u8f49KZ3JpMJFy9eRHR0NFQqbqpmT3l5OZKSknD+/HluEOhmPLeewfPqOTy3nsNzK40gCKioqEDHjh2hVtvPAvHpGRG1Wo3OnTsrPQy/EhMTw18OD+G59QyeV8/hufUcnlvHHM2EiJisSkRERIphIEJERESKYSASILRaLRYuXAitVqv0UAIOz61n8Lx6Ds+t5/Dcup9PJ6sSERFRYOOMCBERESmGgQgREREphoEIERERKYaBCBERESmGgYifevnll3HjjTciMjISsbGxkh4jCAKef/55JCYmIiIiAhkZGTh16pRnB+qHSkpK8OCDDyImJgaxsbF4/PHHUVlZafcxY8aMgUqlavHvySef9NKIfdeyZcvQrVs3hIeH44YbbsDevXvtHv/ZZ5+hb9++CA8Px4ABA/D11197aaT+R865XblyZav3Z3h4uBdH6x927NiB22+/HR07doRKpcKXX37p8DHff/89Bg8eDK1Wi549e2LlypUeH2egYSDip+rr63HPPffgqaeekvyYV199FW+++SaWL1+OPXv2ICoqChMmTEBtba0HR+p/HnzwQRw7dgxbtmzBxo0bsWPHDjzxxBMOHzd9+nQUFhY2/3v11Ve9MFrf9cknn2D27NlYuHAhDh48iLS0NEyYMAGXL1+2evxPP/2EadOm4fHHH8ehQ4cwZcoUTJkyBdnZ2V4eue+Te26Bpk6g5u/Ps2fPenHE/qGqqgppaWlYtmyZpOPz8/MxefJkjB07FllZWXjmmWfw29/+Ft9++62HRxpgBPJrH330kaDT6RweZzKZBL1eL7z22mvNnystLRW0Wq2wevVqD47Qv+Tk5AgAhH379jV/7ptvvhFUKpVQUFBg83GjR48Wnn76aS+M0H8MHTpUmDFjRvPHRqNR6Nixo7B48WKrx997773C5MmTW3zuhhtuEH73u995dJz+SO65lfp3gn4BQFi3bp3dY/785z8L/fv3b/G5++67T5gwYYIHRxZ4OCMSJPLz82EwGJCRkdH8OZ1OhxtuuAGZmZkKjsy3ZGZmIjY2Ftddd13z5zIyMqBWq7Fnzx67j/3vf/+L+Ph4pKamYv78+aiurvb0cH1WfX09Dhw40OL9plarkZGRYfP9lpmZ2eJ4AJgwYQLfnxacObcAUFlZia5duyIpKQl33nknjh075o3hBjS+Z93Dpze9I/cxGAwAgISEhBafT0hIaP4aNZ2nDh06tPhcaGgo4uLi7J6nBx54AF27dkXHjh1x5MgRzJ07Fz///DPWrl3r6SH7pKKiIhiNRqvvtxMnTlh9jMFg4PtTAmfObZ8+ffDhhx9i4MCBKCsrw+uvv44bb7wRx44d48aiLrD1ni0vL0dNTQ0iIiIUGpl/4YyID5k3b16rhDLLf7b+0JB9nj63TzzxBCZMmIABAwbgwQcfxL///W+sW7cOeXl5bvwuiJwzfPhw/OY3v0F6ejpGjx6NtWvXon379nj33XeVHhoRZ0R8yZw5c/DII4/YPaZ79+5OPbderwcAXLp0CYmJic2fv3TpEtLT0516Tn8i9dzq9fpWCX+NjY0oKSlpPodS3HDDDQCA3Nxc9OjRQ/Z4/V18fDxCQkJw6dKlFp+/dOmSzfOo1+tlHR+snDm3ljQaDQYNGoTc3FxPDDFo2HrPxsTEcDZEBgYiPqR9+/Zo3769R547OTkZer0e27Ztaw48ysvLsWfPHlmVN/5K6rkdPnw4SktLceDAAQwZMgQAsH37dphMpubgQoqsrCwAaBH0BZOwsDAMGTIE27Ztw5QpUwAAJpMJ27Ztw8yZM60+Zvjw4di2bRueeeaZ5s9t2bIFw4cP98KI/Ycz59aS0WjE0aNHMWnSJA+ONPANHz68VYk537NOUDpblpxz9uxZ4dChQ8KiRYuENm3aCIcOHRIOHTokVFRUNB/Tp08fYe3atc0fv/LKK0JsbKywfv164ciRI8Kdd94pJCcnCzU1NUp8Cz5r4sSJwqBBg4Q9e/YIO3fuFHr16iVMmzat+esXLlwQ+vTpI+zZs0cQBEHIzc0V/vrXvwr79+8X8vPzhfXr1wvdu3cXRo0apdS34BPWrFkjaLVaYeXKlUJOTo7wxBNPCLGxsYLBYBAEQRB+/etfC/PmzWs+fteuXUJoaKjw+uuvC8ePHxcWLlwoaDQa4ejRo0p9Cz5L7rldtGiR8O233wp5eXnCgQMHhPvvv18IDw8Xjh07ptS34JMqKiqa/5YCEP7+978Lhw4dEs6ePSsIgiDMmzdP+PWvf918/OnTp4XIyEjh2WefFY4fPy4sW7ZMCAkJETZv3qzUt+CXGIj4qYcfflgA0Orfd99913wMAOGjjz5q/thkMgkLFiwQEhISBK1WK9x8883Czz//7P3B+7ji4mJh2rRpQps2bYSYmBjh0UcfbRHg5efntzjX586dE0aNGiXExcUJWq1W6Nmzp/Dss88KZWVlCn0HvuOtt94SunTpIoSFhQlDhw4Vdu/e3fy10aNHCw8//HCL4z/99FOhd+/eQlhYmNC/f39h06ZNXh6x/5Bzbp955pnmYxMSEoRJkyYJBw8eVGDUvu27776z+ndVPJcPP/ywMHr06FaPSU9PF8LCwoTu3bu3+JtL0qgEQRAUmYohIiKioMeqGSIiIlIMAxEiIiJSDAMRIiIiUgwDESIiIlIMAxEiIiJSDAMRIiIiUgwDESIiIlIMAxEiIiJSDAMRIiIiUgwDESIiIlIMAxEiIiJSDAMRIiIiUsz/B+lma6S5lWvuAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "n_train = 250\n", "n_test = 400\n", "sigma_obs = 0.3\n", "\n", "fn = (\n", " lambda x: 20.0 * np.sin(x / 40.0)\n", " + 0.1 * np.power(x, 3.0)\n", " + 0.1 * np.power(1.0 + x, 3.0)\n", " + 1.2 * np.cos(4.0 * x)\n", ")\n", "\n", "x = np.linspace(-1, 1, n_train)\n", "y = fn(x)\n", "y += sigma_obs * np.random.normal(size=(n_train,))\n", "y -= np.mean(y)\n", "y /= np.std(y)\n", "\n", "x_test = np.linspace(-1.3, 1.3, n_test)\n", "\n", "x_true = np.linspace(-1.3, 1.3, 1000)\n", "y_true = fn(x_true)\n", "\n", "plt.scatter(x, y)\n", "plt.plot(x_true, y_true)" ] }, { "cell_type": "code", "execution_count": 120, "id": "675505a3-c052-4df7-8bab-2af8d20184a6", "metadata": {}, "outputs": [], "source": [ "mesh = mesh_1d(x_test, degree=2)" ] }, { "cell_type": "code", "execution_count": 121, "id": "8061e405-f68d-40c0-ac50-1f324de77e70", "metadata": {}, "outputs": [], "source": [ "sigma = 3 * np.std(y)\n", "rho = 0.5 * np.max(x) - np.min(x)\n", "spde = spde2_pcmatern(mesh, alpha=2, prior_range=rho, prior_sigma=sigma)" ] }, { "cell_type": "code", "execution_count": 122, "id": "e1e95e98-2328-459a-8b12-ea058676e68a", "metadata": {}, "outputs": [], "source": [ "A = make_projection_matrix(mesh, loc=x)\n", "Ap = make_projection_matrix(mesh, loc=x_test)\n", "index = spde.make_index(\"ix\")" ] }, { "cell_type": "code", "execution_count": 123, "id": "571efcc8-a2b8-4dbb-ab61-da58262f6ac6", "metadata": {}, "outputs": [], "source": [ "st_est = inla_stack(tag=\"est\", data=dict(y=y), A=A, effects=dict(), index=index)\n", "st_pred = inla_stack(\n", " tag=\"pred\",\n", " data=dict(y=np.repeat(np.nan, x_test.size)),\n", " A=Ap,\n", " effects=dict(),\n", " index=index,\n", ")\n", "st_full = combine_stacks(st_est, st_pred)" ] }, { "cell_type": "code", "execution_count": 124, "id": "607dfe5c-f1e1-41c2-92d9-6f804ca5aa8a", "metadata": {}, "outputs": [], "source": [ "formula = \"y ~ -1 + f(ix, model=spde)\"\n", "data = stack_data(st_full) | dict(spde=spde)\n", "data[\"spde\"] = spde\n", "res = inla(\n", " formula,\n", " data=data,\n", " family=\"gaussian\",\n", " control_predictor=dict(A=stack_A(st_full), compute=True),\n", ")" ] }, { "cell_type": "code", "execution_count": 125, "id": "e491845a-ea43-492e-bb6e-1c4b2b6b84cd", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Time used:\n", " = 1.83, = 0.188, = 0.0102, = 2.03 \n", "Random effects:\n", " Name\t Model\n", " ix SPDE2 model\n", "\n", "Model hyperparameters:\n", " mean sd 0.025quant 0.5quant\n", "Precision for the Gaussian observations 12.42 1.16 10.26 12.38\n", " 0.975quant mode\n", "Precision for the Gaussian observations 14.79 12.31\n", "\n", "Marginal log-Likelihood: -81.72 \n", " is computed \n", "Posterior summaries for the linear predictor and the fitted values are computed\n", "(Posterior marginals needs also 'control.compute=list(return.marginals.predictor=TRUE)')\n" ] }, "execution_count": 125, "metadata": {}, "output_type": "execute_result" } ], "source": [ "res" ] }, { "cell_type": "code", "execution_count": 126, "id": "f70b29e5-f242-4374-b23f-65c005084887", "metadata": {}, "outputs": [], "source": [ "pred_ix = stack_index(st_full, tag=\"pred\")\n", "fitted = pd.DataFrame(res.get_summary(\"fitted.values\"))\n", "fitted = fitted.iloc[pred_ix]\n", "fitted[\"x\"] = x_test" ] }, { "cell_type": "code", "execution_count": 128, "id": "3b8794f5-45a7-4193-9aa9-91f59654205d", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 128, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiIAAAGdCAYAAAAvwBgXAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAA9hAAAPYQGoP6dpAACNIUlEQVR4nO3dd3hb5fXA8e/V9t6xncRJvEji7AQCIWEHCFBW2WUVAi0UfmWV2bJK27ApBcpeLXuWlkDYM4SQZchO7Nix4723rXV/f1zk2LFkDVuWbJ/P8/iBWFf3vpZl6eh9z3uOoqqqihBCCCFECOhCPQAhhBBCjF4SiAghhBAiZCQQEUIIIUTISCAihBBCiJCRQEQIIYQQISOBiBBCCCFCRgIRIYQQQoSMBCJCCCGECBlDqAfQH6fTSXl5OTExMSiKEurhCCGEEMIHqqrS0tLC2LFj0en6n/MI60CkvLycjIyMUA9DCCGEEAEoLS1l/Pjx/R4T1oFITEwMoP0gsbGxIR6NEEIIIXzR3NxMRkZG9/t4f8I6EHEtx8TGxkogIoQQQgwzvqRVSLKqEEIIIUJGAhEhhBBChIwEIkIIIYQIGQlEhBBCCBEyEogIIYQQImQkEBFCCCFEyEggIoQQQoiQkUBECCGEECEjgYgQQgghQkYCESGEEEKETFADkccff5yZM2d2l2hfsGABH374YTAvKYQQQohhJKiByPjx47n77rtZt24da9eu5cgjj+Tkk09m8+bNwbysEEIIIYYJRVVVdSgvmJiYyH333cfSpUu9Htvc3ExcXBxNTU3S9E4IIYQYJvx5/x6y7rsOh4M333yTtrY2FixYMFSXFUIIIUQPbW1QXg7NzZCSAhMmhHY8QQ9ENm7cyIIFC+js7CQ6Opp3332XvLw8t8d2dXXR1dXV/e/m5uZgD08IIYQYFaxW2LQJ8vOhvh46OuCww0IfiAR918zkyZPJz89n9erVXH755Vx44YVs2bLF7bHLli0jLi6u+ysjIyPYwxNCCCFGvMZGWLECPvsMVBVycyEmJtSj0gx5jsjixYvJzs7mySef7HObuxmRjIwMyRERQgghAlRbC598AmVlkJkJJpP2/cJCmDsXDj108K8ZljkiLk6ns1ew0ZPZbMZsNg/xiIQQQoiRqb5eC0IqKiAnB/T63rcP7VSEe0ENRG6++WaOO+44JkyYQEtLC6+88gpffvklH330UTAvK4QQQox67e3w5ZfaTIi7ICRcBDUQqa6u5oILLqCiooK4uDhmzpzJRx99xNFHHx3MywohhBCjmtMJ330HBQWegxBFGQUzIs8++2wwTy+EEEIINzZvhp9+gokTwWj0fFw4BCLSa0YIIYQYQaqr4fvvITYWIiM9H6coQzem/kggIoQQQowQdjusXg0tLZCa6v14pzP4Y/JGAhEhhBBihNi6FbZvD32RMn9IICKEEEKMAI2NsGYNxMWBr5UwJEdECCGEEAOmqntLt/uyJBNOJBARQgghhrnycm2nTHq670mo4bJ9VwIRIYQQYhhzOmHDBq2pnb/dUCQQEUIIIcSAFBVphcvGjfPvfrJ9VwghhBADYrfDjz+CTgcREf7fX7bvCiGEECJgu3ZBcbH/syHhRAIRIYQQYhiy2bTZELMZTCb/76+ni9zdd0DTnkEfmz8kEBFCCCGGoaIiKC2FtLTA7j+182WyKx6C548P6RqNBCJCCCHEMGO3a03tTKbAZkMU1cbc9oe0fxz0Oy3JJEQkEBFCCCGGmd27tdmQ9PTA7j++7nVinSV0GsbAvAsHd3B+kkBECCGEGEacTti0CfT6QGdD7ORUPgBAUdr/gTGA7TaDSAIRIYQQYhjZs0fbKRNobsjY+reJ7tpFh5JIccrFgzq2QEggIoQQQgwTqqp12HU6A6sbguogt/J+AH6MvBK7LnpwBxgACUSEEEKIYaK6GgoLA29sl97wHjGdO7Dq49kUeengDi5AEogIIYQQw8TOndDeDjExAdxZdbLfz7MhRWMux6aLlcqqQgghhPBNczNs2wbJyYHdP7XpI2I7NmPTxVA05reDO7gBkEBECCGEGAaKiqC+HhISArizqnbnhhSPuQSbIcH17ZCTQEQIIYQIc1artmU3Li6w2mNJrd+Q0LYWh2Jh15jfDf4AB0ACESGEECLMlZRAVRWkpAR2/9wKrW5ISfIFWI1jAFCUwRrdwEggIoQQQoQxVdVyQ/R6MBr9v3982zpSWr7EiYHC1P/rdZskqwohhBCiX1VVWkn3QLfs5lQ+CEBZ0pl0mCd0f19mRIQQQgjhVWEhdHRAVJT/943u2EZ64/uoKBSkXt3ndpkREUIIIYRHbW2wfTskJgZ2/5xKrcNuRfyJtEZMHsSRDR4JRIQQQogwVVICDQ2BBSIRXcWMq38TgIK0awd5ZINHAhEhhBAiDDmdWl8Zs1lLVPVXTtU/0OGgOvYomqLmuD1G6ogIIYQQwq3KSigrC2zLrtlWSUbtSwDsTLtukEc2uCQQEUIIIcLQrl3Q1QWRkf7fN6vqMfRqF/VRB1EffbDbYxRFZkSEEEII4UZ7O+zYEVhuiNFez6Sa5wDYmX5tv/t0JRARQgghRB8lJVpfmUACkUnVT2NwttIUMZ3q2GM8Hid1RIQQQgjRh6pqsyGBJKnqHa1kVT8O/LxTxku0IXVEhBBCCNFLdTWUlgaWpDqx9kVMjgZazVmUJ5wy6GMLBglEhBBCiDBSXBxYJVWds4usqkcAKEi7BhTv0ymSIyKEEEKIblartiwTF+f/fcfXvUqErYIO4zj2JJ49+IMLEglEhBBCiDCxZw/U1EBSkn/3U1Q7OVV/B6Aw9f9QdSbv95FkVSGEEEL0VFioBQhGo3/3S2/4D1FdRXQZkihJvsDn+0myqhBCCCEAaGqCoiJITvbzjqqT3IoHACgaczkOvW/JJTIjIoQQQohupaVaMOJvfkhq00fEdm7BpouhOOUSv+474pNVly1bxgEHHEBMTAxjxozhlFNOYfv27cG8pBBCCDHsOJ1akmpkpJ8zFapKbuX9ABSnLMVmSPD5rqNiRuSrr77iiiuu4Pvvv+eTTz7BZrNxzDHH0NbWFszLCiGEEMNKdTWUl/u/LJPU+g0JbWtxKBZ2pV7h93XDYUbEEMyTr1ixote/X3jhBcaMGcO6des49NBDg3lpIYQQYtjYvRs6O/1vcOfKDSlJPh+rcUxA11bV0M6OBDUQ2VdTUxMAiR6K53d1ddHV1dX97+bm5iEZlxBCCBEqgdYOiWtbT0rLlzjRU5j6f0EZ21AYsmRVp9PJ1VdfzcKFC5k+fbrbY5YtW0ZcXFz3V0ZGxlANTwghhAiJsrLAaofkVj6o3T/xDDrMEwO6tqqGfnlmyAKRK664gk2bNvHaa695PObmm2+mqamp+6u0tHSohieEEEKERHGx/7VDoju2k9b4PvBzOfdhbEiWZq688kref/99vv76a8aPH+/xOLPZjNlsHoohCSGGMYdT5YeieqpbOhkTY2F+ZiJ6XZhsARDCD62tWhEzDxkLHuVU/R0FlYr4X9AaMSWgaytK6GdDIMiBiKqq/N///R/vvvsuX375JZmZmcG8nBBiFFixqYI7/7eFiqbO7u+lx1m4/cQ8lkxPD+HIhPDfnj1a7ZCcHN/vE2EtZVzdGwAUpF074DGEOhgJ6tLMFVdcwUsvvcQrr7xCTEwMlZWVVFZW0tHREczLCiFGqBWbKrj8pfW9ghCAyqZOLn9pPSs2VYRoZEL4T1WhoEBbktH58W6cVfkIOuzUxBxKY9S8gK8/KuqIPP744zQ1NXH44YeTnp7e/fX6668H87JCiBHI4VS5839bcPfhzfW9O/+3BYczDOaahfBBfb1WTdWf2iEmWw0Ta/8FQEHadUEa2dAK+tKMEEL4qr/cjx+K6vvMhPSkAhVNnfxQVM+CbD+3HwgRAnv2aDki48b5fp+s6sfRqx00RM6lNuawQRlHqN+qh7SOiBBCeOIt96O6xXMQ0pOvxwkRSq6S7tHRvi+RGBzNTKp5Bvg5NyRc1lYGSJreCSFCzpfcjzExFp/O5etxQoRSVRVUVvpXO2RizbMYHU20WCZTGX/CoIxjVNUREUIId3zN/Zg3MYH0OAuePgMqaDMo8zP93AcpRAiUlGgVVSMifDte5+wgq+qfABSkXQ3KwN++w2VCRQIRIURI+Zr7sW53A7efmAfQJxhx/fvWE6byQ1E97+WXsaqwThJXRViy2WDnToiN9f0+GXUvY7FX027KoCzxjEEbS6hnQ0ByRIQQIeZP7sfJs8fx+Hlz++SSpMVZOGlWOnct3yr1RUTYKy+H2lqYMMG34xXVTk7lwwAUpv4eVfGjBGt/5/05gg91MCKBiBAipPzN/VgyPZ2j89J67a5paLNyxSvr+yzvuHJMHj9vrgQjImwUF2vJqr6WdB9b/zaR1hK6DMmUJp8X1LGFggQiQoiQmp+ZSHqchcqmTrd5IgrajEfP3A+9Tuneoutwqiy653OPOSYKWo7J0XlpUgZehFxHh1bSPSHBxzuoTnIqHwJg15jLcegiB31MoZ4RkRwRIUYZh1NlVWFd2ORR6HWK19yP20/M8xhE+FNfJNx+djH67NkDDQ2+ByKpTSuI7dyKTRdLccolwR1ciMiMiBCjSLj2aVkyPd1j7oe3sfmaY/KvVUVc83o+lc3h9bOL0aWoCPR67csrVSW38gEAiscsxW6IH/TxhMP2XQlEhBglXLU6BjuPYrA64brL/fB2LodTpbaly6fzf7ipqs/3JIdEDKXmZi0/xNfaIUmt35DQthaHYmHXmN8N+njCZfuuBCJCjALeanUEmkfhywyLr4GKvwGNu2v7S3JIxFDas0cLRlJTfTs+t+JBAEqSz8NqHBOUMYV6NgQkEBFiVAhGnxZfZlgAn5aC/F0y8nTtQEiPGjEUXJ12zWbfOu3GtW0gpeULnOgpTP190McWSpKsKsQoMNh9WnyphnrTOxu9lm0H38q799TRqXL7e+6vPRDSo0YEU10dlJX53mk3p1KbDSlLPIMO88SgjEmWZoQQQ2aw+7T4MsPS2G7zeJtrOeTIKalel4zu+O8Wssxp1NYolJXBxqp6qoIQNEiPGhFMZWVap93x470fG925g/TG/wE/l3MPslDPiEggIsQoEEitjv4MdPbAtRzy71XFXgOayuZO7nqtiEi9mcQIC136wQ1C/P3ZhfBXz067vsipfAgFlcq4E2iNmBrcwYUBCUSEGAVctTouf2k9CvQKRnyp1bGvwZo92F3f7tNx37Rt1f6nGaIMg1PeGgL72YXwl6vTri9JqhFdxYyrex2AHenXBXVcihL62RCQHBEhRg1XrY60uN5BRFqcxe/tq64ZloG+dasBvAq22d0v+QQikJ9dCH+VlkJXl2+ddnMrH0KHg+rYo2iKmhf8wRH6YERmRIQYRfyp1dHfdlpvMywqEB9ppKnd5jGpVKfAv78vGeSfcK8og5HoSD1VPQqYpcWaOWf+BCYlRw2o5okQvrLZtGWZuDjvx1qsZWTUvQzAzvTrgzyy8CGBiBCjTM8+LZ74sp3WWzVUwG2g4hLs6uqHWGZw0pw0ojPraegcWLE1IQJVUaF12s3I8H5sduXD6FQbtdGLqI9eEPzBEfrZEJBARIhRa98Zj3kTE1i3u4FPtlTy3MriPse7q0LqbYbFXaAymPLiUyhpbaLVbu3+XrzJwmlZeUyJSmfTRphrTOKEw8Egr3YiBFyddk2m/o8z26qYWPsiMHSzIa7tu6EORuRPU4hRyN2Mh07pf5bCUxXS/mZYlkxPx+mE372y3q/xHZI6gW+qvC/b7Giq44L9ZhOhN9BqtxJnspAdm4ju51fYiRMhP19bmz/44PCpmyBGB1en3fh478dmVT2CXu2kPuoAamMOC/rYXEIdhIAkqwox6ngqIObLUknPKqS+6LKq3PbeFr/HOCEmHrPOe1cwu+rkue3reaXgJ2xOR68gBCAyEsakqrzzbR3/XC4dd8XQKivTOu0metkZbrLXManmOQB2pt8w5BFzqIMRmRERYhTpryKqPyqbOlhVWOcx4dVm06ak311ZT22b/8sybxdtocvp8OlYi95Ag7WTVwo28kN1GRdNnkusyQxAfl0Fb+/aQqO1E74FvpWOu2Lo7NqllXP31mk3s+pxDM42miJmUh179NAMjvCZIZRARIhRxFtFVF/dtXwr9W178zLS4yzc9os85o9Np7QUtm7VkvSKbYFdq9NhJ9JgpN2HrboXT55LeXsLH5TsoKC5nvt+/JZLp8yj3trBs9v6LglJx10xFHzttGuwN5JZ/SQAO9KvD0l0IDMiQoghM1j9VHoGIaAt11z+8npOSZnLOGc6sbEwYaLKnpouaPDv3PEmC7/MzGN64hj+vO5LbTajn2MnxyczNSGF6QljeGrbWqo72vj7plWYPCztSMddMRRKS6GpyXsRs8yapzE6m2m2TKUy/hdDM7gwIzkiQowiwe6n8nnDZkitZVXnZu7Y8CnvFG/16/4mnZ7b5x3BnOR0jDo9p2Xl9Xv8vJT07pyQ1Mho/jBzIXkJKdiczn4Ln/mb6yKEP1QVdu7UkqT767Srd7SQVfUYADvT/wDK0L8lh3o2BCQQEWJUGayKqJ4027t4bMtqvqwo7rWl1ldWp4Oilr1TKLOT0jlqXKbH4z8rK2J9bTlOVWVnUx2bG6o5Ij2TjCgfqkchHXdFcNTWQnm592WZSTXPYXI00GrOpjzh1KEZnBuhDkZkaUaIUcRVEfWyl/zbThuUseAgWylnpm4X05Rixis1jFPqmLCrFTM2dGonoLDQaabZFEE18ZSoqRQ7U9mkZpLvzKaRGJ7fvoFIw6Ze+SSRet9e2mIM0nFXDL7SUmhr67+Imd7ZTnbVIwDsTLsOFO+7xAZbuPSakUBEiFHm6Lw04iONNLYPXs8WX8XTwmL9eg7X5XOobiOxipumd87e/0xUukhUmplEFfPZDj1erwucY/nSOYvPnHNZw2TsP7+ktTvsXscSrViwliXimOp9VwP0X/JeCBe7HbZvh9jY/o+bUPMCZnsN7aYJlCWdOTSD8yDUwYgEIkKMMj8U1Q9pEKLHwRG6fE7Xf82RuvWYlL3bclvUCDark9jozKRITadcTSQ3dRZxEcn8b08JrdYOopROYmhnrFLHBKWabF05s5RCsnUV5OjKydGVcwkfUqfG8J5jIW85DmWLOsnruE6dlMemjQopyTCvn95iDqfKo58X8PzKIho79j5usg1YuFNRoXXbHT/e8zF6Zzs5lQ8B2myIqgxeR2l/yPZdIURIDFVehIUuTtd/zaX65UzUVXd/f4tzIp845/GFYzY/qVk490lV+7ICoB6I1r5+/rS2zvWp7ec4Jp4WDtJt5Sjdeo7Q55OsNHOxYQUXG1awwZnDk/Zf8LFz/z7nd5WAn52UTl0dfP+91pAsJ6fvz7BiUwU3vbPRbeAm24CFO8XF4HCA2ez5mIk1z2KxV9NumkBp8rlDNrZwJYGIECPcvksKydH9vEIOAiN2ztF/xu8N75KsNAPQoEbzpuMw3nYcwnZ1gsdGeP5oJIYVzvmscM5Hb3dwiG4jp+u/5mjdWuboCnjC9HeKnKk87jiJtx2H4vh5TWdu8liiDCacqkpSkkJHB3z9NURFQXqPeGLFpop+c2lkG7DYV0eHtlumv0qqekcbOZV/B7S6IaGaDelJlmaEEEHjrqdMWqw5SDkiKsfq1nKz4RUm6aoAKHGm8IzjeN50HEYHlh5HDi4Her50zuZL52ySaeICw0dcoP+ETF0V9+qe5jf65dxrP4uPnfvzefkuPi/fRZTByAFjxjEjIY2O2kQ+/1xhyRJtp4OrAq33n3jvNmBvHY1B8kxGutJSqK93P7vmMqnmGcz2WtpMk9iTdM7QDc6DUAchIIGIECOWq6fMvq8zVc1dgx4ITDO3cbXzCY7WrwOgRo3j7/bTeN1xeHcC6VCpJY4H7WfyuP0kfqX/nCsM/yFHV85Tpof4wTmZW20XsV2dQJvdxpflxXxZXky8ycKCjjyMn6Wz+GiVtzcW+VWB1pflLndBoeSZjCwFBWA0ek5+1jtaya56GIAd6TeExWwIhD4YkUBEiBGov54yru9ZFCNmg56mAMuwu852nv5TbtO9gUlpw6rqedJxIo/bT6Kd0G6N7cDCs47jecNxOL8xvM8l+g+Yr9vOctMtPOM4noftv+yepWm0dvKhdT2VWzK546c9PpWW78lboThPQaHkmYwcdXVQUgLJyZ6PmVTzNGZ7Ha3mLMqSzhq6wXkQLsmqUtBMiBHIl54ynaqN8/ebxS8nTQ3oGgk084zxfv5ifB6T2ka+OpkTrMt4wH5mUIOQfVcy4iONRJk9779tIZIH7GdyZNcDfOg4AIPi5DLD+3xivoGDdL2XXzZ0FvkdhKTHaUssnvgSFN75vy3SFXiYKy2FlhbP23b1jhZyKv8BaB12VSU85gFCPRsCMiMixIjk686YFlsXTbYuv8+/QLeZvxsfI1VpxI6Jb5Jv4uI9+6H6+NnGrNP73F13X/93ZC4HZSX1yrNwOFUOWvZZnx44PVWQxOW2azjSsZ4/G19gvFLLK8a/8rTjeB6wn4mVwKbJzz5gQr95Ht6CQn/zTET4sdu1Ro8xMZ6Pyax+CpOjnlZzNmWJZwzd4PrhmhEJdTAiMyJCjEC+9pTpsNv4qrzYp2OjDSZA5RL9cl4y/o1UpZFaYxYrp37Ol5Fn+xyEAN1BSCAzw1kpUSzITuLk2eNYkJ2EXqdgMuj426nTfbr/5865HNN1L6/aj0CnqPzWsJz3TH9iP6U0gNHApOTIfm/3NSiUcvPDV3m5VjvE07KMwdHcXUV1R/qNYTMbEi4kEBFiBPKlp0ys0cyK0gLsqhODl2Zb8SYLf523kOUpb/Mn48voFZWSxF/xfd43rLONpbK9JaBxuj6IRZl8L2/tKchaMj2dpQsn+XSOdizcbL+US6zXUavGMlVXynumWzlV943P4/A2Hl9v9/c4EX527QKn03PtkMzqJ3/uKZNLWeLpQzs4H8iMiBBi0Ll6yoDnWYdYk5lmWxdpEdH8KndGv+c7b0IK87cdx7SWd3Ci56fx9/BizK38acNq/rHpez7aUxjwWBUgxmLg3xfPJz7C8/KIgvd8jMV5aX5d+1PnPJZ03cPXjhlEKFYeMj3OXwzPYsJ7nogv4wHvQaGv5xHhqblZ2y3jqcGdwdFEVtWjAOwYe2NIesp4MiqSVb/++mtOPPFExo4di6Io/Oc//wnm5YQQPSyZns7j580lLa73J+14k4XD0ieyp60Zg6LjoslzOSBlPEunzCXe1PfYs9JMnFl6Juld+TSqUZxvvYGzi7N4dvsGGq0DX05QgcrmLgx6HXefNgOFvsGT69+3n5jXbz5GIN2Fa4nj17Yb+bv9lzhVhfMMn/Gm6U7SqPN4H0/jcThVVhXW8V5+GasK63A41X6DQl9/LhG+du+GhgaIj3d/e3bVo5gcjbRYJlOW8MshHZsvVDX0MyJBXahqa2tj1qxZXHzxxfzyl+H3CxBipFsyPZ2j89JYvaueD7/spLPBwrzsGP6y4Svt9oxcxkZpGXazk9KZmZhGYXM9TdZO4kwWotp+4tTyixmjNFLqTOEC200Uqeng9H1nSXyEsVePFk+qWzo5efY4Hj9vbt8ibD7W23C96V/+0nqP1VuzI1KosDbR7tib2OpEx9/tp5PvzOEh4z+ZpdvFe+ZbWWr9A5vUrD7ncDceb3VCBvJzifDkcGhJqtHRoHPzsd5kqyGr6jEAto39U1jNhoSToAYixx13HMcdd1wwLyGE8EKvUxhvSiK5DdInwufl22m1WUmNiOKocb3fZHWKQm5cEi0toK9YxfHNvyJGaWerM4MLrTdRTYLP173yiBwW5iTjVFXOfWa11+NdORKu4CnQCqSe3vRdQUGaPZ1PPlXZbSri/fKtve77pXM2J1r/yrPG+5is28Obpj9zte0KvlEWsH9sJvGGKBYvtHDC/N7j8bVOyEB+LhF+ysq0RFVPDe5yK+/H4GyjIXIulfEnDu3gfBTq2RAIs+27XV1ddHXt3UrY3NwctGvV1MC6ddoTKDsbIiKCdikhQm7bNi2Zzq7v4vPyIgB+MXEyhn0+xrW1wZ49kG38nhNaT8OktLPaOYVLrdfRTJRf18xNjWZBdhIOp0p6nIXKpk63MxQK2sxAzxwJvU4Z0FbW/t70HQ6orlZw5mcSbyrqs7y0R03hdOsdPGr8B4fpf+Jx09/5JtFI48Sj2b1bR+suaJ2qNcoD73VC9u1HI1t0R46dOz0nqUZ0lTCx5jkAto27LXwSMtwIdTASVsmqy5YtIy4urvsrIyMjaNdqbYWNG+HDD+Hdd7Ws51D/MoQIhpoaKCyEMWPgk7JCrE4HGVFxzErsndhZWam1MD8s63tObj0Nk9rGt45pXGi90e8gBPbOcIQqR8L1pt9zmy9o5bcXLIAJGQqHRue5vW8LkVxsu54X7UejQ+Ww+nvpXH8OjVFllJXBt9+C7efVJn/qhIiRo6FBS1JNSXF/+34Vd6NXrdTGHEJtzOFDOjZfKUp4vO+FVSBy880309TU1P1VWhrYvn5f6fWQm6s9oT74ANav16JbIcKBu8THQOzYoc10OMwdfFOxG4ATJ05G6fEJraREe+6fMm81+287DZ2tlaa0BVxi+wOd+Net190uEE+Js2lxlpCUN4+JgUWLIC82nTPS+ybpgtZI73b7RdxhuwCnqnCWsoLZxUvpiNvNli2wdq32Ii51QkanXbu0HTOumbGeoju2kVH3KgBbx94e1rMhEPpgJKyWZsxmM2ZPG7GDRK+HCRO0jolffaV9ypk/333ikRBDZbAapLW0aMsySUnwedku7KqTnNhEpsTvrbxUXg4GA5wwdwPjPzkNrK2QeSjRZ79GwoPfe1xS6Y+7GY5wy5HIyNBmRto/S+f6KWlUOetp6Org3aKttNr3JrK+4FhCrRrHg8Z/cpJ+FUllF7A7azk//BBNQoLUCRmNOjth82Ztp4y7GGNK+V9QcFIZdwKN0QcM+fiGG3m7/VliojZ1vWoVbNgQ6tGI0cyV+LjvdL8r8XHFpgqfz1VUpM34WWJtrKrSZhiPGZ/TPRvS1ARdXbB4biHjvzgdrC0wcRGc8zp6c1T3koqvkqJM/c5weFouCZUZM2DePCgvU5hgSSLBHNErCHF536nNDrWrZhYq+Rxe/gviTPV8+y1MiJA6IaNNURFUV7tflolrW0d64/9QUdg27k9DPzg/hMtETVADkdbWVvLz88nPzwegqKiI/Px8SkpKgnnZgMXFQUKCFozs3Bnq0YjRaDAbpFmt2qe2mBj4vrqELqeD9MiY7tkQm03LCVkwrYKsladAey2kzYRzXgWTVra8e0kl1vtMZWKUkVU3HzWstqLqdHDQQTBt2s9BW4fn5ZOvnbM413oLjWoUY7s28MumE3A2V/PN1wrXHyV1QkYLm03LL4yM1GbU9zW17M8A7Ek8i5YI/wL5UAn10kxQA5G1a9cyZ84c5syZA8C1117LnDlzuO2224J52QFJStKeXCtXam2dhRhKg5n4uHu3FmgkJjn5qqIYgCPHZnbPhuzZA/tNaGTettNQGksgIRPOexsssb3yU+IiTHx9w5Fcs3g/t9dxFSD726kzMBmG3ySr2QyHHqrtnmup6X/5ZIOayxnW22nRpxLbuYVT2o9lc+E2jNVp/OOsvjkwCVFGLl44ibgIk3TXHSF279b+dlJT+96W3Pw5KS1f4lSMbB9785CPbbgKao7I4YcfjhrqUCsA48ZpMyIrV8Jxx4ExsKacQvhtsBIfnU6t0JLBAJuaK2i0dhJjNDMvZSygJdkZFCvHNpyHrmYzRKfC+e9C9Jh+81OeGEFFuRxOtVe+yqJDE9har2JuMdKlei7AVsQEfmW7jSeVWxlr38Wvnaew9Lu7OHHWIr6+/kjW7q7nky2V/Ce/nPo2K8+uLObZlcUB5fiI8OB6rlQ2drJzowWTPhGTaZ8ZLtVB3p5bAShOWUqHedLQDzQA4fAWHVbJquFCUWDSJC0YGTdOW0MWYigMRuKjw6my/Id6VmzpZFySmc/KdgFwSNpEilsaaerqpKnSzPVJd2Cp/AZMMdpMSGKmT4W5vr3xyD4JpwCrCuvCIgnVF+6CLZ0Cvkxa2FUnG60JnK3cyqumv5Ctq+A5w584Z8OfsFgWk50Nz68s9lrcTAwP7p4rcUYLp8flMTtp7+9xfN1rxHVswqaPY0f6jaEYqt9cOSKhDkYkEPHAZNISWNeu1bLrx4wJ9YjEaODqleJP8a+e+rxoNmr/0SkK31bu5oPSHQBcol/OONvLqOhQTn8O0mb4VZirZ1GuwdrhM1Q8BVvegpA4o5l2uw2bqu3xL1FTOct6K6+Z/kKWrpLXTXdxyWojjT+m+FzcTIQ3T8+VJlsnz25bz9Ipc5mdlI7e2c6U8rsA2JH2B2yG4ZOYHOogBGTXTL+Sk7X6Cz/8oPUUECLYBlL8y9NuGwCnqtJs06oWH6Vbxy2GVwC4y3YuK6xa591A8lMGc4fPUOgv2PIkUm/kymkHcv5+s7uDEJc96hjO6rqVEmcKk3RVPGu4FX17ucdzSXGz8OGtTo8vz5W3d23BqapkVT1GhK2CdtMEisf8JrgDD4JQByMSiHiRkaEVhCooCPVIxGgRSPEvX99gpyglPGx8DJ2i8rL9KJ53LOneheNvfspg7vAZKt6CLXfaHTYUFFpsXW5vLyOFs6y3sds5hom6al4x/ZVkmvo9pxQ3C60VmypYdM/nnPP091z1Wj7nPP09i+75nBWbKroDlIc+2e71udJo7aSifjs5lQ8BsHXc7Th1w6deTLhs35WlGS8sFq0Pzdq1WuEz6UkjhoKvxb9cSXQrC2q8vmjG08LTxgeIVjr51jGN2+0XoqJ0f0L3Nz/FnxmUcOmvEmgAUFTRSdZYz49PBUmca/sjr5v+TLaugn+b/sbZ1ltpItrt8cW17QGNQwxcf3lQl720nvhII43tvneXnlf7wM+N7eZRnnDa4A52iIR6RkQCER+kp2uzIlu2SOKqGDreGqS5y83wRIeTh42PkaGrodiZyu9sV2Hv8edf3dLJL2aO9Ss/ZTiWNg+0ummM0UJUeyLxJkufJnkue9QUzrXewhumu5iqK+Vfprs5z3oLLUT2Ofbvn+5gclp0WObQjGS+zOL5E4TkKHuY1/oWAFvG/2XIpxisDgetdittNittdittNhutdisddhtdDjtdDgdWpwOrw0GX047d6URVVZxoy7UGRcf+xqmAmzr1Q0gCER/o9Vq+SH4+5OS47y0gxFDy9KnOk6sMb3OY/ic6VBOX2a6heZ9P6mNiLN35KZe/tB4Fep3bXX7KcCxtPj8zkcQoI/Vtvr3ZuIKvXx2dyOefKSxJzeO10vUejy9W0znXeguvm/7MLN0unjPdywXWm+ig72MgSatDL5ClOc9U7jC9hg4nFfG/oD7m4EE6714ddhs1HW1Ud7ZR19lOo7WThq7On//bQbvd96DJk4iIUiQQGSaSk2H7dti0CRYuDPVoxGjmb8LlUbp1XGV4F4CbbZewTZ3Qfdu+sxyu/BRfaoUMdIdPKOh1CqfOHsezK4t9vs/tJ+YxfZpCWyt8+20652fO5X9lWzzOjOxUx3O+9RZeNf2FA3Q7eMb4ABfbrqcLU/cx4bhsNRoM5uzckboNLFLW41SMbB335wGdq8Nuo6ytmT1tzZS3tVDd2UZ1R5vHvKSe9IpClMFElNFEtMFIlNFEpMGIWW/ApNNj0usx6/SY9AYMig6doqBTFH6sq2RdbTlOnLI0M1woiraFd/NmmDJFq8AqRCj486luolLJQ8bHAXjBfgz/cS7qvs3TLhxf81P8nUEJF4vz0nwKRJKiTPz11Ondwdf++2u9eX76KZ3b5qTxbVUR7xRvdXvfzeokLrTeyL9Ny1io38w/eZjf2q7ptRwGQ79stW8Rt3Cv9zLYBmt2LsWk537TG2CHwjFX0GbJ9vm+XQ47xS2NFLU0UNraRFlbM3VdHR6PjzWaGRMRRZIlkgSThXhzBPFmCwkm7b8RekOvTtq+qu5oA7QlmlCTQMQPCQlarsjmzVpJaCFCwdc3LzNWHjc+TKzSzjpnLn+1n9fr9v4qonrLT3HxZwYlXHibyYG9fXN6lqw3GLTZ0KYmKNmtEOOl/84GNZeLrdfzL9PdHKXfwL08xXW2y1B7bFYcymWr4VbvJRh8+d33x6zXc+mU/Tmu498klpfQaUxjZ/of+r1Ps7WTguZ6djU3sKu5gbK2Zpxurp5ojmBcVCzjImNIjYxmTEQ0KZZIIgzBKe2t+zl4UZFk1WHFNSuydavWJEtmRUQo+PrmdYvhZfJ0u6lVY/md9SpsPf7cE6OMfHX9EYPSG8bXGZRw4ctMjqe+OdHR2oeQ99+Hmhbvv4cf1Klcbruap4wP8kv9tzSpUdxpvwAFZUiXrXypmDsagpH+fve+6HI4iHFUM7nyPgC2jrsThz6m1zE2p4NdzQ1sa6xha2MtZW3Nfc6TYLKQGZvIxOg4xkfFMi4qliijqc9xweR6rqsBhWSDSwIRPyUkQFUVbNsmuSIiNHz5VHesbg0XGj4B4A+2y6ii9xtefZuNdbsbBi0/wdcZlHAxkJmctDRYtAhaViQSqTfS7ug/YfAL5xz+YPstD5v+yUWGj6hXY3jU8cshW7byp2JuuAaPg8nT795sUOiye39TPqj6bxicbdRHzWdP4pkAtFi7+Km+io31VexsqsPq7F0Bc3xULFmxCWTFJJIVm0CCOfR1ILpnRFRVZkSGo5QUbSvvtGkQHx/q0YjRxtunurHUcq/xSQCetJ/Al87Zbs+zsqAmrGcugm0gMzmTJ8OKTZW0V/q2a+E95yLibW3caXyR64xvsXBKHgdNP2GgP4JPhmO9l2Dr+bvv2aDQm3nKdma1/w8VhVVpd/FlxW5+rKuksLm+199hrNHMlPhkpiakMDk+mRhj/8t4oaD0WJoJNQlEApCYqM2IbN8OBx4Y6tGIkc5dgqHrU93Nb22hoXPvm4weB09YniCOdvKdWdxvP8vjeR/9opC315eNqhyBfQU6k+NUVV7bscWv+7zJLzg2IY6DG/7BgYV30fzdRGIPPt3va/trONZ7GQp6nUJTh9Vtg0K3x+Pgz8YXAPhQdyxXbq5Dpa779oyoOGYlpTEtcQzjImMCSiAdSq6FR3f5KkNNApEAKIq2nXfzZpg+HaKiQj0iMVL1l2A4KymdS8akUaerp05t5L3d27ja8DYz2UKLGsn/2f6vV16IO6MtR2Cw/FBUT2Wz/6Xiv0+8ijRDJ1k1TxH98W/pjI3DMv3oXscN9s6W4VjvZSj4uw3+1/oVTNPtplGN4tb2k1GBzJgEZielMTspjURL38J14UyRpZnhLylJ20FTWAgzZ4Z6NGIk8pZgeMnUuUR0pjMnN4nnt5dwkG4LVxjeA+CVmJsorU31eo3RmCMwGAKdPWiydbE54x6MtnoyGt9C//Z5vLTlEcpjZ5IQaaK8qYP39lkmiI8wctHCSVx5ZG5Avx9fdwnNm5gQ0M80XPmzDX4stVxr0CqoLrP/ioykyVw5aSrJwyz46EmHa2km9DMi0vQuQDodxMbCxo1g9b60KIRfvCUYqsCr27cwdpxKeVszBbUFPGB8HB0q31lOYlltjs/Xko6w/nE4VWpbvBeacifOZAFFx7/ibudr5xyMaicnbr6Kz7/6gr9+sJXnVxb3yVVo7LDx0Kc7mfeXTwLqZtxfR2eX+jYbh933Rdh1Sw4Wh1Plmx21Ph9/h/FFopQu1jj34w3HYRyaPmlYByEQXjkiEogMQEqKtoNm9+5Qj0SMNL58Wmt1dlJuq2d5yQ5uN77IOKWOUjWVSxtPDuiaoy1HIBCurq13LXdfyMybDruN/LoKnt6+kd9Yf88a537EKe38y3Q3GUpVv/dtbLdx2UvrefjTHR5b13viqaNzT66ZtpEejLh+h//8yreW6sfo1nCMfh02Vc8ttkuIM0WSHRs+1YID5ZpcU5GlmWHNaNSKHG3eDNnZ2iyJEIPB16BgV3M96Y3/4zTTNzhUhautl9FGYFsDR1uOgL/87e/jzr92/ojx50+inZhZav0Dr//cJO8l4zJOt95ODf0vkTz06c7u//enINmS6ekcOSWVg5Z95naHyEhcpts336ahzcoVr/j+O4yigzuMLwLwlOMEdqrjWZqV1731dThTwmhpRgKRAUpLg5ISqKiAceNCPRoxUvgaFOyp28pzxmcAeMJxIuvUyX5fKxx7woQbfxMbe4o3WTglcwrfVpRQ0FxPz0WdZqK5wHoTb5nuZKKumn+b7uYs66007dOU0BN/k43X7W7od5vqSNrK6y7R298iZtca3mKsUk+JM4V/685h6ZQ5zE4aGUnd4ZSsKp/hBygiAux2LXFViMHiSjDs73NXlMHA760PkKi0stk5kb/b/d8KGs49YcJJIF1bZ0ZM5Iq8AzkvdxaqCoePzSTKTbnuGhI4z3YzVWo8U3SlPG+6l0h8u5br/ePO/23xaZlmtGzldc1e7fs78+f9dq6yg4v0KwBYO+5ebt7/uBEThEDPZNXQk0BkECQnw86d0NgY6pGIkcKXBMPTlU85Up+PFSNX267wulXXnbQ4i2zd9UEgb8yJkWZe3P4jj25ezYs78nlm2zocHj56lqqpnG+9mUY1irm6Ap4wPoQJ34ql+ZNsPBq28g5k9srFjJX7jE+iU1T+qx6BPv3UEbEc05MSRjkiEogMgvh4rRFWUVGoRyJGEk8JhgowQaniBp22dn2P7Sx2quN9Pu81i3N5+OzZvHrpQXx745EShPjAnzdmBYiPNPJl3Q5aHb0DmE6H3eP9dqgZ/Np6I22qmUP1G3nY+Ch6HB6P35cvwZK3mTYFLe8kHJbpHE6VVYV1fifmBjJ7ta9rDG+RraugSo3nT13nUNg88naU7W16F/o5EckRGQSKom3ldZV9Nw1t7yIxgrlKUT/9n3o++LGSn7qKASf3Gp8iSulilSOP5xxLfDrXaOu0Oph87doaSCO1nvLVHC61Xcfzxns5Tr+GZeoz3GD/DZ7nxfbyJVjypeFfOCzTDaRT8ECXleYoO7lUvxyAW2xLaSaa7Y01NFk7iTNZyI5NHBGzI66lGWeop0OQGZFB49rKW1oa6pGIkaa8TMFWnkixsxKAc/WfcZBuK22qmevtv+nVVt6TaxbnyuzHAPiyVAbaUtc1i3NpbPe+rKLzcKbvnNO51n4VDlXhTMNX/NHwMv2FN/7OYniaaQuXZTpP+R2+bi8eyLKSa0lGr6i87VjEZ855AHy0p5AXd+Tzj03fc/vaz8mvG/5bnHuWoA91LCIzIoPEaNS27+7YoW3lFWIw2GywZg0Ut9XTbOtkvFLDzYZXALjHfjZ71DH93j8pysRfT50e8jeXkcBT19bEKCOnzh7H4rw05mcm8v5P5T6d74ysaXy8p4AGaydRBiOHpU8i2RJFvNlCduzxvL87iZPrbuFSwwc0qVE86ji1zzlcbyW3njDVr7LwA2n4F0y+dAq+47+bibEYqW3t6jNuh1PF6VSJjzDS2OFbjk1P1xteJ0dXTrUaz59tF7g9ptHaybPb1rN0ytxhnbzaPSMiSzMjS0oKFBdDba2WwCrEQG3fruUeVRsrAZW/GZ4hSulitXMK/3Yc7fX+fzphqgQhg8iXN3BfP5GnRkZz1YwFPLTxO5qsXWxpqOHK6VmY9QacqsrtNVP40XE+txn/zR+Mb9JEFP92HNPrHGlxFk6alc5dy7f6vYyxb8M/V05GKAMTXzoFVzZ3ce4zq7u/5/pZgT5BojdGRYde0dHptHOI7icuMXwIwE22S7xuoX571xZmJqYN22WacEpWlUBkEMXGQnm5FoxIICIGwuFU+WpLPR983kmVrZ1vGos5S/8lh+o30qkaudF2qU9LMmlxgRU3E55569jrSz5JlMFIZkwCBp2OK6YdyN83rqK4tZGnt67jt3n7U9zSSKO1k+c4jjiljasM73CX8QWa1Ujecy7qPs9Z+2fw8Gc7PfYj8nWpZSA5GYMpkPyOyqZOLntpfUDXu2DybGYmprFy9xqur30CgBftR/O5c67X+zZaO9nZVMfk+OH5Yq+EUbKq5IgMsvh42LZN+s+IwLlKUF/80ve8VZ7PN407SKOOPxpeAuA++5kUq/2/OYTT7ofRxpd8kja7jTvXfUF+XQXpkTFcnjcfk07P9qZantm2jrrO9u5jH7KfxvP2YwF4wPgER+nWdd/2/HfuW9i7+hHd8d/NXnebDDQnYzAFkt/hy9vovr8Hg6Lj4slacTIdcK39QVKURiqUXLak3MNZ033r1fTc9vXDNl+ku46IGvocEQlEBllSElRXS9Kq8Ky/bYnu3xRUlhmfIVbpYL0zh+cdx/V7/nDa/TBa+dLbxZVrkF9XwaSYeH4zdX+MOh1bGmr4orxnLQCFP9vP523HIgyKk38aH+Yw3Y8ANHnJg6hs7uLRzz33VPGWkwG+F0sbDL4U8guECkz9eebCrNdz69zDmZM8FoCJtc+S1vQhdkxs3e8Fbjw/m1MW+DbL0W63df8OhxuZERnBXEmrBb71UxKjjGu245ynv+eq1/I55+nvWXTP56zYVIHV7uSWdzf1eVn4pe4bjtD/SJdq5Hrbb3F6+bNNiDLy2K/mSG5IiC2Zns5X1x9BYlT/+/nf3rUFp6oyOT6ZK6YdiEVvoLy9pdeuGhUdN9h+yweO+ZgVO08ZH2SRbqNP43jo0x0eZzV8yckYys7Mvu5OCsS2Rq3b7tnZM0i0aEuWse0bmVb6R+328Xdy4EnTSUz0PyBy/Q6HE9fPJoHICJWcrOWJ1I+8GjhiAPqbAr/spfXMuevjPn1AUmjkduO/APi7/TQKVfcNjXq+YNa32bhr+dYR30V1OPDW2wW0mRFXwazs2ESOzchBoe9uBgd6fm+7ko8d8zArNp4x3s8C3WafxnHtGz+ysqC2z8zGp1sqfbr/UJZ892U2KRAqMDspjVijmbU1ZZTUFzGv8Hz0aicVMccw6azLiIrSju0ZEPmi5+9wuOhZ0CzUMZQEIkEQGwvNzVozPCHAtynwtq6+VTRvNf6bOKWdjc5JPOU4weP5PSUrSjASWr6+gTdZtePy6yp4r3ibx8+odgxcYbuKTx1zsCg2njPex4HKVq/nb7c6OPeZ1d2zb6A9J9/NL/NpfMEq+e5pmXLJ9HS+vfFIXr30IB4+ezYvLz2QtNiBLdlE6I3sam7gkc2r+deO9exXsJRoaxH1unGYz36S6Jjeb4eugCg+om9/IHdcv8PhYm8gEnoSiASBokBMDGzdqjXEE6Obw6nywsoiv8tOH6b7kZP0q3CoCjfZLsWB3uf7hmJ9X/Tl6xu40mnBqaq8vWuL12NtGPid7Wq+cMwiQrHynOle9le2+XSdngHqD0X11Ld5r7WRFGUKStJzf8uUsHd30smzx7EwN5k7ThrYkk2Hw0azTet9/Fv9+xytX0eXauCCjsv5oaHL7X2WTE/nsXO976ABiDMNr/48rscxHJaUJBAJkuRkrdJqhXwgHdVcL7Z3Lff+qbUnC13cZXgOgOcdS9isZvp97aFe3xd9+dLbJSnCQnRHIhsr62n08VO1FSOX2a7ha8cMopQuXjDdyzxlu9f79QxQK5t9u9bJs8cOetJzIDt1PC3Z+DI0k25vEL9At5nrDa8DcLv912xWs/oN2A/KSvKaLxL/c+n34WRv/RNZmhmxzGZwOmHXrlCPRASDLw25PL3Y+uIqwztM0NVQpibxoP2MAY11uLd0H876S750/fuuU/M4cL5CaY1/v6cuTPzGdi0rHdOIVjr5l+lun3JGXAFqfav7WYB9HZ2X1ud7ruf/u+v38Ow3u3h3g++N6QayU2ffJZtXLz2Ivxzf/4xFpMGI1akte05SKvin8WH0isqb9kN5zXGE14DdlwTa07Lyhl1hM0Uqq44OiYlQWAj77093EpQY/nwp/jSQVuSTlRIu0X8AwG22X9OOhWijiUPTJvJB6U6/zzecW7qPBJ5Kw6f1eM7Y7bCxysJHa/07dydmltr+wFM8yKH6jbxgvJfLbFfzhXOO1/s2VptIjbFQ3eK58Jq7WjTunv89j/dWBM2fnTruCsf1LCi3ezfs2AUnxM1lZceWPjNKCnBY+iQ+LN1JLK08a7yfBKWVfGc2f7JfTM/Qor+A3dPvMCXKwoHGPPaLGH471HrmiIR6RkQCkSCKj9e28e7ZA5Mnh3o0YjC4Zjm8VbL8flddQDMhCk6WGZ/BqDj40HEAnznnEWUwctf+R9HYoPCNvpQWh2/nVdDe7KSoWeh5Kw1vMMBFv0jk6Y0WGrp8f94YFIVO1cyltut4lEc4Wr+OJ40PcpXtSj50HtjvfUu3R3BkfB6vtvjeidfT89+lwoeKrr7O0PV3nNMJmzfDypVaP6Zjp6VzDGkUNtdT0FTHitICnKgcnDqB7NhEDNj5p/FhsnUVlKlJXGq9ji56b6v2FrC7+x3uPzGRTz9R2LoVcnN9+rHChmzfHSX0ejCZYOfO0EecYuB8nVL+4KcKrng5sJLTv9J/zlxdAS1qBHfYLgTg7JwZqE6F/LJ6Fk3Spsm9TQJLUbPw0zP5ckF2Up/fS4RF4c+n+rZl1KjTXrrtqopeUZgQm8bLiQ/ynuNgTIqDR43/4DTd1x7vH2+ycOjURCbo0vlF/Fzizd478fo6y6cCN7290e12YfB9hs7Tca2t8MUX8MknWt2mSZO0DQI6RaGus707CAFYWVXCv3ds4G7TiyzSb6ZNNbPUej01xHefz58qxPv+Do0GhdmzwWLRdkoOJ90zImHw5iQzIkGWnKxVWa2v16quiuHL1ynl370SWBCSQgM3Gl4D4H77mXQY01maPQ2A29Z+TqujE1q0YxWld3CrU6Dna35aCPqEiIE7cXY6CnO59d0tvWZG4oxmFqZNICUiirifEyNLWht5p2grRS0N7GyuZ2czbDJfS7vdzDmGL3jA9ARJtiaecvyCfUPXg1MnYDQojB8PCQnpjN+dhjqhnnE5neSOtzBvYgLrdjfwXn5Z9+yNt+d/T40dNs59ZrXbpRpvvXj2nclzONXuWQhnm4X23YmU7VGYMKH3kvfamjJeLvipz/l+o/6L0/Wf4VAVfm+7km3qhF7XgoEF7OnpMHUqrFunlW4YLpQwqiMigUiQxcRAWZkWjEggMrwFO+nzNuO/iVXa+dGZRX7sudw1dT4/1Vfy7La+gY0r6Fi6cBKL89K63zjCqaW7CMwvZmtLAP/6sJ51mztJjbMwNyOxTzLkpJgErpmxgF0tDXxRXsRPdZXUdlm5mUtoJYJLDR9wi/FV0pV67rKf36sib0pEZPf/R0XBlMkKRUVJ2Athh72Ca9/4ok8O1PHT+yateuOu+Z4r+fPyl7wvCbnLR4nWWTgjJ4+oqJ6zNU5ecROEXKJfzhWG/wJwh/035JsXQof7PJ2BmD5d65Td2KgtyQ8HrmTV0M+HDFEg8thjj3HfffdRWVnJrFmzeOSRR5g/f/5QXDosxMRoT9IZM7TlGjE8BTPp83BdPifqv8ehKtxiu4Q9rS3k11XwbpHnbb8K8MGmSm45Ic9rR1gxvJhNCktPTGJeBnz7LVRVap+896UoCtmxiWTHJtLQ1cH7u7fzQ00Zf7WfR6WayK3Gl7jI8BFpSj1X267ozovYt+aFXg/Z2fDZ9gre29E38K1s6uTZlcV+/xwq2vP0zv9t4ei8tO7g2JcEXk/5KK3OTp7fsR6dbi6zk9JRVZXntq/H5nT2Ou50/Vf8yfgyAHfbzubfjsP499mzMBiUQQ/YU1K0WZE1a4ZPIBJOlVWDHoi8/vrrXHvttTzxxBMceOCB/P3vf+fYY49l+/btjBkzJtiXDwvJyVBeDpWVMM59hW4xDPjS3t0XkSY97da9VVQj6OQuw/MAPOc4js3qJLB18fz2Df2ex9vuAjG86XQwdy5ERsJXX2ltIyZO1Jbl3EkwR3Bu7iy2N9bSZOviWcfxVKqJPGj8J8fp15Cs/I3LrNdgM44hMyaBnU11NFk7u5d6AL5qcV9QzfV833cJ0Bc9n6euJR5XIPDV9UewbncDlU0d1LdZSYw2Exdhwmp3csd/+89Hea1gI9Pix/De7m38VF/V67bT9V9xr+EpAJ60n8ATjhMBqO/o4uTZwXkRzsvTilgOl1mR7mTVUEchDEEg8uCDD3LppZdy0UUXAfDEE0+wfPlynnvuOW666aZgXz4smM1gtWovJBKIDF/eppR9/XN++vz9KSixce+nm2lTu7jK8C4Zuhr2qMk8ZD/d73FJnZCRbcoULRj54gttF15mprbLxh2donB69rTu5bzlzoOotcbxtOkBDtDt4D3zrfzWei03rLb1mkGIN1k4ODXDa0G1gRTp/WRLJde+kd9nyeekWen898eKXt+Psxhp6uy/6mub3cYNqz/GrvaeCTlb/zl/MzyLTlF5yX4Uy+y/wvW2G8xZzZQU7Xe1du3wCETCaUYkqLtmrFYr69atY/HixXsvqNOxePFiVq1a1ef4rq4umpube32NFAkJ2otIp7xnDGueqjumxVn456/meK2imR5nYUZaEqbqdI5LnM0UpYRL9MuBvTVD/CV1Qka+CRPghBO0GZGdO6G93fOxs5PSWTplLvE/L7+sVqdyqvVOitV0xiu1vGW6gyV82+s+jdZOn2vUXLxwEukBNKR7bmVxn2TXiqZOnvy6b/sDb0GIiysIWZiaQbzJwnn6T7jb+Aw6ReV5+7HdtUL82RkzEFOnajk3w+Gty5Ws6vRy3FAI6oxIbW0tDoeD1NTUXt9PTU1l27a+vRGWLVvGnXfeGcwhhUxiIhQVaUs0WVmhHo0YiP5qQuh0Sr9JeLf9Io8N6xX27IHI+A6WtTyDQXHygWM+nzt962nR85xSJ2T0SE6G446Dld+pfLi2HoepkwkpFnLi+iayzk5KZ2aiVldDW345iE0RZ7Dnp1NZpGzgEdOjTLPv5j77mX71MAKt0uofT8jjh6J6Kps6qG3t4tEvCmnq8Bw8BLKk449N9VU8kvg5hzZpS5xP2U/gbz/PhAzlVvbUVK2eyE8/hf8OGl2PSiKhFla7Zm6++Wauvfba7n83NzeTkZERwhENHoNBW9stKpJAZCTwlBzqLQkvx5LOBz/C2LGQ2vgOc36uGXKn7YJ+r+drwSkxsn29q4I7fujRJ6YaYg3aLpLZSb2zWXWKQm7c3ufojqY6Hu26jusNr3O54X9cZvgfc3U7uNp6BeUk+3R916zCvs//jMRILn9JWw5yt2QZzCDEiJ3beJxDm74B4En1DJbZT8H1V+L62zs6L41VhXVB31mWlwfbtkFbW3hX1HbFrs4wWJoJaiCSnJyMXq+nqqp3IlFVVRVpaX23gpnNZsxmczCHFFJJSVog0tKi7aQRI5OnGZOWZoX//lcrwpQSUcmhhfcBcJ/9TKpwP6uRHmfh1hPyuGu5590FYvjpWRvD1zdFT7tImu2dPLttPWeOm8vCCenoPCy4N1k7caLjHvs5bHJmcrfxaebrtvOB+WZutP2Gj5wHeB332fvl4XQo6Pe5Rn8B+PHT0wLaceOLJJp41PgIC/RbsKs67nJeSlfa1ZzWaiZrWidzp2qP7SdbKll0z+f9tmUYLOnpWh7Pzp3aTqRwpeuxiKwlrIbuA01QAxGTycS8efP47LPPOOWUUwBwOp189tlnXHnllcG8dFiKi9PyRMrKtKQmMXLt+4nRZtO2YVZXa1O304pvwuRsId+ZxUuOo92eQ4HuF8pjp3suDy6GF196Fe3Ll6qm/6vYSGezkckJSaSOUdj3M13PLbvLnQex0ZrJP4yPMFu3iydND/GK/Qj+Zj+XViJxx6TTs20bfBoJBx+svZ715CkA/6GofkCBiB4Fh5uffLZSwD9Nf2esUk+bauYK2+/50jkHSjdw1by5XH78OHQ639syDBZFgWnTtECkq4s+v4dw0XM5L5gzVr4Ieon3a6+9lqeffpoXX3yRrVu3cvnll9PW1ta9i2Y00em0T8MFBVLyfTRRVS2TfutWrRx1WsvHjGt4F7uq4xbbJZwwcWp3YqFLaoylTxGo/sqDi+HBU0dm15viik0Vbu/nS1XTDqeN/7as5qmKz/myoIKdO6GhYe9rTXZsYq/nWYmayhnWO3jCrm1t/ZXhCz4238AROvfbxq1OBx+2rOe/6yv4z3+05QeHo/cx7p6n8zMTiY809jt2T/RK3yBEwclS/XLeMN3JWKWeQmc6J1n/ogUhP3ujYAsq6oA6/Q5ERoaWYFxZOainHVRKr0AktG9IQc8ROeuss6ipqeG2226jsrKS2bNns2LFij4JrKNFUpLWBK+hQUtgFSPfjh1aoaP0dIg0tDF9t5YH9ZzjOFLSDueY8TksHpfNzsZ6tpd0csB0C0tPkhmPkcbbm6K7wl8u/mzRbrZ18n7jeq6cPZeulnS2b9eWgseMUTgtK69XpV4bBu62n8NXzpksMzzDJF0Vz5vu4z3HwSyznUMlffOgVnVtYXpbGh9+qFBQADNnwvjxeFwSWrGxksZ233bB7MuhqkQajBw4ZjxflBcxjhoeMD3BQTqt0N8KxwH8wfbbPrM4rrolrv/3JFi1ePR6bVakqAjsds/brUOp5zNsxAciAFdeeeWoXIpxJyZG2zmzZ48EIqNBZSV88402PRsXB/vtuZsoWyl71GRe0J/HtRO1tsw6RSHemsSCsXD6kfRZgxfDn6+9ity9KfqzRdsV1LxdtIXll6VRWqKwZQvs3g3RpPOrCXP5oHJLr5ohPyqzWWK9m2sNb7FU/wEn67/jaN06nrCfyFOOE+hk7/pCo7WTrth6MsxJFBTArl1afaTsbO05bjJp3XHb26GqxskNH27qd7wGRUeU0USTmxomE6Pjta3IRj1HtL3EqZ3PEaN00Kaa+av9PF5xHImn3AZ/grdg1OKZNEnbRVNdrSWnh5ueSzODPSPkrzCM00Y2RdEyqXfs0PoTePoUIYa/lhb4+mutW2h2NsS2bySr6jFAqxly8pT9Meu1P0G7XZslW7w4/Lf9icD4+mbn7jh/q/q6gprt9fXMn5FIS0Q9bbs76ai3ML4+jV/FpNEeVY9q0SqrOlWVRzev5m/2c/mvYwG3G//FAbodXGt8i3MMn/OE/URecxzRXSK+ydpJbpz2vO7ogKoq7dO/omivaaoKOzsr+KJlI51q/7MhdtWJzW4nJzaRwub6Xj/f7tYGvvnxEW7SvUC2UgYKrHHux3W2yylR+59V9yd4C0YtHrNZa+vx8ceQlhZ+r/UK4ZMjIoFICCQlaX+4NTVaxCxGHldyakkJ5OQAqoPpu69Ch4Pljvl0pJxMXsLeFgfl5dr09tSpoRuzCC5f3+zcHddfVd/+uKtmOibawq+m5BFXl05LPYzLAEuESrzJQqO1k01qFmdYb+cE3WpuNr7CeKWWO40vcqXhPzxtP57XHEf0SnyNiNByIlwcDtjYUMGH233vQt3utFPQXN/jOyqH6DZyteFt5um0Qmt1agx328/hLcehqP2kN+5bX8efTr+DLStLm/luaAi/pqc9S8+EekYkzGK00SEyUpu2LC0N9UjEYHA4VVYV1vFefhmrCuuwO1R++AE2b9amZ/V6mFjzHEnt62hRI3hEfxmnZeZ137+jQ8uu339/sEiR1BHLNavRX+ZPYpSRyuZOVhXW9Xlz8FTVtz/uqpnWtHby8Nr1JMyoYObMn/tgVSi9npOgsNx5EEd13c8fbRezR00mRWniFuOr/GC+kmMrr6Om9D0KGqv75BcoOpW3i9z3rPEmljYu0H/ECtNN/Nt0N/N0O+lUjTxrP44juh7gTcfhXoMQ2FtfxxXA9bzN07HBEBOj7ZCsqQnK6Qek59JMqPvNyIxIiMTGatu7Zs8Oz0Qm4Rt3WzGTIy0cZMpjQUY6ZjOYrRXst+cOAO63n83x047sXpIBLV9o8mQpdDfS+TKrUd9m45rX8wH3W3pdW2S/L6zjilfW0xhANVNXDslDX23h6+vTmDBBYdUqiKxP59c5c3mjeCPtdu28XZh42bGY1x2Hc6r+Wy7Wr2CqroRpre8yrfVdaqri+IKDUFNPIi79eGyGRHY21XntWdNTGnUcrv+Rxbp1LNJtwqJo1+5QTbziOIon7L+ghgTfztXjMXPVaumyO7l68X68+kPJ3kJwDF0tHlel1XCrH9V7aUYCkVEpKUmrJ1JVJY3whitP9Qlq2zt5v309qWlzmU06WcXXYVFbyXdmUz/2MvaP2fui2tCgzZDNm6fNnIiRzVPhL3c81bnQ6xQOyk7iooWZPPTpjj7386WaqSuHZO3uehZMSSI5WUuq3rkzndtmpPF1bQFfVhR1ByR2DLzlOJw3HYcxRyngHP3nLNGvIUVp4kw+guqPoBpqjNkYu8ZwrSGVUnUMTWoUzUThVBWMip0IrKQp9YxV6shWypmlKyRNaeg1tu3O8bziOIp3HYtoxntp0iuPyCY3NaZXfR13HxDSYs1csziXSclRQ1qLJyVF+5CxZUt4BSLhVEdEApEQMZu1PIKSEglEwlV/1S99KTD19q4tHKKuI6dlOXZVx7MRN7A4Y/Le8zu0QPSQQ7RkNm/XFCNDz8JflU0d3LV8K/Vt1j7HedrS6+5Ntid/qpm6EmOTk2HJEi2R/scfFQ4Zn8uxGTl8VLo3INGe6wob1Fw22HP5o30pC3SbOUa3lgN128jVlZFiK2SJrpAlfiz621Udm9RMPnHM41PnXLarGfhT5XNhTkqvXUaePiBUNXfx90938vh5cwd1q643iqLNeG7dGn4FznQoOFFlRmQ0i4/XlmfmzdO2vInw4a36pS8FprqsjUzZfSMArygnsWjKab0+hZSVaQmqM2f6dk0xcrgKf60qrHMbhLjsu6XX05usyzWLc7nyyFyfq5n2TIyNiIAjjtDeKNesgfqoSj4o7Tvj4mLDwNfOWXztnAVAIs3M1BWSrZSTrVSQrtQRq7QTRxsKKlYMdGGkWk2gTE1mj5rCj84sNquT6Aig67S7RNOB1GoJpvHjtQ+c1dW9E3tDTfl5+izUyaqjMxBRVRLW30Wy41RgRsiGkZio7e0vL9eSGkV48KUkdJfde/Psqw1vk0YNe9RkbJPvIcGwt7pka6u2ZXf+fG1pZqjLUIvw4M+WXm+zcArw2ppSrjwy1+t2X0+7RYxGWLRISzj93Yf+JZzWE8uXzjl8yRzvBw+Qp0TTgdRqCSaDQWuG9+GH2kxouCzD6hQFhxr6GZHRuWvm+38Sv/EBzmg/kdj2H0M2DKNR22+/e3fIhiD24WtJ6OTo/udX85Riluo/BOCLlD8zJmZvk0enU0tQnTlTWzsOVRlqEXr+bOn15012ILtF9HrQp9XT6hz8Il+BUvYZZlqcxW1wPpBaLcGWmanlBtbXez92qLgSVh3eP1cF1egMRGafS2fyPCLUBhbsOIm4Nve9FYZCQgIUFmpbOEXo+fpij4rHrZg6nCwzPoNBcfK+40AeqE4jv07rIeJUVb4rqKM6ogw1pQ6nqvr1BiNGFm9behW059n8zES/32Q9bff19CbeU21b+AQhsLdfztKFk3j10oP49sYj3Y5/ILVagi0yUqsTFE6BiGupWLbvhkJEPJVHv4vpjdMZ6/iBg3aezPe579IUNW/Ih+IKRMrKfi58JULK1xf72rYu/nhcHle+1rdo0wX6j5ml20WzGsmdtgtoRGvTftS4TNZUV9Bs067x5kvam8zx09P6nGMgYxPDR39beveduQjkTdZTR1xv+RGBvFFHGYy02X3rKXN4+iTy6yp7bfONUsygd3o8hwJ8sKmSW07wXPcj0CWpoZKdDRs2QHNzeFRQdj2KsjQTIqopjrcj36E+6iBMjiYW7DiF+La1Qz4OvV77Kioa8ksLN/x5sU93pDMnIrPX99Op4w+GNwC4x352r/oHn5UVdQchLpVNnT63SA/FpzgRfL7OXPgze9JTIJ2b52cmkhbr+/aOo8ZlcnaO7/l2M5PSuHP/I7kw4yCOjp7NpdkHcdsxs/sNZHyZGQx1ATNvXFt5q6tDcvk+XB14JRAJIZsSw/e5b1MXvRCjs5mDdpxCfNu6IR+HK2m1pWXILy324euL/ThTIt+vVim09W7bfqfxBaKVTtY69/u5IVf/XH/+OsXzhkVPbzBi5FgyPZ1vbzySVy89iIfPnu12+WEo32T1OoVz5k/w+fh1NRXMTEzjoslzvG68jTdZSDcksnOHwjhjEpceM45rz0vCktDl07X2nRnct7Lx0XlpAS9JDYXJk7Wcly7fftyg0uEKREI7jtG5NNODQx/N6pw3mV94Fskt33Dgzl/y3X7LaYmcPmRjiIuDggJteWbKlCG7rOihZ/2Osw+YwN8/3eG2+qUK3Hj0VH5YrbCrpY5m+94XxWN1azhGvw6bqucW29J+S1Hvy/VC4G16XoxcrpmL/ngqiBaMKqGTkr0XE3NptHZS2FzP3OSxKCg810+fmUOi8qipVpgxA+bO1WYJILD8jv62vH9745FhWZNn/HitG29Njfb/oRQuMyKjPhABcOij+CH7NQ7a+UsS21azYOcprJz8IW2W3CG5vk6n7aDZtUsCkVBw92IWH6lttW1s7ztV/Of3tzLfoJCS6oSfi0JG084dxhcBeNLxC3ao/hcLuHjhJD7cVBn0NxgxvAWa9+Evf5cCm37O95iTnM5SZS5v79rSKwckRm9hUVQeh2WnM28eTJzYuyOtv/kdw3XLe8+tvE5naLvyup4yod6RJ4HIz7SZkTc4eMeJxHX8xIIdJ7Fy8od0mCcNyfUTE7UmeI2NWqEzMTQ8vZg1tds81muo7+xkBes5PmFvoPoHwxukK/UUOVN5xH5qQGM5Oi+NP56QF5af4kR48WX2ZKBcgYG3wn3dY7JasNm0N9oZ8enk5KWxtbae6uZOInXac3n2LIXsbO2DV5/7+5G4G66Fy3w1adLerbzJyaEbh0J4zIiM6hyRfdkN8Xyf+y4tlslE2MpZsONkLNbyIbl2bKyWSV1WNiSXE/hWv6M/KytLMOv0zFYKuED/CQB/tC+lC//K5PbMAQkksVCIYOiZk+JNvMlCmj6RPXu0ZebiYmioV8iKTuLsg8fxf2clcdovFaZMcR+EgPb3GBdh4qKFk0iI6v03tG9+x3Df8h4dreWK1NWFdhw6RXJEwpLVmMyq3PdYuH0JUdZiDtp5Mt9NXoHVENxPH4qi7TPfuVObttu3gI8YfL6Uae9Pk60LA3aWmZ5Bp6i87VjEd07/coskB0SEsyXT03nivLnc9M5Gt8uUrmfs3WfmcUimQkuL1kNLp9NKxsfEeA48enK3PJoYZeTU2eNYnJfWZ2YwnAuX+So7G/LztSrL0dGhGYNs3w1jXaZ0Vu33XzqM44jp3MH8nWegd7QF/bqJiVBREfooebQYjBeppfoPmaoroV6N5q+28/y+f7hk8gvhyZLp6az709Fcs3g/4iN6RxU9n79RUVrzxowMra9KYqLvQcjlL63v86Ggoc3GcyuLaeqw9gnSw7lwma9SU7UlmlBu5ZWCZmGuwzyR73PfZeH2Y0loX8e8XReyJudVVMWHv6wARUdrSzNlZaFdNxwtBvoiNV6p5mrD2wD8zX4u9fheoSg+wshj587loCxZfhGh42u3Z71O4arFuVx5ZM6g5jAFmusR7oXLfKEo2uaE7du1WSRfgrZBH4OrxLsEIuGrNWIyP+S8wYIdJ5Ha/Amziq8kf9LjoARvIikqCnbsgBkzQptNPRp4ezHrn8oywzNEKFa+c+TxluNQv+7d2GFDpygShIghs2/Q0dBm5a7l/nV77i9J1tegpqdAm9T5k9gaziZM0GaRamq0Lb1Dbe+MyNBfuycJRLxoiJ7P2uwXOaDgHDLqX6PLOIat4+8K2vWSkqCyEqqqIF1m64PK24tZf3+bZ+i/4hD9JjpVI7fYl+K5HJln4bx+LUYWdzkY7gS69bW/eh79nWcguR5DWVMlWIxGmDYNPvpIe70f6txA1/Xssn03/FXHHctPE//B7N1XkFP1D7qMaexKvSIo14qI0Cru7dkjgchQ6O/F7JrD8vhhDXxUvYUWx97bxtDArYaXAHjAfgbFamC/qHBevxYjh6ct6u4EsvXVl3oenuqeDDTXY6hqqgRTZqaWT9PQoP13KLmWZiRHZJgoTT4Ps72aqWV3Mm3PLXQZx1CWeEZQrhUbqy3PzJ4dmnXD0cbdi9mc8Yl88rFCuh3ump/GrpZ6vq4oJr+ugr8anydWaSffmcVzjuP8vl7P9etAprOF8FV/ORieeFoO8ff8rqDmpnc2csd/t1DZ3He25Oi8tAHnegxFTZVgio2F3FxYu3boAxHX0owUNBtGClKvwWyrIqv6CWYXX06XIYnaWO/9RPyVlKQVN6us1DLQRfDt+2K2bp0WDE6cqN1W3NJAfl0lJ+hWc7R+LVZVzw223+JA79d1eq5ff7KlMqDpbCF8NZAt6r4sm/iS46Ft++299bfnbMlIyPUYqNxc+OknaG/XyjgMFSloNhwpCpvHL6Ms4ZfoVBv7F15AbPvGQb+MyQQOB5SUDPqphQ/Ky2HNGi0gNJthTc0e/rt7Owk0c6fxBQD+6Ti53zLuSxdO4prF+5EW677xFuB2y6LrBXrFpt7N9IQIxEDykHxZNgn0/K63PdcSUDg3qRsKaWnah56qqqG9rhQ0G64UHfmTnsBsqya59VvmF5zJt1M+pdM0blAvk5CgVSicNw8skkowKHxZBmlrV/nXinqK6juISrQSWWXitUIt2LzN+G+SlWa2OTN4zH6K22vsO6PhbrsjwKJ7Ph+25anF8BFIHpI/W18HkufUcwloJOR6DIROB1Onaq/5drtWJn8ouJJVZWlmGHLqzKzNfomF248lpnM7BxacwcrJH2LXxw3aNRISoKhI+3SelTVopx21fMnq/3BjBX98ewv1nT8f07z3/kfoNnCqfiUOVeEG22+w7fOnc+UROSzMSe7z4tlzyccVCK0sqAloy6IQ/vJ3i7q/yyED2wKvcc2qDPdcj4GaMAHGjNEKnA3VVl5dmCSrytJMgGyGBFbnvEWnIZXYjs3sX3gBitM6aOc3GLRotaho0E45anmq3NhzGWTFpgouf3n93iCkh2ja+ZvxWQCecRzPT2p2n2NyU6P77Q2zYlMFi+75nHOe/p5Hvyj0adyyvVcMVM9+Mb7MLfi7HOLv+d2R3WMas1lr79HcPHR1PZTuZNWhuZ4nMiMyAB3mCfyQ8zoH7ziBlJYvmVVyFfkT/zlom8GTkrRApKVF69kg/OdLVv8d/92MU/X8O7vV8FJ3Z92H7Ke7Paa/F1N/tk/6ek4hfOVpi3p6nIVbT5hKQpR5QMshHrfAx5rptDs9drIeDtVPh1p2tpYo39iozYoHm2smItTJqhKIDFBT1BzWZb3A/IKzyKh7hXbTRHaMvWlQzh0Xp60ZlpVppYCF/3zJ6q9s7vJ4+1G6dZxl+BKnqnCD7bd0Yu51u7cX00C2T8oLtBhswc7B8HT+T7ZUjvodMf6Ii4P99tO28g5FIKIosmtmxKiOO4afJjwIwOSKZYyvfXlQzqvTaXVECgpCX4J3uBrI8kYCzdxtfAaApx3Hs0btHQ368mLq7/ZJeYEWweLKwTh59rh+lxEH8/yu2ZLRvCPGX7m52jJNe3vwr6ULk0BEZkQGSUnKRURaS8itfJBZu39Ppyl9UGqMJCdrNUXq67WlGuGfgSxv3GV8gRSliR3OcTxo71u8Li7CyEULJ3F0XprHc/gbCA2n8tRC+GK074jxV3q61pW3sDD4GxX2VlYN7nW8kUBkEG0beysR1lLG17/J/oUXsHLyCloipw/onDEx2tLMnj0SiATCl6z+aJ0Zg0Gh0bo3aDhR9x2/0H+PXdVxre1yujABYNbpMZt0NHfaaOyw8dCnO3ltTanH4MHXQMjTrhshRoLRviPGH66uvDt2BL8rb/f2XVmaGUEUHT9OfIza6EUYnS0cWHAGFmvZgE8bHa21inY4BmGMo4wvWf1n5kzjtKy87n+PoYG7jM8D8KjjFDapez+WdDkdNHe6rxLprgiZKxDydG0FLWnwmqP3C8p0uRBi+Jk4UZsZCXaBM51UVh2ZXDVGWiyTibCVc2DBGRgcTQM6Z3KyVu69snKQBjnKeFqnjtZZuDB7LnNS0pkcl0y00QSo3G18mniljZ+cmTz6c+GyKMxEe/ho0rNK5L6FgfoLhCQfRAjhjtEI06dDays4g7i1NlySVWVpJghcNUYWbVvcXWNkde5bqEpgc2wWC1itUFwM4wa3gOuw5k/DuJ7r1Bu2dVK4ycLscYnEx2nHv1W0hVablUvN33Ckko8NEx+m3MevIvPoqLcQn6DydMFqj2PprwjZSGhXLoQn0rgxOLKytOX4ujpISQnONbqTVaWOyMi0b42Rmbuv4seJjwVcYyQxEXbuhLlzISJicMc6HPlSKXVfep1ChjmJn/bAlESIj9M+CSwv2c4P1XvIVfZwg+4FUGHH+NvJSz2a+npojwZLZhkUeB+Xp+RUSdgTI1Egf4fCN1FRMG0afPNN8AIR16tPqGdEZGkmiFw1RlR0TKh7mf0q7gn4XImJWmS8Z88gDjAMOJwqqwrreC+/jFWFdT71PPClUqo7DQ3w5ZfQ1aU1mcqvq+BPaz7l4z2FmLHysPExjGonhZaF7BrzO6xWrdzyvHkwZZJvSaf9JacGe/ukEEMp0L9D4bucHG3DQmNjcM4v23dHCVeNkVklVzO5YhkdpvGUJp/n93n0eq3s+86d2pNzkIq3hlQgn6Z8qZTqrmFcRwd89ZWWZ5ObqwUhz25b3337DYbXydPtplaN5ezGczmpvoqo+nQmT4bZs0Fv6H/3jRQhE6NJoH+Hwj9JSVqBs/XrIT5+8M+vjPRk1b/+9a8cfPDBREZGEh+MR3AYKUm5iJ1p1wIwc/fvGdP0SUDnSU6G3bu1miLDXaCfpnyplFrR1Mn3hXXd37PbtenNnTt/3pevqLy9a0v37Yfr8llq+BCAP9h+Sw3xvFmwhegYlYMOApNJkk6F6MnXv8MfikbAi1WITZ6sFThraxv8c4fLjEjQAhGr1coZZ5zB5ZdfHqxLDCvbxt5KaeLZ6HAwb9cFxLet8/scsbFaFnVx8eCPbyh5+zQF2qcpq93ZZ9nG1wJhV7yiBTOqCmvWwE8/ad0t9QaVr8qLumuGJNPEfcYnAHjefixfOucA0GzvJDqznjFj9p5TqkQKofH171AaNw5cejpkZgZn1+TeXTODf25/BG1p5s477wTghRdeCNYlhhdFx4+THsVsr2FM82fMLziDlZM/ps2S49dpYmNh2zaYOTO4hW6CyddPUwct+4z6tr0djdPjLJx9QIZP12jssHHZS+s5afIkjNVpzBqbyI6OSt7evKU7CFFw8oDxcVKUZrY6J3C3/Zxe54hI7DtGSToVwvdCfdK4ceAURUta3blTy28zm73fx1eumQhVckT26urqoqtrbwOy5ubmEI5m8KmKkbVZ/+LgHb8gvn0DB+08lW+nfEKX0XOJ8H2lpGgzIqWlwS//Gyy+fkrqGYSAtmzz0Kc7iY80euzoua//bi8GivmwxUi7vXchsiv1/+Ew/U90qCZ+b7uyu3qqy74zHy5SJVKMdt4qFkvO1ODKyNBmdMvKtPLvg8U1IyKVVXtYtmwZcXFx3V8ZGb59+h1OHPpoVue8Sas5i0hrCQfuPN2vgmdGo9YMb+fOIA4yyIprA1vsdCXBBWLfIGShbiPXGN4G4E+2i9mpju91e7q8iArhkeRMDS29HmbM0Eq+22zej/eVLkx6zfgViNx0000oitLv17Zt2wIezM0330xTU1P3V2lpacDnCmdWYwqrc9+hy5BCXMdG9i88D53Tcyv6faWkQFGRtp13uHE4VV5ZXRLw/VWgsd3G1Yv3Iz4isLWpVOp52PgYOkXlVfsRvO08tM8xZx8wQV5EheiH5EwNrUmTYPx4qBjEXdHdvWZCnCTi19LMddddx69//et+j8kawHqB2WzGPJgLYGGs3ZzJ6ty3OHj7CaS0fM3s4stYn/ksKN5jw9hYKC/XgpFAGuGFshLio5/vpKrF96DLk6YOK4+dO5dzn/Fc7dQdA3YeMT1CstLMZudE7rBf6Pa4ScmRAx6jECOd5EwNHaNRyw1cvnzwmuF1d9/1aaE7ePwKRFJSUkgJVom3UagpcjZrsv/NgQVnMq7hHbqMY9g8/m6fioTEx8PmzVo/Aosf+WChrIS4YlMFD306OGtKz60sZv+JCV476+7rBsPrzNdtp1mN4He2q/rkhbhIkp0QvpGcqaGTlQVjx2rN8MaP9368NyN++25JSQn5+fmUlJTgcDjIz88nPz+f1tbWYF1yWKqNPZINk7Tto1nVT7BfxTKf7peSAjU1/m3lDWUlRNeW3cGiAHct38qtJ+R5PdblZN23/MawHIDrbb9lt9o3SdjVDVfyQ4QQ4cZk0mZF2tq0+kgDtbeg2cDPNRBBC0Ruu+025syZw+23305raytz5sxhzpw5rF27NliXHLbKE09nY8b9AEyuuIfsyoe93kev13rObNrk/gm5b+l0q93pU+2OYK0Vetuy21N8pPc5R9cW353VrcRYvB8/UynkHuPTADxmP4mPnPM9HitJdkKIcJWTo82KDEZdEdfLXKhnRIK2ffeFF16QGiJ+KB5zKQZnC1PL7iSv7Dbs+ih2p1zS731SU7XeM6WlWsEbF3fLL4lRRurbPKdb99c91mUguSX+FDa6+5czWFNUz7Mri70e+9CnO7wek0IDT5kexKLY+MQxl/vtZ3o89urF+0mSnRAibJnNMGsWfPjhwHNFXEszUkdEdCtIuxaDo5XcygeYWXIdDl00e5LO9ni8K69382aYOFHb1utaftn3adVfENKTp4BhoLklvuZcXPNzIBAXYfIpEPHGjJWnTA+RpjSwwzmOa2y/Q+1nItBTkqq0OhdChIucHBg3TpsVGUiVi3DpNSOBSJjZNvZW9M42sqqfYHbx5dh1kVQmnOTx+PR0KCzUZkbGjfdcOt1XyVFmVhXW9XrD/WRLpdvgxpVb4stWPW8FkADSYs1ceWSOz8d7p/I347PM0RXQqEZxqe06Wul/N4y7gElanQshwonZDHPmaDtorFYtdyQQ3TMigzi2QEggEm4Uhc3jl2FwtDKh7iXmFV3MGuUlquOXuD08MlIrRvPTT7DH6nseRp/LAnGRRq5780cqm/eeIy3WTKfdOeAum64CSJe/tB6F3k98173uOGla9zn6O95XVxve5jT9N9hVHVfYfu82ObXnGNxVgvQ0w+RPECaEEIMtO1ubCS8vD7zaquu1Vyqrir4UHT9O/AdlCaehU20csOs8UhuXezx87FgoKIDNuwIPQlyFwnoGIQCVzV00tvuWW+KNvwWQPB0frbOwODW332udqf+Cqw3vAHCr/SJWOmd4PNZTJUhfm/OFuhiQEGL0MRph9mxts0JngL0FFckRCZ2WThvba9qos+uJ79Bh0Okw6vQYddr/630oKhZ0ip4NmU8BMK7hbfYvvIB1Wc+7XaYxW1TKbPU0bGkJ6FKpP8969BdweONrMqq/BZCWTE9n/rg03vqqnp92dmJ2WliQk4jRCGsbSrsb2PV0uC6fvxmeBeDr2N/yv/qjgb19a3RK7+1qaR6WWfxpdS51FIQQQy0zE3JztZYfOf71TwV6Ls1IIDLkfixt4rJ3fq7I6eaDvEVvINJgJMpgIsqo/TfWZCbBHEGC2UKiOZIEs4Vog6k7ogwGVTGwIfMpVEXP+Po3mLfr16zPepaKhFO7j8mvq+DtXVvcviF7c+UROSzMScapqn5XKN3XvrkV/SV3+loAqaEBduyAzZsVGmuTODBdK+TmclpWHs9uW9/rPjOUXTxmfBiD4uTHqJNpyrmHvwLfbq9HtXRy9CILR8xIYN3uBq+BkLQ6F0KEM71emxXZvRuam7Wq2/4IlzoiozIQ0elgTLSFljYnqs6J3enErjq7b+902Ol02Knv6uj3PEadjhRLFGmRMaRHRpMaEU1aZDQpligMusGZVVEVAxsmPYGKnoz6V5m7ayn5mXbKEs8gv66izxuxP3JTo1mQncR7+WUBn8NdbkWgyZ1OJzQ2aoXadu/WStg3N2tl7KdM6VtwNi9+DBlRsZS2aV2aJysl/Mt0N1FKF7ssB1O63zOgKFRXQk5sEscdp3WwBHwKhKTVuRAi3I0bB9OmwZo1EBPjU2Hubq5jZWkmBA7OTub1c4/i3Xe1aS3Qti/ZnU6sTgftdittNhttdittdhttNivN1i7quzpo6OqgvquDZlsXNqeT8vYWytt7L4noFYWxkbFMjIljYnQ8E6LjSYuM7p4G85uiJ3/SY6iKngl1LzGn6FIM9gZu3Z09oMehuLYdCPyN1F1uhb/Jna2tWhOnigpt509TE7S3a5F+UpK2K8jdw9Zus/LI5tXsaWvGoOi4cKyRaxruI9rZSkPkPLbv9zqqzkRDg3a+Y47ZG4T4SlqdCyHCnaJodUV27YLqaq2+lK9c3Xfdb0cYOqMyEHFHpyiY9HpMej3RRhNE9H+8zemgoauT6o5WKttbqehopbK9har2VrqcDkrbmihta+JbtE6zJp2eiTHx5MYmkRuXxKSYeP9mTRQ9P058BIfOQmbNM8wsvZ6LnL/kIU6jbyNu37y2poQrj8zpfsP1d8fNvrkV3pI7e+6wqa9T2LpVW3ppbAQUlUZjPTZDJ2PGWMiMS/QYuH1buZs3d23u3vueRiW/qf0z0Uo9TREzWJ37NnZ9LK2tUFsLixZpMyr+8mWnj1RhFUKEWkICzJsHn34KiYm+FzmTgmbDnFGnZ0xEFGMiopieuDcEVVWVuq4OSlobKWlpYndrIyWtTVidDnY21bGzqQ5KtWWdrJhEcuOSmJqQQkZUrPd8E0XHpoz7sRpSmFyxjKsM75BIM7fbf40zgA1QPRMtT5qVzpNfF/l0vyuPyGFBVhIoUNvaxarCOuZnJvqc3Hn7S0U015jRdVmYmZ5Ie2Il7xT3znOJN1k4LSuP2Um9l3I+2VPAf3dv7/73eKWGV41/JVWpZ6dzHP9JegJ9m5PatjLa6yycfHAi8+YpPk1Xustrce3c2XepyVOCqxBChMLUqVpNqZISbWuvL1zvOU4vxwWbBCKDTFEUki2RJFsimZs8FtCWfao6WilsqmdHUx07m+totVnZ3lTL9qZa3i/ZTpzJzPSEVKYnjmG/uGRMer2nC7Bj7E3ssVs4ovoOzjd8SpLSzLW2y+nE7Pd4q1s6WbGpgqd8DEIAtlQ08fb6PX1yQI6f7rlOR08vbdna/f8ftxtpt/fdrdNo7eTZbetZOmUus5PScaoqn+wp5P2SvUFIllLOy6a/ka7UU+xM5VzrLdQU7EZld/cx339n4Y5k7wGDt7wWaXUuhAhnJhMceKDWmbehQZsl8cb18TXUMyKKGuoR9KO5uZm4uDiampqI9Tcd2IuiInrliAwlVVWp7GhlR1MdOxpr2dZYi9Xp6L7dqNMxOS6ZucnpzEhMw2LoGy86VZWVa+/iL8rfMSkOfnRmcYn1Omrw4dnXw8tLD+QPb/0YcCE0l0CLjnkTb7Jw8eS5vFu8haKWxu7vT1V28y/TMlKUZgqcYznXegtV9M3VcIUK/RUe85TX4st9hRAinHz/PXzzjTYr4m2J5uuKYt7ctZlZiWm8d8O8QR2HP+/fMiMSAoqikB4ZQ3pkDIelT8L287LNpvpqNtVX0WDtZFNDNZsaqjHqNjI9IZV5KWPJS0jBqNNmSnSKQkzWpZy3PZInTA8xS7eL/5pv5RLrH9isTgIgymCkzc1sA+xNtERhwEEI7A1CFEWr9DpYGq2dPLjxOwAMig676mR/ZRvPmu4nTmlns3Mi51tvph73T/SehcfcVX/1J69FZkCEEOFuzhyt2mpxsfcP2q5kVakjIjDq9OQljCEvYQxnZE2jrL2FH+sqWVdTTk1nGxvqKthQV4FFb2BOUjoL0jKYFB2v5U9MPpcLdqXzd/Uv5OjKedN0J3eoVxKTdSmAx+29KnDNYXlUN3cN6s8SrPm1+WPGMz1hDLU7H+U+4xOYFTvrnLlcZL2BZqK83t9T4TEpWiaEGEnMZli4EOrrtd2I6f1M5nbniEgdEdGToiiMj4plfFQsx2fkUtrWzLqaMtbXVtBo7WRVdSmrqktJj4xhYWoGB4wZz8z9f8VHDQdDxRXkdK7kXuUBdrW1smX8X1g6ZW6fgmcxeguLIvOozU+nmroQ/rS+OTt7BgvGjCdx910cYnoUgBWOA7ja9ju/8mI+2VLZJ5iQomVCiJEmNVULRj7+WCuJEBfn/jjXJK/MiAiPFEVhQnQcE6LjOHnSVHY117OqqpQNdRVUtLfwVtEW3tu9jdlJ6RySPpFtef9FLf8LuZUPkVXzJAlta7FkvcDM/Y+ksLmeJmsncSYL2bGJKCh0dICpNZEYvYUWR3i+0cabLMTprJjWn8ghyrcAPGU/gbvt5/i9U+i9/HL+eELv7bZStEwIMRJNmaKVRli5UssViXTTeNxVWTXUqaISiAwTOkUhJy6JnLgkTrNPY01NGd9VllDe3sKamjLW1JQxMTqew8deypKs+czbfRkJ7es4dOshbJzwILrE0/qcMzISIiMVztT3LZUeLn6R0MXpu08iV1eGTdVzp/0CXnIcHdC56tqsfZZYpGiZEGIkUhQ44ACtoOP69VpBx32Dkb3bd6X7rvBTpMHIYemTuGn2IVw382AOHDMeg6Jjd2sjL+7I57ICE3+N/zd1EbMxORqZV3Qxc3ddjNHuvkPu7KR0lk6ZS7wpdJ/6owy907vjTWYeGrebWxovIldXRqWawNnWPwUchLjsu8TiKloGfcvCSdEyIcRwZjBoBR1nz9bqi7Ts0xc1JzaRE5PncPS4ADrmDSLZvhuC7bvB0NjVyfsl2/mxrpJOhx0As+LgbzGfcar13+hw0GlMY9P4ZVrTPDcVvpyqSmFzPT/VVfJlRfGQjNtVuGxmYlr38lG6roWTG+4gvelDAFY6pvF725XU4WGh0w+vXnqQ26TTQPvjCCFEuLNatW2969dDRASkpWk910BrrTFuHJzUt7H7gMj23VHGUwfeLlXPdc3H8KKSzaOWp5hgK2X/oouoqX2RjRPup83SOwrTKQrZsYn8a0d+UMaZHhnDqZOmYtDpeuWruMoM58YmMq7+DaaX3ojJ0YAdAw/afsmTjhNx4KHAm4+8LbFI0TIhxEhlMmkzIykp8MMPsH271qk3Ohps7is8DCkJRIY5Xzrw/qRmc3THXVym/x+/M/6XlJYvOWzzAopTLqYg/Q9YjWO6jy1sru8T0AyGhWkTODNrusf+MXFtG5heeiOJbasBaIqYyXvJd/PPnY1+XyvQvjB6nSJbdIUQI5JOp5WBHz8edu6EggKtu3ldHUyaFNqxSSAyjDlVlbd3bfHpWIdi4WHHabzrXMQdhhc5Up9Pds2TjKt5ka+ifsVy01lU2C0U96heOhgi9AbOzJ7O/inj3N4e3bmD3Ip7GVf/Fgoqdl0kO9P+QGHa74lVDcQWfU6z3Xtg5JrxuPWEPO5aLn1hhBDCnZgYmDtXyxvZtAkqK0M9IglEhjV/Zi/sqtbWqERN5WLb9Sx0bOJ6w+vM1u3i2LbnOKz137zjWMTzjuOoYbzfY9ErCjFGM4nmCMZGxpAcEcmE6PheSy/dVJWEtjVkVj/B2IZ3UH6ev9iTcDqbx93JusYIKndUo7NaWJyUxztVfbvf9tRzxmPJ9HSOnS5LLEII0R+dTttF42un3mCSQGQYawp4CUVhpXMGK63TOUa3lqsM7zJNV8yvDF/wK8MXrFen8I59AcsdB9LgoXR6TxfkzuKAMZ6DF1cSrKOjlLldXzCv7W3iOzZ2314W8wt+iLiBT5vG8NX6LbQ69/5cqaqF3x6ayetr99DY7n4xc98ZD1liEUII7ywWLX/E4fB+bDBJIDKMxQ14u63Cx84D+Ni6P/dk6zm07WXSGpczV9nGXOM27jC8yAY1h28dM/jOOY0t6kTaiOhzlnhz7++5Ao+2zkaiW9ZgbPiSRaxnf2UHOkWb17ApZnZFnMH36m9pMc+kTK1geVPfXJfqlk6e9NIZ+NYTpsqyixBC+Mls1mZErNbQjkMCkWEsOzaReJNlEJJLFbboZxKZfRwWazljG94hsfpl0q1bOEDZwQG6HVzD2wAUOVMpVtOoUBOpJgF0kSzu2IXSCUZ7E01tZXQ0b+dUdQ8ZSg0GxdmrWs16Zw7LHQfxluNQjkk4hDMXjGXsOJXTXnCf6+Jtb7kC3LV8K8dOT5flFyGE8IPFIkszYoB0isJpWYNTFdU1u9JpGsuu1CvZlXol5s7dVO9+g4Smj5mjK2CsUk+mropMqnrfuWyfkyl0J25UqfH84JzCaudUPnPMpYK9SyZvV2zgCLOCvcNEZXNgwZQ0pRNCiMC4ZkRCTQKRYc5VFdVdHRGX/hI9QSsqlh3bt75Gl2UicZOvZ33defxx1xZ0XdVM1e1mvFLLJH0z82MhxWxAp9pQVVhZ30K1w0y1mkChOpZdznSqSKBvzVKNU4XfvbKeo6ak+P+D70Oa0gkhhH9MJglExCCZnZTOzMQ0Piot4IPSHX1u97a8cVpWnsf6Hj3P37Nx3tjYRMoVhfKfj9nZVMc/Kr4PaPyfb6sJ6H49SVM6IYTwn8kU6hFIIDJiOJ3wdVlJv8fsOzPiKq8+O8l9oqfDoVXdczhAVRXS9UmMiwK9Hhx2UHvkftR3BD4jMZAeA9KUTgghAmc0hr66qgQiI4DDAV9vre+17dUdFTgvbyrtDWbosDDWmEhMl0JtrXa7zQadnXuflHq99iTV6/e2pnE6tS+HQ/uvqmq3ddmCPyMRaMVUIYQQ7hkMYLeHeAyhvbwYKFWFXbvAHN8JTd6PP2CmmWMmj6OsDCoqtKp6bW3abbGxkJkJSUlaD4KoqL3JTDqddi2HQ/uy210zJT8HKUoiXz9vobqlMygNpa9ZvB+vrSmRiqlCCDGIFCX0eSISiAxze/ZAYiJMnWHhtd3ej0+OMpNfXkd1aydjxlr45cJEUBUUZW83xsAo3HlyHpe/1H8VVHd0ihbQeLpPepyFK4/M4cojc6RiqhBCDKKpU6Fs352PQ0wCkWGsqUmbmVi0CDKzEkn/2EJlk/sZCQWIizRy3Zs/9toqO5it7pdMT+fx8+Zy5/9693qJjzDQ2NF37s8VQlx6SCZPfV3k09KLbNEVQojBk5OjfYWSBCLDlN0O5eWwcCFkZ4OiKNx+ovsZCde/tRLpvbOSKps6ufyl9Tx+3txBC0aOzuvb6+WTLZV9ApSeSytzJiT0e7sQQoiRSVFVNRhL+oOiubmZuLg4mpqaiI313vPEH0VF8O67kJs7qKcdMkVFkJYGJ54IET0qrK/YVNH3DT3WTKfd6bFXi2vnybc3HhnwUofDqbpdNun5/eQoMyhQ29rldmnF0zmEEEIML/68f8uMyDDU2qr994ADegch4H5GwqmqnPvMao/nc1UnfWFlEckxZr+DAHfBT3qchZNmpfPfHyv6fP/2E/PcLrFIszohhBh9JBAZZlRVSyyaMwcmTXJ/zL5v6O/l+5aJdNfyrd3/72vuyIpNFVz+0vo+eSkVTe6b1fm6FCSzI0IIMTpIIDLM1NZCfLwWiPRTDLWXQKqO+hIwOJwqd/5vi187ZFS0paA7/7eFo/PS3AYXnmZYJF9ECCFGngFt2BRDy+HQApHZs7Utu76an5lIepzFQ8cX91zBxZ3/24LD6T7U+KGovlew4M+5XY3q9uWaYdn3vK7AaMWmCr+vJ4QQInwFLRApLi5m6dKlZGZmEhERQXZ2NrfffjtWqzVYlxzxqqogPV3b9+0PvU7bUQOe2s+511/AAANvNLfv/fubYfElMBJCCDH8BC0Q2bZtG06nkyeffJLNmzfz0EMP8cQTT3DLLbcE65Ijmt0OLS3abEhUlP/3d9X4SIvzf5nGU8Ax0EZz+97f2wyLt8BICCHE8BO0HJElS5awZMmS7n9nZWWxfft2Hn/8ce6///5gXXbEqqyEceNgv/0CP4drR82jnxfw0Kd9u/R64ingcC35eCqi5omnRnW+zrAMdCZGCCFE+BjSHJGmpiYS+0lu6Orqorm5udeX0JrQtbdrsyFm88DP99qa/rv0uihoSaL7BgwOp8qqwjre/6mcsw+Y0H2sr+cE943qfJ1hGehMjBBCiPAxZLtmCgoKeOSRR/qdDVm2bBl33nnnUA1p2Kiq0mZDsrMHfi5/E0z3DRjc7WiJj9Q6JvUsmOapjkh/1VK9zbB4mkkRQggxfPkdiNx0003cc889/R6zdetWpkyZ0v3vsrIylixZwhlnnMGll17q8X4333wz1157bfe/m5ubycjI8HeII4rdrnXHPfxwMJkGfj5flzUUBR47p/fWXU81Q5rabajANYtzmZQc1avuxw1LpvpcD8SVVOupTD24n0kRQggxfPkdiFx33XX8+te/7veYrKys7v8vLy/niCOO4OCDD+app57q935msxnzYKw9jCCu3JAeD+mA+LqsoaqQELU38vG2o0UBXltT2qdMvL/VUj01zpO+M0IIMTL5HYikpKSQkpLi07FlZWUcccQRzJs3j+effx7dwPrMjzoOhzYbcuihg5MbAtryR3yEkcYO931neuo5e+LPjpaBlmn31DhPZkKEEGLkCVqOSFlZGYcffjgTJ07k/vvvp6ampvu2tLS0YF12RKmu1hrbDUZuiItep3DRwkk89OlOr8f2nD0Z6h0t0ndGCCFGh6AFIp988gkFBQUUFBQwfvz4XreFccPfsOF0QlMTHHggWAZ5k8iVR+by/HfFXrvx9kwKlR0tQgghgiFoayW//vWvUVXV7Zfwrq4OkpIgJ2fwz63XKdz9yxlut9x6Sgr1Vibe01ZfIYQQoj+StBGGVFXrKTNtWmBVVH3hSgpN36fSalqcxW2ju/7KxMuOFiGEEIGS7rthqLFR67A7kCqqvvA3KVR2tAghhBhsEoiEoZoa2H9/SEgI/rUC2V4rO1qEEEIMFglEwkxrq5acOnlyqEfimexoEUIIMVgkRyTMVFZqxctSU0M9EiGEECL4JBAJI11dWmn1vDztv0IIIcRIJ4FIGKmqgowM2KfsihBCCDFiSSASJux26OjQtuzq9aEejRBCCDE0JBAJE7W1Wjn3SZNCPRIhhBBi6EggEgZUVasdkpc3+OXchRBCiHAmgUgYaGjQCpgNZnM7IYQQYjiQQCQM1NRoVVTj4kI9EiGEEGJoSSASYm1t2nJMsMu5CyGEEOFIApEQq6rSElTT0kI9EiGEEGLoSSASQjYbOBwwdaoUMBNCCDE6SSASQtXVkJ4OEyaEeiRCCCFEaEggEiJOJ7S0aAXMjMZQj0YIIYQIDQlEQqS+HpKStAZ3QgghxGglgUiI1NXB5MkQHR3qkQghhBChI4FICLS2QmQk5OSEeiRCCCFEaEkgEgJVVZCZCWPGhHokQgghRGhJIDLErFatt8yUKbJlVwghhJBAZIhVV8O4cZCREeqRCCGEEKEngcgQcji0ku55eWAwhHo0QgghROhJIDKEXFt2MzNDPRIhhBAiPEggMoTq67Vy7lFRoR6JEEIIER4kEBkizc1aAJKdHeqRCCGEEOFDApEhUl2tBSEpKaEeiRBCCBE+JBAZAp2d2lbdKVNCPRIhhBAivEggMgSqqrQOu+PGhXokQgghRHiRQCTI7HZtRiQvD/T6UI9GCCGECC8SiARZbS2kpcGkSaEeiRBCCBF+JBAJIlWFhgZtNsRiCfVohBBCiPAjgUgQ1ddDYqJs2RVCCCE8kUAkiGprYb/9IC4u1CMRQgghwpMEIkHS0gKRkTB5cqhHIoQQQoQvCUSCxFXAbMyYUI9ECCGECF8SiARBZ6f23ylTtEJmQgghhHBPApEgqKyEiRMhIyPUIxFCCCHCmwQig8xmA6sVpk0DnTy6QgghRL/krXKQVVXB2LGQmRnqkQghhBDhL6iByEknncSECROwWCykp6dz/vnnU15eHsxLhpTDAW1tMHMmGI2hHo0QQggR/oIaiBxxxBG88cYbbN++nbfffpvCwkJOP/30YF4ypGpqtF0yWVmhHokQQggxPBiCefJrrrmm+/8nTpzITTfdxCmnnILNZsM4wqYMnE5obITFiyEiItSjEUIIIYaHoAYiPdXX1/Pyyy9z8MEHewxCurq66Orq6v53c3PzUA1vwOrqICkJcnJCPRIhhBBi+Ah6suqNN95IVFQUSUlJlJSU8N5773k8dtmyZcTFxXV/ZQyT/a+qqgUi06dDTEyoRyOEEEIMH34HIjfddBOKovT7tW3btu7jr7/+ejZs2MDHH3+MXq/nggsuQFVVt+e++eabaWpq6v4qLS0N/CcbQvX1kJAg5dyFEEIIfymqp6jAg5qaGurq6vo9JisrC5PJ1Of7e/bsISMjg++++44FCxZ4vVZzczNxcXE0NTURGxvrzzC9KiqCd9+F3NyBnUdVYft2OOQQOOigwRmbEEIIMZz58/7td45ISkoKKSkpAQ3M6XQC9MoDGe5csyFTpoR6JEIIIcTwE7Rk1dWrV7NmzRoWLVpEQkIChYWF3HrrrWRnZ/s0GzIcqKq2ZXfRIoiPD/VohBBCiOEnaMmqkZGRvPPOOxx11FFMnjyZpUuXMnPmTL766ivMZnOwLjukXLMhU6eGeiRCCCHE8BS0GZEZM2bw+eefB+v0Ied0arMhhxwisyFCCCFEoKTXTIBqa7W6IXl5oR6JEEIIMXxJIBIAh0Nblpk1CwZ5M48QQggxqkggEoCqKkhPl9wQIYQQYqAkEPGTzQYtLTBnDkRGhno0QgghxPAmgYif9uyBSZNgv/1CPRIhhBBi+JNAxA/t7VrtkLlzYYQ1DxZCCCFCQgIRP5SWahVUJ00K9UiEEEKIkUECER/V1mo7ZObOBZ08akIIIcSgkLdUH9hsWiAybx4E2GZHCCGEEG5IIOKD0lLIzITp00M9EiGEEGJkkUDEi4YGLTH1wANhhLTIEUIIIcKGBCL9sNm04mVz50JGRqhHI4QQQow8Eoh4oKpQVAS5uVrxMiGEEEIMPglEPCgrg4QEWLRIlmSEEEKIYJFAxI26Oq2x3aGHQnJyqEcjhBBCjFwSiOyjqUlLUD34YMjJCfVohBBCiJFNApEeGhu15NSDDoLZs0M9GiGEEGLkM4R6AOGishLa2mDhQpg/X6qnCiGEEENh1AcibW1QXg4xMXD00ZCXB4oS6lEJIYQQo8OoDkR0Oq10+9SpWvn2MWNCPSIhhBBidBm1gUhKChx1lLYrJi1NZkGEEEKIUBi1gUh0NMyYEepRCCGEEKObpGQKIYQQImQkEBFCCCFEyEggIoQQQoiQkUBECCGEECEjgYgQQgghQkYCESGEEEKEjAQiQgghhAgZCUSEEEIIETISiAghhBAiZCQQEUIIIUTISCAihBBCiJCRQEQIIYQQISOBiBBCCCFCJqy776qqCkBzc3OIRyKEEEIIX7net13v4/0J60CkpaUFgIyMjBCPRAghhBD+amlpIS4urt9jFNWXcCVEnE4n5eXlxMTEoChKqIcT1pqbm8nIyKC0tJTY2NhQD2dEkcc2OORxDR55bINHHlvfqKpKS0sLY8eORafrPwskrGdEdDod48ePD/UwhpXY2Fj54wgSeWyDQx7X4JHHNnjksfXO20yIiySrCiGEECJkJBARQgghRMhIIDJCmM1mbr/9dsxmc6iHMuLIYxsc8rgGjzy2wSOP7eAL62RVIYQQQoxsMiMihBBCiJCRQEQIIYQQISOBiBBCCCFCRgIRIYQQQoSMBCLD1F//+lcOPvhgIiMjiY+P9+k+qqpy2223kZ6eTkREBIsXL2bnzp3BHegwVF9fz7nnnktsbCzx8fEsXbqU1tbWfu9z+OGHoyhKr6/LLrtsiEYcvh577DEmTZqExWLhwAMP5Icffuj3+DfffJMpU6ZgsViYMWMGH3zwwRCNdPjx57F94YUX+jw/LRbLEI52ePj666858cQTGTt2LIqi8J///Mfrfb788kvmzp2L2WwmJyeHF154IejjHGkkEBmmrFYrZ5xxBpdffrnP97n33nv5xz/+wRNPPMHq1auJiori2GOPpbOzM4gjHX7OPfdcNm/ezCeffML777/P119/zW9+8xuv97v00kupqKjo/rr33nuHYLTh6/XXX+faa6/l9ttvZ/369cyaNYtjjz2W6upqt8d/9913nHPOOSxdupQNGzZwyimncMopp7Bp06YhHnn48/exBa0SaM/n5+7du4dwxMNDW1sbs2bN4rHHHvPp+KKiIk444QSOOOII8vPzufrqq7nkkkv46KOPgjzSEUYVw9rzzz+vxsXFeT3O6XSqaWlp6n333df9vcbGRtVsNquvvvpqEEc4vGzZskUF1DVr1nR/78MPP1QVRVHLyso83u+www5Tr7rqqiEY4fAxf/589Yorruj+t8PhUMeOHasuW7bM7fFnnnmmesIJJ/T63oEHHqj+9re/Deo4hyN/H1tfXyfEXoD67rvv9nvMDTfcoE6bNq3X98466yz12GOPDeLIRh6ZERklioqKqKysZPHixd3fi4uL48ADD2TVqlUhHFl4WbVqFfHx8ey///7d31u8eDE6nY7Vq1f3e9+XX36Z5ORkpk+fzs0330x7e3uwhxu2rFYr69at6/V80+l0LF682OPzbdWqVb2OBzj22GPl+bmPQB5bgNbWViZOnEhGRgYnn3wymzdvHorhjmjynB0cYd30TgyeyspKAFJTU3t9PzU1tfs2oT1OY8aM6fU9g8FAYmJiv4/Tr371KyZOnMjYsWP56aefuPHGG9m+fTvvvPNOsIcclmpra3E4HG6fb9u2bXN7n8rKSnl++iCQx3by5Mk899xzzJw5k6amJu6//34OPvhgNm/eLI1FB8DTc7a5uZmOjg4iIiJCNLLhRWZEwshNN93UJ6Fs3y9PLzSif8F+bH/zm99w7LHHMmPGDM4991z+9a9/8e6771JYWDiIP4UQgVmwYAEXXHABs2fP5rDDDuOdd94hJSWFJ598MtRDE0JmRMLJddddx69//et+j8nKygro3GlpaQBUVVWRnp7e/f2qqipmz54d0DmHE18f27S0tD4Jf3a7nfr6+u7H0BcHHnggAAUFBWRnZ/s93uEuOTkZvV5PVVVVr+9XVVV5fBzT0tL8On60CuSx3ZfRaGTOnDkUFBQEY4ijhqfnbGxsrMyG+EECkTCSkpJCSkpKUM6dmZlJWloan332WXfg0dzczOrVq/3aeTNc+frYLliwgMbGRtatW8e8efMA+Pzzz3E6nd3BhS/y8/MBegV9o4nJZGLevHl89tlnnHLKKQA4nU4+++wzrrzySrf3WbBgAZ999hlXX3119/c++eQTFixYMAQjHj4CeWz35XA42LhxI8cff3wQRzryLViwoM8Wc3nOBiDU2bIiMLt371Y3bNig3nnnnWp0dLS6YcMGdcOGDWpLS0v3MZMnT1bfeeed7n/ffffdanx8vPree++pP/30k3ryySermZmZakdHRyh+hLC1ZMkSdc6cOerq1avVb7/9Vs3NzVXPOeec7tv37NmjTp48WV29erWqqqpaUFCg/vnPf1bXrl2rFhUVqe+9956alZWlHnrooaH6EcLCa6+9pprNZvWFF15Qt2zZov7mN79R4+Pj1crKSlVVVfX8889Xb7rppu7jV65cqRoMBvX+++9Xt27dqt5+++2q0WhUN27cGKofIWz5+9jeeeed6kcffaQWFhaq69atU88++2zVYrGomzdvDtWPEJZaWlq6X0sB9cEHH1Q3bNig7t69W1VVVb3pppvU888/v/v4Xbt2qZGRker111+vbt26VX3sscdUvV6vrlixIlQ/wrAkgcgwdeGFF6pAn68vvvii+xhAff7557v/7XQ61VtvvVVNTU1VzWazetRRR6nbt28f+sGHubq6OvWcc85Ro6Oj1djYWPWiiy7qFeAVFRX1eqxLSkrUQw89VE1MTFTNZrOak5OjXn/99WpTU1OIfoLw8cgjj6gTJkxQTSaTOn/+fPX777/vvu2www5TL7zwwl7Hv/HGG+p+++2nmkwmddq0aery5cuHeMTDhz+P7dVXX919bGpqqnr88cer69evD8Gow9sXX3zh9nXV9VheeOGF6mGHHdbnPrNnz1ZNJpOalZXV6zVX+EZRVVUNyVSMEEIIIUY92TUjhBBCiJCRQEQIIYQQISOBiBBCCCFCRgIRIYQQQoSMBCJCCCGECBkJRIQQQggRMhKICCGEECJkJBARQgghRMhIICKEEEKIkJFARAghhBAhI4GIEEIIIUJGAhEhhBBChMz/AxImruKEBjH0AAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.plot(fitted[\"x\"], fitted[\"mode\"])\n", "plt.fill_between(\n", " fitted[\"x\"], fitted[\"0.025quant\"], fitted[\"0.975quant\"], alpha=0.3, color=\"b\"\n", ")\n", "plt.scatter(x, y)\n", "plt.plot(x_true, y_true)" ] }, { "cell_type": "code", "execution_count": null, "id": "42d674ea-467e-470b-a091-941ddb5b7cb8", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "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.4" } }, "nbformat": 4, "nbformat_minor": 5 }