{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Scripting landmark warps\n", "\n", "Landmark warping can be scripted using the **mwarp1d.warp_landmark** function.\n", "\n", "Let’s first run an example script then consider the details:" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZMAAAEJCAYAAABR4cpEAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOydd3xV5fnAv08GhJWEJIQNicywguwdhiBCFUfdOKu01aptbfvTOutoa+tsXXXbFrV1oFZxRENAZEaUvXcA2QRCQsh4fn+898ZLyDg3OeeO5Hw/n/M5957znvd9bnLufc77PktUFRcXFxcXl7oQEWwBXFxcXFzCH1eZuLi4uLjUGVeZuLi4uLjUGVeZuLi4uLjUGVeZuLi4uLjUmahgCxAskpKSNCUlJdhiuLi4uIQV33zzzQFVbVXxeINVJikpKeTk5ARbDBcXF5ewQkS2V3bcXeZycXFxcakzrjJxcXFxcakzrjJxcXFxcakzDdZm4uLi0vAoLi4mNzeXEydOBFuUkCcmJoYOHToQHR1tqb2rTFxcXBoMubm5tGjRgpSUFEQk2OKELKrKwYMHyc3NJTU11dI1IbPMJSKviMg+EVlVxXkRkb+JyCYRWSEiA3zOXSMiGz3bNYGT2sWl4TBzJqSkQESE2c+c6cfFW2fC+ynwRoTZb/XnYvs4ceIEiYmJriKpAREhMTHRrxlcyCgT4DVgcjXnzwG6ebYZwHMAIpIA3AcMBYYA94lIS0cldXFpYMycCTNmwPbtoGr2M2ZYVChbZ8KSGVCwHVCzXzIjaArFVSTW8PfvFDLLXKo6T0RSqmkyDfinmpz5i0QkXkTaAmOBTFU9BCAimRil9KazEjtEQS7sng3tpkDTDjU237MH/vlPOH7cvI+KgosvhrQ0h+V0aRDk58Nrr8Hdd0NBwannCgrgrrvgyiuruPj4TtjzCXz7OyitcHFpASy/C1Krutgl3AilmUlNtAd2+rzP9Ryr6vhpiMgMEckRkZz9+/c7Jmit2TkLZveDJT+FD1JhwVVw+LtKm+7fD7/9LXTpAnfcAQ89ZLb77oM+feDqq2Hz5gDL71JvKCiAxx6D1FS45RbIy6u83Y4dFQ4c3QAr7jX38QedzL1cXMXFBRUvbjjk5uYybdo0unXrRpcuXbjttts4efLkae12797Nj3/84xr7mzJlCkeOHKmVLPfffz+PPvpora71JZyUSWVzLq3m+OkHVV9Q1UGqOqhVq9OyAQSP0hOw5Ofw1YXQvAtMmAPdb4bcWfDpQNg375TmmZlwxhnw+ONmFrJpE5SVmW3fPvj1r+Gdd6BHD3j++SB9JpewZf166NkTfvMbGDAAFi2Czp0rb9u8OWhpKWx+GT4bCh/1gNUPQ6MEOPOvMHUNNO1U+cWNk5z7ECGMqnLhhRdy/vnns3HjRjZs2EB+fj533XXXKe1KSkpo164d77zzTo19zp49m/j4eKdEtkQ4KZNcoKPP+w7A7mqOhw8rH4BNz0Pa72Di19B6LAx8EqbtgCYdIOcWKCsBICsLzjvPKJPVq+H1183sxEurVvDXv8KWLTB5Mvz85/DCC8H5WC7hx8aNMG4cFBVBdjZ89hkMHQoPPwxNm57aNioK2jVfw+5XBsDiG6C0EM58FKbthLOyIe03EJcG6X+EyAoXI1B0qMqZd30mKyuLmJgYrrvuOgAiIyN54okneOWVV3j22We5+OKLOffcc5k0aRLbtm2jT58+ABQUFHDJJZfQr18/Lr30UoYOHVqeEiolJYUDBw6wbds20tLSuPHGG+nduzeTJk2isLAQgBdffJHBgweTnp7ORRddREHFdcs6EjI2Ewt8CPxCRN7CGNvzVHWPiHwG/NHH6D4JuDNYQvpN4V5Y/xR0vhzOfOTUc40TYMBjMP9i2Pwic3f/nB/9CLp2hS++MIqjKtq0gXffhQsvhJ/+1Hzxr7/e2Y/iEt5s3mwUSXExzJljlku9eO0id91llrY6dSzl4Wuf48oet7D9QCde3f1frr3nx0hEJQsFXrvI8rvM0lbTTtDrDlj9IHx1MZzzDUTHOv8BK/DLX8J3Nuuy/v3hySerb7N69WoGDhx4yrHY2Fg6depESUkJCxcuZMWKFSQkJLBt27byNs8++ywtW7ZkxYoVrFq1iv79+1fa/8aNG3nzzTd58cUXueSSS3j33XeZPn06F154ITfeeCMAd999Ny+//DK33HJLnT6vLyGjTETkTYwxPUlEcjEeWtEAqvo8MBuYAmwCCoDrPOcOiciDwFJPVw94jfFhwZo/Q1kR9L2/8vMdL4LW4yhZdjdX3XoJqamJfPll9YrES+PGRqGcfz7ccAO0bAkXXGCr9C71hIMHYfx4KCw8XZF4ufJKj1LZOxcWTofC79He9/P0q7/j0SebsFPg3nurGCD1ytON7fG94ctxsPhGGPkWNBAvK1Wt1FPKe3zixIkkJCScdn7+/PncdtttAPTp04d+/fpV2n9qamq5ohk4cGC5Qlq1ahV33303R44cIT8/n7PPPtumT2QIGWWiqpfXcF6Bm6s49wrwihNyOUpBLmx8DlKvgdjulbcRoTj9KeTTM/nDRfdyzr3PkJxsfYiYGJg1C8aMMa6cI0fi1/UuDYNbb4Xdu2HBAqjiNwq0DFbeD6seghZdYdJCJHEQf3kcdu+DBx6AqVOhwkN31SSPhn4PwfI7ITkDut9k06exRk0zCKfo3bs377777inHjh49ys6dO4mMjKRZs2aVXmd+AmumcePG5a8jIyPLl7muvfZa3n//fdLT03nttdfIzs6u3QeognCymdQ/Vj0ElEGfe6pt9ufn+vJs5k1cO+p52jRd7/cwTZoY986jR+EXv6idqC71l/ffhzfegHvugcGDq2hUUgDzL4FVD0Lq1TB5GSQOAsyE4umnoXVruPZaqMQpqWp6/Q7aTobvfgcnDtT1o4QFEyZMoKCggH/+858AlJaWcvvtt3PttdfStKJhyodRo0bx3//+F4A1a9awcuVKv8Y9duwYbdu2pbi4mJl+RZxaw1UmwSJ/i/GA6TIDmqdU2WzFCnjwQVgbcZeZGm95rVbD9e4N998Pb79tNhcXgEOH4Gc/M2v9d1ZlaSzca5ajdr4HAx6HYa9CdPNTmrRsaRw9Vq0yLuqWkQjTZ2khrP1LrT9HOCEizJo1i7fffptu3brRvXt3YmJi+OMf/1jtdTfddBP79++nX79+PPLII/Tr14+4uDjL4z744IMMHTqUiRMn0rNnz7p+jNNR1Qa5DRw4UIPKst+ovhmtenxXlU1OnlQ980zV5GTVAwdUdc4U1VkdVctKazVkcbHqwIGqrVqp7ttXS7ld6hVXXqkaFaX67bdVNMjfofp+qupbTVR3zKqxv6uvVo2MVP3mGz8F+Xq6GaNgj58X+seaNWsc7d9JSkpKtLCwUFVVN23apJ07d9aioiJHx6zs7wXkaCW/qe7MJBiowo53oPVZ0LRdlc3+8Q/49lt47jlITARSpkPBztPiTqwSFWWWu44cMRHNLg2br74y6VB+/3szMzmNE/sg6yw4eRAmZEPH82vs88knjU3u5pvNbW6ZvvdB2UnjkOJSKQUFBYwaNYr09HQuuOACnnvuORo1ahRsscpxlUkwOLwMjm+DTlVHthYWwh//aAzn5R5YHaZBVAvY+q9aD92njzHEv/IK+HgdujRA7r3XuJDfcUclJ08ehqxJ5uEl42NIGmKpz5YtzXLqokXw6ad+CNOiq3FE2fg8FOzy48KGQ4sWLcjJyWH58uWsWLGCc845J9ginYKrTILBjndAIo1yqIJ//MPk3XrgAR+Pyaim0Oki2PkOlBTWevg774TISD/Xtl3qFXPmmKDE3//eOGicQkkBZE+Fo2thzPuQPMqvvq+91mQVvvdeP2cnfe4BLYXV1dsOXEITV5kEmvIlrgnQOLHSJsePw5/+ZPz+MzIqnEyZDsVHYdf/ai1C+/bG6Praa27+roaIqvmhb98ePDFsP1BWAl9fBgcXw8g3oe0kv/tv1Mh4huXkwEcf+XFh8xQ44zrjmHLysN/jugQXV5kEmiMrIH9TtUtczz1ncmz94Q+VnEweC03a12mpC8zSRqNGxlPMpWHxxRcwf76JZo+J8TmhCktvMg8qg56GjhfWeoyrrjJpfvyenXS/yQTxbv13rcd2CQ6uMgk0O94x7pAdKjdm5ufDI4/ApEkwqrLVhYhISLkC9nwKJ2qf+bhNG5O361//gg0bat2NS5jhnZV06lRJep1VD8HmF6HXndDt53UaJzrajPPddyaOxTIt+0PCQCOHX1rIJdi4yiSQqMLOt83sIqbyfCgvvQQHDlQxK/HS+XLQEtj9cZ3E+b//MylX/tIw3PtdMHaSRYvMrMQnUBpy/wcr74WUqyD9YVvGuuIK6N7dOJL4RZcb4chKOLi05rZhyK9+9Sue9Am/P/vss7nhhhvK399+++08/vjjAZOnefPmNTeygKtMAkneaji6vsolrrIyePZZGD4chg2rpp+W6SZ99945dRInORmmTzfRz4fdJeoGwTPPGDfzq6/2OViwCxZfZ2YFQ1+0LUdWVJRJ05KTA0uW+HFhyuUmy/Dml2yRI9QYMWIECxYsAKCsrIwDBw6wevXq8vMLFixg5MiRNfajqpSVlTkmp79YViYicoaI/F5EnvG87yEivZ0TrR6S+z4g0KHybItffmlSgN9caQYyHyQCWo+DvVl1Xgq46Sbjhvzaa3XqxiUM2LXLLDn95Cc+tpKyUlh4lampM/ItiGxcbR/+ctVVpubJs8/6cVF0LHS+FLa/CcX5tsoTCowcObJcmaxevZo+ffrQokULDh8+TFFREWvXriUtLY0JEyYwYMAA+vbtywcffABQnmL+pptuYsCAAezcuZPmzZtz++23M2DAACZMmIC38N/mzZuZPHkyAwcOZPTo0axbtw6ArVu3Mnz4cAYPHsw991SfyskfLCV6FJGJwHvAHExm35uBJOBuTG12FyvsnQPx/aBJm0pPP/OMyQZsobCaUSY73ob8zcZHv5b07w8jRpgv+223QYQ7V623vPCCmf3+7Gc+B9f82dyXw16F2B62jxkbaxTKK6/Ao49CktV6WF1uhC2vwo7/QJef2C4XAN/80v56Ki37m1pE1dCuXTuioqLYsWMHCxYsYPjw4ezatYuFCxcSFxdHv379aNq0KbNmzSI2NpYDBw4wbNgwzjvvPADWr1/Pq6++yrMeDX38+HEGDBjAY489xgMPPMAf/vAHnn76aWbMmMHzzz9Pt27dWLx4MTfddBNZWVncdttt/PznP+fqq6/mmWeese2jW/3p+DNwsaqeB5R6ji0DBtgmSX2ntAgOLDCFryphxw743/9MqvjGVh4OW483+71ZdRbt5ptNtcbMzDp35RKinDxplMmUKaYULwAHFsHK+4wNLvUax8a++WZTbOsVf/J6Jw2DuF6w6UXH5Aom3tmJV5kMHz68/P2IESNQVX7/+9/Tr18/zjrrLHbt2sXevXsB6Ny5M8N81sEjIiK49NJLAZg+fTrz588nPz+fBQsWcPHFF9O/f39++tOfsmfPHgC+/vprLr/cJGm/6qqrbPtMVlPQd1FVbzyrAqhqoYhE2yZJfefgUrOUkFwxcMTwj3+Y/U9/arG/Ft2hSTv4Pgu6zqiTaBddBL/6lZmd2FziwCVEeP99+P57s6wJmOWjBdOhaQcY/JyjtUR69zbxUs8/D7ffbgJma0QEUq812YTzt0Lz1Bov8ZsaZhBO4rWbrFy5kj59+tCxY0cee+wxYmNjuf7665k5cyb79+/nm2++ITo6mpSUFE6cOAFQZYp6LyJCWVkZ8fHxfFdF9a/K6qnUFaszk50ickq5HBFJB7bZLlF9ZV+22SePOe1UURG8+CL86EdV19o+DREzO7HBbtK4sZkRffQRbN9ep65cQpRnnjEzksmTPQeW/dpkrh7+L2hkPfNsbbnpJti61c8UK15HlZ3vOSJTMBk5ciQfffQRCQkJREZGkpCQwJEjR1i4cCHDhw8nLy+P5ORkoqOjmTNnDtur+WKWlZWV14l/4403GDVqFLGxsaSmpvK2J0W4qrJ8+fLysd966y0AW1PRW1UmfwPeE5HpQKSIXAT8G3jCLkFEZLKIrBeRTSJyWrYgEXlCRL7zbBtE5IjPuVKfcx/aJZOt7Jtr7CWVRL2/9x7s32/B8F6R1uOhaL/xEqsj3hmRWy++/rFmDcybZ+KKIiKA3A898SS/MwWqAsAFF0DbtiYg1zLNU6HlmbDz3Zrbhhl9+/Ytt4X4HouLiyMpKYkrr7ySnJwcBg0axMyZM6tNGd+sWbPyUsBZWVnc6yl3OXPmTF5++WXS09Pp3bt3uRH/qaee4plnnmHw4MHk5eXZ96EqSyVc2QbMAFYC+cAq4Aar11roOxLYDJwBNAKWA72qaX8L8IrP+3x/xwxoCvqSIpNee+ktlZ4+5xzVTp1US/3NLJ+/TXUmquueqruMdZHDJaS54w6TFn7vXlUt3Kf6TivV2f3NfRksOayy8iFzjx/PtUWGcE5BXxXNmjVzrG9HUtCr6guq2ldVm6tqH1W10wl8CLBJVbeo6kngLaDqLIhwOfCmjeM7y6EcU/ynEuP7vn3w+ecmwMtvT6pmnaH5GbYY4cHU996xA77+2pbuXEKAsjITRzRpkqdc87e/geIjZnkrMrDpy6+8EkpLwVMs0BodLzL7XH/C6F2CgaWfLxFpV9VmkxztgZ0+73M9xyqTpTOQCvj+gsaISI6ILBKRKosuiMgMT7scry92QPDaS1qdbi/5z3/MF2z69Fr23Xo87M028QJ15PzzoVkz+LebFqneMH++eUCYPh34/kvY+k9I+z+I71PjtXbTpw+kp/t5f8X1hNi0ernUZRf5+aERi2P1WTgX82Nf2WYHlbkWVGVVvgx4R1V9fz07qeog4ArgSRHpUtmFntnVIFUd1KpV5elMHGHvXIjrAzGnO9n/+9/mC9a7tuGfrcdDcR4c/rZuMmIUyfnnmyfHoqI6d+cSAvz73+b/Ou1HJ2DJz6B5V+j9+6DJc+WVsHixCc61TMeLjM2xDrnofFE355cl/P07WVUmqRh7hncbDXwGXOvXaFWTC3T0ed8B2F1F28uosMSlqrs9+y1ANnCmTXLVnbJiOPB1pS7BGzeaNBO1npXAD95h++1Zm5o+3VRi/OQTW7pzCSJFRfD228b43WzbH0226iHPQVTFAiaB4/LLjSPiG2/4cVGni0DLIPeDOo8fExPDwYMHXYVSA6rKwYMHiTklrXT1WIozUdWKfmnbReQazFJT3XKhG5YC3UQkFdiFURhXVGwkIj2AlsBCn2MtgQJVLRKRJGAkEDqpCw99AyXHK7WXzJxpvlie+KHa0bS9SUl/0J/kR1Vz1llmbf3f/zazFJfwZfZs82Aw47J1JtI9ZTq0OSuoMnXoAGPHmvvr3nsthrfEpxvb4M53oesNNbevdvwO5ObmEtBl7jAlJiaGDh06WG5vNWixMgoAq1ER1aKqJSLyC8xsJxLjqbVaRB7AeA543X0vB97SUx8r0oB/iEgZZqb1Z1VdY4dctrBvrtlXiC9RNcpk3DhTpKhOJA01xYxsICoKLrvMBFEeOQLx8bZ06xIEZs6E5GRlZMxtUNAUBjwWbJEAM/v9yU9g6VIYYqUasIhZ6lr/JBQfg+gWtR47Ojqa1FQHAiBdLBvgr6iw3Qj8D7DnFwxQ1dmq2l1Vu6jqw55j9/ooElT1flW9o8J1CzxeZume/ct2yWQLBxabdeqY5FMOL1liUphceaUNYyQOMTm6ThywoTMjU1ERvOvaPMOWI0dMep6Hf/4BEXs/h75/OO0eDBYXXWQCZf0yxLebYpaM65gp28U5rNpMHq6w3Qpsxz6bSf3lUA4kDj7t8H/+YyodXnSRDWMkDvWMZU/9h8GDoWtXI6NLePLBBxChhUzv9SuI620qGIYIcXEm28PbbxvXZUskjYCo5qYonEtIYkmZqGpqha2vql6nqrucFjCsKdwLBTshYdAph1Vh1ixjn4izI5NFwkBAzCzIBkTgwgthzhy3zkm4MmsWPHjFo8SUbINBf4eI0Eqjd+GFJlfYokUWL4hsZDwXd3/iVmAMUdyE405y6BuzTzxVmSxfDtu2GS8bW4huYZ4+bTLCg/myl5TAx3Ur5ugSBI4fhzVLd3DrhD9Bp4tNuYIQY+pUU9p31iw/Lmo3GY5vg2P++BW7BIoqDfAikknVsR7lqOokWyWqTxxaCojJL+TDe++ZaHdPeQJ7SBpqooRVbckAO3gwtGtnZK2T67JLwPn0U7hn2l1ERiqc+ddgi1MpcXEwYYK5v/7yF4u3bFtPSus9n0Jsd0flc/Gf6mYm84GvLWwuVXEwB2J7nuZ9MmsWjBrlSW9hF4lDoOigMcTbQESEcQ3+9FMoKLClS5cA8V1WDleN+jfS81cm5U6IcsEFsGULrFxp8YLmZ5jSC7tdu0koUuXMRFX/EEhB6h2qxvje9tSJ28aNsGoVPGFbvmUPXiP8wSV1qrzoy4UXmhonn3/uxpyECyeLlLNb/YajJ1sR2/e05NshxbRppurjrFnQr5/Fi9pONhmPS09ApPWAOhfn8ctmIiJNRKSjiHTybk4JFvYU7oYT359mfPeuEdv+4xzXGyKb2maEBxgzBlq2NEsRLuHBqs//x6juc9nW4n5TSz2Ead3alIz26/5qN9kkTd33lWNyudQOq3EmZ4jIAkz6+W3AVp/NpTIO5Zh9JcpkwABISbF5vIgo49VloxE+OhrOPdfEKxQX29ati1OUFdNu7+9Yv6cn3c+5MdjSWOLCC2HFCrPcZYnkDIho7LoIhyBWZyZPY5I6pgPHgH7A+8BPHJIr/DmYAxIJLdPLD+3ebVwhbfPiqkjSUJPwsfSkbV1eeKEJgJs717YuXRyibNMrtGm6nve3P0JM09ByBa4K73fBsldXVFOTTcJVJiGHVWUyFFMMaxWAqq4Gfgr81inBwp5DOWbpKapp+SFPoTPnlEniECgrgiPLbety0iRo2tRPF06XwFNSSPG3D/D1hhF0GnZusKWxTGqqyZrt1/3V9mzIWwMFbphbKGFVmZQBhZ7X+SISDxwCXJtJZXiN7xWWuD7+GM44A3r1cmjcRE+iI+8Smw00aWIUykcfubFiIc2m52lcupv73n2Ic6bU3TU8kEybBgsXwsGDFi/wxs3sc6fLoYRVZbIak40XTD6uJzB14V2bSWUc3w5FB05Jo1JYCFlZJljLhjCQymnaCaLj4bB9MxOAKVNMgaU1oZM+08WX4nxY/ScWbZtASeK4sEvOOWWKSavy2WcWL4hPh+g4UxTOJWSwqkxuxcxEwCxttQcGYZa6XCpSifE9O9solClTHBxXBFr2gyMrbO32nHPMfvZsW7t1sYv1T0HRfm575WFn7y+HGDQIkpL8uL8iIo3dxFvB1CUksJqba4WqrvS83qKqk1R1mKourOnaBsmhHJMLKb5v+aHZs82SUcbpNbLsJT7dKBO1mkGvZjp0MHEArjIJQU4ehrV/ZXvJuSzZPDQslUlkJEyebAJkS61Wn04ea9KquHaTkMGqa/AXntTzbpSQFQ4v98R9NAaMrWH2bJM+oonTRe5apptiXPlWfS2tMXWqqSeel2drty51Ze3jUJzHk1kP0qlTHco/B5mpU43NZKnVxNfeYnOu3SRksLrM9TXwELBHRJ4TkdNzqtcREZksIutFZJOInBa6KyLXish+EfnOs93gc+4aEdno2a6xWza/ObIC4n8I6d2wwfjRB+Sp0TuuA3aTkhLIzLS1W5e6UHQQ1j9Jafsf89K76UyZ4qA9zmEmTTIpfCzPfsvtJm59k1DB6jLXfap6BnAR0AyYIyKrRORXdgghIpHAM8A5QC/gchGpzOfpP6ra37O95Lk2AbgP4748BLjPU8o3OJw4YKLffZSJN/Ou1/bgKHF9QCJsdQ8GGDbMVF10l7pCiLWPQclxlp64j/z8AD2sOERCAgwf7keWaq/dxDXChwx+pVNR1SxVvRpTrnc78KhNcgwBNnnsMSeBt4BpFq89G8hU1UOqehjIBCbbJJf/eI3fPspk9mzjDmx71HtlRDUxyfBsnplERcHZZ8Mnn/hR0MjFOU7shw1/g86X8vbnfWjUCMaPD7ZQdWPKFFi2DPbssXhB63GQvwkKch2Vy8Ua/ubmSheRJ4G1QF/gTzbJ0R4TYe8l13OsIheJyAoReUdEOvp5bWCooEyOHYN58wL81BifbvvMBMxn+P57+PZb27t28Ze1fzU5qvrcx8cfw9ix0KxZsIWqG97vyKdWg9uTx5r9XtduEgpYNcDfIiLLgIVAa2A60FlV77ZJjspWeiuGyP0PSFHVfsAXwOt+XGsaiswQkRwRydm/f3+tha2WIysgpjU0aQ3Al1+avFZTpzozXKW0TDexLieP2Nrt5MlmTd5d6goyhXthw9PQ+XI2H+jJ+vUBvr8cIj3d1NCxbjfpZ+KqXBfhkMDqzORq4EWgnaperqqfq9oaD50LdPR53wHY7dtAVQ+qapHn7YvAQKvX+vTxgqoOUtVBrVq1skXw06hgfP/kE2jRAkaOrOYau4n35AM7YrVQhDWSk01MwCef2Nqti7+se8ykzelzb/lTfEDscQ4jYj7H558bZ48ace0mIYVVA/xgVX1OVe191P2BpUA3EUkVkUbAZcCHvg1EpK3P2/MwS20AnwGTRKSlx/A+yXMs8JSVQN7qU5RJZiaMG2cy8AaMls54dIHxulmyxHURDhonD8PG56DTJRDbncxMY4vrak8Jm6AzaRIcPeqHi3Byhsdu4sabBJuQqAGvqiXALzBKYC3wX1VdLSIPiIi3uO2tIrJaRJZjIvKv9Vx7CHgQo5CWAg94jgWeY5tM0R6PMtm8GbZuNV+QgNKkPTRKcMRuMnGiCSyb43pkBocNz0BJPvS6g+Ji83+YODF8XYIrMmGC+SyWXdBbjTL7/W7R12ATEsoEQFVnq2p3Ve2iqg97jt2rqh96Xt+pqr1VNV1Vx6nqOp9rX1HVrp7t1WB9hvIfb8/M4PPPzduJEwMsh4ixmzgwMxk+3Bh63XiTIFByHNY/Ce2mQMt0liwxT/EBf1hxkMREGDjwh+EDwk4AACAASURBVO9OjSScCZFNXGUSAoSMMqkXHFlhapjEpgHmB7dTJ+jWLQiyxKdD3ioos5qfwhqNGhnPIVeZBIFNL5lAxV53AuZ/IBL+LsEVmTjR1P05etRC44hoU7L6gKtMgo2rTOzk8AqI7QmRjSkpMVmCg7YE0TLduI4e22h71xMnmlr227bZ3rVLVZSehHWPQqvRkGyWdjIzjUNEQkKQZbMZ71JqdrbFC1qNhMPfmezJLkHDsjIRkVEi8oKI/M/zfqCIjHFOtDDEx5MrJ8cYqYO2BFHu0eWM3QTc2UlA2TbTBOf1NrOSvDxYvDgIS6gBYMQIU5DNut1kJGgpHFzsqFwu1WM1zuQK4APgBOBVIAo84JBc4cfJI1Cwo7xM7+efmxnJhAlBkicuzZNWZZXtXaelmXgAV5kECC2DtY+YB4S2JrnDnDnm6b0+KpPGjU12bct2k6ThgLh2kyBjdWZyF3C2qt6KqboIsAoI0xylDlAh8j0zEwYMMAbFoBAZA827GldlmxExM64vv/QjZbhL7cn9EI6uh153lK+ZZmYaR4jhw4Msm0NMnGgSpO7YYaFxo3iI7+MqkyBjVZm0U1VvLVhvsGIJEGm/SGHK4R+UybFjxoAYdC+buN6OKBMwX/ZDh9zUKo6jCmv+BM1SodOPyw9nZpqn98aNgyibg3i/O5Znv0kj4cBC2x1OXKxjVZlsFpERFY6NANbbLE/4cmQFNGoJTdqRnW0ieIO+BBHX2wR0lZ6wveuzzjJ7d6nLYfbNhYNLoNdvISIKgO3bjQNE0O8vB+nVy8+l1FYjoeSY8WB0CQpWlclDwAcicjcQLSK3A2/i2kx+IG+1J/27kJlpimCNqKh+A018H7PefnRdzW39JDnZ5FJylYnDrHkEYpIh9dryQ96/eX1WJiLmgeWLLyxmqXaDF4OO1XQq7wNXYGqGbAfGA9erqpulCcxSRN4aMxPAuASPHh0CSxAeeTjizFLXWWfBggWmtr2LAxz+DvZ8Cj1+aUoLeMjKgjZtzNN7fWbCBFN9caWVFHPNOkOTdrB/vuNyuVSOZddgVc1U1XNVtY+qTlVV95nUy4nvofgIxPVi715YvTpEAsladAeJcsxuMm4cFBXBwoWOdO+y9lGIag7dfl5+SNUok3Hj6k8KlaoYN87ss7IsNBYxS13uzCRoVKlMRKSdlS2QwoYs3h/ruN7lOatCQplENoIW3RxTJqNHQ2SkxS+7i38c3wnb34IuNxpvJQ/r1sHevSFyfzlMx44me4Tl+ytppHHPd4tlBYWoas7lUkVdEA/iOe96dOWtMfu4XmRlQVwcnHlmcEUqJ643HHbG5So2FgYPdpWJI2z4m9n3vO2Uw96/dUNQJmA+5xtvGIeWqOp+rQBaeYyUBxad4vnmEhiqW+ZKBc6oZvOed8lbbbL0xrRmzhzjslnjjR8o4vtA/hYoKXCk+/HjTbrwY8cc6b5hUnwUNr0AHX9sbAE+ZGVB586Qmhok2QLM+PHm3lq2zELj+HSIaGyUiUvAqVKZqOp2K1sghQ1Z8tZAXC927BQ2bfphrTckiOsNqCMeXWC+7CUlMN+1e9rHppeMQkm7/ZTDZWUm8n38+PpvL/EydqzZW5r9RjaChIEm3sQl4PiTm+tSEflURFZ59pc5KVjYoOpxC+4VWvYSL+UeXc74348YYTIJu0tdNlFWAuufMhUEEwefcmr5cjh8OMTuL4dJToY+ffyxmwyHQ9+YxJguAcVqbq7/A/6GKT71hGf/pOd4w+bEXlP9Lq43WVmQlGRu/pChRVeTptshI7w3nsZVJjax4x1jRO55+2mnvA8rITXzDQDjx5uZb1FRzW1JGmZKGjuQ4NSleqzOTG4GzlHVe1T1ZVW9BzgHUx3RFkRksoisF5FNInJHJed/LSJrRGSFiHwpIp19zpWKyHee7cOK1zqKx/iusb3KXTYjQimxf0Q0tOjhmDIB82X/9luTXsWlDqjCuseNB177H512OisLevSA9u2DIFsQGT/exDIttpIUOGmY2btLXQHH6s9eU2BFhWMrgSaVtPUbEYkEnsEoqF7A5SJSMSTrW2CQqvYD3gH+4nOuUFX7e7bzCCSeH+mth3qTmxuiSxDxfRxXJqowd65jQzQMDiyAQ0tNkKKc+tUsLjZ/35C8vxwmI8M8oFma/TbtYDbXCB9wrCqTfwK/qnDsl8DrNskxBNikqltU9STwFjDNt4GqzlFVr0vSIqCDTWPXjbw1EB1P5vw2QIh+2eN6w/FtjhUPGjzYZLB1l7rqyLrHTX63M6457dQ330B+fojeXw4TH28ycFu+vxKHucokCFQXtJgpIp+LyOdAP+BhEdkmInNFZBvwMJBukxztgZ0+73M9x6riJ4BvKpcYEckRkUUicn5VF4nIDE+7nP3799dNYi9H10B8b+bMEdq1C1KJ3prwGuGPrnWk+0aNYNQoV5nUifytkPs+dP0ZRDU77bTXXpKREWC5QoRx40wm7gIrHu5Jw+D4Vijc67hcLj9Q3cxkPvC1Z/sKozxeBbI8+z96jttBZY6OlQZMish0YBDwV5/DnVR1ECZ/2JMi0qWya1X1BVUdpKqDWrVqVVeZyz25NLYX2dnGjTEkXTYd9ugC89nXrIF9+xwbon6z/m9ABHS/udLT2dnGscOO2zYcGTvWLPVZSt2T5CnyctCdnQSSKkPrVPUPAZQjF+jo874DsLtiIxE5C1OoK0NVy307VHW3Z79FRLKBM4HNTgoMQNF+KDrIviKTk8vrEx9yNO8CEY0cm5nAD5993jz4sRt87B8n82DzS9D5Mmh6+oS8uNh4M11/fRBkCxFGjTJ2k+xsC9VLW55pHE8OLIIO02po7GIX/sSZRIpITxEZLSJjvJtNciwFuolIqog0Ai4DTvHKEpEzgX8A56nqPp/jLUWksed1EjASWGOTXNXjMWovXW+e/ENWmUREmqSPec79WQYONHaT7GzHhqi/bH4ZSvKhZ0WzpCEnxyzvhOz9FQBiY809Zun+imoC8f1dj64AYynph4gMAN4DOmGWn7x5uUqBRnUVQlVLROQXwGeYXF+vqOpqEXkAyFHVDzHLWs2Bt8WsJe3weG6lAf8QkTKMcvyzqgZImZhhZi/oRbt20LVrQEatHXG94FBOze1qSXS0eXp0lYmflJXChr9Dq9GQMKDSJt6/6Ri7Ht3ClLFj4amnjGJt2rSGxknDzWyvrKS8qJiLs1idmTwJzALigKNALGaWcK1dgqjqbFXtrqpdVPVhz7F7PYoEVT1LVVtXdAFW1QWq2ldV0z37l+2SqUby1qDRccz6tF3o2ku8xKYZI2+Jc8VHxo416fft8m1oEOz60HjaVTErAeMS3JDtJV7GjoWTJ40hvkaShkFpARyxUgzFxQ6sKpO+wB2qegwQVc0HfkdDr7SYt4YTjdL4/nsJfS+buDRA4dgGx4bw/g3mzXNsiPrH+iehWQq0rzw8ymsvCfn7KwD42k1qJGmo2R9c4qRILj5YVSbFPq/zRCTZc6yN/SKFEUfXsv1IGhAG69mxRk7ynDPCDxpklh/cpS6LHPoW9s2D7rcYu1YlfPMNHD8eBvdXAIiNNfEmlu6vZqnQOAkOWgmbd7EDq8rkG8BbcTob+BcmsLBiVHzD4eRhOLGXZRvTaNs2RONLfIntbqKqHfTocu0mfrL+KRNT0qVqNy3XXnIqY8eatCo1looWgcQh7swkgFhVJjcA3sxpvwZ2AEXAdU4IFRZ4nvA/XZgW+vYSgMgY87TmoDIB82Vftcq1m9RI4few/U0447pTKilWJDsbevc22XNd/LSbJA41TjLFR50WywWLykRVd6nqTs/rg6p6o6pepqrrnRUvhPH8KH+9Ki18liBi0xxd5oJT401cqmHj81B2ErrfWmUTr70kbO6vAOCX3SRxKKBw0DkvRpcfqC6dylCf1yOq2gIjZgiSt5YSbcy2/SnhYxyNSzMG+LISx4bw2k3cpI/VUFoEm56HdlMgtur10WXLjL0kbO6vAOAtiW1NmXjqwbh2k4BQnQP2F0ALz+uq6ug13BrwR9ey61gPkltH0r17sIWxSGyaeRrO31rtj1hdiI6GkSN/yCXlUgk73jZ1cHrcVm0z7w+mq0xOZexY+Pvfjd2kSXV5yxsnmHT+rt0kIFRXtreFz+uIKraGqUgAzVvLd1vSyMgIA3uJlziPR5fDdpOMDGM3OXjQ0WHCE1VjeI9NgzYTq206dy6kpbn2kopkZBi7yRIrOiJxqJmZaKWp/lxspEabiYhEicjHIhITCIHCgpJCOL6NZZvTwuupMTZwygTgK7vSgNYnDiwymQh63FLtU0hJiRtfUhWjR5s/naWl1MShULgHCnIdl6uhU6MyUdUSYCDg3EJ7uHFsPYKydleYKZNGcdCkneNG+MGDISbGtZtUyoa/QXQcpFxVbbPvvoNjx1xlUhnx8ZCeblWZDDF7d6nLcay6Bv8LG0v0hj2eH+PvC9JISwuyLP4Sm+b4zKRxYxg2zFUmp1Gwy9R47/ITiG5ebVPv386NL6mcMWNMOvqTJ2to2DLdZMx2jfCOY1WZDAD+4qnP/oW3aJancFbD4+haSssiaNuje/jYS7zEedyDHV5DzsgwT9d5eY4OE15sfB60FLrX/Fw2b55JHNquXQDkCkMyMowBPqcmr9/IxiYlvatMHMeqMpmHKY71L0xBrK99tgbH8T1r2bLvDEaObhxsUfwnNg1KjkHhLkeHycgw+mp+VX6ADY3SE7DpH9D+XGieWm3TsjJjb3KXuKrGO2OzbDc5mOOoS7yLxRT0AS6UFfIUH1hr7CVXB1uSWhDnk6OraQfHhhk2zLgJz50LU6c6Nkz4sP0/pphaDe7AACtXwuHDrjKpjqQkkxlg7ly4884aGicOMbaqvNVm2cvFEUKlOFb4UFZCc93A1oNp9O0bbGFqQYA8upo0gSFDXLsJ8IM7cFxvaD2uxuauvcQaY8bA118bz7dqKc8gvNRxmRoylpSJpzjWZkwFw2zPNgcT2NiwyN9CVEQxGpdGhGVVHELEtIboeDi6zvGhMjJM1tv8fMeHCm0OLIDD30KPWy0FJc2bB507m82lajIyzL317bc1NGzeBRoluHYThwmZ4lgiMllE1nuM/HdUcr6xiPzHc36xiKT4nLvTc3y9iJxtl0yVcWibeaJv1SXc3Lg8iPxghHeYjAwoLYUFCxwfKrRZ/zdo1BJSrqyxqapRJu4SV814/0Y1zn7dDMIBISSKY4lIJPAMcA7QC7hcRHpVaPYT4LCqdgWeAB7xXNsLUzO+NzAZeNbTnyPsXG1+hHsN7enUEM4TAPdggBEjIDKygS91FeTCzneNO3BUsxqbr11rMi67yqRm2rSB7t39iDfJWwUlxx2Xq6ESKsWxhgCbVHWLqp7E1EqZVqHNNOB1z+t3gAliisFPA95S1SJV3Qps8vTnCIXfr2X3kXb0HRjn1BDOE9vT5IY6edjRYZo3N4kfG7Qy2fgcoNDtZkvNvX8rV5lYIyPDeL6VltbQMHEIaBkcWhYQuUKVNfOXsvAvF7NjzRbb+w6V4ljtgZ0+73M9xypt44nKzwMSLV4LgIjMEJEcEcnZX8uCGyfLmrItfxRRlvzgQpQ456suehkzxuRQqrGYUX2kpNDHHTjF0iXz5pnYkjPOcFa0+sKYMSaWaWVNpd7LI+Ebtt1kz/K5DO/wDs3iqg+arQ21LY61HXuLY1VmlawYVVdVGyvXmoOqL6jqIFUd1KpVKz9FNIy5/TlG/O4/tbo2ZCj36AqMEb642GIxo/rG9reg6KAld2Aw9pK5cwmv5KFBxrLdJKaVKQ7XwO0m0XlLyD2SQmJ7+7OH1rY41gybi2PlAh193ncAdlfVRkSiMM4Ahyxe6+JLsxSIaBwQu4m3mFGDW+pSNbENcX0geaylSzZtgj173CUuf+jYEVJT/bCbNGBlUlICKS0Ws7fYGSuAVdfgTSLyexGpdPnIBpYC3UQkVUQaYQzqH1Zo8yFwjef1j4EsVVXP8cs83l6pQDeg4d4xVoiIhNgeAVnmiouD/v0boDLZPx8Of2fZHRhce0ltycgwy4NlZTU0TBwCx7dD4d6AyBVqrFz6PZ0SdxCRPLTmxrXA6jLXw8AkYKuIfCIiF4tItF1CeGwgvwA+A9YC/1XV1SLygIic52n2MpAoIpswS213eK5dDfwXEwPzKXCzqtZkjnOJ7RmQmQmYde1Fi6CoKCDDhQZ+uAN7mTfP1C7p0cNBueohGRmmds7amm7n8uDFhvmsuS3HfO5OZwZxZqKqr6rqWCANyAH+CuwRkafsEkRVZ6tqd1XtoqoPe47dq6ofel6fUNWLVbWrqg5R1S0+1z7sua6Hqn5il0z1mtg0U3GxxHnLeEYGnDgBSxtKAPLxnZA7C7rcAFFNLV82d65RvK69xD8s5+lqeSZIZINVJif3LKGkLJLErgMc6d+vGG5V3ayq9wDDgMW4aenDl7g0QOHYRseHGj3a7BvMUle5O/BNli/Ztg127HCXuGpDaip06GDh/opqCvF9G6QyKS2FVpFL2H28r18POP7gb26uaSLyPrANiAd+6ohULs4ToBxdAImJ0LdvA1EmJYWw+QVof55ld2Bw7SV1QcT83ebOtVBZwWuE15oMLPWLlSvKGNBpCUXNHQvBs2yAfxzjIfU0xqbRT1VHqupLjknm4iyx3QEJiBEezFLEggXGTbhes/0Nv9yBvcydCwkJJhOui/9kZMDevbBhQw0NE4dC8ZGAzMhDieVfbyS+WR6J3Z0xvoP1mUlHTB6uzqp6p6rW9C9zCXUiY0xdjQAZ4TMy4PhxWFafA5BVjeE9vh8k+zfFmDfPLAeGZfLQEMA7o5s3r4aGDbSMb95mE6yZ0C3IMxOP4fsT1QY2N6zvBChHF/hZzChc2TcPjqzwyx0YYNcu2LzZXeKqC926mVxdNd5fsWkQ1RwONJxIeFVoVriEwpLmPyxvO4D7HNSQiUuDoxugzHlP6tatoWfPeq5M1j9lUp13vsKvy1x7Sd0RMQ8sNdpNIiIhcXCDSquyZg30bb+EIxGDzOd3CFeZNGRi06CsCI5vC8hwY8aYMr41JuULR/K3wa4PoOsMiGri16Vz50JsLKS7RQDrREYG5ObC1q01NEwcCkeWm1LKDYD5c0/Qv/N3NGnv3BIXuMqkYRNAjy6AsWPh6FH47ruADBdYNj4LiF/uwF7mzjX2kkjnHhobBGPHmn2Ns9/EoVBWDIdqqqpVP9ixYjmNooqJ6+Kc8R2se3NVWkBUREbbK45LQInz1GQJkEeXdxknOzsgwwWOkuOw6UXocAE061hzex/27IH163/4IXSpPWlp0KqVhfurPBK+/i91qVJuH5Kk0JiZfFTF8Q/sEsQlCDRqCTFtAjYzadfOFDOqd8pk6z+Nu2nPX/p9qfcp2lUmdUfE/B2zs2uwmzRpC007NAhlsm4d9G6ziOPa3nxmB7GqTE5zTRGRFoDr3RXuxPWCvDUBG27sWOO+WW/sJlpm3IETBkHSCL8vz8429pL+/e0XrSGSkWEyCWzbVkPDxKENwqMrOxuGdlmMJji7xAU1KBMR2SgiG4AmIrLBd8MUpJrjuIQuzuJVJjWGDttDvbOb7Pnc1IXpcVutkmplZxt7SVgXWwshvDO8Gme/iUPh+FY4UbsieeHCsgX76NJ6C806D3N8rJpmJg8Bf8SU6H3YZ3sQuBiTKt4lnInrBSXHoHBXQIazXMwoXFj/lFkq7HSJ35d+/72xl7guwfbRqxckJVm4vxqA3UQVTuzy2kucVybVPg+p6usAIrJOVRtirbz6T2wvs89b4/iaKhi7Sbdu5snx1792fDhnyVsHez6Ffg9CZCO/L3ftJfbjazeploSBJoPwgcXQ/kcBkCzwrF8PPZIWU6aRRCQMdHw8qxHwi0QkRkT6isgI381pAV0cJs5HmQSIemM32fA3U7Gya+3ynWZnQ4sWcOaZ9orV0MnIgO3ba7CbRDUzVTDr8cwkOxuGdV3EyWb9HMsU7ItV1+DzgD2YOvDzfbavnBPNJSDEtILGSQFXJnl5sHx5wIa0n6JDsOV1U/wqplWtunDtJc5g2W6SNBQOLq23GYTnZpcytOsSGrdzfokLrHtzPQb8AWiuqhE+W53DrEQkQUQyPcb+TBFpWUmb/iKyUERWi8gKEbnU59xrIrJVRL7zbK5fjL/E9YKjgVMm9SLeZNMLUFpQK3dgMPaSdevcJS4n8NpNLBnh62kGYVXYvXYdLWKOBcReAtaVSWtVfVJVCxyQ4Q7gS1XtBnzpeV+RAuBqVe0NTAaeFJF4n/O/VdX+nq2++AkFjtjAenS1b2/sJmFrhC89CRv+Dm0mmmJLtcCb3dZVJvYTEfFDfZNqSfQY4Q/UP3Pwhg3QtaXncyU57xYM1pXJ5yLilHqbBrzuef06cH7FBqq6QVU3el7vBvYBtVtbcDmduF5w8jCc2BuwIceONV/2sLSb7HgbCndDz1/VugvXXuIsY8cam0m1dpO4NIiOgwMLAyNUAMnOhqFdF1Ma2RJadAvImFaVyTbgQxH5u4j83nezQYbWqroHwLNPrq6xiAwBGgGbfQ4/7Fn+ekJEGldz7QwRyRGRnP3767d/uV8EwQg/bpyxm4RdvIkqrHscYntC27Nr3U1Wlkl86dpLnGHcOLOfU10knER4ghfrnzLJyoLRaYuISB5qPmcAsDrKQGA10AeY6LOdZeViEflCRFZVsk3zR1gRaQv8C7jOp7bKnUBPYDCQAPxfVder6guqOkhVB7Vq5U5syonzlPcLsBEezE0fVuz/Cg4vgx6/rPWXdNcu47bp/cFzsZ9evSA52cL9lTQc8lZB8bGAyBUIyspgyYJj9GizCkkMzBIX1BBn4kVV63Tbq2qVSkdE9opIW1Xd41EW+6poFwt8DNztG/PindUARSLyKvCbusjaIIlpbfJ0BdAI37atScyXlQW//W3Ahq076x6HxomQelWtu/A+LY8fb5NMLqchYpR1VpaZTFaZnCBpuPHmOrgE2kwIqIxOsXo1pMTmECEKATK+gx8p6EUk0hNbcqnnfVMR8a9wQ+V8CFzjeX0NlSSPFJFGwCzgn6r6doVzbT17wdhbVtkgU8NCJOA5usD8mH71VRjVhT+6EXI/hK4/q5PfflYWtGzp1i9xmvHjYffuGurCe43T9WipKysLhnf1fJ5EZzMF+2I1zqQL5kd6NvCy5/Ak4EUbZPgzMFFENmKWzv7sGXOQiLzkaXMJMAa4thIX4JkishJYCSRhUsC4+EtscJTJ8eOwdGlAh60965+AiGjo/os6dZOVZZ6a3XrvzuKd+VW71NUo3tT1qUceXVlZcFb/BeZzNU4I2LhWb+e/A29hbBLe58hsoM71TFT1oKpOUNVunv0hz/EcVb3B8/rfqhrt4/5b7gKsquNVta+q9lHV6aqaX1eZGiRxvaBof0AT32VkmElRWNhNig7CltcgZTo0aVPrbrZuNdHZ7hKX83TpAh071mCEB7PUdXBRwFzjnaS0FObNK2PIGQuhVWATlFhVJkOAhz1GbwVQ1SNAfLVXuYQPXo+uANU2AUhMNEs9YaFMNj4HpYXQs24Jxbyf1TW+O4+IUdpz5hijdJUkDTcPC/UgePHbb6F1kw00jz5Uq5IIdcGqMjlKBcUhIu2AwAUmuDhLENyDwXzZFyyAE6Fcjrv0BGx4GtpOhvjedeoqKwtatzbOBy7OM348HDgAq6qzpCYNN/t6YDfJyoIR3ReYNyGqTN4DXhGRDgAikgg8iVn6cqkPNGkPUS3gSGD9F8aPh6IiWBjK3+NtM01AZ1rdHAVVzZd9/PhalT5xqQXeGWC1s996FLyYlQVThyyARgkQ2z2gY1tVJvcA+cAOzAxlH1CEqXXiUh8Qgfg+xuc+gIweDZGRIbzUpWXGHTg+HVrXzdCxbp3JyeXaSwJHx44mdU+191c9CV48edJ4R47svsDMtgIUrOjFagr6QlW9AhOdPgRoo6pXqWooL064+Et8XziyMqCGyNhYGDw4hJXJro/N0l/ab+o8nfB+RleZBJbx403qnpKSahrVg+DFJUugsRyiTdO1ATe+gx9xJgCqesDjZeXmIqmPxPWFk4egcE/NbW1k/HjzRTgWat9jVVjzJ2iWAp3rXlT0yy+hUydITa27aC7WGT/elIr+5ptqGvkGL4YpX34Jw7t5kzuGkDLxpHXfUtMWSGFdHMabAffIioAOe9ZZ5qkx5FLS759vlj563g4RdUuiVVJiZiYTJ7r2kkDjtVFlZlbTKGkYILD/60CJZTuZmXDRmAWmgmTi4ICPX93M5G6MreQe4FWgCfAGpgb8m0AM8IrTAroEkHJlsjKgw44YAU2b1vBlDwZrHjGFw7pcX+eucnJMYsuJE22Qy8UvkpJMduZq769GcRDfz+ReC0OOHoVFiyCj9wJo2d9UkgwwVT5uqepM72sRmQecq6o5PsfexXh0uRHn9YXGCcarK8DKpHFjk0H3888DOmz1HF4Buz829d1tKHmamWmejifUj/RPYcfEifD442YptUWLKhq1GgVbX4OykjrPRANNdjagJaS0WAxJPwmKDFZtJv2BisnCV3iOu9Qn4vtCXmCVCZgv+/r1sHNnwIeunLV/MU933W6ypbvPPzdPx0lJtnTn4icTJ5occNUWzEoeDSXH4XC41UUw99fQ7iuI1IKg2EvAujJZD1SsBPRLoLoUai7hSHxfyFtrns4CyKRJZh8SS135W2H7W9Blhi25jY4dM0sQ3s/oEnhGjoQmTWq4v1qNMvv98wMik51kZsJVkz32niB4coF1ZXIz8H8isk1E5orINkx53Zsdk8wlOMT3hbKigKeW6N3bpKUPCWWy9q/GRz/tdlu6y842BnjXXhI8YmLMUmq191fT9tAsNezsJjt2mMzI4/t+BU07QbNOQZHDapzJEuAM4C5MTZG7gC6quthB2VyCQZCM/DV0ugAAIABJREFU8CLGq+uLL2rIo+Q0hd/D5lcg9Rrz42IDmZnmqXjkSFu6c6klEyfC2rWQm1tNo1ajzMwkjJI+GgWppDafB8ljgiaH5TgTVT2qqjNV9S+efZ6TgrkEidg041oYYGUC5st+4ECQS/mufxK0GNJ+Z1uXmZnmqbhxlQWlXQKBd2b4xRfVNEoeDSf2hVXSx8xMGNVvI9Ele0NfmYjhchF5TERe8N2cFtAlwEQ2hhbdAx5rAmZmAkFc6jp5BDY8Cx0vhthutnS5c6dJo+LaS4JP374myWa1XoOtPFU1wsRuUlZmlOP1P/J4FoS6MgGeA54GOgLRFbY6ISIJIpIpIhs9+5ZVtCv1KYz1oc/xVBFZ7Ln+P56qjC51Ib5fUGYmbduaL3zQXIQ3PAMlx6D3HbZ16VWMrr0k+IiY/0O1S6mxPUxsUZjYTb79Fg4ehLG95kFMsnkQDBJWlcnFwBBVvURVr/PdbJDhDuBLVe0GfOl5XxmFPoWxzvM5/gjwhOf6w0BwnKzrE/F94fjWoOQpmjgR5s+HgoIAD1xyHNY/BW3PMUFfNpGZCW3aQJ8+tnXpUgcmToT9+2H58ioaiBi7yb7wmJl4H1Y6NZkHrcYENb2CVWVSgMkY7ATTgNc9r1/H1HG3hKfu+3jgndpc71IFXiN83uqADz15ssl+WmN1PLvZ8IypNNnnbtu6LC2Fzz6Ds892U6iECt7lxk8/raZRq1GQv8k4Y4Q4n3wCUzK2E3liR1CXuMC6MvkLcK/nx9tuWqvqHgDPPrmKdjEikiMii0TEqzASgSOq6g2KyAXsccFpyATJowuMobppU/j44wAOWnzMBCm2nWyrj/6iRXD4MEyZYluXLnWkTRsYOLCG+6vcbhLaS12HD8PXX8N1U+eZA8kZQZXHas6AW4HOwC0iss/3hKrWuEgnIl8AlRXOvsvi+ACdVHW3iJwBZInISkwFyIpU6dMnIjOAGQCdOgXHFzssaNYZopoHxQjfuLExxM+ebbwzA/JEv+FpU7a17x9s7Xb2bFOrxTW+hxZTpsDDD8OhQ5BQWUxqwpnm/t+bDZ0uDrR4lsnMNLPf0T3mQlG8qUcURKwqkzrl31LVs6o6JyJ7RaStqu4RkbaYwluV9bHbs98iItnAmcC7QLyIRHlmJx2A3dXI8QLwAsCgQYPCx5E80EgEtEwPWlqJqVPhww9NTECvXg4PVnzUBCm2+xEkDbG169mzTWxJfHzNbV0Cx9Sp8OCDxtHjssoqC0REmyWjvaFaZMcwe7ZRhsnMMy7NAS6GVRGrQYuvV7XZIMOHwDWe19cAH1RsICItRaSx53USMBJYo6oKzAF+XN31LrWg5UA4/C2UlQZ86HPOMfvZswMw2Lqn4ORh6He/rd3u2mXiZdwlrtBj0CCTI63a+6v1eDi6DgqqfDYNKmVlxl5y2bQ9SP7GoNtLwI+gRREZ5Ykt+Z/n/UARseMT/BmYKCIbgYme94jIIBF5ydMmDcgRkeUY5fFnVV3jOfd/wK9FZBPGhvKyDTK5JAwwHk5BCN7q2NG4CDuuTE4egXWPQYdpkDDQ1q69Bl5XmYQekZHG0eOTT6pxEfaWaN4baE8QayxbBvv2wWUTPHadVmGiTETkCswM4gTglVqBB+oqgKoeVNUJqtrNsz/kOZ6jqjd4Xi9Q1b6qmu7Zv+xz/RZVHaKqXVX1YlUtqqtMLvzw43qouvJ0zjFliqlnnedknoW1j0Jxnu22EjAG3o4dXZfgUGXKFJNtYenSKhq0TIdGLUN2qevjj409cVDHbJPdOuHMYItkeWZyFzBJVW8FvLp8FdDbEalcgk9sT4hsAoeXBWX4KVNMcsRqU1/UhRP7TOqUTpeaHw4bOXnSGEenTHFdgkOVSZMgIqKa2a9EQOtxIatMZs+GoUOhSd4XkDzW2HmCjFVl0s6nMJbXcF0CRNovkktIEBEF8elBm5kMHw5xcQ4uda15BEoLoe/9tnc9fz7k57tLXKFMYiIMG2bBbnJ8mylJEELs22dmVJeft90sQ7ep0r8poFhVJptFpKID/ghMnROX+krCAGOE18Cn8Y2ONsF+XhdhWynYZYIUU6+GuJ42d25kbtTIraoY6kydasop791bRYPW48w+xOwmn31mvhPThnqm7W1CI1ePVWXyEPCBiNwNRIvI7Zg68HW2mbiEMAkDjevssc1BGX7qVPj+e/OFt5VVDwFl0Oc+mzs2X/IPP4SxY6FZ4Mtwu/jB1Klm/9FHVTSITYOY1iG31PXhhyb4slOjLyCmDcQ57T9vDauuwe8DVwBDge2YFCbXq+onDsrmEmwSBph9kJa6fvQj43kza5aNnR7bBJtfgi43QPMUGzs2rFkDGzfCBRfY3rWLzfTrBykp1dxfImapa29WyNQ3KSw0XmgXnF+G7P3SLHGFiGHOn3ommap6rqr2UdWpqhoKNfFcnCSuN0Q0CpoRPiEBMjJsVibf3WHS7Pe5x8ZOf2DWLPPdnjbNke5dbETEKP3MTFNauVJaj4fCPXA0NFb0v/gCjh+H6eeuNLnkQsReAv7FmXQSkTtE5GnPPsU5sVxCgohok44+SDMTMF/2detMNHyd2f817HzXFL5q0taGDk/nvfeMYbetM9272MwFFxjvuyoN8W088SbfO+VW6B/vvWccU4Z08tpLwkyZiMhkjLF9KhDn2a/1HHepzyQMhEPLgjbNP9+T0rPOsxNVWHY7NGlnW233imzbZupLuEtc4cOIEZCcXM391fwMaNEN9gR/Rb+kxNhLzj0XovZ/YWw6NpWWtgOrM5O/Aj9R1dGqepWqjgauBx5zTjSXkCBhABQfMfVNgkCHDjBkiA3KZMfbcHAx9HvIBHk5wPvvm72rTMKHyEizJPnxx3DiRBWN2p5j7CYlhQGVrSLz5pnklBedXwT75oXUrASsK5MU4K0Kx/4DuKl36zvlkfDBsZuA+XHOyYEdta2oU1pkbCXx/Yw7sEPMmmXSwHTt6tgQLg5wwQUmLujLL6to0G4KlJ6AfdmBFOs0Zs2CmBg4e9BCKC0IW2WSDYytcCwDmGunMC4hSFwfYzsJot3kwgvN3vvk7zfrHjczqwGPQYQzcbb79pn0L+6sJPwYPx5iY/+/vfMOs6q6+vD7Y0bKgEgH6SqgQUVUREU+RBApYiFWRD/UqAm2JCYx1hiNGE3TfLYEiYLSREHB2FCUGCVgKIomKgwgMNKRpoCUWd8f64xexzv1tin7fZ7znHvO3ufstefcuevstddeq5jRb/OTPRrE6nREHo1Pfr7LN2AA1NnyGigLmvfOmDzxKK0yyQWekzRW0p2SxgLPA4sl3VKwpU7MQMbIqgUHHOkmogzRqZOHoi+XqWtHnq8raT0kpW9y06f7tExQJpWPWrV8zcm0aT4v8R2yakPzvq5MMjR3OG+eR6IeMgT47O/Q9CTYr35GZCmK0iqTrsAC3KzVK9ovwHOK9Iu2ijXmCiSPpj1h41zI35MxEYYMcZvxhg1lvHDhL4B8H5WkkKlTfc3CUckN8xVIE0OGeODHt4tK/d5qEHyxDLYvTqtcBUyd6vM7Z/Zd4UnrWp2ZETmKo7SLFk8pxdYn1cIGMkSznm6j/XxhxkQ4/3wf6j/zTBkuWvcPWDHJXYHrHZQy2TZs8LUK551XYdaPBcrIoEEesWDixCIqHBgl2cmAqSs/HyZN8gykDb6Mluu3OiPtcpREWdaZZEnqIemC6DhHUp3UiRaoMFSAnNhdung493HjSnlB/l6Yfz3ktIXOv0ypbJMnu3nk4otT2kwghdSt627okyfDV/GSWNRr72FLMqBM3nkHVqyIvl9502H/TlC/xGzpaae060wOwUPOv8Q3yadOAx5LkVyBikSdFlCvQ0aVCfg/07/+BcuWlaLy4ofcHHDMnyA7J6VyjR/vXlxduqS0mUCKufhi2LLFw5XEpeUgd8nd80Va5Ro/HnJy4OzB292jrHXFM3FB6UcmD+KuwY2AAsP5LOB/EhVAUiNJr0laEu0bxqlziqT3YrZdks6OysZIWh5T1jVRmQJxaNYTNryd0RhFQ4f6fvz4Eip+uRIW3eb//G2+n1KZli51BTdsWEqbCaSBU0/1BYxFjn5bDoL83bCuKB/i5LN7t4+Wzj4b6m2b4e1XQBMXlF6ZdAdGmlk+UT4TM9sCNEiCDDcBM82sIzAzOv4WZvammXU1s654kMkdwIyYKr8oKDez95IgU6AwTf8HvtrkebEzRNu2Hqtr3LhidJoZzLvW990eTvkkRoFiu+iilDYTSAPZ2XDhhR5FeMuWOBWanAT7HQCrkhksrnhefhk2b45MXJ9N9+yPTQpnA6kYlFaZbKOQ4pDUEigqE0BZOAsYG30eC5xdQv1zgZfNbEcS2g6UlqY9fZ9hU9ewYbB4McwvatlL3vPw2QvQ5c6URAWOxcyVyckne4reQOVn2DCfM5kyJU5hVk1ofbZ/x/alJzv4uHHQtCn067vP52tanu6J6yogpVUmU4HHJbUGkNQYeIDvroovD83NbA1AtG9WQv0L8VwqsYyUtEjS/ZJqFXWhpKskzZM0b0OZfUyrOft3hNrNYH1RvpPp4dxzPfFUXFPEnm0w7zrPEHnoj1Muy7x5rtjCxHvV4bjjoGPHYkypbc+HPVthbeqDpm/dCi+84KOl7K1z4KuNFdbEBaVXJrcDXwAr8RHKeuAr4J7SXCzpdUkfxtnKFKhb0oHAkcCrMadvBg4DjsPndIp03TGzUWbWzcy6NW3atCxNByQ3dWV4ZNKwoS8wmzQpzgKzedfDrjXQfVRacmKPH++K7dxzU95UIE1IPjqZNQvy8uJUaHGqm5pWTE65LFOm+Chp2DAgb5p/pw/sn/J2y0tp15nsNLOLgKb4/EmLKOBjUaHRCl9/apQHpfA2DVgXKYkCZbG+mFudDzxnZl+vnjOzNeZ8BTwRyRdIBU17ek7sHfH+y9LHJZd4qtVvhQ1fOQWWj4XDb4Umqf8K7NrlymTwYGiQjJnDQIVh2DA3YY4ZE6cwq6ZHU/hsmsfrSiFPPOGjpO7H5ft6qRb9oOYBKW0zEcqyzqQT0AnYZmbJtBFNB4ZHn4cD04qpO5RCJq4YRSR8vuXDJMoWiKVZ5LyXYVPX4MHQsiU8+mh0YsdqePcqaHRcypJeFebZZ33F9I9+lJbmAmmkQwf37PrrX4sIr9L2fDeprpkRpzA5LFrkq/F/+EPQhn/AjlXQ/pKUtZcMSlQmks6StBL4CHgHz2OyssA1NwncC/STtAQPy3Jv1G43SaNj5GgPtOG7wSXHS/oA+ABoguerD6SCBkdBdr2Mm7r22w+uugpeeQVyl+TDnMtg307o8VRazFsADz/sMcP69k1Lc4E0c801buaKmx++RR+o2QhWps7U9cgjHiH4ssuA5eMge/8Ku76kgGKViaRjgGdwN9z+QOdoPwOYLOnYRAUws01m1tfMOkb7z6Pz88zsiph6n5pZq8g9Ofb6PmZ2ZGQ2u9jM0ruiqDpRI9vnTdbMyHhO7CuvdFfO/0z9P1g7w2Nv1T80LW0vWABz5sCIEVCj1GP7QGVi8GD30Hv44TiFNfbz9Ut501KS42TrVncwGToUGtXfCauehbbnpHzxbaKU9K/wE+A+M7vCzF43s4+j/RXAb6PyQHWi9RnwRW7Gc2K3bAk3XLqAgQfeyN4WZ0KH9NmbHnnEVyRfemnamgykmexsNzG9/jp8Eu+r3vZ82PsFrHkl6W0/+aTneb/mGnxtyZ5tFd7EBSUrkx7AqCLKHovKA9WJloN9/9n0zMqx5wvuOPVC1m9rxuQVj6ctwuLmzTBhgk/Shon3qs0VV7hJ9eu5uVianwK1m8OyJ5Lappm/rBx/PBx7LLD8KajTCpqdnNR2UkFJyqSxma2KV2BmeUDj5IsUqNDUbQMNu/rCwEwy/zrq7Mvl9pfG8aeHG6fN6jZmDOzcCVdfnZ72ApmjeXN3+x4zxkcK36JGNhxyhecW+XJF0tp84w34+OPo+7Vrg4982g9LWVK3ZFKSMimpPATcro60OhM2zoZdGzPTfu4oWDYGHXEbxw3uzfz5MHt26pvdu9dt6D16QNcQAa5acM01Pocxdmycwg5X+Yg4tyjjTdn585+hSRNPucCKSWD74KCKb+KCkpVFrdhMioU3oGY6hAxUMFqdAZafmTSmn/0d/j3C80sc8SuGD/fgfL/+deqbHjfOAzveeGPq2wpUDHr08O23v/W1Rd+ibls3+y4dDft2J9zW/Pm+4v2666B2rXxY8ig0PBoaHJHwvdNBScpkDt9kUoy3zUmpdIGKSaNjoE7L9M+bbJwLb5/v/2A9J0ONbOrWhZtu8onSt95KXdN79sBdd7kd+8yK7aEZSCKSP/e8PBg9Ok6FjiNg13pYNTXhtu64wyM8/OQn+EvTto/gez9P+L7pQpZhF89M0a1bN5s3b16mxai8vPtD+HQCnLPR88Snmm2L4bUo73W/2VCn+ddFO3fCwQfDYYfBm2+mpvnRo90d+cUXPStfoPpgBr17w5IlPjKtE5sS0PLhhY6Q0xpOLbwErvTMnQsnnAD33AM33wy81tMjTZyRW+ECO0qab2bdCp8PXvKB8tHqTHeNXDcr9W19uRLeOBUQ9H7lW4oE/J/7lls8nlIqlMnu3fCb37iHzcCByb9/oGJTMDpZs8ZXxX+7sIa7pa9/C7aUP/jGHXf4XMm11wIb3vHtsJ9VOEVSHEGZBMpH8z6QVSf1pq6da2FmX/e17zMD6neMW+3KK6FVK/jVr5K/nvLxx2HlSv9BCTneqycnnwx9+vjcyXc8uw6+DGrUgo//WK57z54Nr77qc3H77w/893dQqzEccnnCcqeToEwC5SO7DrQaDCsmwp7tqWnjq03wRj+PBNz7ZXdJLoLateHWWz2e0bTioruVka1bfVRy0knQr1/y7huofNx5J6xfD3/4Q6GC2k2g0zWw/EnY8kGZ7pmfD7/8pTuRXH01sPW//oLW8VrIrps02dNBUCaB8nPYz2D3Zsh9LPn33p4LM06E7Uug1/PQ9MQSL/nBDzwP+4gR8PnnyRHjhhtg7Vq4//4wKqnu9OzpuUVGjoQPCuuMw2+F7PqwsMgMGHF56CF/AbrvPqhbF/hwpI/4O12bNLnTRVAmgfLT5HhfCfzxH5ObeW792zDjBNj9OfR53XNIlIKaNX2B2caNkUdMgrzyipu4brzRkyYFAg8+6B5Xl17qHn5fU6sRHHErrHkZ1pYuR3xurnsinn46DB8OrH4VVkyAw27w0U4lIyiTQGJ0vgl2roZPi0pNVwbMIHc0vNHXo7KeNgea9SzTLY4+2ifjn3rKffbLy9atPg/TuXN61rAEKgdNmnh4lQUL4Pe/L1TY6Vqo2w4W3uheXsWQn+8j6Zo1fVJfe7d7GoX6h8ERt6WuAykkKJNAYrTo5+s+Pvod5O8r/312b/Y1JO9e6ZGJT/sX7N+hXLe69VY48kgP1LdpU9mvN3Pz1urVnqCoVho8nwOVh+9/Hy64wF8yFi2KKciqDV1GwuYFHlOrGB56yNdFPfCAO47w/q2es+T40X6fSkhQJoHEkHx0su0Tzz5XHla/Ci8dBXnPQ9f73GurVvnDvhWYuzZtgv79YcuWsl1/551u3vrlL6F7yNsZiMODD0LjxjBggK8/+Zr2Q6FJD4/SsDH+mu7Jk+GnP/Uw98OH427Aix/ykU3Tk9IifyoIyiSQOG3OgXqHwHs3ly1e17Yl8I8zYdYAn3Q8bTZ0vtF99xPkmGNg6lR/c+zf381WpeHuu12ZXHaZfw4E4tG0qUdd2LMHTjnFFzMC/t3t9ZxHiPjHGb7YNoYpU+CiizxEy8SJoM0L4Z/nemiWo+5Jf0eSSMaViaTzJP1HUr6k76yqjKk3QNInknIl3RRz/iBJcyUtkfS0pBAvLN3UyILj/+bRU2cN8DUhAMvHw/PtYUIN3y8f77bktW94dsSXDvdFj11/B4MWQePkznKffjo884zbtwcOhA2Fkk2PHw/t23uCq3btPELs7bd7jvnHHguJrwLFc/jhMHOmR2A45RSP9gtA7WZwSpTn5M0BsHMtZvD00+4N1r07vPQS1Nv+Grzey5Nt9X4J9quXsb4kBTPL6AZ8DzgUmAV0K6JOFrAUOBgPLvk+0DkqmwxcGH3+CzCiNO0ee+yxFkgyeS+YTcg2e62XWe7jZpNyzMbzzTYh22xyQ/88ub7Z3B+Z7ViTcrGmTDHLzjarW9fsllvMNm0yGzfOLCfHzGdIvtlOPNFs796UixSoQixYYNaokVlWltnll5stXx4VbJhjNqmO7ZlQ3yb+4mfWpvEKO+nEr2z7infNPrjb/x9e7GL2ZV4mxS8zwDyL85taYWJzSZoF/NzMvhMwS9KJwK/NrH90fHNUdC+wAWhhZnsL1yuOEJsrRXw6CWZf5MN9izMhn1UHTnjCw7Fk1/lueYr45BOfMH36affn37XLQ8oXpm1bWJG89BSBasK6dXDvve7plZ8PHTr4dOJBDRcxrOtvOe/4Z6hRA5S1H8qPwg+3OM0DltY8ILPCl5HKHpurFRCbpCsvOtcY2GJmewudj4ukqyTNkzRvQ2GbRyA5tL8Qej4dX5EA7NsF7S5IqyIBOPRQt1EvWuS5teMpEoBVcVPBBQLF07y5L2xduhSuv95NYJ07Q50Du7Cuw0T2DlpGjcNvRB1HuAI5exX0ebXSKZLiSEsUMUmvAy3iFN1qZqVxAYq39tiKOR8XMxtFlIa4W7duFWNIVhVpex7ktIMdcV7xc9qmX54YjjgCRo2CGTPij0DaZla8QCWnVas44VYAaAuNK/cEe0mkZWRiZqea2RFxttL6kuYBbWKOWwOrgY1AA0nZhc4HMs1RIyEr59vnsnL8fAVg5EjIKSReTo6fDwQCZaeymLn+DXSMPLdqAhcC06PJoDeBc6N6w4EkhvkLlJuDhkH3UT5CQb7vPsrPVwCGDfMRSrt2bttu186Ph1UM8QKBSkfGJ+AlDQEeBJoCW4D3zKy/pJbAaDMbFNUbBDyAe3Y9bmYjo/MHA5OARsBC4GIzKzFQVJiADwQCgbJT1AR8xpVJpgjKJBAIBMpOZffmCgQCgUAFJiiTQCAQCCRMUCaBQCAQSJigTAKBQCCQMNV2Al7SBqC8gTOa4GtcqhOhz9WD0OeqT6L9bWdmTQufrLbKJBEkzYvnzVCVCX2uHoQ+V31S1d9g5goEAoFAwgRlEggEAoGECcqkfIzKtAAZIPS5ehD6XPVJSX/DnEkgEAgEEiaMTAKBQCCQMEGZBAKBQCBhgjIpI5IGSPpEUq6kmzItT7KR1EbSm5I+kvQfST+OzjeS9JqkJdG+YaZlTTaSsiQtlPT36PggSXOjPj8dpT+oMkhqIOlZSR9Hz/vEqv6cJf00+l5/KGmipNpV7TlLelzSekkfxpyL+1zl/F/0e7ZI0jHlbTcokzIgKQt4GBgIdAaGSuqcWamSzl7gZ2b2PeAE4JqojzcBM82sIzAzOq5q/Bj4KOb4PuD+qM+bgR9kRKrU8WfgFTM7DDgK73uVfc6SWgHXA93M7Ag8ncWFVL3nPAYYUOhcUc91INAx2q4CHi1vo0GZlI3uQK6ZLTOz3XgelbMyLFNSMbM1ZrYg+rwd/4FphfdzbFRtLHB2ZiRMDZJaA6cDo6NjAX2AZ6MqVarPkuoDvYC/AZjZbjPbQhV/zniq8jpRdtYcYA1V7Dmb2VvA54VOF/VczwKeNGcOnrn2wPK0G5RJ2WgFrIo5zovOVUkktQeOBuYCzc1sDbjCAZplTrKU8ABwI5AfHTcGtpjZ3ui4qj3rg4ENwBORaW+0pLpU4edsZp8BfwBW4kpkKzCfqv2cCyjquSbtNy0ok7KhOOeqpG+1pHrAFOAnZrYt0/KkEkmDgfVmNj/2dJyqVelZZwPHAI+a2dHAl1Qhk1Y8onmCs4CDgJZAXdzMU5iq9JxLImnf86BMykYe0CbmuDWwOkOypAxJ++GKZLyZTY1OrysY/kb79ZmSLwWcBJwp6VPcdNkHH6k0iMwhUPWedR6QZ2Zzo+NnceVSlZ/zqcByM9tgZnuAqUAPqvZzLqCo55q037SgTMrGv4GOkfdHTXzybnqGZUoq0VzB34CPzOxPMUXTgeHR5+HAtHTLlirM7GYza21m7fFn+oaZDQPeBM6NqlW1Pq8FVkk6NDrVF/gvVfg54+atEyTlRN/zgj5X2eccQ1HPdTrwv5FX1wnA1gJzWFkJK+DLiKRB+FtrFvC4mY3MsEhJRVJP4J/AB3wzf3ALPm8yGWiL/1OeZ2aFJ/kqPZJ6Az83s8GSDsZHKo2AhcDFZvZVJuVLJpK64g4HNYFlwGX4C2aVfc6S7gQuwL0WFwJX4HMEVeY5S5oI9MZDza8D7gCeJ85zjZTqQ7j31w7gMjObV652gzIJBAKBQKIEM1cgEAgEEiYok0AgEAgkTFAmgUAgEEiYoEwCgUAgkDBBmQQCgUAgYYIyCQSKQFJbSV9IaplpWeIRyXZionWSjaT+kv6ZwPWHR5GMayVTrkBqCcokkHEkzZL0laTtkrZKWibpKUnHZlIuM1tpZvXMrEKuiI5k+xf4+hhJe4urkw6idQv342sbCs4NkbRK0jpJIwrV/4uiNAcFmNl/8PUe16ZD5kByCMokUFH4jZntb2YHAKcAK4A5koZkWK5A2TgNXwT5JnydtuEx4Aw86vbdkppHZb2BI4EH49znceA6SeE3qpIQHlSgwmFmK8zsNuBJ4MHobRdJjSU9KWmNpLWSxkpqVHCdpE8l3SZP7vWFpA8kdZE0NEr+szWKjpsdc80T0Vvzdkn/lXRRTFl7SRaFp0fSryXNlHSPPPnQ+mhFdVxirr9C0mJJWyRNk9Qspk5Jfbpe0vJIvs8k3RNTZpJ6Rma4l4GsqN9fSBoeWyfmmnMkvR/9Ld6PVdaSLo3+TtdLypO0WdJfI4VQWs4GXrdvVkM3BvaY2XtmtgJYArSdxr3sAAAEBUlEQVSTVAdfeX2FmeXHuc9bQAugaxnaDmSQoEwCFZlJeKiLgvhR44GGeGKy7+HhIp4qdM1w4Oqo3vvAc/hI5yj8LfhM4PyY+m/jP1gNgLuAMSo+4VkvPBxFS/xt+xZJJ5XQj/+NrmuLh6gZF1NWZJ8kdQLuBQab2f7A4cSJBReZ4QYC+yKzVj0zG1u4XjR3Mh6PDtwYD5MzUdLxMdXaAc2BQ4DjgPPweGWl5Rg83lUBG4B9ko6VdBDQHsgF7gYmmNlH370FROFMlkT3C1QCskuuEghkjLxo3zh6++4PdDKzzQCSbgA+lnRgTHC6UQU/UJImAMOAE8zsS+BLSbPwH8kJAGb2t5j2Jkn6OR7XKPYHMZbFZvaX6PNcSe8B3YB3iunHnVFgRST9AlgSM6lfZJ/w+FECDpe0IkpeNaeYdkriMmCKmb0cHb8o6Tngcjz2GsBO4Fdmtg/IlTQz6t/4UrbREPg6ZYGZmaSheFbHmsAIoANwMtBL0gO4wlgKXB8lZCtgGx4vK1AJCCOTQEWmdbTfxDdhspfHlC+N9rEhtGMjnu7A39Y3FDq3P4CkGpLukvRJZPbZgo9gmhYjU+GIql8W3K8YPo3zuTUl9MnMluHK8EpgtaS3JZ1WQlvF0QYP6BjLUr7991sfKZICStO/WDYD9WNPmNk/zaynmXUHXgRG4X26BGhoZr3wv2vhfCr1+W7GwEAFJSiTQEXmAuAz4BO+yQbXPqb84GgfmymuLAzFo8aeg/+oNcBNY/ESBiVC+zif8yhFn8xsqpn1w81fk4FpknLitBFv3qEwq/DEULEcTPn/fvFYiJvsiuI24EUzW4gr7tnR+bfwrJ4AyFM8dIzuF6gEBGUSqHBIahNNbF8K/DjKT70amAH8UVIDeda8PwIvlzf/Av7muxe369eQdDn+A5dsbpfUXJ53/T5gppmtLqlPkg6VNCBSHnvwNLNGfMWxFp+AL6wsYhkDnCNfB5IlaSDwfeCJpPXUQ533jVcgqQue6fCu6NRSoF/kEDEQn0spoBcePj0ok0pCUCaBisLtkcfSNvwttQPQw8ymxNS5GNgOfBxtW/DJ7fIyFp8ryMVHQJ3xXC7JZlx031X4vMHFMWXF9akmvl5jTXT+euAcM9tVuAEzWww8ArwbeY1dEqfObNxB4Q+4Oep3eO6OROZhCvMqsFfu9vs1kUfYaGBETK6Qv+J93IQr8XtiLrkceLAIT69ABSTkMwkEUoSk9vh8SBszyyu+dtVB0gDglmgupDzXd8a98LpU5iRV1Y2gTAKBFFFdlUmgehLMXIFAIBBImDAyCQQCgUDChJFJIBAIBBImKJNAIBAIJExQJoFAIBBImKBMAoFAIJAwQZkEAoFAIGH+HzgMzqJoYkTyAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "import numpy as np\n", "from matplotlib import pyplot as plt\n", "import mwarp1d\n", "\n", "\n", "#define a plotting function (for re-use below)\n", "def plot_warp_results(y, yw, x0, x1):\n", " plt.figure()\n", " ax = plt.axes()\n", " c0,c1 = 'blue', 'orange'\n", " ax.plot(y, color=c0, label='Original')\n", " ax.plot(yw, color=c1, label='Warped')\n", " [ax.plot(xx, y[xx], 'o', color=c0) for xx in x0]\n", " [ax.plot(xx, yw[xx], 'o', color=c1) for xx in x1]\n", " ax.legend()\n", " ax.set_xlabel('Domain position (%)', size=13)\n", " ax.set_ylabel('Dependent variable value', size=13)\n", " plt.show()\n", "\n", "\n", "#define landmarks:\n", "Q = 101 #domain size\n", "x0 = [38, 63] #initial landmark location(s)\n", "x1 = [25, 68] #final landmark location(s)\n", "\n", "\n", "#apply warp:\n", "y = np.sin( np.linspace(0, 4*np.pi, Q) ) #an arbitary 1D observation\n", "yw = mwarp1d.warp_landmark(y, x0, x1) #warped 1D observation\n", "\n", "\n", "#plot:\n", "plot_warp_results(y, yw, x0, x1)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The input arguments to warp_landmark include:\n", "\n", "* **y** — the 1D data to be warped\n", "* **x0** — initial landmark positions\n", "* **x1** — final landmark positions\n", "\n", "Let’s consider what happens when only a single landmark is specified:" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZMAAAEJCAYAAABR4cpEAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOydd3xW5dn4v1dCIKwkhB3CiAxlK0tZJkwRUBRH3VprtdXWDvv+aqvVqrWvra9V39dtXa2zDhABR9gioESRLXtkMMIKJJBAkuv3x/08EELGM854nuR8P5/zOc9zxn1fT3LOuc59X0tUFQ8PDw8Pj3CIcVsADw8PD4/ox1MmHh4eHh5h4ykTDw8PD4+w8ZSJh4eHh0fYeMrEw8PDwyNsGrgtgFu0atVKu3Tp4rYYHh4eHlHFt99+u09VW1feXm+VSZcuXcjKynJbDA8PD4+oQkR2VLXdm+by8PDw8AgbT5l4eHh4eISNp0w8PDw8PMKm3tpMPDw86h8nTpwgJyeH4uJit0WJeOLj40lNTSUuLi6g4z1l4uHhUW/IycmhefPmdOnSBRFxW5yIRVXZv38/OTk5pKWlBXROxExzicirIrJXRNZUs19E5H9FZLOIrBKRARX23Swim3zLzc5J7eHhEU0UFxfTsmVLT5HUgojQsmXLoEZwEaNMgNeBCTXsvxjo7ltuB54HEJFk4EHgfGAI8KCItLBVUg8Pj6jFUySBEezfKWKmuVR1kYh0qeGQKcC/1OTMXyYiSSLSHsgAMlX1AICIZGKU0jv2SmwTR3MgbzakTIQmqbUevmsX/OtfUFRkvjdoAFddBT172iynR2RQ8APkzYSyYtAyiGkI7cZB8kCw4KFZWAivvw5795rvIjB2LIwcGXbTHnWMiFEmAdAByK7wPce3rbrtZyAit2NGNXTq1MkeKcMhexp8/RM4fhCkAXS+BnreAy3OPePQ/Hz4+9/h2Wfh2LFTzw1VeOghuP56ePBB6NrV4d/gYT9aDrkzYeP/we45Z+5f+Udo0hE6Xgm9fg+N2wbdxdGj8Pzz8NhjsG/f6dfXww/DuHHwyCNw/vlh/pZ6Sk5ODnfddRfr1q2jvLycyZMn8/jjj9OwYcPTjsvLy+Puu+/mgw8+qLG9iRMn8vbbb5OUlBS0LH/+859p1qwZv/vd74I+tyKRNM1VG1W9ZmkN28/cqPqSqg5S1UGtW5+RDcA9yorhm5/Dl1OhWVcYMx963AU50+CzgbB30WmHZ2bCWWfBP/5hRiGbN0N5uVn27oXf/hY++ADOPhteeMGl3+RhDyX7YcEkWDQFDv8A/R+Fy/PgmuNwbRlcsR8ueAOSB8CmZ2BWL9j2b6MFAmTDBjjnHPjd72DAAFi27NT1VVQE//M/sGIFXHAB3H13UE17YIzbU6dO5bLLLmPTpk1s3LiRwsJC7rvvvtOOKy0tJSUlpVZFAjB79uyQFImlqGrELEAXYE01+14Erq3wfQPQHrgWeLG646pbBg4cqBHDij+ovoXqd/9PtbTk1Pbi/arTOqnO6qdadkJVVefOVY2PV+3XT3X9+uqb3LVLddIkVVB98UWb5fdwhvxlqtM6qr7TUHXDMyeviWo5tF7182Hm2pp3seqx/Fq72LhRtX171TZtVBcsqP64I0dUf/lLc3394heq5eVB/haXWLdundsi6Jw5c3TkyJGnbSsoKNDk5GR99tln9corr9TJkyfrqFGjdNu2bdq7d29VVS0qKtKrrrpK+/btq1dffbUOGTJEly9frqqqnTt31vz8fN22bZuec845etttt2mvXr103LhxevToUVVVfemll3TQoEHar18/nTp1qhYVFamq6oMPPqiPP/54lbJW9fcCsrSKZ2o0TXPNAH4hIu9ijO0FqrpLRD4H/lrB6D4e+INbQgbNsT2w4WnofC2c97fT9zVKhgFPwOKrYMvLLMz7OZMnQ7duMGcO1DS4atcOPvwQpk6FO+4wtpRbb7X3p3jYyI7/wNIboHEHGL/E2ERqI/EcGLsINj0LK/4fzB0FY+ZCfJsqD9+yBUaNghMnYP586NOn+qabNYOnn4aGDeGJJyA2Fp580hIzjWP8+tfw/ffWtnnuufDUUzUfs3btWgYOPP3/l5CQQKdOnSgtLWXp0qWsWrWK5ORktm/ffvKY5557jhYtWrBq1SrWrFnDueeeOf0NsGnTJt555x1efvllrr76aj788ENuuOEGpk6dyk9/+lMA7r//fl555RV++ctfhvV7KxIxykRE3sEY01uJSA7GQysOQFVfAGYDE4HNwFHgx759B0TkEWC5r6mH1WeMjwrWPQblJdD3z1Xv73gFtB1F6Xf3c+PdV5OW1pK5c2tWJH4aNTIK5bLL4LbboEULuPxyS6X3cIK8z2DJ9dBqKFw43bxkBEpMLJx9NyT2hoWXwJwMo1Aatz/tsP37YfRoY3+rTZH4EYHHH4fSUqNYkpPhgQeC+2n1EVWt0lPKv33cuHEkJ5/5P168eDG/+tWvAOjTpw/9+vWrsv20tLSTimbgwIEnFdKaNWu4//77OXToEIWFhVx00UUW/SJDxCgTVb22lv0K3FXNvleBV+2Qy1aO5sCm5yHtZkjoUfUxIpzo/zTy2Xk8dMUDXPzAs7Sp+sWySuLjYdo0uPBCuP12GD6coM73cJn8pfDlFZDUB9I/gYaJobXTbgxkfAoLJ8GcdDNiadzu5O6774a8PFiyBKp5RlWJiBmR5Ocbw/ykSTAwgEFTJFDbCMIuevfuzYcffnjatsOHD5OdnU1sbCxNmzat8jwN0DjVqFGjk59jY2M5duwYALfccgvTp0+nf//+vP766yxYsCC0H1AN0WSAr3us+QtQDn3+VONhjz3fl+cy7+SWES/QrsmGoLtp3Ni4dx4+DL/4RWiierhAwXrz8G+cAhmfha5I/LRNh1Gfw9FcWHw1lJ8AYPp0ePtt+NOfYPDg4JsVgWeegbZt4ZZb4Pjx8MSs64wZM4ajR4/yr3/9C4CysjLuuecebrnlFpo0aVLteSNGjOA///kPAOvWrWP16tVB9XvkyBHat2/PiRMneOutt0L/AdXgKRO3KNwKW16BrrdDsy7VHrZqlXHBXB9znxkab309pO5694Y//xnef98sHhFO6TFjK4tpCKO/CMm9t0paD4chL0H+l7Dyjxw4AD/7mZnr/0MYlsYWLeCll2DNGvjLX6wRta4iIkybNo3333+f7t2706NHD+Lj4/nrX/9a43l33nkn+fn59OvXj7/97W/069ePxMTAXzAeeeQRzj//fMaNG8c555wT7s84k6qs8vVhcd2b67vfqb4Tp1qUW+0hx4+rnnee8azZt09V50803jzlZSF1eeKE6sCBqq1bq+7dG6LcHs6w/BfGCyvvc3va/+Yu1bfQJ3/7oTZooLpihTXN3nSTamys6rffWtOe1USCN1eolJaW6rFjx1RVdfPmzdq5c2ctKSmp5azwCMabyxuZuIEq7PwA2o6FJinVHvbii8af//nnoWVLoMsNcDT7jLiTQGnQwEx3HToE998fmugeDpA7EzY+A2f/BtqPt6ePAU9wuOH53NrrFv7nwU1U4xgUNE89ZWxyd93lxZ9YzdGjRxkxYgT9+/fn8ssv5/nnnz8jyNFNPGXiBge/g6Lt0OnKag85dgz++ldjOD/pgZU6BRo0N0FoIdKnjzHEv/oqVPA69IgUju2GZT+GpP5w7n/b109sI257/X2UGH5x/s8se/K3aGGmU5ctg88+s6RJDx/NmzcnKyuLlStXsmrVKi6++GK3RToNT5m4wc4PQGKNcqiGF180ebcefriC736DJtDpCsj+wMyph8gf/mDiAry57Qgk65dQWgjD34bYRrUfHyLz58P7szqyouyvxObPgx3WpbK75Rbo0sW4CXujk/qDp0yc5uQU1xho1LLKQ4qK4L//2/j9p6dX2tnlBjhxGHI/CVmEDh2M0fX1102gmkeEsHuueVHo9UdI7GVbN6rmQd+hA1xw4x2QPBi++y0cP2RJ+w0bGs+wrCyYOdOSJj2iAE+ZOM2hVVC4ucYpruefNzm2Hnqoip1tMkwUdBhTXQD33mtu+kceCasZD6soPwHf3g1N06DXf9na1Zw5sHgx3HcfxDeOhSEvQEk+rLTOkHbjjSbJqDc6qT94ysRpdn4AEgOpl1W5u7AQ/vY3GD8eRoyo4oCYWOhyHez6DIrzQxajXTv4+c/h3/+GjRtDbsbDKjY+BwXrYOCTEBtvWzf+UUmnThXS6yQPgO53wabnYH+WJf3ExZl+vv/exLF41H08ZeIkqpD9vhldxFedD+Wf/zQpv6sclfjpfC1oKeTNCkuc3//epFz5+9/DasYjXIr3wuoHof1F0OFSW7tasMAYx++7z/zvT9LvEWjUyqSvt4jrroMePYwjiccpfvOb3/BUhfD7iy66iNtuu+3k93vuuYd//OMfjsnTrFkzS9rxlImTFKyFwxuqneIqL4fnnoOhQ01672pp0d/c+HvmhyVOmzZwww0m+vngwbCa8giHVQ9CaREMfNr2TInPPmvczG+6qdKOhomm9snuTNj7pSV9NWhg0rRkZcE331jSZJ1g2LBhLFmyBIDy8nL27dvH2rVrT+5fsmQJw4cPr7UdVaW8vNw2OYMlYGUiImeJyB9F5Fnf97NFpLd9otVBcqYDAqlVZ1ucOxc2bTI++jUiMdB2FOyZF/aE9J13Gjfk118PqxmPUCnaAVtfgW63Q8LZtnaVm2umnH7yE5Oz7Qy6/xzi28GqP1lm6LjxRpNh+LnnLGmuTjB8+PCTymTt2rX06dOH5s2bc/DgQUpKSli/fj09e/ZkzJgxDBgwgL59+/Lxxx8DsH37dnr27Mmdd97JgAEDyM7OplmzZtxzzz0MGDCAMWPGkJ9vpr+3bNnChAkTGDhwICNHjuSHH34AYNu2bQwdOpTBgwfzpz/VnMopGAJK9Cgi44CPgPmYzL53Aa2A+zG12T0CYc98SOp3WoK9ijz7rMkGfGX1tvlTtB0FO9+Hwi3QvFvIIp17LgwbZm72X/0KYryxqrOseRQQ6G1/1YSXXjKj35/9rJoDGjQxcnz7K/Oi0m5M2H0mJBiF8uqrpqhWq1ZhN2kd3/4aDlqcg77FuTCw5gySKSkpNGjQgJ07d7JkyRKGDh1Kbm4uS5cuJTExkX79+tGkSROmTZtGQkIC+/bt44ILLuDSS80U6IYNG3jttdd4zqehi4qKGDBgAE888QQPP/wwDz30EM888wy33347L7zwAt27d+frr7/mzjvvZN68efzqV7/i5z//OTfddBPPPvusZT890EfHY8BVqnopUObb9h0wwDJJ6jplJbBvCbTNqHL3zp3wyScmVXyjQMIL2o426z3zwhbtrrtMtcbMzLCb8giGwm2w9TXo+lNokmprV8ePG2UycSKkpdVwYLfbjSyrrHPDuusuKCkxCsXD4B+d+JXJ0KFDT34fNmwYqsof//hH+vXrx9ixY8nNzWXPnj0AdO7cmQsqzIPHxMTwox/9CIAbbriBxYsXU1hYyJIlS7jqqqs499xzueOOO9i1axcAX331Fddea5K033jjjZb9pkBT0HdVVX88qwKo6jERibNMkrrO/uWmPG+byoEjhhdfNOs77giwveY9TDbZ3fPMAyAMrrgCfvMbMzqxuMSBR02sfdQErzowKpk+HXbvNtOaNRIbD73vg+U/h12fQ8qEsPvu3dvES73wAtxzjwmYjQhqGUHYid9usnr1avr06UPHjh154oknSEhI4NZbb+Wtt94iPz+fb7/9lri4OLp06UJxcTFAtSnq/YgI5eXlJCUl8X011b+qqqcSLoGOTLJF5LRyOSLSH9huuUR1lb0LzLrNhWfsKimBl1+GyZOhc+cA2xMxoxML7CaNGpkR0cyZsGNHWE15BErhVpMButvt0KSD7d09+6wZkUwIRDecdasZnaz7W+3HBsidd8K2bV6KFT/Dhw9n5syZJCcnExsbS3JyMocOHWLp0qUMHTqUgoIC2rRpQ1xcHPPnz2dHDTdmeXn5yTrxb7/9NiNGjCAhIYG0tDTe96UIV1VWrlx5su93330XwNJU9IEqk/8FPhKRG4BYEbkCeBN40ipBRGSCiGwQkc0icm8V+58Uke99y0YROVRhX1mFfTOskslS9i409pIqot4/+sgUF6rV8F6ZtqNNsFnB2tqPrQX/iOill8JuyiMQ1j4G0gB6nXGpW866dbBokYkrCsgmFtsQzv6VeQE68J0lMlx+ObRvbwJyPaBv374nbSEVtyUmJtKqVSuuv/56srKyGDRoEG+99VaNKeObNm16shTwvHnzeMBX7vKtt97ilVdeoX///vTu3fukEf/pp5/m2WefZfDgwRQUFFj3o6pKJVzVAtwOrAYKgTXAbYGeG0DbscAW4CygIbAS6FXD8b8EXq3wvTDYPh1NQV9aovpuY9Xlv6xy98UXq3bqpFoWbGb5wu0mTfkPT4cvYzhyeATHsXzVd+NVv77dke7uvdekhd+zJ4iTSg6pvtdcdfF17sphMdGcgr46mjZtalvbtqSgV9WXVLWvqjZT1T6q+k/LNBoMATar6lZVPQ68C1SfBRGuBazLTGc3B7Kg7FiVxve9e+GLL0yAV9CeVE07Q7OzLDHCA1x/vXEE+OorS5rzqI4tLxv7WY+7be+qvNzEEY0fH2S55oaJ0O2nsPM9KNppiSzXXw9lZeArFuhRxwjo8SUiKdUtFsnRAciu8D3Ht60qWToDaUDFJ2i8iGSJyDIRqTpPiTn3dt9xWX5fbEfw20tan2kvee89c4PdcEOIbbcdDXsWQHlZrYfWxmWXQdOm8OabYTflUR3lJ2Djs9BuLCTZH6a1eLF5QQjp+jr7V2a94X8tkaVPH+jf37u+rKawsNBtEYDAbSY5mId9VYsVVOVaUJ1V+RrgA1Wt+PTspKqDgOuAp0Ska1Un+kZXg1R1UOvWVaczsYU9CyGxD8Sf6WT/5pvmBusd6nOl7Wg4UQAHV4QnI0aRXHaZeXMsKQm7OY+qyP4IjuWeelDbzJtvmv/rlJrG+dXRtBN0ugo2vwTHrZlbv/56+PprE5zrFuplngyIYP9OgSqTNIw9w7+MBD4Hbgmqt+rJATpW+J4K5FVz7DVUmuJS1TzfeiuwADjPIrnCp/wE7PuqSpfgTZtMmomQRyVwyjss35q5qRtuMJUYP/3UkuY8KrPhaWjWFVIm2t5VSQm8/74xftfiTVo959wDpUdgyyuWyHTttcYR8e23LWkuaOLj49m/f7+nUGpBVdm/fz/xVaZKqJqA4kxUtbJf2g4RuRkz1RReLnTDcqC7iKQBuRiFcV3lg0TkbKAFsLTCthbAUVUtEZFWwHAgclIXHvjW5F2qwl7y1lvmxvLFD4VGkw4mJf1+a5IfjR1r5tbffNOMUjwsZP9y2LcUBjxlUuLYzOzZ5sXg+uvDaKTlIGg9HDa/AOf8Omy5U1MhI8NcXw88YHsqsir6TyUnJwdHp7mjlPj4eFJTAw+mDTRosSqOAoFGRdSIqpaKyC8wo51YjKfWWhF5GOM54Hf3vRZ4V09/regJvCgi5ZiR1mOqus4KuSxh70KzrhRfomqUyahRpkhRWLQ6H/Z/HWYjhgYN4JprTBDloUOQlGRJsx4AG/7PlF3u+mNHunvrLfNiMHZsmA11+xksvdGkA7IgxcoNN5j8YMuXw5AhYTcXFHFxcaTVmALAI1QCNcBfV2n5KfAJYM0TDFDV2araQ1W7quqjvm0PVFAkqOqfVfXeSuct8XmZ9fetrRmPW8W+r6FZN4g/3ZXmm29MCpOw3hr9tBxicnQV77OgMSNTSQl8+KElzXkAHD9oyg+k3QBxCbZ3d+iQSc9zzTXmBSEsOl0JDZNh0wuWyHbFFSZQ1jPE1y0CHbM+Wmm5G9iBdTaTusuBLGg5+IzN771nKh1ecYUFfbQ839fXcgsag8GDoVs3I6OHRWx707gDd/2pI919/LHJx3XdGZPFIRAbD2f92GS9PrYr7OYSE022h/ffN67LHnWDgJSJqqZVWvqq6o9VNdduAaOaY3vgaDYkDzptsypMm2amHxITLegneSAgZhRkASIwdSrMn+/VObEEVRNbkjwQkp3xDZk2zdgnLJtG6naHKchmkSF+6lSTK2zZMkua84gAvITjdnLgW7NueboyWbkStm83XjaWENccEntbZoQHc7OXlsKs8Io5eoD5vxxa7diopKgIPv/cXF+WGbgTupvYmM0vWRLTNGmSKe07bZoFsnlEBNUqExHJFJEvalucFDbqOLAcEGhx+tvoRx+ZaPdLrazQ2up8OPCNZWnDBw+GlBQjq0eYbHkZYptAl3Dc9gLns8+guNjClxU/3X5mRtp5s8NuKjERxowx15fnpVs3qGlkshj4KoDFozr2Z0HCOWbkUIFp02DEiCDTW9RGyyFQst8Y4i0gJsa4Bn/2GRw9akmT9ZMTR2DHu9D5GkcM72Cur5YtYeRIixtOvRTi25rKkBZw+eWwdSusXm1Jcx4uU62fh6o+5KQgdQ5VY3xvP/60zZs2wZo18KRl+ZZ9+I3w+78Jq/JiRaZONTVOvvjCizkJmR3vmDijbs5McR0/bkoJTJ1qgRdXZWLiIO1G+OEpKM6H+PCySEyZYqo+TpsG/fpZJKOHawRlMxGRxiLSUUQ6+Re7BIt6juVB8e4zjO/+OWLLH86Jvc1UikVGeIALL4QWLbyprrDY8qpJpeNX9jYzfz4UFNgwxeUn7WZjiN8efgh727amZLR3fdUNAo0zOUtElmDSz28HtlVYPKriQJZZV6FMBgyALl0s7i+mgfEWstAIHxcHl1xi4hVOnLCs2frD4Y0mmPSsmx0L9Z42zaROGTfOpg6S+pjrbNvrljQ3dSqsWmWmuzyim0BHJs9gkjr2B44A/YDpwE9skiv62Z9lSrK26H9yU16ecYW07a2x1fkm4WPZccuanDrVBMAtXGhZk/WH7W+a9COdrQj2qJ3ychNfMnEiBJFSKXjSboGD38PBlWE35b8XPK+u6CdQZXI+phjWGgBVXQvcAfyXXYJFPQeyzNRTgyYnN/kKndmnTFoOgfISOBT+Te5n/Hho0sS72YNGy2Hbv6HtGGhiVaWGmlm2zMRu2HZ9+elyrbGfbH0j7KbS0kzWbO/6in4CVSblwDHf50IRSQIOAJ7NpCr8xvdKU1yzZsFZZ0GvXjb129IXoeafYrOAxo2NQpk503PhDIr8r6BoO6Td5FiXM2dCbCxcfLHNHTVqCR0uMSOv8vDnP6dMgaVLYf9+C2TzcI1AlclaTDZeMPm4nsTUhfdsJlVRtANK9p2WRuXYMZg3zwRr2TZ93qQTxCVZMv1QkYkTTYGldZGTPjPy2fZvaNAUOto9TDjF7NnG5dyR5Jxpt0BJPuSFX6tg4kQzRff55+GL5eEegSqTuzEjETBTWx2AQZipLo/KVGF8X7DAKJSJdpaxEIEW/eDQKkub9b/pzg4/Vq1+UFYMO/8DqVONQnGA3FyTWcHW66siKROgUWujNMNk0CBo1cq7vqKdQHNzrVLV1b7PW1V1vKpeoKpLazu3XnIgy8wpJ/U9uWn2bDNllH5mjSxrSepvlIlal0EvNdXEAXg3e4DkfmKqX6bd6FiX/mJmjimTmDjodDXkzYQTh8NqKjYWJkwwAbJl4Wdq8XCJQF2D5/hSz9vpI1J3OLjSF/fRCDC2htmzTfqIxo1t7rtFfxMkV2itr+WkSaaeeIE11VvrNtvehMYppqSyQ8yaBZ06hVH+ORS6XGdGYdnTw25q0iRjM1luTeJrDxcIdJrrK+AvwC4ReV5EzsypHiYiMkFENojIZhG5t4r9t4hIvoh871tuq7DvZhHZ5Ftutlq2oDm0CpJOhfRu3Gj86B15a/T3a4PdpLQUMjMtbbbucfwg7PrUpE+JiXWky5ISmDPH/I8crVzYaig07Wyi/MNk/HiTwscb/UYvgU5zPaiqZwFXAE2B+SKyRkR+Y4UQIhILPAtcDPQCrhWRqnye3lPVc33LP33nJgMPYtyXhwAP+kr5ukPxPhP9XkGZ+DPv2u5lAybaWmIsdQ8GuOACY9j1bvZayJ5uPJw6X+NYl4sXQ2Ghg1NcfkTM79ydadKrhEFyMgwd6mWpjmaCSqeiqvNU9SZMud4dwP9YJMcQYLPPHnMceBeYEuC5FwGZqnpAVQ8CmcAEi+QKHr/xu4IymT3buANbHvVeFQ0aQ/Melo9MGjSAiy4yc/NeQaMa2PEuNDvrDLdwO5k92xRaG+3crNopOl8LWgbZH4Td1MSJ8N13sCv8+lseLhBsbq7+IvIUsB7oC/y3RXJ0wETY+8nxbavMFSKySkQ+EJGOQZ7rDJWUyZEjsGiRw2+NSf0tH5mA+Q27d8OKFZY3XTcozoc9c6HTjxydb5o1CzIyTBoVx0nqB4m9YHv4U13+e+Szz8JuysMFAjXA/1JEvgOWAm2BG4DOqnq/RXJUdedVDpH7BOiiqv2AOYA//DaQc82BIreLSJaIZOXnhzcsr5ZDq0ya7sZtAZg71+S1mjTJnu6qpEV/E+ty/JClzU6YYJ6R3lRXNWR/aN7SO//IsS63bIENGxy+vioiYkYn+V9C0c6wmurf39TQ8a6v6CTQkclNwMtAiqpeq6pfqFoaD50DdKzwPRXIq3iAqu5X1RLf15eBgYGeW6GNl1R1kKoOat06vPTZ1VLJ+P7pp9C8OQwfXsM5VpPkywd2yNpCEW3amJiAT8OPU6ub7HjX1K9Jci6fuv8t3hF7XHX47UM73gurGRHzO774wjh7eEQXgRrgB6vq86pq7avuKZYD3UUkTUQaAtcAMyoeICLtK3y9FDPVBvA5MF5EWvgM7+N925ynvBQK1p72MMnMhFGjTAZex2hhj0cXGK+bb77xXITP4Gge7F3k+BRXZqaxxXWzpoRNaDTvZmxEO98Pu6nx4+HwYc9FOBqJiBrwqloK/AKjBNYD/1HVtSLysIj4i9veLSJrRWQlJiL/Ft+5B4BHMAppOfCwb5vzHNls/O59ymTLFti2zdwgjtK4AzRMtsVuMm6cCSybP9/ypqObne8D6ugU14kT5v8wbpzDLsFV0ekqU6a6aEdYzYwZY36L54IefUSEMuL+AU4AACAASURBVAFQ1dmq2kNVu6rqo75tD6jqDN/nP6hqb1Xtr6qjVPWHCue+qqrdfMtrbv2Gkw9v38jgiy/MV9tqS1SHiLGb2DAyGTrUGHq9m70SO9/zGaN7OtblN9+Yt3jHX1aqouMVZr3zw7CaadkSBg48de94RA8Ro0zqBIdWmRomCeaBkplpopK7d3dBlqT+ULAGyq3NT9GwofEc8pRJBY7mwL6l5u3cQTIzzXuDKy7BlWneFVqca4mL8LhxJp3+4fCytHg4jKdMrOTgKmOAjW1EaanJEuzaFESL/lB2DI5ssrzpceNMLfvt2y1vOjrJ9hXj8L+dO0RmpnGISE52tNvq6XilUapHc8Nqxj+VumCBNWJ5OEPAykRERojISyLyie/7QBG50D7RopAKnlxZWcZI7doUxEmPLnvsJuCNTk6S/aGJtXBwiqugAL7+2oUp1JrodKVZZ4dX1H3YMFOQzbu+ootA40yuAz4GigG/AlHgYZvkij6OH4KjO0+W6f3iCzMiGTPGJXkSe/rSqqyxvOmePU08gHezA8V7TYyFw6OS+fPN23tEKZOEs02C0zCnuho1Mtm1PbtJdBHoyOQ+4CJVvRtTdRFgDeBkjtLIplLke2YmDBhgDIquEBsPzboZV2WLETEjrrlzvZTh5Ew36f5dmOJq2tQ4REQUHa+EvV/Csd1hNTNunEmQujO8OEgPBwlUmaSoqr8WrD9YsRRwJi1qNHDwlDI5csQYEF33sknsbYsyAXOzHzjgpVZh54fQrKujgYpglEl6unmLjyg6XQmoUbJh4L93vNFv9BCoMtkiIsMqbRsGbLBYnujl0Cpo2AIap7BggYngdX0KIrE3FPpiXyxm7Fizrtc3e8kB2DPPPEAd9LLYscM4QLh+fVVFYm+TaHRneFNdvXp5U6nRRqDK5C/AxyJyPxAnIvcA7+DZTE5RsNaX/l3IzDRFsIZVVr9Ok9THTMEc/qH2Y4OkTRuTS6le3+y5M0BLXZnigghVJiKm7v3ehaa2SxjNjB1r6rR4Waqjg0DTqUwHrsPUDNkBjAZuVVUvSxOYUooF68xbGcYleOTICJiC8MnDIXumusaOhSVLTG37esnOD6FJJ0fTzYO5vtq1M2/vEUnq5UbJ5oZXnGTMGFN9cbW1KeY8bCJg12BVzVTVS1S1j6pOUtX6/E56OsW74cQhSOzFnj2wdm2kBJL1AGlgm91k1ChT5W/pUluaj2xOFJqiUB0vd3SKS9Uok1GjIiCFSnW0HAyN24dtNxk1yqznzbNAJg/bqVaZiEhKIIuTwkYs/od1Yu+TOasiQpnENoTm3W1TJiNHQmxsPb3Zd30O5SWQepmj3f7wA+zZEyHXV3VIDHSYArs+g9LQh60dO5rsEfXy+opCahqZ5GCKTlW3+Pd7FKwz68RezJsHiYlw3nnuinQSGz26EhJg8OB6erPnfGySabYe4Wi3/r91RCsTMEq2tMgUCwuD0aNh4UIvJX00UJMySQPOqmHx7/coWGseLPFtmT/fuGw2aOC2UD6S+kDhVig9akvzo0ebdOFHjtjSfGRSfgLyZkKHSyDG2X/0vHnQuTOkpTnabfC0HQVxCWFPdY0eba6t776zSC4P26hWmajqjkAWJ4WNWArWQWIvdmYLmzefmuuNCBJ7A2qLRxeYm720FBYvtqX5yGTvl8ZTKXWKo92Wl5vI99GjI9he4ie2IaRMgpwZYSUbzcgw63o5+o0ygsnN9SMR+UxE1vjW19gpWNSg6nML7hVZ9hI/Jz26rE+rAsb9uWHDenaz53xsMgy0dzYqdeVKOHgwwq6vmki9DEryYd+SkJto0wb69Kln11eUEmhurt8D/4spPvWkb/2Ub3v9pniPeUtN7M28edCqlbn4I4bm3SAmzja7iT+ept7c7OqL7m43Dho0dbRr/8tKRI18ayJlAsQ0PJVVOURGjzYj35KS2o/1cI9ARyZ3ARer6p9U9RVV/RNwMaY6oiWIyAQR2SAim0Xk3ir2/1ZE1onIKhGZKyKdK+wrE5HvfcuMyufais/4rgm9TrpsxkRSYv+YOGh+tm3KBMzNvmKFSa9S5zm00iT0dNiLC4zCPvts6NDB8a5DIy4B2o7xBXdq7cdXw+jRJpbp668tlM3DcgJ97DUBVlXathpobIUQIhILPItRUL2Aa0WkckjWCmCQqvYDPgD+XmHfMVU917dcipP4HtLbDvQmJydCpyCS+tiuTFSN102dJ3u6z/X1Eke7PXHC/H0j8vqqidRLoXALHF4fchPp6eYFrd6MfqOUQJXJv4DfVNr2a+ANi+QYAmxW1a2qehx4FzjNuqmq81XV75K0DEi1qO/wKFgHcUlkLm4HROjNntgbirabQDsbGDzYZLCtFzd7znRoNQziWzva7bffQmFhhF5fNdFhslnnhD5hkJRkMnDXi+sriqkpaDFTRL4QkS+AfsCjIrJdRBaKyHbgUaC/RXJ04PSYlRzftur4CVAxlUu8iGSJyDIRqXb+QURu9x2XlZ+fH57Efg6vg6TezJ8vpKS4VKK3NvxG+DDeDmuiYUMYMaIe3OxFO800VwdnB79wyl6Snu541+HRJBWSB5qprjAYNcpk4j5qj4e7hwXUNDJZDHzlW77EKI/XgHm+9V99262gKkfHKidZReQGYBDweIXNnVR1ECZ/2FMi0rWqc1X1JVUdpKqDWre24M3S58mlCb1YsMC4MUaky6bNHl1gfvu6dbB3r21duE/uJ2ad6rwyWbDAOHZYcdk6TodLYd8yOLYn5CYyMsxUX71M3RMlVBtxpaoPOShHDtCxwvdUIK/yQSIyFlOoK11VT/p2qGqeb71VRBYA5wFb7BQYMG6PJfvZW2Jycvl94iOOZl2NV41NIxM49dsXLYIrr7StG3fJ/cSkp0k429FuT5ww3ky33upot9aReimsfhDyZkHX0H7EiBHGbrJggYvVSz1qJJg4k1gROUdERorIhf7FIjmWA91FJE1EGgLXAKeNi0XkPOBF4FJV3VthewsRaeT73AoYDqyzSK6a8Rm1l28wb/4Rq0xiYk3SxwL7/iwDBxq7yYIFtnXhLieOwJ75rkxxZWWZ6Z2Ivb5qI6k/NOl4amQXAgkJ5hqrs9dXHSCgXBAiMgD4COiEmX4S37oMaBiuEKpaKiK/AD7HVG98VVXXisjDQJaqzsBMazUD3hczl7TT57nVE3hRRMoxyvExVXVImZhuZi/pRUoKdOvmSK+hkdgLDmTVflyIxMWZt8c6e7Pv+gLKjzvuxQWn/qYXWvXq5jQiRglvfc0kfmwQmhNoRgY8/bRRrE2aWCuiR/gEOjJ5CpgGJAKHgQTMKOEWqwRR1dmq2kNVu6rqo75tD/gUCao6VlXbVnYBVtUlqtpXVfv71q9YJVOtFKxD4xKZ9llK5NpL/CT0hMJtYWVxrY2MDJN+3yrfhogi9xNTSbP1cMe7Xrgwiu0lflIvhbKjpjJliGRkwPHjxhDvEXkEqkz6Aveq6hFAVLUQ+H/U90qLBesobtiT3bsl8r1sEnsCCkc22taF/2+waJFtXbhDeZmZ70+Z6HhiR7+9JOKvr9pokw4NmptUNCFS0W7iEXkEqkxOVPhcICJtfNvaWS9SFHF4PTsO9QSiYD47wchJgX1G+EGDzPRDnbvZ9y+Dkn2uTHF9+y0UFUXB9VUbsY1MepW8maaUdAgkJJh4kzp3fdURAlUm3wL+itMLgH9jAgsrR8XXH44fhOI9fLepJ+3bR2h8SUUSepjIbRs9uuqs3ST3E1Oxsv0Ex7uOentJRVImw7FdcHBFyE1kZJi0KvW2VHQEE6gyuQ1Y6fv8W2AnUAL82A6hogLfG/5nS3tGvr0ETJbbpmm2KhMwN/uaNXXMbpIzw0zTNEx0vOsFC6B3b5M9N+pJmWheaHJC9+ry7CaRS0DKRFVzVTXb93m/qv5UVa9R1Q32ihfB+B7KX63pGT1TEAk9bZ3mgtPjTeoER3x5pfxpQRzEby+JmuurNuJbQauhYbkIe3aTyKWmdCrnV/g8rLrFGTEjkIL1lGojtud3iR7jaGJPY4Avt68Gqt9uUmeSPubNMmsX7CXffWfsJVFzfQVCymQ4+B0czQ3pdH9JbE+ZRB41jUzmVPi8uJrFqnQq0cfh9eQeOZs2bWPp0cNtYQIkoaeJlSjcZlsXcXEwfPipXFJRT+5MSDgHmleZocdW/A/MOqVM/ErZr6RDICPDTHN5dpPIoqayvc0rfI6pZol1RszIQwvW8/3WnqSnR4G9xE+iz6PLZrtJerqxm+zfb2s39nPiMOxd4MqoBMzormfPOmIv8ZPYy9juwrCbpKcbu8k331gol0fY1GozEZEGIjJLROKdECgqKD0GRdv5bkvP6HprTHBOmQB8Ge3j1l2ZUH7CFXtJaWkdiS+pjIj5e+6ZA6WhpQAeOdI0U2emUusItSoTVS0FBgL2TbRHG0c2ICjrc6NMmTRMhMYpthvhBw+G+Pg6cLPnzYS4JFO/xGG+/x6OHKmDygTMSK+sOORo+KQk6N+/DlxfdYxAXYP/jYUleqMe38N499Ge9OzpsizBktDT9pFJo0ZwwQVRfrNrOeTOgpSLHY96h1N/uzoRX1KZNhdCg2ZheXVdeKFJR3/8uIVyeYRFoMpkAPB3X332Of6iWb7CWfWPw+spK4+h/dk9osde4ifR5x4cRk3uQEhPN2/XBQW2dmMf+5ebEgMu2UsWLTKJQ1NSXOneXmIbQfuLjHNDiNdheroxwGfZl7vUI0gCVSaLMMWx/o3x4PqqwlLvKNq1nq17z2L4yEZuixI8CT2h9AgcC801M1DS081zYvFiW7uxj9xPQGLNQ89hysuNvalOTnH56TAZjuXBwe9DOt0/Yovq0W8dI6Dxu8OFsiKeE/vWG3vJTW5LEgKJFXJ0NUm1rZsLLjBuwgsXwqRJtnVjH7kzTYbgRsmOd716NRw8WMeVSfuLATF/5+Tzgj69VSuTGWDhQvjDH6wXzyN4IqU4VvRQXkoz3ci2/T3p29dtYULAIY+uxo1hyJAofXMsyja13lOc9+KCOm4v8dO4LbQcYpwcQuTCC+Grr4znm4f7BKRMfMWxtmAqGC7wLfM5PbCxflC4lQYxJ9DEnsQErIojiPi2xkPp8A+2d5WebrLeFhba3pW1nIx6d0eZLFoEnTubpU7TYTLs/ybk2vDp6ebaWhF63kgPC4mY4lgiMkFENviM/PdWsb+RiLzn2/+1iHSpsO8Pvu0bRMTWSe4D280bfeuu0ebG5UPklBHeZtLToawMliyxvStryZ0Jzc4yke8Oo2qUSZ2e4vLjV9Z5s0M63f83isrRbx0kIopjiUgs8CxwMdALuFZEelU67CfAQVXtBjwJ/M13bi9MzfjewATgOV97tpC91jyEe53v/IPGMhxwDwYYNgxiY6PsZi89CnvmmikuF1z11q83GZfrhTJJ6g+NO4ScWqVdO+jRI8qurzpMpBTHGgJsVtWtqnocUytlSqVjpgBv+D5/AIwRUwx+CvCuqpao6jZgs689Wzi2ez15h1LoO9D5dOSWkXAOFO8xNVlspFkzk/gxqm72PfNMQJ1LU1z+v1W9UCb+aPhdn0NZaAEj6enG862szGLZ6ihrF69g2eNXsHPdFsvbjpTiWB2A7Arfc3zbqjzGF5VfALQM8FwAROR2EckSkaz8EAtuHC9vwvbCETRwPo7NOhLtr7ro58ILTQ6lqEnKlzvTBNS1ccf6vWiRiS056yxXuneeDpOhtBDyQ6tZcOGFJpZp9WqL5aqjHFj1MYPbT6dpkvUvw6EWx9qBtcWxqppPqBzNVN0xgZxrNqq+pKqDVHVQ69atgxTRcOE9zzPs/70X0rkRw0mPLmeM8CdOREkxI1WjTNqPN4F1LnS/cCHRlTw0XNqONoXbckPz6vLsJsHRtnQmK3OH0jKlleVth1oc63aLi2PlAB0rfE8F8qo7RkQaYJwBDgR4rkdFmnaBmEaO2E38xYyi4mY/tNIEc7rkErx5M+zaVU+muPw0aGIUSu4nIUXDd+wIaWlRcn25TOmRPHq0+padZfZc34G6Bm8WkT+KSJXTRxawHOguImki0hBjUJ9R6ZgZwM2+z1cC81RVfduv8Xl7pQHdAS85dU3ExELC2Y5McyUmwrnnRsnNnuszBKdMdKX7emUvqUiHyVC4FQ6H9m6anm6mB8vLLZarjpH9tfGaiz/LRWWCSaUyHtgmIp+KyFUiEmeVED4byC+Az4H1wH9Uda2IPCwil/oOewVoKSKbMVNt9/rOXQv8BxMD8xlwl6p65rjaSDjHkZEJmHntZcugpMSR7kInd6YJpGvc1pXuFy0ytUvOPtuV7t0jxZciIcQAxvR0UztnvTOXc9RSumMm2/M7c256b1vaD3Sa6zVVzQB6AlnA48AuEXnaKkFUdbaq9lDVrqr6qG/bA6o6w/e5WFWvUtVuqjpEVbdWOPdR33lnq+qnVslUp0noaSoultpvGU9Ph+JiWL7c9q5Cp3gv7P/61IPNBRYuNIq33thL/DTtBEn9QrabeHm6AqCsmI5xmSzeNpl27e25wIKK4VbVLar6J+AC4Gu8tPTRS2JPQOHIJtu7GjnSrCP6Zs+bDahrWYK3b4edO+vhFJefDpMhf3FI7uppaZCaGuHXl8uU7VpAfIOj7I+3zx4YbG6uKSIyHdgOJAF32CWYh804lKMLoGVL6Ns3wm/23JkmgK7Fua50X2/tJX5SJoOWQd7nQZ8qYv5uCxfaXlkhajmweiZFxU1o3SfDtj4CNcD/A+Mh9QzGptFPVYer6j9tk8zDXhJ6AOKIER7MVMSSJcZNOOIoO24C5zq4E/UO5kGYnGwy4dZLWg6BRq3Cspvs2QMbN1osV11AlUb7ZjJn7VhGpttXfT3QkUlHTB6uzqr6B1X1/mXRTmw8NEtzzAifng5FRfDdd450Fxx7F5rAOZei3sEY30eOJDqTh1pBTKzxosv7FMqDTwPsH9EtCi32sW5TsJaE2B0sy55Mx461Hx4qgRrgr1LVT1XVc76rSziUowsi3EiaO9Mo17aj3ek+F7ZsqcdTXH46TIbjB2Bf8BGu3bubXF0ReX25jOaY8sjFyfa6vNfX9yAPMEb4wxuh3H5P6rZt4ZxzIvBmVzUBc23HmgA6F6j39hI/7caDNAhpqkvEvLB4dpMzObb5E7K2DqT/BXaFCRo8ZVKfSegJ5SVQtN2R7i680JTxjaikfIfXQ9E2V6e4Fi6EhATo3981ESKDhokmJ1ruJyGdnp4OOTmwbZvFckUzxXtpXLSMT1ZcYnuxNU+Z1Gcc9OgCyMiAw4fh+9DKftuDP7bBZWUycqRJ11/v6TAZCtaZGKggycgw64gb/bpJ3mxElKy8yaSl2dtVoN5cVeo0ERlprTgejpLoq8nikEeXfxpnwQJHuguM3JnQ4jxoYu8UQHXs2gUbNpx6ENZ7/HE+IYxOevaE1q0j7PpyGc2dya5DKbTsPsB2R8VARybVTWJ+bJUgHi7QsAXEt3NsZJKSYooZRczNXrIf9n3l+qgEPGVykubdzIg5p3JqvtoRMX/HBQs8uwkAZSWU537OjG8nk5Fhv8t7oMrkDElEpDngeXdFO4m9zLSCQ2RkGPfNiLCb5H0KWu5a1DuYB19CgkmG6eGjwyXGXft4QdCnpqebTALbt1svVtSxdyGx5YV8suISR5w7alQmIrJJRDYCjUVkY8UFU5Bqvv0ietiKX5k49CoXUXaT3BnQuD0kD3RNhAULjL0kqoutWU2HS0BLTSBpkPhHeBEz+nWT3E8oKW3MhoNjHCm2VtvI5C/AXzEleh+tsDwCXIVJFe8RzST2gtIjpo6HA0RMMaOy45D3ma/Wuzt+KLt3G3tJvXcJrkyrodCoZUh2k169oFWrCLi+3EYVzf2EBT+M5YIRjR1J7FDj+5CqvgEgIj+oajTUyvMIloReZl2wDpqk2t5dSooJMFuwAH77W9u7q578RUaJujjF5dlLqiEmFtpPhLxZJho+JvBhW0W7Sb2mYA1StIMPlt1HxvXOdBloBPwyEYkXkb4iMqziYreAHjaTWEGZOERE2E1yZkBsY2g3xjURFiyA5s3hvPNcEyFySb3UFw2/NOhT09Nhx456bjfxjepmrZjk2Mg3UNfgS4FdmDrwiyssX9onmocjxLc2CfYcViYFBbBypWNdno4/6r2de1Hv4NlLaqT9eIiJC2mqy7ObADkfs/ngYGKaptC1qzNdBjpZ/ATwENBMVWMqLGGHWYlIsohk+oz9mSLSoopjzhWRpSKyVkRWiciPKux7XUS2icj3vsXziwmWxF5w2Dll4nq8ScEaE/Xv4hTX7t3www/eFFe1xCVAmwzjJBEkfrtJvVUmx3bB/m/4YOmlZGQ4lwg7UGXSVlWfUtWjNshwLzBXVbsDc33fK3MUuElVewMTgKdEJKnC/v9S1XN9SyT4CUUXCc56dHXoYOwmrhlJ/W+7LmcJBk+Z1EiHS01d+CBrw8fEnKpvUi/xXd9vLZri6PUVqDL5QkQusEmGKcAbvs9vAJdVPkBVN6rqJt/nPGAv0Nomeeofib1MhbviPY51mZFhbnZX7CY5MyB5sHELdgnPXhIAqZeadU7wsdEZGcZmUi/tJjkfc7g8jTXZfRz1FAxUmWwHZojI/4nIHysuFsjQVlV3AfjWbWo6WESGAA2BLRU2P+qb/npSRBrVcO7tIpIlIln5+fkWiF5HcMEIP2qUsZs4Hm9ybDfs/8bVKS6AefNM4kvPXlIDTTuZVDchKJNRo8x6fn2LhDtRCLvnsnjbpaSmCt26Odd1oMpkILAW6AOMq7CMDeRkEZkjImuqWKYEI6yItAf+Dfy4Qm2VPwDnAIOBZOD31Z2vqi+p6iBVHdS6tTewOUmir7yfw0Z4MA9VR8n9BFBIDerSs1aEXBNf4n/gedRA6hTj0XUsuFFzr17Qpo0L15fb7P4Cykt4+dMpjB7tbOHQgN6LVDWsy15Vq1U6IrJHRNqr6i6fsthbzXEJwCzg/ooxL/5RDVAiIq8BvwtH1npJfFuTp8tBI3z79iYx37x58F//5Vi35i23aRok9XWw09Pxvy2PdqcWV3SRehms/rOpcdL1JwGfJmKU9bx5xhToUjVm58mZQWlsC2Z+M4KXb3e264BDf0Uk1hdb8iPf9yYi0tgCGWYAN/s+30wVySNFpCEwDfiXqr5faV9731ow9pY1FshUvxBxPEcXmIfpl186WBf+xBHYPcc8oFx8usybBy1aePVLAiKpHzTtDNnTgz519GjIy6tHdeHLSyFvJpuLJlJaFuf4yDfQOJOumIf0bOAV3+bxwMsWyPAYME5ENmGmzh7z9TlIRP7pO+Zq4ELglipcgN8SkdXAaqAVJgWMR7AkuKNMiopg+XKHOtz1uSkG5uIUFxhlMmpUPa73Hgwi0GEK7JkDpUVBneof+dWbqa59S6BkPx9/O4WuXaFzZ2e7D/Ry/j/gXYxNwv8euQAIu56Jqu5X1TGq2t23PuDbnqWqt/k+v6mqcRXcf0+6AKvqaFXtq6p9VPUGVS0MV6Z6SWIvKMmHYuccE9LTzbPCsZs9Z7rJ+dR6uEMdnsm2bSY625viCoLUKVBWDLu+COq0rl2hY8d6ZITP+RiNacgzH17kyvUVqDIZAjzqM3orgKoeApJqPMsjevB7dDlU2wSgZUsz1eOIMik/AbmzjBdXELmerMb/Wz3jexC0GQlxSUF7dYkYpT1/PpTX9WIZqpA9jcPxY8jZkxDRyuQwlRSHiKQAzgUmeNiLC+7BYG72JUuguNjmjvYughOHzJSJi8ybB23bGucDjwCJiYMOk4wnXnlpUKeOHg379sGaum5JPbQKiraxNGcq4M7LSqDK5CPgVRFJBRCRlsBTmKkvj7pA4w7QoDkccvauGz0aSkpgafD5/IIjZ7pJ7Nh+vM0dVY+qUSZOu2zWCVIvN4kf84NLB+h/qNZ5u0n2RyAxvDbnUnr3Ni8sThOoMvkTUAjsxIxQ9gIlmFonHnUBEUjqY/JWOcjIkRAba/PNrmqmSNqPdzWx4w8/mJxcnr0kBFImmJeB7I+COq1jR5O6p84rk5xplLcczszMNq5dX4GmoD+mqtdhotOHAO1U9UZVtXtywsNJkvrCodWOFtBOSIDBg22+2Q98C0ezI2KKCzxlEhINmkL7CUaZaHAGkNGjTeqe0uBmyKKHI1vg0Gq2l17O0aPuXV9BOSeq6j6fl5WXi6QuktjXTCUc21X7sRYyejR88w0cOWJTB9kfgsSeyvXkEnPnQqdOkJbmqhjRS8epcCzPpMMJgtGjTanob7+1SS63yZkGwIwVlyNi0vS4QbXKxJfWfWtti5PCetiMPyr80CpHux071rw12pIyXNUok7ajjFuwS5SWmpHJuHGevSRkOkwGaRD0VJffRpWZaZNcbpM9DVqcywefdmHQIEhOdkeMmkYm92NsJX8CXgMaA29jasC/A8QDr9otoIeDnFQmqx3tdtgwaNLEppu9YA0c2QQdr7Sh8cDJyjKJLceNc1WM6KZhkqmMmf1RUFOxrVqZ7Mx1Upkc2w37llLceirLlrl7fVXrcK+qb/k/i8gi4BJVzaqw7UOMR5cXcV5XaJRsvLocViaNGpmh+RfBxaQFxs4PADEpVFwkM9O8HY9xr0pw3aDjVPjmDnONtugX8GnjxsE//mGmUps3t1E+p8n5GFCWZl9OWZm7yiRQm8m5QOVk4at82z3qEkl9ocBZZQLmJtiwAbKzLW44+0MT9NbYBV/JCnzxhXk7btXKVTGinw5TAAl6qmvcOJMDrs4VzMr+AJp358O5vWnSBIYOdU+UQJXJBuA3lbb9GqgvKdTqD0l9oWB90MFh4TLeF/5h6VTE4Q1QsBY6XmFho8Fz5AgsW3bqN3qEQeO25uUgJzhlMnw4NG5cx6a6ivfBnvnQ6SoyM4WMDDPKd4tAlcldwO9FZLuILBSR7ZjyunfZnJQAgAAAIABJREFUJpmHOyT1NckQj2xytNvevU1aektv9uwPzbrjVAsbDZ4FC4wB3rOXWETqVDPNFUQ53/h4M5Vap5RJznTQMnbFXcnGje5fX4HGmXwDnAXch6kpch/QVVW/tlE2DzdwyQgvYry65syxMI/Szg+h5QXQJNWiBkMjM9O8FQ93L79k3aKTz5li5/s1H1eJceNg/XrIybFBJjfY+T4068rspcbaEBXKBEBVD6vqW6r6d9+6wE7BPFwioaeJyXBYmYC5Gfbts6iUb+FWOPgddHJ3iguMMrnwQnenIOoUTTpA6xGw8z9BneZ/2M6ZY4NMTlOyH/bMPTnFlZJiqku6SaD1TERErhWRJ0TkpYqL3QJ6OExsI2jew/FYEzAjE7BoKmKH70HjsktwdrZJo+LZSyym09Xmhacg8CzXffuanFW2eA06Tc7HoGWUp17JnDmREb8U6MjkeeAZoCMQV2kJCxFJFpFMEdnkW7eo5riyCoWxZlTYniYiX/vOf89XldEjHJL6uTIyad/e3PCW3Ow73oVWQ6FZFwsaCx2/YnR7CqLO0fEKQIIanYiY/4OlU6lusfN9aJrGiu0D2L8/Mq6vQJXJVcAQVb1aVX9ccbFAhnuBuaraHZjr+14VxyoUxqqYF+NvwJO+8w8CgReK9qiapL5QtM2UuXWYceNg8WI4ejSMRgp+gEMrodOPLJMrVDIzoV076NPHbUnqGE1SjFdXCFNd+fmwcqVNcjnB8YOm/HSnq8icY4Yj/lG9mwSqTI5iMgbbwRTgDd/nNzB13APCV/d9NPBBKOd7VIPfCF+w1vGuJ0yA48fDrI638z1AoNNVVokVEmVl8PnncNFF7k9B1Ek6XW3q7xwK/Dr1Tzd+9plNMjlBzsegpdDpSj791MQvuZFyvjKBKpO/Aw/4Ht5W01ZVdwH41m2qOS5eRLJEZJmI+BVGS+CQqvqDInKADjbIWL9wyaMLjKG6SROYNSvEBlTNFFebdPP26iLLlsHBgzBxoqti1F1CmOpq1w4GDgzj+ooEdrwHTbtwUAbx1VeRc30FqkzuBn4PHBSRjRWXQE4WkTkisqaKJZic4J1UdRBwHfCUiHQFqlJu1SbtEZHbfQopKz/fS3xcLU07Q4NmrhjhGzUyQ/bZs0PMhH9oFRz+ATq7P8U1e7ap1eIZ322icTvz0rDzP0FdLBMnmmJsBw7YKJtdFO+F3ZnQ+Voy5whlZZGjTAIthh1W/i1VrXZGT0T2iEh7Vd0lIu0xhbeqaiPPt94qIguA84APgSQRaeAbnaQCeTXI8RLwEsCgQYOcK9oRbUgMtOgPB63w0Q2eSZNgxgwTExC0u+OO94xrs8tR72CUyfDhkJRU+7EeIdL5R7D858ZG1iKw7E6TJsEjjxhHj2uusVk+q9n5PmgZdLmO2U+ZDMHnn++2UIZAgxbfqG6xQIYZwM2+zzcDH1c+QERaiEgj3+dWwHBgnaoqMB+4sqbzPUKgxUA4uALKyxzv+uKLzXr27CBP9E9xtRsL8a0tlysYcnNNvEykvDXWWTpeadLSb3+r9mN9DBpkcqQFfX1FAtvfhqS+lCf04dNPjY0xNtZtoQwBBy2KyAhfbMknvu8DRcSKMiyPAeNEZBMwzvcdERkkIv/0HdMTyBKRlRjl8ZiqrvPt+z3wWxHZjLGhvGKBTB7JA6C0yPG0KmBKrfbtG8LNvv8b44UWAV5cfgOvp0xsJr4VpFwM298J+MUnNtY8hD/9NMpchAu3w74l0Pk6vvsO9u6NrOsr0KDF6zAjiGLAr0AUeDhcAVR1v6qOUdXuvvUB3/YsVb3N93mJqvZV1f6+9SsVzt+qqkNUtZuqXqWqJeHK5AEkDzTrA+6Up5s4Eb780tQACZht/4bYeNdzcYEx8Hbs6LkEO0KX6+FYLuwNPCXwxIkm28Ly5TbKZTU73jXrztcwa5bxELzoIndFqkigI5P7gPGqejfg1+VrgN62SOXhPgnnQGxjk5LEBSZONMkRA059UXYcdrxjUpQ3TLRVtto4ftzEl0yc6LkEO0KHS4zDSBBTXePHQ0xMlE117XgbWg2DZl2YPdvYSiKppEGgyiSlQmEsv+G6FIiQ2ToPy4lpAEn9XRuZDB0KiYlB3Oy7PjX169NuslWuQFi8GAoLI2sKok7ToIlxuMj+AMqKAzqlZUu44IIoUiaHVpuly3Xs3WtGVJF2fQWqTLaIyLBK24Zh6px41FWSBxgjvDo/sRwXZ4bwAbsIb/sXxLeB9u774c6eDQ0belUVHaXL9XDiMOQGHkAyaZIpp7xnj41yWcX2t42XYqer+Pxzc09MmuS2UKcTqDL5C/CxiNwPxInIPZg68GHbTDwimOSB5gY9ssWV7idNgt27zQ1fIyUHIPcT6HydGVG5iKpxa87IgKZNXRWlftF2NMS3g+1vBnyK/2E8c6ZNMllFeZl5WWp/EcS3YcYME3x5boTVuQ3UNXg6JljwfGAHJoXJrar6qY2yebhN8gCzdmmqa/Jk43kzbVotB+78D5SfiIgprnXrYNMmuPxytyWpZ8TEQudrIW+WebkIgH79oEuXAK4vt9n9BRzLg64/4dgx44V22WXG5hNJBFPPJFNVL1HVPqo6SVXrUs0yj6pI7A0xDV0zwicnQ3p6ADf7tn8ZWQMMWrOTadOM0X1KMLkdPKzhrJvNS8X2twM6XMQo/cxMU1o5YtnyKjRqBSmTmTMHiooi82UlmDiTTiJyr4g841t3sU8sj4ggJs6ko3dpZALmpvnhBxMNXyVHNsO+pWZUEgGuUx99ZAy77du7LUk9pEV/aDEAtr4a8CmXX2687yLWEF+8D3I/hi43QmxDPvrIOKZkZLgt2JkEGmcyAWNsnwQk+tbrfds96jLJA+HAdyEmygqfy3wpPasdnWx5xaR/6XK9YzJVx/btsGJFZL411hu63mqcRg6sCOjwYcOgTZsInura8bYZbXX9MaWlxh53ySXGwSPSCHRk8jjwE1Udqao3qupI4FbgCftE84gIkgfAiUMmstwFUlNhyJBqbvbyE7D1NUiZZEq5usz06WbtKRMX6XIdxDQKeHQSG2umJGfNguLAvIqdQ9W8LCUPgqS+LFpkklNG6vUVqDLpArxbadt7QCdLpfGIPE5GwrtjNwFz82Rlwc7KFXVyP4HiPdDtdlfkqsy0aSYNTLdubktSj2nYwmRA2P5WwDEnl19u4oLmzrVZtmA5uMJkwe56K2Cur/j4yIp6r0igymQBkFFpWzoQeP4Cj+gksY+xnbhoN5nqy47if/M/yeaXoXEHaO/+bOvevSb9S6S+NdYrut5qqhFmV75gqmb0aEhIiMCpri2vmPRAna+lvNzIN2FC5LqcB6pMNgPTROQNEXlIRN4ApgMbReSP/sU+MT1cI7YRJPaF/V+7JkKPHiYV/Wk3e9EO2PW5eXC4HFsCZi5b1VMmEUHb0aYmT4BTXY0amZiTjz82KXwighNHTK65jldBw//f3pmHR1VeDfx3IIIE9GMLIPsqFQVFIohFFJcilqoUNwSLKFoVty4uxbpW3D77uaAoqAjKJhVcWisuCEWlUkFBtOybRFAWSdiXkPP9cW7CEGaSSWYmkwzn9zzz3Jn3vvfe8+ZO5tz3nPOeU5O5cy0TdXn+fkWrTE4CvsTMWt2D7ZdYTZFzg1c5qELsJISMbrBpjvkokkSfPjBrltXvBuypDaDVNUmTKZSpU23NwoknJlsSB6kELQdZnfTt0fn6+vSxxI+ffppg2aJl1WuQuw2OvQmw71flyrb2qrwS7aLFHlG8zkq0sE6SqNcN9u+MOkImEVx6qaUL/9vfgLxci70/pqc9gSaZjRttrcIll5SL6GQH7CFDKsHS56Lqfv75Zj6aODHBckWDKix91hzvdTuTlweTJlkF0tq1ky1cZEqyzqSyiJwmIpcFn9NFpFriRHPKDRmn23bjJ0kToUMHS+c+bhyw7p+Wcrz1tUmTJ5TJk808MmBAsiVxCkhvbMkfV7wE+7YX2716dQtDnzwZ9iS7iMWGmbB1ERw7BIDPPoM1a8r/9yvadSatsJTz/+RA8alfAC8mSC6nPFGtAdRonVRlAvbP9O9/w675T9uPRaNfJVWefMaPtyiuDh2SLYlzEG1vhX05sPq1qLoPGADZ2ZauJKksfQ6q1C4o8jZ+PKSnH1hzVV6JdmYyHAsNrg3kG85nAqfHKoCI1BaRD0VkWbCtFaZPDxGZH/LaLSIXBfvGiMiqkH3Jz6mRitTrBhs/TdriRYB+/eCEJguptvVjaDPEosySzIoVpuD6J3/NpFOYul0ttH3JM1F9b885xxYwjos+V2T82ZkFWW9Bq8GQVo29e222dNFFUKNGEuWKgmiVSWdgmKrmEdQzUdVsoGYcZLgLmK6qbYDpweeDUNUZqnqSqp6EJZncCXwQ0uX2/P2qOj8OMjmFyTgd9myGrYuTJkLTpvDoVU+za181tFX5MHGND+oxXXFFcuVwwiBis5Oti80ZXwxpaXD55ZZFODu7DOQLx7KRVvKhzfWAzZK2bCn/Ji6IXplspZDiEJGGQDwqAVwIjA3ejwWKm8xdDLynqjvjcG0nWjK62TaZpq7dm+jZdjyvzrqSed/USZ4cAaqmTM44w0r0OuWQppfCkfVhydNRde/f33wmU6YkWK5w5O6AZSPMfFujBWCzpIwMOPfcJMhTQqJVJlOB0SLSGEBE6gBPceiq+NJQX1XXAwTbesX0vxyrpRLKMBH5WkSeFJGqkQ4UketEZK6IzN1YEGPqRMVRbaz41IYkxk6uGEWa7GbEx7cm1xQRMHcuLF1aMZ4aD1sqV4XW11tq+q3F1/I75RRo0+bAjLNMWfGyVQttdycAOTnw97/bbCkt+UupiiVaZXIPsB34DpuhbAD2AA9Hc7CIfCQi34R5lShRt4gcA7QH3g9p/hPwM+AUzKdzZ6TjVXWUqmaqamZGRkZJLu2ImKkrWTOTvH3mmGxwLq06tmPSpOQvMBs/3hLuXXxxcuVwiuHYG6FyNfjvo8V2FbHZycyZkJWVeNEKyNsHi56w/7EMK2o7ZYrNkiqKPy7adSa7VPUKIAPznzQIEj5GlfxGVc8J6qAUfr0N/BgoiXxlsaGIU10KvKmqBavnVHW9GnuAVwL5nESQ0Q12rDYnYVmzZpIVCGp7K1deaaVWk5k2fPduUya9e0PNeHgOncRxZD3L37bqNdi+utju/fubCXPMmIRLdoDVE2HnWmh3wGX8yis2S+pcQX7RSrLO5FjgWGCrqsbTRvQOMDB4PxB4u4i+/Shk4gpRRIL5W76Jo2xOKPWC4L2yNnXl5cI3f7HaKg170bs3NGwIzz9ftmKE8sYbtmL6+uuTJ4NTAo673WqoRzE7ad3aIrtGjiyj2a/mwaLHoWZ7aNgLgK+/ttX4v/1txVkIW6wyEZELReQ7YBHwGVbH5Lv80Nw48Chwrogsw9KyPBpcN1NEXgqRoznQhEOTS44XkYXAQqAuVq/eSQQ1T4S0GmVv6lo9AbYtg/b3g1TiiCPguutg2jRYvrxsRcnnuecsZ9jZZyfn+k4JSW8ELa+2kgVRzKyHDDEzV5nUh//+Xcj5Fo67s0BzjBhhGYIHDSqD68eJIpWJiJwM/A0Lw+0JtAu2HwCTRaRTrAKo6mZVPVtV2wTbn4L2uao6OKTfalVtFIQnhx5/lqq2D8xmA1S1+OWuTumolGY23fUflN16k/xZSa2ToPGB55drrzWn5AsvlI0YoXz5JXz+OdxwQ/mrw+0UQbs7QffDf/+32K69e1uE3nPRZWMpParw7UOWFqiZLVLMybEorn79ynf6lMIU969wG/CYqg5W1Y9UdXGwHQw8Eux3Dica/wq2L48qMiYurB5v1zvhvoPm+w0bWnK+0aNhZxkHiY8YYSuSr7qqbK/rxEiN5tDiSlgxCnb9UGTXtDQzMX30ESxJ5Fc9603Y/B844d6C7Nevvmp13ocMSeB1E0BxyuQ0YFSEfS8G+53DiYZB2tLv30n8tQpmJR2h8aGBf0OG2IKu119PvCj5bNkCEyaYk9Yd7xWQ4++279XC+4vtOngwHHFEAn1zebmwYCgcfRy0+A1gE5URI6BLF+gUs92nbClOmdRR1bXhdqhqFpD8lWNO2VK9iZmcvv974q+1cjRsXxH4Sg71QnbvDscfb6aIsrK6jRkDu3bBjTeWzfWcOHNUa2hzI6x4EbKLjtWpX9/CvseMsZlC3Fk5xmb4Jz5cMCv5+GNYvLhifr+KUybF7a8gcQZOXGl0AWyaDbs3Je4ae36yp7aMbhETOorYP928eTB7duJEySc31xTXaafBSZ4BruLS/l5IOxq+ur3YrkOGmA9j7Nhiu5aM3F02O6rb9aBZ99NPQ926VnKholGcsqgaWkmx8AuoUhZCOuWMRr+ycMZ1CVzo8fW9Vno189kiYyMHDrTkfPffnzhR8hk3zhI73nFH4q/lJJCqdeCEe2D9NFj3fpFdTzvNXo88YmuL4sbS4VZG4aRHC77f8+bZivebb7ZIroqGaBH2ARGZSZDYMRKq2iPOMpUJmZmZOnfu3GSLUTHRPHiriT1Vnf5G/M+/ZT5M62TmiMzhxXZ/8kn4/e/hX/8y01ci2LcP2ra16Jovvqg4sf9OBPbvgXePt5XxveZDpcoRu06fbutOhg+Hm26Kw7V3ZsE/2kG97nDmgdjj3r1thr16tdWkL6+IyDxVzSzcXuTMRFXPLK7CYuJEdsotUgka9bYa7PvjXElIFebeZPUcOjwY1SHXXw8NGsB998VXlFDGjoVVq+DBB12RpASVq0LHxyHnG1j6TJFdzzrLHlIeftj8ZTEz92bQXMg8cN05c+Ddd+H228u3IikKj5J3SkejCyB3O/w4M77nXTkGNn5m0/8qh5S2CUu1ajB0qOVTmjEjvuIA7N0Lf/mLRdj06hX/8ztJonEfM9kuGApbl0bsJmIPEevX26r4mFj7ptUraX8f1GhZ0HzffeYricvMJ0m4MnFKR/2zzEQQzxDhbcth3i02/W9ZsqW/114LjRrBvffGP7Jr9Gj47juflaQcItB5pH2PP78K8vZH7HrGGTZDeeSRGCK79m21WUnNDvCz3xc0z54N779vvrijjirlucsBrkyc0pFWzUxdaybCvm2xn2//XvisH0gadB1nprQScOSRcPfdls/o7aKyu5WQnByblfz85xWjpoRTQqodA52Gw6Z/w5Ini+z6wAOwYQM88UQprzX/LktW2vnFgiqheXlw550WRFIRw4FDcWXilJ6f/cEirpa/GPu5vr4HfpoLXV6ytSyl4JprrA77DTfATz/FLhKYY/+HH8zJ77OSFKX5FRaeu+DPkL0wYrdu3ay2yLBhsDByt/CsngDLnoe2t0HdA2mAn33WHoAeewyqVy+l/OUEVyZO6anbBer3gMV/jc0Rv26aZU1tdS007Vvq01SpYgvMNm2C2+KQ6GfaNDNx3XGHFU1yUhQROOUFqFob/nUB7I5cBWP4cKhVy1Lp7NsXsdvB/PQlzLnGzLcdHytoXr4c7roLfvlLC3Gv6LgycWKjXTB1X13K0nSb5sCnF5sduVPRZoZo6NjRnPGvvWYx+6UlJ8f8MO3alc0aFifJVGsA3d+B3T/CrD6wP/yikrp1Lb3Kl1/C/xafLxJ2b7TzVa0L3f52kHnrmmvsAWjkyNSY9boycWKjwbmWO2vR40U6MMOS/Q3M7AVHNoAe0yAtPvP8u++G9u0tUd/mzSU/XtXMW+vWWYGiqhELQTspRZ1M6PqqZXf4/JqIkRy//jVcdpk9ZHz9dRHn25sNn/SBPRug+1tWpCvg2Wdh1ix46ikLHEkFXJk4sSFis5OtS+D7Eni+ty2HGb+wSJqzPjRHaJzIN3dt3gw9e0J2dsmOf+ABM2/deWfFqXLnxImmF0OHh2DNBPjPtRYYEobhw6FOHTjvPFi2LEyHHWvhw26WEbjrq1D7QNbGyZPhd7+zRYqpYN7Kx5WJEztN+kKNVjD/T9Hl61r7Fkw7xfwsPT6AGi3iLtLJJ8PUqfbk2LOnma2i4aGHTJkMGmTvncOQ44daduEVL8OMnrDn0OltRoalp9+3D3r0sDQ7BWxZAB+camV4z5wGTS8p2DVlClxxhaVomTgxNcxbBahqUl/AJcC3QB6QWUS/84AlwHLgrpD2FsAcYBnwOlAlmut26tRJnTjyw0zViVVV3+ukujfH2laOU32zmep4se2KMapzb1Mdj+p7marbViZcrLfeUk1LU+3aVXXDhoP3jRun2qyZqohq06aqffuqguqVV6rm5iZcNKe8s/I11YlVVN9urbp+umpe3iFdFixQrV1btUkT1SULf1Kd9wc75s3GqlsWFvTLy1OdNOnAd3Hr1rIcSHwB5mq43+hwjWX5Ao4D2gIzIykToDKwAmiJJZdcALQL9k0GLg/evwDcEM11XZkkgKy/q05IU/2wu+ry0aqT0k1xFLzEtl/crJq7u8zEmjLF/omrV1cdOlR182ZTJOnp9h8Q+ura1RWJE8KG2apTG9r39p8d7QFpx/eqeftt/74duuTTT3RYvwd188haun+c6NYPB6nuXF9wipkzVU8/3b5fXbqo5uQkaSxxIpIyKTLRY1kSJJX8o6oekn1RRLoC96tqz+Dzn4JdjwIbgQaqmlu4X1F4oscEsXoSzL7CFh1qGId81XrQ98cyF2vJEnOYvv66xfPv3m0p5QvTtCmsWVPm4jnlmf27YdU4WPx/sHWRtUmaRYDt+sHybAGLss9jwF8fZ+Ha9rRubSasvXstBPiYYywwZPDgih/QESnRY1oyhCkFjYDQIl1ZQBesOFe2quaGtEeMjRCR64DrAJo2bZoYSQ93ml9uGVg/jVCQYc/GspUnoG1bs1HffTc88wy8GGGd5dqwpeCcw5rKR0LrwdDqatgwyxTKjrWW/bd6E6hzKtTtwnFH1uOdC22Ba+gDyZAhFllYrVryhlAWlIkyEZGPgAZhdt2tqtGEAIVzU2kR7WFR1VEEZYgzMzPLx5QsFWl6CaQ3g51hHvHTk6vETzgBRo2CDz4IPwPxZwwnIlIJ6p9prwg0ahRDupUKTplEc6nqOap6QphXtLGkWUBojo3GwDpgE1BTRNIKtTvJ5sRhUDn94LbK6dZeDhg2DNILiZeebu2O45ScihIa/AXQRkRaiEgV4HLgncAZNAO4OOg3EIhjmj+n1LToD51H2QwFsW3nUdZeDujf32YozZqZbbtZM/vcv3yI5zgVjqQ74EWkDzAcyACygfmq2lNEGgIvqer5Qb/zgaewyK7RqjosaG8JTAJqA18BA1S12ERR7oB3HMcpOZEc8ElXJsnClYnjOE7JKVXZXsdxHMeJBlcmjuM4Tsy4MnEcx3FixpWJ4ziOEzOHrQNeRDYCpU2cURdb43I44WM+PPAxpz6xjreZqmYUbjxslUksiMjccNEMqYyP+fDAx5z6JGq8buZyHMdxYsaVieM4jhMzrkxKx6hkC5AEfMyHBz7m1Cch43WfieM4jhMzPjNxHMdxYsaVieM4jhMzrkxKiIicJyJLRGS5iNyVbHnijYg0EZEZIrJIRL4VkVuD9toi8qGILAu2tZIta7wRkcoi8pWI/CP43EJE5gRjfj0of5AyiEhNEXlDRBYH97trqt9nEfld8L3+RkQmisiRqXafRWS0iGwQkW9C2sLeVzGeCX7PvhaRk0t7XVcmJUBEKgPPAb2AdkA/EWmXXKniTi7wB1U9DjgVGBKM8S5guqq2AaYHn1ONW4FFIZ8fA54MxrwFuCYpUiWOp4Fpqvoz4ERs7Cl7n0WkEXALkKmqJ2DlLC4n9e7zGOC8Qm2R7msvoE3wug54vrQXdWVSMjoDy1V1paruxeqoXJhkmeKKqq5X1S+D99uwH5hG2DjHBt3GAhclR8LEICKNgV8CLwWfBTgLeCPoklJjFpGjge7AywCquldVs0nx+4yVKq8WVGdNB9aTYvdZVWcBPxVqjnRfLwReVeNzrHLtMaW5riuTktEIWBvyOStoS0lEpDnQEZgD1FfV9WAKB6iXPMkSwlPAHUBe8LkOkK2qucHnVLvXLYGNwCuBae8lEalOCt9nVf0eeAL4DlMiOcA8Uvs+5xPpvsbtN82VScmQMG0pGVstIjWAKcBtqro12fIkEhHpDWxQ1XmhzWG6ptK9TgNOBp5X1Y7ADlLIpBWOwE9wIdACaAhUx8w8hUml+1wccfueuzIpGVlAk5DPjYF1SZIlYYjIEZgiGa+qU4PmH/Onv8F2Q7LkSwA/By4QkdWY6fIsbKZSMzCHQOrd6ywgS1XnBJ/fwJRLKt/nc4BVqrpRVfcBU4HTSO37nE+k+xq33zRXJiXjC6BNEP1RBXPevZNkmeJK4Ct4GVikqv8XsusdYGDwfiDwdlnLlihU9U+q2lhVm2P39GNV7Q/MAC4OuqXamH8A1opI26DpbOC/pPB9xsxbp4pIevA9zx9zyt7nECLd13eA3wRRXacCOfnmsJLiK+BLiIicjz21VgZGq+qwJIsUV0SkG/AJsJAD/oOhmN9kMtAU+6e8RFULO/kqPCJyJvBHVe0tIi2xmUpt4CtggKruSaZ88URETsICDqoAK4FB2ANmyt5nEXkAuAyLWvwKGIz5CFLmPovIROBMLNX8j8B9wFuEua+BUn0Wi/7aCQxS1bmluq4rE8dxHCdW3MzlOI7jxIwrE8dxHCdmXJk4juM4MePKxHEcx4kZVyaO4zhOzLgycZwIiEhTEdkuIg2TLUs4Atm6xton3ohITxH5JIbjjw8yGVeNp1xOYnFl4iQdEZkpIntEZJuI5IjIShF5TUQ6JVMuVf1OVWuoarlcER3I9m+w9TEikltUn7IgWLfwJLa2Ib+tj4isFZEfReSGQv1fkKDMQT6q+i223uOmspDZiQ+uTJzywl9U9ShV/R+gB7AG+FxE+iRZLqdk/AJbBDkDCso2vAj8Csu6/ZCI1A/2nQm0B4YPiRi8AAAE1UlEQVSHOc9o4GYR8d+oCoLfKKfcoaprVPXPwKvA8OBpFxGpIyKvish6EflBRMaKSO3840RktYj8Way413YRWSgiHUSkX1D8JyfIjpsWcswrwVPzNhH5r4hcEbKvuYhokJ4eEblfRKaLyMNixYc2BCuqwxJy/GARWSoi2SLytojUC+lT3JhuEZFVgXzfi8jDIftURLoFZrj3gMrBuLeLyMDQPiHH9BWRBcHfYkGoshaRq4K/0y0ikiUiW0RkZKAQouUi4CM9sBq6DrBPVeer6hpgGdBMRKphK68Hq2pemPPMAhoAJ5Xg2k4ScWXilGcmYaku8vNHjQdqYYXJjsPSRbxW6JiBwI1BvwXAm9hM50TsKfgC4NKQ/p9iP1g1gQeBMVJ0wbPuWDqKhtjT9lAR+Xkx4/hNcFxTLEXNuJB9EcckIscCjwK9VfUo4HjC5IILzHC9gP2BWauGqo4t3C/wnYzHsgPXwdLkTBSRLiHdmgH1gVbAKcAlWL6yaDkZy3eVz0Zgv4h0EpEWQHNgOfAQMEFVFx16CgjSmSwLzudUANKK7+I4SSMr2NYJnr57Aseq6hYAEfk9sFhEjglJTjcq/wdKRCYA/YFTVXUHsENEZmI/khMAVPXlkOtNEpE/YnmNQn8QQ1mqqi8E7+eIyHwgE/isiHE8ECRWRERuB5aFOPUjjgnLHyXA8SKyJihe9XkR1ymOQcAUVX0v+PyuiLwJXI3lXgPYBdyrqvuB5SIyPRjf+CivUQsoKFmgqioi/bCqjlWAG4DWwBlAdxF5ClMYK4BbgoJs+WzF8mU5FQCfmTjlmcbBdjMH0mSvCtm/ItiGptAOzXi6E3ta31io7SgAEakkIg+KyJLA7JONzWAyipCpcEbVHfnnK4LVYd43ppgxqepKTBleC6wTkU9F5BfFXKsommAJHUNZwcF/vw2BIsknmvGFsgU4OrRBVT9R1W6q2hl4FxiFjelKoJaqdsf+roXrqRzNoRUDnXKKKxOnPHMZ8D2whAPV4JqH7G8ZbEMrxZWEfljW2L7Yj1pNzDQWrmBQLDQP8z6LKMakqlNV9VzM/DUZeFtE0sNcI5zfoTBrscJQobSk9H+/cHyFmewi8WfgXVX9ClPcs4P2WVhVTwDESjy0Cc7nVABcmTjlDhFpEji2rwJuDepTrwM+AP4qIjXFqub9FXivtPUXsCffXMyuX0lErsZ+4OLNPSJSX6zu+mPAdFVdV9yYRKStiJwXKI99WJlZJbzi+AFzwBdWFqGMAfqKrQOpLCK9gF8Dr8RtpJbq/OxwO0SkA1bp8MGgaQVwbhAQ0QvzpeTTHUuf7sqkguDKxCkv3BNELG3FnlJbA6ep6pSQPgOAbcDi4JWNObdLy1jMV7AcmwG1w2q5xJtxwXnXYn6DASH7ihpTFWy9xvqg/Ragr6ruLnwBVV0KjAD+E0SNXRmmz2wsQOEJzBz1OFa7IxY/TGHeB3LFwn4LCCLCXgJuCKkVMhIb42ZMiT8ccsjVwPAIkV5OOcTrmThOghCR5pg/pImqZhXdO3UQkfOAoYEvpDTHt8Oi8DpU5CJVhxuuTBwnQRyuysQ5PHEzl+M4jhMzPjNxHMdxYsZnJo7jOE7MuDJxHMdxYsaVieM4jhMzrkwcx3GcmHFl4jiO48TM/wP8O/MmzTZLugAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "x0 = [38]\n", "x1 = [25]\n", "yw = mwarp1d.warp_landmark(y, x0, x1)\n", "plot_warp_results(y, yw, x0, x1)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Note that using a single landmark causes the entire domain to be warped.\n", "\n", "To restrict warping effects to a local area, add more landmarks, like this:" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZMAAAEJCAYAAABR4cpEAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOydd3xV5f34358MsncIEEJIIGHvJTvIEqFVcdZtrdpWu/3299U62trab5fVDrXFulqxbtQqjigkiAwJInskkAAhQMIKhOzc5/fHcy/GkHGTnHNH8rxfr/M6957znOf53OTc+znP81milMJgMBgMhs4Q4G0BDAaDweD/GGViMBgMhk5jlInBYDAYOo1RJgaDwWDoNEaZGAwGg6HTBHlbAG+RmJio0tLSvC2GwWAw+BUbN248ppTq2fR4t1UmaWlp5OXleVsMg8Fg8CtEZH9zx80yl8FgMBg6jVEmBoPBYOg0RpkYDAaDodN0W5uJwWDoftTV1VFcXEx1dbW3RfF5QkNDSUlJITg42K32RpkYDIZuQ3FxMVFRUaSlpSEi3hbHZ1FKcfz4cYqLi0lPT3frGp9Z5hKRZ0SkVES2tXBeROQvIlIgIltEZFyjczeLSL5zu9lzUnuf1UuXUvx4Go6lARQ/nsbqpUu9LZKhi9IV7rXq6moSEhKMImkDESEhIaFdMzifUSbAc8CCVs5fDGQ6tzuAJwFEJB74OXABMAn4uYjE2Sqpj7B66VLG1t5BStx+AkSRErefsbV3+OWX3ODbdKV7zSgS92jv38lnlIlSahVwopUmlwL/Upp1QKyI9AEuArKVUieUUieBbFpXSl2GtFP3ERFS+ZVjESGV9D95n5ckMnQ1zh7dyydPPUr/E/c0e6+lnTL3mkHjM8rEDfoCBxu9L3Yea+n4eYjIHSKSJyJ5ZWVltgnqCY4dKiM59kCz5/rGHuCmm2DvXg8LZega1J6idvMfOPr8OCI+zmBGxE/oG3+o2abJMc3fg4bWKS4u5tJLLyUzM5OBAwfywx/+kNra2vPalZSUcOWVV7bZ38KFCzl16lSHZPnFL37BH//4xw5d2xh/UibNzblUK8fPP6jUEqXUBKXUhJ49z8sG4DdsfDebkA8GUHIyudnzxSf78dprMHgw/P3vHhbO4L9Ul8Hm+2h4oz89tv8/Cvf3YMmGP7IptYiSU6nNXnLoZF+UwxTYaw9KKS6//HIuu+wy8vPz2bNnDxUVFdx331dnefX19SQnJ/Paa6+12efy5cuJjY21S2S38CdlUgz0a/Q+BShp5XiXZNN7KxhWdgklpwewpeEeztaEf+X82ZowDtRfzL59sGABfPe7sGSJl4Q1+AfKAXuegLcHoLb/H8s3zWfOHz+nJmsddzx6N2On96co9uFm7rVw9h/rx6pHf2AUSjtYsWIFoaGhfPOb3wQgMDCQRx99lGeeeYYnnniCq666iq9//evMnz+foqIiRowYAUBlZSVXX301o0aN4pprruGCCy44lxIqLS2NY8eOUVRUxNChQ7n99tsZPnw48+fPp6qqCoCnnnqKiRMnMnr0aK644goqKyubF7CD+JNr8NvA90TkJbSxvVwpdVhEPgB+08joPh+411tC2skXH+Yy+MjXKC7PIOGqjxjctyerl8aRdvI+kmMOUFKeSlHlDKYPeAOS/srrrwdz+eXw7W9DUBDcequ3P4HB5zhTAOtvg9JcKmPms/CBx9h+cCgrV4LzNwyA6ddfz+qlfPVei3mY+thNzOr1CLmPBTLzR48iAf5j3P7Rj+CLL6ztc8wYeOyx1tts376d8ePHf+VYdHQ0qamp1NfXs3btWrZs2UJ8fDxFRUXn2jzxxBPExcWxZcsWtm3bxpgxY5rtPz8/n//85z889dRTXH311bz++uvccMMNXH755dx+++0A3H///Tz99NN8//vf79TnbYzPKBMR+Q8wC0gUkWK0h1YwgFLq78ByYCFQAFQC33SeOyEivwI2OLt6SCnVmiHfL9nz2RYyihdRciad2Cs/JrGvXqabfv31wPWAnpKlFL8Fq16Aw9mE9F3I66/DZZfBbbdBXBwsXuy9z2DwMUo+gE8uh4BgKoY/zfBLvklFhZynSFycd68BynEduY/Vk9X7z+Q8Hs+s7z/oyU/glyilmvWUch2fN28e8fHx551fvXo1P/zhDwEYMWIEo0aNarb/9PT0c4pm/Pjx5xTStm3buP/++zl16hQVFRVcdNFFFn0ijc8oE6XUtW2cV8BdLZx7BnjGDrl8gbqaOhyf3kxlWCTRl31Mz5Sklhv3uRh6xEPRv6HvQkJDYdkymDkT7rgDpk2DpFYuN3QTDi6DT6+BmOGQ9V++fXsKJSWwZg208BvVLBIgzPzRo3z6SBnTez/Ezk8XMXTa+LYv9AHamkHYxfDhw3n99de/cuz06dMcPHiQwMBAIiIimr1O/wS2TUhIyLnXgYGB55a5brnlFt58801Gjx7Nc889R05OTsc+QAv4k82k2/LpU79lSK8vKEz8O0n9e7feOLAH9L8Git+EutMAhIXBc8/B6dPwve/ZL6/Bxyl8AVZfBXHjYc5K3vwwhRdfhAcegIkT29+dBAjDb/4bZRW9CNpwC7XV53slGb5kzpw5VFZW8q9//QuAhoYG7r77bm655RbCw8NbvG769Om88sorAOzYsYOtW7e2a9wzZ87Qp08f6urqWGpDfJBRJj7Ons+2MDX6V3x66DouuPIy9y5KuxEaquHgG+cODR8Ov/gFvPqq3gzdlKIXYe1NkJQFs7M5URHLd76j1/rv7YSlMTYpjuI+S8hM2saaJb+2Tt4uiIiwbNkyXn31VTIzMxk0aBChoaH85je/afW6O++8k7KyMkaNGsXvfvc7Ro0aRUxMjNvj/upXv+KCCy5g3rx5DBkypLMf4zzE3alTV2PChAnK14tj1dXUUfDEBSSGHyLwkh3E90lw70Kl4L+ZENEf5nx87nB9PUyeDAcOwPbt4Mfe0YaOUPIe5F4CPafBrPcgKIwbboCXX4YNG7RC6Syr/3Azk3svJX/gZwydOq7tCzzMzp07GTp0qLfF6BANDQ3U1dURGhrK3r17mTNnDnv27KFHjx62jdnc30tENiqlJjRta2YmPsyaf/2Dob02sTf+SfcVCYAIpN0AR1dCZfG5w0FBernr1Cm4/37r5TX4MGVr4ZMrIHYEzHwLgsL45BNYuhR+9jNrFAnAiJsf41hFEvVr7zLuwhZTWVnJ9OnTGT16NIsXL+bJJ5+0VZG0F6NMfJSqiioG1/+GL0pmcsEVHXDBSrseUHDwza8cHjFCG+KfeQYaeR0aujKn90DuIgjrC7Pehx56aeTBB6F3b7jnHuuGik2KY0/ILxjZZx1577xvXccGoqKiyMvLY/PmzWzZsoWLL77Y2yJ9BaNMfJTPlv6D3jGHYdRDHfPdj8qAsD5w/LPzTt17LwQGwq/N0nbXp65Cu/9KAMz+AMJ6AbByJeTk6FlJWJi1Q06+7hYOnkgjouBBMzvpRhhl4oOcLT/LMP6Pz0tmM2Z+Vsc6EYGESXB8/Xmn+vaF73xHL3mZ/F1dGKV0QOLpnTDtJYgccO7wgw/q+8AZw2YpPUJ7UBT1AMN657HhrXesH8Dgkxhl4oNsePFJekaVEjTml53rKGESnNkDtSfPO3XPPdCjB/zqV50bwuDD7P4zHHgZRv0aes89d/ijj2D1arjvPggNtWfoyd+4kf3HBxJdZGYn3QWjTHyMilMVjAj8HXmH5jNq7vTOdZYwSe+Pn++11ru3ztv173/Dnj2dG8bggxxbD5v+B1Iug2FfGkVcs5LUVHvT6wSHBHMw9kGG9PqC9W+82fYFBr/HKBMfY+PL/yQx8hihEzs5KwGId3rvNWM3Afjf/4WQEPj97zs/lMGHaKiGdbdAWDJMfk4veTrJyYF16/SspFGgtC1MvuY6Co8NIvpA6/ET3Y0f//jHPNYo/P6iiy7itttuO/f+7rvv5k9/+pPH5ImMjLSkH6NMfAhHg4PU2ifYengKI2ZN7nyHPWIhenCLyiQpCW64AV58EU6evxJm8Fe2/gJO74JJT53z3HLx+OOQkAA33WS/GEE9gjgQ/gOG9c5j+6rm78HuyNSpU1mzZg0ADoeDY8eOsX379nPn16xZw7Rp09rsRymFw+GwTc724rYyEZEBIvIzEXnc+X6wiAy3T7Tux6b3PyY9IZ/TvZtNQdYx4p1G+BaCU++8E6qqtDHe0AU4vgF2/gEGfguSv5rI79AhePNN+Na37LOVNGXs4hs5Ux3JyfVPeGZAP2DatGnnlMn27dsZMWIEUVFRnDx5kpqamnOBgnPmzGHcuHGMHDmSt956C+Bcivk777yTcePGcfDgQSIjI7n77rsZN24cc+bMwVX4b+/evSxYsIDx48czY8YMdu3aBUBhYSFTpkxh4sSJPPDAA5Z9LrcSPYrIPOANYCU6s+9dQCJwP7o2u8EC6nc8Tll0Tybc2HZlNbdJmKSTPlYWQ0S/806PGQNTp8ITT8APfwgBZq7qvzTU6OWt0D4w9pHzTi9ZAg6H9uTzFNEJ0aw6eSOTkp7heMkfSUhO9NzgbbHxR3DS4hz0cWNgfOsZJJOTkwkKCuLAgQOsWbOGKVOmcOjQIdauXUtMTAyjRo0iPDycZcuWER0dzbFjx5g8eTKXXHIJALt37+bZZ5/liSe0gj579izjxo3jkUce4aGHHuKXv/wlf/vb37jjjjv4+9//TmZmJuvXr+fOO+9kxYoV/PCHP+S73/0uN910E48//rhlH93dn47fAlcppS4BGpzHPgd8L1+Cn3JozwEm9Pkv26tvIyTcwsXsc0b4lpcZ7roLCgogO9u6YQ1eYMfvoHwHTFpy3vJWba1WJgsXQnq6Z8VKnnUXocE1bH2zyyb2bjeu2YlLmUyZMuXc+6lTp6KU4mc/+xmjRo1i7ty5HDp0iKNHjwLQv39/Jk/+chk8ICCAa665BoAbbriB1atXU1FRwZo1a7jqqqsYM2YM3/72tzl8+DAAn376Kddeq5O033jjjZZ9JndT0A9USrnCWRWAUqpKRIItk6Sbk//eP+idABkLvm1tx3GjISBYK5PUK5ptcsUV8OMf69mJxSUODJ7iTAFs/w2kXgN9F553+s034cgRvazpaTLGD+eLlVkMCP07DXV3Exgc6HkhmqONGYSduOwmW7duZcSIEfTr149HHnmE6Ohobr31VpYuXUpZWRkbN24kODiYtLQ0qqurAVpMUe9CRHA4HMTGxvJFC9W/mqun0lncnZkcFJGvlMsRkdFAkeUSdUNqKmsYHvoUeYe/Rsrg/tZ2HhgCsWNanZmEhOjiWe+8A/v3Wzu8wQMoBXnfg4AeMK55L6DHH9czkgULPCybk8qUO0mNL+Tzd02KFdAzk3feeYf4+HgCAwOJj4/n1KlTrF27lilTplBeXk5SUhLBwcGsXLmS/a18MR0Ox7k68S+++CLTp08nOjqa9PR0XnWmCFdKsXnz5nNjv/TSSwCWpqJ3V5n8BXhDRG4AAkXkCuAF4FGrBBGRBSKyW0QKROS8bEEi8qiIfOHc9ojIqUbnGhqde9sqmTzFxjffoGdUGUFDLTS8NyZhEpzIA0dDi02+7ZwQmXrxfsjB1+DwBzDqVxCefN7pHTtg1SodV+Qtm9jExYs5Ut4Hlf+kdwTwMUaOHHnOFtL4WExMDImJiVx//fXk5eUxYcIEli5d2mrK+IiIiHOlgFesWMGDD+pql0uXLuXpp59m9OjRDB8+/JwR/89//jOPP/44EydOpLy83LoPpZRyawPuALYCFcA24DZ3r3Wj70BgLzAA6AFsBoa10v77wDON3le0d8zx48crX+Gz31+sDv4tVTXUN9gzwN7nlVqKUie3tdrs4ouVSk1VqsEmMQw2UHtaqTeSlVo+RqmGumab3HOPUoGBSh096mHZmrDy0XtU3b8CVelB7wmyY8cOr41tFxEREbb13dzfC8hTzfymuv2copRaopQaqZSKVEqNUEr90zKNBpOAAqXUPqVULfAScGkr7a8F/mPh+F6jrLiUsb0/pKDuOgICbXpsdMMID3D99brWyaef2iOGwQa2/x9UlcDEv0PA+SZQh0PHEc2f7/1yzX2nX09QYAM73nvFu4IYbMGtXy8RSW5ps0iOvsDBRu+Lnceak6U/kA6saHQ4VETyRGSdiLRYjlBE7nC2y3P5YnubHe+9TFBgAykzbrBvkOhBEBgK5dtabXbZZRARAS+8YJ8oBgs5exB2P6rLDSRe0GyT1av1A8INNt5e7pI5YQS7j44m4Yy5waykoqLC2yIA7ttMitE/9s1tVtCca0FL2eG+AbymlGpsAEhVuvLXdcBjIjKwuQuds6sJSqkJPX2kzGDCmRfYfXQ0GeNtjP+UAIgarKOiWyEiQiuUV16Bmhr7xDFYxJb7tfF99MMtNnnhBf1/vbS1eb4HORx6PSN6r6doa77XZFDdtLpse2nv38ldZZKOtme4thnAB8At7RqtZYqBxhF1KUBJC22/QZMlLqVUiXO/D8gBxlokl60Ubc1nRO/POBzqgcfG6CFtKhPQT7CnTsF779kvkqETnPwCCv8Ng3+gyzM3Q00NvPoqLF6sFYovMHj+tTgcQlHOi14ZPzQ0lOPHjxuF0gZKKY4fP05oO1IluBVnopRq6pe2X0RuRi81/dt9EVtkA5ApIunAIbTCuK5pIxEZDMQBaxsdiwMqlVI1IpIITAP8InVhUc5SUuOEwRdda/9g0UPgwCs6CWBgyzfI3Ll6bf2FF/QsxeCDKAWbfgo94mD4z1pstny5fjC4/noPytYGfQamsGnZLNJCXkA5HuxY4bdOkJKSQnFxMb6yzO3LhIaGkpKS4nZ7d4MWm6MSsCQoQilVLyLfQ892AtGeWttF5CG054DL3fda4CX11ceKocA/RMSBnmn9Vim1wwq57EQ5FGks5YsjFzJuQLPmIWuJHgIoOJMPsSNbbBYUBN/4BvzjH/qHKDbWftEM7eRINhz5CMY9qpN5tsDSpfrBYO7cFpt4hYqeNzA2+FtsX72B4TMneXTs4OBg0j2dAqCb4G5urqazhAj0D/v5Zfw6iFJqObC8ybEHm7z/RTPXrQFa/nX0UbZ/8hkjEgo4WHevZwaMcfqpn97VqjIB/ST7l7/A66/rpIAGH0Ip2PKAXtrK/G6LzU6dgv/+V+fhCurMI6MNjFp0BdXv3smxDS+Ah5WJwT7cvc2aWvgqgDx0okdDBzi28WVqEnowalHzKU4sJ2qQ3pe3bTeZOBEyMuDll40y8TlK3tMu3pOW6OwGLfDWWzof13XnLRZ7n5jEGNaWfo0hUa/iaHjMPpd4g0dx67+olEpvso1USn1TKXXIbgG7IsqhGNhjGVtK5xKTGNP2BVYQFK6fZt0wwovA5ZfDypWmzolPoRRs/TlEpMGAW1ptumwZpKTAJB998Hf0vZxe0UfYnrvO26IYLMI8EniBPZ9tpl98EdWJiz07sJseXaCVSX09vPuuzTIZ3KfkXZ0WZ8T9OnlnC5w9Cx98oL24bMjnZwkj5i+itj6Y45uXeVsUg0W0uMwlItm0HOtxDqXUfEsl6gYc/uwNMuIDGDLvEs8OHD0ESj8B5dCxJ60wcSIkJ8Mbb/hGwFu3RyldQTFyAKS3Xibx/fehulorE18lJjGGDUfnkB72Bsrxe497dRmspzWbyWqPSdHN6OtYxtYj0xmT4uH8FtFDoKESKg81WyirMQEB2jX42WehshLCwz0ko6F5Dr0DJzbCBc+0OisBvcSVkAAzZnhItg5SlbCY/qHfZk/eVgZNGuVtcQydpEVlopT6pScF6S4Ubc0nM2kbuWcsS7jsPtGNPLraUCagl7qeeAI+/NDEnHgVpWD7ryEiHdJbL2ZUW6tLCVx+ue95cTVl6NxLcaz6DiXrlxll0gVol81ERMJEpJ+IpLo2uwTrqhR9oteIB2Z54de5sTJxg5kzIS5OL3UZvMjRj7UH1/B7mk3m2JiVK6G83LeXuFz0TO3FtqNT6dNgbrCugLuJHgeIyBq0S3ARUNhoM7SDhKpl7Dw6jpQhaZ4fPLQXBMe4rUyCg+HrX9fxCnV1NstmaJltD0NYX0i/uc2my5bp1Cnz5nlALgs4EXE5g5O2cGDHPm+LYugk7s5M/oZO6jgaOAOMAt4ETBRCOzhSWMLIPus4Guylx0aRdnl0gV4uOXUKcnNtlMvQMmVroDQHhv5Pq3EloNPNv/WWrvPejpRKXmVglv4u7FtlvLr8HXeVyQXoYljbAJRS24FvAz+1S7CuyJ4VutJZyhQvrkG0U5nMn6+N78vMd907bH8YQhIh4/Y2m65bp+u8+8MSl4t+Q9PZfXQ08WfNDebvuKtMHECV83WFiMQCJwBjM2kHoSfeZf+JAQwcO8x7QkQP0cWU6k671TwsTCuUd97RdmCDBzmxCUqWw+AfQVDbaX/feQcCA+Hiiz0gm4UcDryU4b3WcuLwcW+LYugE7iqT7ehsvKDzcT2KrgtvbCZuUlVRxYjEFRTVLPKuT/05I/xuty9ZuFAXWNrh8+kzuxg7/wBBUTDoLreaL18O06f7X3LOxNELCQxwsHPlB94WxdAJ3FUmP0DPREAvbfUFJqCXugxusG1FDuEhVURkLvSuIO306IIvn3SXL2+9ncFCzh7UJQMy7mg1M7CLQ4dg82at+P2NoVMncKwiEQ6ZG8yfcTc31xal1Fbn631KqflKqclKqbVtXWvQVOYvp7ImjOGzs7wrSNRAkMB2zUxSUmDUKKNMPMqev+j94O+71dxVzMwflUlgcCC7yhcwJOZ9Guoa2r7A4JO46xr8kYhcJyJ+4iPiWyiHIi1kOduOzSEsMsy7wgQEQ3gqVLRvhXLRIl1PvLzcJrkMX1J3BgqWQL8rW6yi2JR334XUVBhuY/VnOwlIWURC5HF2frrB26IYOoi7y1yfAr8GDovIkyIy0WpBRGSBiOwWkQIRuaeZ87eISJmIfOHcbmt07mYRyXdubTvje5jCrXvoH7+P6ngfeWyMTIeK9vn1L1yoEz9mZ9skk+FL9j6tHSSG/MSt5jU18NFH+n/kq4kd22LohfNpcARwbIuZ/vor7i5z/VwpNQC4Al0Ya6WIbBORH1shhIgEAo8DFwPDgGtFpDmXp5eVUmOc2z+d18YDP0e7L08Cfu4s5eszHFirU+8OmOYjbjaRA+Bs+2Ymkydrw65Z6rIZRz3s/jP0nA6J7uWPX70aKir8c4nLRVzveLYfnUKvBpOm2l9pVzoVpdQKpdRN6HK9+4E/WiTHJKDAaY+pBV4CLnXz2ouAbKXUCaXUSSAbWGCRXJYQXbGcgrJh3ol6b47IdKg+CvVn3b4kKAguukivzTscNsrW3Sl+E84WuT0rAa3ge/SA2bPtE8sTnAhdyNBen3O06LC3RTF0gPbm5hotIo8BO9Glcv/PIjn6oiPsXRQ7jzXlChHZIiKviYgrU6G713qFMyfPMCJpFcUOH3psjBig9+20myxcqIPiNm2yQSaDZvef9cyxr/vlCd59F2bN0mlU/Jk+4/V3ZE/u+16WxNAR3DXAf19EPgfWAr2AG4D+SimryvY2t9LbNETuv0CaUmoU8BHwfDuu1Q1F7hCRPBHJKysr67Cw7WHHxx/TI6iO2GGLPDKeW0R2TJksWKDX5M1Sl02c2ARlqyHzLggIdOuSvXth927tIOHvDJo0msPlyQSVmhvMH3F3ZnIT8BSQrJS6Vin1oVKWxkMXA41zoqcAJY0bKKWOK6VqnG+fAsa7e22jPpYopSYopSb07NnTEsHboqboPU5XRTF81rS2G3uKyHS9b6cRPikJJkz40g3VYDF7/gqB4TDwVrcved/5EO9vUe/NIQFCwdmLGRb/IfW19d4Wx9BO3DXAT1RKPamUOmWTHBuATBFJF5EewDeAtxs3EJE+jd5egl5qA/gAmC8icU7D+3znMZ8gLTSbnScuJDik9YJGHiUkEYIi222EB51a5bPPjIuw5VQfg6IXdRVFN4IUXWRnQ1oaZGTYJ5onCUyZT0zYaeMi7If4RA14pVQ98D20EtgJvKKU2i4iD4mIa/H4ByKyXUQ2oyPyb3FeewL4FVohbQAech7zOvu37yU1vpDqOB+rbCzSIfdg0KnNGxp03QyDhex9Chw1MOh7bl9SV6f/D/Pm+a9LcFOGZM3B4RCObzM+6P6Gz9RiU0otB5Y3OfZgo9f3Ave2cO0zwDO2CtgB9q/7kP5hkDrJB4tLRA6AM3vbfdmUKdrQm51tqi9ahqMe8p+EXnMg1v2ow88+g9On9WyxqxDfJ4EdpeOJVx8CD7bZ3uA7+MTMpKvS40Q2xSdTSRuR6W1RzifCOTNpp+mrRw/tOWSCFy2k+C2oPOh26hQX2dl6RuLvLsFNKQ2cx7CkdZw+7l5ma4NvYJSJTdTX1jM0bgWFVfO8myW4JSIHQEMlVJe2+9J58yA/H4qKrBerW5L/uE6bkvy1dl2Wna0dIuLjbZLLS8QMnkdQYAM7c3O8LYqhHbitTERkuogsEZH/Ot+PF5GZ9onm3+xak0dMeDmBKT66BuFyD+6AEd5VEtbMTiygfCccXQmZ33XbHRi0A8T69f5Tnrc9DJs5lbM14VQXmhvMn3A3zuQ64C2gGnApEAU8ZJNcfs+xrR/icAhDsuZ4W5Tm6aB7MMDQoZCcbJSJJeQ/CQE9YID77sCgDe8NDV1TmYSEh7DjWBapwR96WxRDO3B3ZnIfcJFS6gfoqosA2wA/zVFqP3G12ewqHUd8nwRvi9I8EWl6387ARdDr9PPnw8cf6x80Qwepq4DC5yH1aghtX9xTdrZ2hJgyxSbZvMzZ6HmkJ+7h0J4D3hbF4CbuKpNkpVSe87XLYlsPuD8v70acOXmGYT3XURroo0tcAEHhENq7QzMT0E/EJ06Y1Cqdomipzg6ceWe7L83OhqwsCAmxQS4fIGWC/u7sW2Omv/6Cu8pkr4hMbXJsKuB+haVuxM6cHIKD6okZ7ONrEJEDOqxM5s7Ve7PU1UGUgvwnIG4MJE5u16X792sHiK64xOVi4NhhHC5PJrDM3GD+grvK5NfAWyJyPxAsIncD/8HYTJqlal82lTVhDJvZVP/6GB1IRe8iKQlGj1EAj3MAACAASURBVDbKpMMcWwOntuhZSTsjDl1/866sTCRA2Fsxl8ExH+FoMGmq/QF306m8CVyHrhmyH5gN3KqUMlmamiE5aAU7js0gJNzH1yAi03V8g6OuQ5fPnQtr1kBVlcVydQfyn4TgaEi7rt2XrlgBvXvDsOYq/nQles8hIfI4+XlbvS2JwQ3cdg1WSmUrpb6ulBqhlFqklDLPpM1QduAomT23UxHhB5FkkQNAOeBsx4ycF16oq/ytXWuxXF2dmhNw4DVIuxGC2pc3XimtTC68sOukUGmJAVMuBODwFyu8LInBHVpUJiKS7M7mSWH9gfxPddKqpJF+oEwiOu4eDDBjBgQG6h83Qzso/LfOw5Vxe7sv3bULjh7telHvzZGc0Y/C45mEnzY3mD/QWm6uYlqoC+JEnOeNR1cj6g+toDw2hkEXjPW2KG1zrq5Jx5RJdDRMnGiUSbtQCvYugYRJEDe63Ze7/tbdQZkAHKyZzZj4F6mvrSeoh8+kEjQ0Q2vLXOnAgFY213lDI1JDVrLrZJZ/3PhhyRAQ3GEjPOgftQ0b4MwZC+XqyhxbC+U7YGD7ZyWglUn//pCebrFcPkpQymyiw86we93n3hbF0AYtKhOl1H53Nk8K6+sc2nOAtIQCqqMv9LYo7hEQCOH94GzH/42zZ0N9PaxebaFcXZmCJbqWTP9vtPtSh0NHvs+e3fXtJS4yp88CoGybmf76Ou3JzXWNiLwvItuc+/Z/G7o4heu0vaTPWD9ag4jo3yllMnWqziRslrrcoPYUHHhFe3AFR7b78s2b4eTJ7rPEBdAzJYn80hFEnTU3mK/jbm6u/wX+gi4+9ahz/5jzuMHF0RUcq0gkY/wIb0viPp1UJmFhWqEYZeIGRUuhoQoy7ujQ5a6CZBf6ycTXKkocsxmauJqaypq2Gxu8hrszk7uAi5VSDyilnlZKPQBcjK6OaAkiskBEdotIgYjc08z5n4jIDhHZIiIfi0j/RucaROQL5/Z202s9gXIoBkSsIL/8QgIC/Sizf3h/qDoMDbUd7mL2bJ1W5YRP1Lf0UZSCgqd0xHv8+A51sWIFDB4MfftaLJuPE9p/NuEhVez6dL23RTG0gru/euHAlibHtgJhVgghIoHA42gFNQy4VkSahmRtAiYopUYBrwG/b3SuSik1xrldghfYv72A5Nhi6uL9bA0ioj+goKq4w13Mnq1/K3NzrROry3Hyczi1ucOG97o6/fftTktcLobMzKLBEcDJXWb668u4q0z+Bfy4ybEfAc9bJMckoEAptU8pVQu8BFzauIFSaqVSqtL5dh2QYtHYlnBgg77R+030s297RKred2Kpa+JEncHWLHW1QsE/ITC0QxHvABs3QkVF91QmMT1j2V06jrgac4P5Mi36r4pINl/GmQQA3xORu9DpVPoDfYBPLJKjL3Cw0ftidOqWlvgW0DiVS6iI5KEzGf/Wmf7lPETkDuAOgNTU1E4J3JSg4ys5HJnsmyV6WyPCuVrYCWXSowdMn26USYvUV8L+F6HfldAjtkNduOwlWVkWyuVHlMqFTE16jMrTlYRHh3tbHEMztBYM0dTZ0yrF0RzNOTo2GzApIjcAE4DGX6tUpVSJiAwAVojIVqXU3vM6VGoJsARgwoQJ7St+3grKociIyiH/zBz6+GKJ3tYI76f3nVAmoOvC33svlJbqJJCGRhx4TaeaH3hbh7vIyYERI6Bn+8qedBkiBsyix9k/sG31WsYt9NGCc92cFpWJUuqXHpSjGOjX6H0KUNK0kYjMRRfqylJKnXPtUEqVOPf7RCQHGAucp0zsonDLbgZEH2VXyCxPDWkdgSEQ1qfD+blczJql96tWwZVXdl6sLsW+pyEyA5I6VuW6rk7H8dzavmKMXYrB06fT8F4ApwtyAKNMfJH2xJkEisgQEZkhIjNdm0VybAAyRSRdRHoA3wC+4pUlImOBfwCXKKVKGx2PE5EQ5+tEYBqwwyK53KJ4Yw4AqRNmeXJY6whP7fTMZPx4bTfJybFGpC7D6T1QugoybutwpGFeHlRWfqmwuyPRCdHsKh1PXF2Ot0UxtIBbOT9EZBzwBpCKXn5y5eVqAHp0VgilVL2IfA/4AJ3r6xml1HYReQjIU0q9DfwBiAReFf2lPOD03BoK/ENEHGjl+FullEeVSdCJXA5HJtN/eIYnh7WOiP5wonPpKoKDtd3EKJMm7HsGJBDSb+5wF66/6UyrHt38lDKZxZSEPxu7iY/i7szkMWAZEAOcBqLRs4RbrBJEKbVcKTVIKTVQKfWw89iDTkWCUmquUqpXUxdgpdQapdRIpdRo5/5pq2RyS26HYmBUDvsqZiH+Zi9xEdEfKg/odPSdYNYs2L4dysqsEcvvcdTBvucgeRGE9e5wN7m53dte4iJiwCxCgmvZ/ek6b4tiaAZ3lclI4B6l1BlAlFIVwP/DVFqkcOseekUfoSHBj91sIvqDoxaqj3aqG5en0apVFsjUFSh5X/9NB36rw1247CXd1YurMYOnT6fBEUB5fo63RTE0g7vKpHEpvnIRSXIe6/jjVhfBZS/p56/2EtA2E+i0EX7CBAgPN0td59j3NIT2guSLO9zFxo1w9mz3tpe4iE6I1vEmtTneFsXQDO4qk42Aq+J0DvBvdGBh06j4bkfg8VyOlPfxv/iSxlgQawLGbvIVqo7AoXe0rSQguMPdGHvJVymVWQzpuZ6qClMr2tdwV5ncBmx2vv4JcACoAb5ph1D+gnIoBkbmsNef7SVgmTIB/QS9bZuxm1D0AqgGGNC5r0hODgwfbmJ3XESka7vJrtXGbuJruKVMlFKHlFIHna+PK6VuV0p9Qym1217xfJuibfn0jjlMQ8Isb4vSOXrEQHCMZcoEurndRCnY+zQkToWYIR3uxmUvMUtcXzLIZTfZk+NtUQxNaC2dygVKqfXO11NbaqeUWmOHYP7Awbwc0kMhZXwXsI52MhW9C5fdJDcXrrjCArn8kWPr4PQuuOCfnerm88+1vcQY378kJjGGHaVjiXXkeFsUQxNaizP5CIhyvm6pjl63rgEfeDyXo1G9SR85yNuidJ7wVO0e3EmCg2HatC9zSXVL9j0DQRGQenWnunHZS4wy+SqlMospPf9KVUUVYZGWJC43WEBrZXujGr0OaGHrtopE1y/JZe+ZLP+2l7iwaGYC+sdv2zY4ftyS7vyL+rOw/yVIvQqCo9pu3wq5uTB0qLGXNCU8LUvHm6z5zNuiGBrRps1ERIJE5F0RCfWEQP7CgZ376BN7iPr4LvLYGNEf6sqhtrzTXbmepD+xMzWor3LwDaiv6LThvb7exJe0xODpM3A4hFO7TQEdX6JNZaKUqgfGo9O7G5wcyNM3cvLYLvJtt9Cja+JECA3tpsWy9j0LkQOh54xOdfPFF3DmjFEmzRHTM5Y9ZaOJrumON5jv4q5r8L+xsERvV0BKczlWkcjAMUO9LYo1WFAky0VICEye3A2VSUUhHF0JA27pcFJHF66/nYkvaZ6jaiZDEtZSW93xctMGa3FXmYwDfu+sz/6RiHzo2uwUzpdJC8+loHxm17CXwJczEwuM8KCfqL/4Aso7v2rmP+x7HhBIv6nTXa1aBRkZkJzcebG6IiH9sggPqWL3mjxvi2Jw4q4yWQU8jJ6hfAJ82mjrdhTv3k9K3H5qYrrQGkRoLwjoYakRXim97t8tUA4ofB56z/lyltdBHA5tbzJLXC2TOU1P2Y7v7G7TX9/FrRT0Hi6U5fMUbVhFSgD0Ht2Fvu0S4KxrYs3MZPJk7SacmwuLFlnSpW9Tmgtni2D0w53uautWOHnSKJPWSEhOJL9sOJG1ucC93hbHgJvKBHRxLCAT6EmjMrtKqW4X66yO5HIqNpbMCSO9LYq1RHS+SJaLsDCYNKkb2U32PquzCKQs7nRXxl7iHofrZzI24d/U19YT1MPtnzKDTbi1zOUsjrUXXcEwx7mtRAc2djv6heSy++QMAgLdLlTpH0RYE7joIitLZ72tqLCsS9+k7jQcfA36XwNBnQ+iW7UK+vfXm6FlgpKziAqtYM/6Td4WxYAPFccSkQUisttp5L+nmfMhIvKy8/x6EUlrdO5e5/HdInKRVTI1x5HCEtISCqiK6oJrEOH9oeowNFjjIZOVBQ0NsKarJ9w58Co0VHU6tgS0nWnVKrPE5Q4ZU/UfqXRbd5n++jY+URzLuYT2OHAxMAy4VkSGNWn2LeCkUioDeBT4nfPaYeia8cOBBcATzv5sYd86vaqXNKILftsjnFWZqw5Z0t3UqRAY2A2WuvY9C9FDIOGCTne1c6fOuGyUSdsk9e9N4bFBhFd09RvMP/CV4liTgAKl1D6lVC26VsqlTdpcCjzvfP0aMEd0MfhLgZeUUjVKqUKgwNmfLdSX5HK6KopBk8bYNYT3sDDWBCAyUid+7NLK5PQeKPtUz0o6GVsCX/6tjDJxj4O1WQyO+4SGugZvi+IX7Fidx9o/XMmBHfss79tXimP1BQ42el/sPNZsG2dUfjmQ4Oa1AIjIHSKSJyJ5ZZ0ouLH9xIKuafCzqOJiY2bOhM8+g6quWsto33PaEy79Rku6W7VKx5YMGGBJd12ewN4ziQkvp2DjVm+L4heUbs5mSt/XiYiNtrzvjhbH2o+1xbGae6RTbrZx51p9UKklSqkJSqkJPXv2bKeImpl3P8mUn77SoWt9nvB+em+xEb6uDtZ1xVpGjgYo/Bf0WQBhfTrdnVJ6ZpKVZckkp1uQPllP4Y5s6crTX+uIrMwlv2w4CcmJlvfd0eJYd1hcHKsY6NfofQpQ0lIbEQlCOwOccPNagzsEhUFokqUzk+nTISCgiy51HflI25csMLwDFBTA4cNmias9JGf048CJdELKu+INZi31tfUMTviUw/X2+Jy76xpcICI/E5Fml48sYAOQKSLpItIDbVB/u0mbt4Gbna+vBFYopZTz+Dec3l7p6FgYk5u6o4RbF2sCEBMDY8Z0UWWy71noEQ99v25Jd8Ze0jH2V2WRGbMKR4PD26L4NHvWbyIqtIKgZHtuMHeXuR4G5gOFIvKeiFwlIsFWCeG0gXwP+ADYCbyilNouIg+JyCXOZk8DCSJSgF5qu8d57XbgFXQMzPvAXUopY43rKBH9LV3mAm03WbcOamos7da71J6E4jch7ToIDLGky1WrdO2SwYMt6a7bIL2ySIg8zr4vdnpbFJ/G5ULtcqm2GneXuZ5VSs0ChgJ5wB+AwyLyZ6sEUUotV0oNUkoNVEo97Dz2oFLqbefraqXUVUqpDKXUJKXUvkbXPuy8brBS6j2rZOqWuFKqqGbNTh0iKwuqq2HDBsu69D5F/wFHjWVLXKBnJjNnGntJe0mdoJdtDm3qitNf6wivyKXw2CCS+lvlhPtV2hXCrZTaq5R6AJgMrMekpe96RKRCQyXUnrCsyxnO0h5daqlr33MQOwrixlrSXVERHDhglrg6Qr8h6ZScSiH4ZFe6wayloa6BwXGfcLDWvhvMbWUiIoEicqmIvAkUAbHAt+0SzOAlLCyS5SIhAUaO7ELK5NR2OLHBstgSMPaSziABQuHZLDKiclEO62bUXYmCjVuJCS8nsLd9Cd/cNcD/Ce0h9Te0TWOUUmqaUuqftklm8A4R1seagF6+WbNGuwn7PfueBQmCtOst6zI3F+LjYfhwy7rsVjQkZpEUfZTCrXu8LYpP4nKddrlS24G7M5N+6Dxc/ZVS9yqlzH+sq+IKXLTYCJ+VBWfPwuefW9qt53HUQdG/oe/XILRjsUrNsWqVXg4M6GK5Qz1Fv/H6R7L4826XxNwtQspXceBEOskZ/dpu3EHcNcBfpZR6TyllfO+6OiGJEBhmy8wEusBSV8n7UF1qqeH90CHYu9cscXWGtBGZHD3dm8Bj/n6DWY9yKDJiVrG/yt4bzDwHGb6KiKV1TVz06gVDhnQBZbLvWR3YmXyxZV0ae0nnkQBh75mZpEcYu0lT9m7aQWLkMaSXUSYGTxNubV0TFzNn6jK+Df4aBVRdCof+C2k3QoBlYVbk5kJ0NIwebVmX3ZK6uCySY4s5uKvQ26L4FC6XaZcLtV0YZWI4nwjryvc2ZtYsOH0avvjC8q49Q+ELoOph4K2Wdpubq+0lgbYVTuge9B03C4D9G/x9+mstwSdzOXSqH/2GpNs6jrveXM2qNBGZYa04Bp8gvD9UH4EGa0PWXcs4OTmWdusZlIJ9T0PCZIhpWmqn4xw+DLt3a0Vr6BwDxwyl7ExPpDTH26L4DMqhyIzOofDsLCTA3mhYd2cm77Rw/C2rBDH4EC734MqDrbdrJ8nJMGiQnyqT4xugfIctsxIwysQKJEAoOD2LtPAcYzdxsm/zLnpGlaKSZtk+lrvK5DyVJiJRgPHu6orYFGsC+kdz1So/tJvsexoCw3WddwvJydH2kjFdsNaaN6iNyyIl7gDFu4u8LYpPcOjzHABSJ9jv3dGqMhGRfBHZA4SJyJ7GG7og1UrbJTR4HptiTcBP7Sb1lToXV+pVEGxtUaGcHG0vCeqCtda8Qd+xswAo+izHq3L4CkEncig5lULqUPurrbV1C/8aPSt5Ep052IUDOAKssEkugzcJTwHEcvdg+NJukpsL48db3r09HHwd6s9YvsR15Ii2l3zrW5Z2260ZOHYYxzYmIuW5WFe7zz9RDkVGVC57zswj2WZ7CbShTJRSzwOIyC6lVFeslWdojsAQCOttizJJTobMTP1E/pOfWN69Pex9GiIzoKe1/ibGXmI9EiDkl2u7SXencMtuBkQfZXfoLI+M524E/DoRCRWRkSIytfFmt4AGLxGRZosyAT+zm5zOh9JcGGhdUkcXOTkQFQVjrUk8bHBSG5tFStx+incVeVsUr1K8MQfwjL0E3HcNvgQ4jK4Dv7rR9ol9ohm8is3KpLwcNm+2pXtr2fcMSACk32J518ZeYg/Jxm4CaHvJ4VN9SR020CPjuevN9QjwSyBSKRXQaOt0mJWIxItIttPYny0icc20GSMia0Vku4hsEZFrGp17TkQKReQL52b8YqwgIk0b4B3WTx/8Jt7EUa/rliQvgvBkS7s+cgR27TJLXHYwcOwwjlUkwtEcb4viNbS9JId9HogvceGuMumllHpMKVVpgwz3AB8rpTKBj53vm1IJ3KSUGg4sAB4TkdhG53+qlBrj3PzJT8h3ieivM+RWH7a86759td3E5/N0lSzXwZsDrbeQr3ImtzXKxHoCAgPIL8+if5iv32D2Ubh1D0nRR3EkzvLYmO4qkw9FZLJNMlwKPO98/TxwWdMGSqk9Sql85+sSoBSwLv+34Xwi0vS+osiW7mfN0srEp+0me/8Job0heaHlXRt7ib3Uxs6iX3xRt7WbuOwl/TxkLwH3lUkR8LaI/FVEftZ4s0CGXkqpwwDOfVJrjUVkEtAD2Nvo8MPO5a9HRSSklWvvEJE8EckrKyuzQPQuTGSa3p8tsqX7Cy/UdhOfjTepLNEzkwE3W5rU0cWKFTrxpbGX2EPfcRcCULi+e4bCBR1fQcmpFPoPz/DYmO4qk/HAdmAEMK/RNtedi0XkIxHZ1sx2aXuEFZE+wL+BbzaqrXIvMASYCMQD/9vS9UqpJUqpCUqpCT17molNq7gCF21SJq7lnRW+GqlU+DyoBhhg/RLXoUM6vuTCCy3v2uBk4NhhlJ1JQkp99QazD0eDg8zolew7O9tj9hJoO2gRAKVUp257pVSLSkdEjopIH6XUYaeyKG2hXTTwLnB/45gX16wGqBGRZ4H/6YysBidBYXqJxyZl0qcPDB2qlclPf2rLEB1HOXRsSVIWRGda3v1K58Py7NmWd21wIgFC/ukLGRi5AuVQHv1R9TYFn29nUFQZu8M8+7Tidgp6EQl0xpZc43wfLiJhFsjwNnCz8/XNNJM8UkR6AMuAfymlXm1yro9zL2h7yzYLZDKAtpvYZDMB/WP6ySc+WBf+aA5U7LXF8A5agcbFmfoldlOfMJs+MSXdri784U16NpZ2gQ8qExEZiP6RXg487Tw8H3jKAhl+C8wTkXz00tlvnWNOEJF/OttcDcwEbmnGBXipiGwFtgKJ6BQwBiuITLNtZgJamZw9Cxs22DZExyhYAsGx0O9KW7pfsUIvcZl67/bSf5Ke+hXnda+lrtDyFew/PpCUwf09Oq67t/NfgZfQNgnXc2QO0On8Ekqp40qpOUqpTOf+hPN4nlLqNufrF5RSwY3cf8+5ACulZiulRiqlRiilblBKVXRWJoMTG2NNQMebiPiY3aS6DIqXQfpNeqnPYgoLYf9+s8TlCVKHDeTQqX4En+g+RviGugaGxOWyv8bzN5i7ymQS8LDT6K0AlFKngNhWrzL4NzbGmgAkJOilHp9SJoX/AkctZNxuS/euz2qM7/YjAULh2dlkxqzE0dA9qmXsXr+JmPBygvr6rjI5TRPFISLJwFHLJTL4DjbHmoB+Ql+zBqqrbRvCfZSCvU9B4hSIHWHLECtWQK9e2vnA4AF6zyYx8hgFG7uHKbV0i35ayZzm+acVd5XJG8AzIpICICIJwGPopS9DV8WlTGzK0QVamdTUwNq1tg3hPmWr4fRuyLjDlu6V0spk9mzLc0YaWmDAFP2j6jJKd3Uiz64gv2w4PVN7eXxsd5XJA0AFcAA9QykFaoDf2CSXwReIcBrwbDTCz5gBgYE+stRVsEQXv0q9ypbud+3SObmMvcRzJGf0o/B4JqGnfeEGs5fa6lqGJn5CSb13bjB3U9BXKaWuQ0enTwJ6K6VuVEr5wuKEwS6CwiC0l63KJDoaJk70AWVScxwOvgZp10NQhC1DuD6jUSae5WDNbIbE51JfW+9tUWxl1+rPiAipJLS/DysTF0qpY04vK5OLpLsQkWarMgH94/rZZ3DmjK3DtM6+56ChGjK+Y9sQH38MqamQnm7bEIZmCEqZTUzYaXat3ehtUWzlxM6PcTiEQdNnemX8FpWJM637vrY2Twpr8AI2By4CzJ0L9fVeTEmvHJD/JPScDnGjbBmivl7PTObNM/YSTzN45mwcDuHY1mxvi2Ir8bXZ7CydQFzveK+M39rM5H60reQB4FkgDHgRXQv+P0Ao8IzdAhq8TGQaVO7XP7g2MXUqhIdDtre+64ezdcR75p22DZGXpxNbzptn2xCGFkhITmR32Vhiq7uuMjl9/DTDktZRFui9G6zF3FxKqaWu1yKyCvi6Uiqv0bHX0R5dJuK8KxORpmNNqg5DeF9bhggJ0Rl0P/zQlu7bJv8JCE2CflfYNkR2tp6RzJlj2xCGVjgq85ia9CfOnDxDVFyUt8WxnJ25OVwQ2EDsEO8pE3dtJmOApsnCtziPG7oyHvDoAv3Evns3HDxo6zDnc3Y/lLwDA2+HwB62DfPhh7p2SWKibUMYWiF60Dx6BNWxM6drFsyqLvyQszXhDJ0xxWsyuKtMdgM/bnLsR0D3yqDWHfFA4CLA/Pl67/GlroIlem9TbAlox4J16778jAbPMyxrGpU1YVQVds2lrtTgbHYcm0VIeIvlnGzHXWVyF/C/IlIkIrkiUoQur3uXbZIZfAMPzUyGD9dp6T2qTBpqdDXFvl+HiFTbhsnJ0QZ4Yy/xHqERoew4NpOUwK6nTA7tOUB64h4qo717g7kbZ/IZMAC4D11T5D5goFJqvY2yGXyBoHBtT7BZmYhor66PPgKHp9IoHXgFqkttNbyDVpBhYTBtmq3DGNqgImoeA3vu5PDeYm+LYin71mgF2XeCHygTAKXUaaXUUqXU7537cjsFM/gQEQOgwn4v8Hnz4NgxD5by3f1XiB4Cve39EmZnaweDEO+tQBiA5HH6/7z304+8LIm1BJZlc7g8mYFjh3lVDnfrmYiIXCsij4jIksab3QIafICoDDhTYPswc531OD2y1HVsPZzYAIO+Z2vgx8GDOo2KsZd4n8wJIyk93Qs56i23QetxNDgYHPMReyvmeb2apLszkyeBvwH9gOAmW6cQkXgRyRaRfOc+roV2DY0KY73d6Hi6iKx3Xv+ysyqjwUqiMqDyoI4Qt5E+fWDkSA+5CO/+CwRF6bolNuJSjMZe4n0kQNhzZh6Doj/qMinpd6/bRELkcSTZ+zeYu8rkKmCSUupqpdQ3G28WyHAP8LFSKhP42Pm+OaoaFca6pNHx3wGPOq8/CdhTa7U7E5UJKI8tda1eDZWVNg5SdQQOvgoDb4Vge2MOsrOhd28YYU9Ge0N76T2PnlFl7Plss7clsYSjm/XTyqAZc70sifvKpBKdMdgOLgWed75+Hl3H3S2cdd9nA6915HqDm0Rm6L0HlroWLIDaWlhpZ3G8gn/oQMxMe50RGxrggw/gootMChVfYdBMvd545PP3vSyJNcRWvsfOo2O9knK+Ke4qk98DDzp/vK2ml1LqMIBzn9RCu1ARyRORdSLiUhgJwCmllCsdaDFgT5h2dybKpUzybR9q5kydWuXdd20aoKEW8v8OfS6G6EybBtGsWwcnT8LChbYOY2gHSf17s+PIeOKr7brBPMep0pOM6PUpRwN84wZrMZ1KE34A9Ae+LyKljU8opQa1dbGIfAT0bubUfW6OD5CqlCoRkQHAChHZiq4A2RTVihx3AHcApKbaF1fQ5QiJhx7xHpmZhIRoQ/zy5bqYlOWPL/tfguojMPgHFnd8PsuX61otxvjuW5QGLWRG/MOcPHLCa0kRrWDHimymBjaQMNK/lEmn8m8ppVpc0BORoyLSRyl1WET6oAtvNddHiXO/T0RygLHA60CsiAQ5ZycpQEkrciwBlgBMmDChRaVjaIaoDKiwX5kALFoEb78NO3fCMCu9HZWCXX+CmGHQ5yILO26e5ct1bElsbNttDZ4jcfQiAg/9ip0rP2Tqtd/wtjgdxnFwOSdi4xk24wJviwK4H7T4fEubBTK8DdzsfH0z8FbTBiISJyIhzteJwDRgh1JKASuBK1u73mABkRkeWeYCuPhivV++3OKOj66EU5thyE9sN2IcOqTjZcwSl+8xdOoEjlUkmD6SqwAAHgxJREFU4ii2+gbzHI4GB4Oj32PnqQUEBgd6WxygHUGLIjLdGVvyX+f78SJiRRWW3wLzRCQfmOd8j4hMEJF/OtsMBfJEZDNaefxWKbXDee5/gZ+ISAHahvK0BTIZmhKVCWcP6BQkNtOvn3YRtlyZ7PqTjuZPu97ijs/nfad91ygT3yMwOJBd5QsYEvOe37oI71r7OT2jSiHZd24wd4MWr0PPIKoBlwJRwEOdFUApdVwpNUcplencn3Aez1NK3eZ8vUYpNVIpNdq5f7rR9fuUUpOUUhlKqauUUvb/2nVHojLQ7sGFHhlu4UL45BNdA8QSyndCybvagysw1KJOW+bdd7VSNC7BvklA34UkRh5jx+oN3halQ5RueheHQxhyof3Lte7i7szkPmC+UuoHgEuVbwOG2yKVwffwoEcXaGVSX69zdVnC7se0Esn8rkUdtkxtrY4vWbjQuAT7KkMunE+DI4Bjm/1zqSuxbjnbj15AQrLv1DRwV5kkNyqM5TJc1wO+sVhnsJ8opxuth4zwU6ZATIxFS13VpVD4Lx3tHtrTgg5bZ/VqqKgwS1y+THyfBLYfnUxSvf8pk7LiUoYlbeB4iG/dYO4qk70iMrXJsanoOieG7kCPeAiO9Yh7MEBwsA72c7kId4pdj2lbz5CfWCJbWyxfDj16mKqKvs6J0EUM651H2YGj3halXezO+YCAAEWvcYu8LcpXcFeZ/Bp4S0TuB4JF5G50HfhO20wMfoKIM+GjZ5a5QLsIHzmi66d3mNpTkP84pF4J0YMtk60llNJuzbNmQUSE7cMZOkGfCfrHeNeKd7wsSfsIPPw2R0/3ZvAFvlXo1l3X4DeB64ALgP3oFCa3KqXes1E2g68RlemxmQnA176mg/6WLetEJ3v+BnWnYfjPLJOrNXbsgPx8WLzYI8MZOsGgiaM4eCKNkLLO3GCepaqiipEJ77H77GUEBLrtjOsR2lPPJFsp9XWl1Ail1CKlVNcrWWZonagMqNyvU5J4gPh4yMrqhDKpq9CG9+RFEOeZp7hly/Qk7tJLPTKcoRNIgLC3djGjk7I5c/KMt8Vxi63ZHxEZepaIwb73tNKeOJNUEblHRP7m3KfZJ5bBJ4nKBOWwvepiYxYv1vVAdu7swMUFS6DmOAxvT9aezvHGGzB5sk6nb/B94kYtJiS4lq3v+4chvqbgDcorYxg5d5a3RTkPd+NMFqCN7YuAGOd+p/O4obsQ6Vn3YIDLnCk92z07aaiBXX+EXhdCzymWy9UcRUWwaZNZ4vInRmRNpexMEhT7/lJXfW09w2PfZuvJr9Mj1PfKNrk7M/kD8C2l1Ayl1I1KqRnArcAj9olm8DmiPJeK3kVKCkya1AFlsvefUHXYo7OSN9/Ue6NM/IfA4EB2nrmUkQnvUn3W3uJvnWXrx6uIjzhBUJpv3mDuKpM04KUmx14GTOrd7kRIIgTHeCzWxMXixdqj64C7FXUaqmH7b6DndOg121bZGrNsmU4Dk5HhsSENFhA2aDFRoRVszf7Y26K0yukdy6iqDWXkfN+Jem+Mu8okB5jV5FgWkGulMAYfR0TbTU57Nrzo8sv13vXk3yYFT0FVCYz8pcdC0EtLdfoXMyvxP0bNnU15VTTVBb671OVocDAobBmbyxYQEeObPufuKpMCYJmIPC8ivxSR54E3gT0i8jPXZp+YBp8hdgSc2urRIQcN0qno3Vrqqq+CHf8HSTO1vcRDvP22jjExysT/CAkPYduJRQyNfov62vq2L/ACOz/No0/sIep7++4N5q4yGQN8jl7Wmuncf46uKTLPuXm/CLHBfmJH6eJS1c2WnbGNxYth1SooK2ujYcESbSvx4KwEtBdXWhqMHu2xIQ0WEpC6mMTIY2xbudrbojRL2edvUN8QyPB5X/O2KC3ibtDihW5snlucNniP2FF67+HZydVXg8MBr77aSqP6KtjxW0iaBb1meUgyreCys+Gqq0xiR39l1IKFVFRHcHrrf7wtynk4GhxkBL/EF0fm+nRlyPbEmQSKyFQRucb5PlxEwuwTzeCTxI7Uew8rk1GjdDr3F174/+2dd3hVVdaH3x+EYmih96rIANJDEbuiwMioM+IoxcGCfT7rjN2xO/rNzKczthEVRI0FFUccxXFAUBQLIAhYAAGRUAOGCKGEkPX9sU/0Gm/qbcnNfp8nz7lnn33OXjs7Ob+721olZFrxd9dr6h1fLz/TpjkPx+PHx7VYTxSp16gen20/jd5p09i3u3JFsVg25wPaNV7HvtaV+w+srPtMDsa5nH+TH4NPnQQ8HiO7PJWVui2gbkvYsTTuRY8fDx9+CGvWhLm4N8vNlbT9FbQ4Kq52ZWS4VVy9e8e1WE+UqdNtPGmpO1jyZuXyEpWzJIPcfan0GXVaok0pkbL2TB7ELQ1uAuwP0uYCEf/XSmoi6b+SVgXHxmHyHCdpScjPXkmnBdeekrQ25Frl8n6WjKT1ToiYjBnjjhkZYS4uvxPyc6HvfXG1afVqJ3DjYh+80RNj+o4YRtbOFhSsLan7G1/y9ubRq9E0lmw7jfpp9RNtTomUVUwGAXebWQFBPBMz2wGkRcGG64HZZtYVmB2c/wQzm2Nmfc2sL87J5G7g7ZAsfyy8bmZLomCTpyTSekPO51AQ35UvHTo4X13PPlvELf33q2DVo3DwRGjUPa42FQrb2LFxLdYTA1Jqp/BF7ln0a/FvcrJ2JNocABa/MZPG9bKp3a1yD3FB2cXke4oIh6Q2QDQCAZwKTA0+TwVK68uNBmaa2e4olO2pCGm93cbAOO6EL2TcOFi5EhYtCkn87AaoWQd63RZXW8ycmBxzjAvR66n6NEsfR91a+1j6xiuJNgWAgjXPkrWzOf1GnJhoU0qlrGIyHZgsqR2ApKbAA/x8V3xFaGlmmwCCY4tS8p+Fi6USyt2Slkq6X1Kd4m6UdKGkhZIWZpW6xtRTLD9Mwsd/qGv0aBd46oeJ+K3zYP0r0P1aOKhVXG1ZuNAJm594Tx56HDmQtdu7Un9buLHU+JKzLYd+LV7ni11nkVI7JdHmlEpZxeQWYBfwLa6HshXYB9xTlpslzZK0PMxPuRx1S2oN9AL+E5J8A/ALYCBuTue64u43s0lmlm5m6c2bxz58a9LSqDuoZtxXdAE0buyCZr3wAuTn7YcFl0JqB+h+TdxtychwwjZ6dNyL9sQI1RDrGEefVnPZtDozobYsfeMV6tbaR9P0qjEhV9Z9JnvMbCzQHDd/0ipw+Fgmz2hmNiyIg1L05zVgSyAShWJR0m643wKvmlnhIgDMbJM59gFTAvs8saRmXRe1MAE9E4Czz4YtW2DFjAcgZzmkPwgp8XUxsXevE5NRoyAtGjOHnkpD5+PGUaOGsWLmUwm1o+G2Kazd3pWeR1WNV1p59pkcChwKfG9m0RwjmgFMCD5PAF4rIe8YigxxhQiRcPMty6Nom6c4ErSiC9wLPL37eg7efZtbCtzulLjb8PLLsG0bXHxx3Iv2xJiOhx3Cog3DOLTGYwlzr7Lyk6X0af0+62pdhGpUjZ2wpYqJpFMlfQt8CXyAi2PybeHS3ChwL3CipFU4tyz3BuWmS3oixI5OQHt+7lwyQ9IyYBnQDBev3hNr0nq7IFl5OXEvulYtePrKKyg4YHzT7B9xLx/g4Yedz7ATTkhI8Z4Yk9/5MtqkZbJoRmLiw29+7xH25NWlz2nnJqT8ilCimEjqD7yEW4Y7HOgRHN8GpkkaEKkBZrbdzE4ws67B8bsgfaGZTQzJ942ZtQ2WJ4fef7yZ9QqGzcab2a5IbfKUgUK3KjkJ6AhmzqB7/Ve567U/8dCUTnEv/tNP4aOP4JJLoEblCsPtiRIDThnFhh3tSVn7cNzLztmWQ/8mz7Jw25hK7T6lKKX9K1wJ3GdmE81slpl9FRwnAn8OrnuqI4la0ZWXDQsuhrTerK1zNZMnw+44LxJ/5BFITYVzzolvuZ74kVI7hVUFFzGg7SzWfBbfkAtLXn2a+nVzaXb4ZXEtN1JKE5OhwKRirj0eXPdUR1Lbu0BZ8V7Rtegq57F4yBQuvrQ22dnw4ovxKz47G557zu138RPvyU3PUyaSl1+L9bMfjVuZVmC02/MIyzcPpvsREQ/8xJXSxKSpma0Pd8HMMoGm0TfJUyWQoHEf2L4wfmVueAPWToUeN0CT/hx9NPTs6eYvfrIjPoY89RTs2QOXXhqf8jyJo3mHlizYMpq+DZ8iNyc3LmUufusdDm72FTuaV70/sNLEpLTrVWOZgSc2tDgWshe5oadYk7cDPrkIGvWEw24GnJ5deqnbDT9/fuxNyM93wjV0KPT1HuCqBQ36X0aj1BwWvjS19MxRYP/nf2fbrmak/+a3cSkvmpQmFnVCIykW/QFqx8NITyWl1TCwAtgyN7blmMFH5zr38kOmONcpARMmQIsWcNttsTUB3K771avh2mtjX5anctDruKEs3TSUrnl/Zm9umbbVVZgvP1jE4Lavs3zf/1C3Xt2YlhULShOTj/gxkmK4n49iap2nctNsCKTUh82zYlvOF/dB5r+g31+g6cCfXKpXD66/HmbNcpEYY8X+/XDHHTBgAJwS/20tngShGiK/+x20Scvk4+eeKP2GCNg5/1aycxvTf2zVXNdUopiY2bGlRViMl6GeSkiNWtDiGNj839iVsXkWLL0JOpwJ3cL/k118MbRqBbfeGjszpk6FtWudoPhoitWLfiOOZ8nGo+mWfw97du2JSRnL3/2YQW3f4LO8P9KwacOYlBFr/Cp5T2S0GgY7V0Huuug/O3cdfDAGGv4CBj9R7Fv8oIPgxhth7lyYMyf6ZuTlwZ13wuDBMHJk9J/vqdyohqD3HbRqtIlPMh6LSRl7P7mVbbuaMWDs72Py/HjgxcQTGa2GuePm2dF97t5tMGc4FOyHo6ZDrZIDA11wAbRtC3/6U/RXdk2eDN9+63sl1Zm+Jx3DpxuPpwd/jvrKrmXvzCe97X9YfuBaGjRuENVnxxMvJp7IaNQT6raK7lDX/p0w95euZ3LM686pZCnUrQs33QTvvw+vleTdrZzk5LheyRFHwImVP6SEJ4ak9L2d5g22suDpv0btmQUHCihYfB1ZO1swcGzVWw4cihcTT2RIrneyebZb2RUpB/bBe7+G7E/hiGnliud+/vkuDvsll8B330VuCsDVV8PmzXD//b5XUt3pPexI5m84i6GN7mblguhs1p335EP0af0+Kw66j3qN4uv5Otp4MfFETqthsC8r8t3w+Xtg3m9gy2wY/CS0+1W5bq9d220q3LYNrozCgpi33nJDXNdeCwMHlp7fk/x0G/cgOXsac+D9c9i/b3/pN5TAuuVfM7DW9Xyy4WSOGD+h9BsqOV5MPJHzw7xJBEuEC4e2Ns6EQY9Bl4r9c/Xr5ybjn3kGXn+94ubk5Lh5mB494rOHxVM1aNqmGaubPEr3lp/ywRN/qfBzCg4UkP3W+ew/UJv2pz9WZdzMl4QXE0/kpLaFht1hQwXf3vu+g3dOhKx5cPgzcMiFEZlz003QqxdcdBFs317++83c8NbGjTBlCtQpNhC0pzoy5IzfMH/DmQxtcBsrP6mYo9N5Tz5E3zbvsaz2A7Tu0jbKFiYGLyae6HDw+bD1XdhQzvgPmTPgzcMgezEc+TJ0jjxEaeFw1/btMHw47NhRvvtvv90Nb113HQyqGkHuPHHm0LEPkr27KQ0XjeCbZavKde/8F6ZxZOpVfLJhVFIMbxXixcQTHbpdDo16wMLL3dxHaezZAh+MhfdOhTrN4aQPoX204q1B//4wfTosXeoEJaeMMbzuusuJybnnus8eTziatW3OzoGzSKm5nzrzjmPd56vLdN+H015h0P6xLN8ylB4XPJ8Uw1uFJFxMJJ0h6XNJBZLSS8g3QtIKSV9Luj4kvbOkjyWtkvSiJO8vLBHUqAXpD0PuWvjiXpe2NgP+1Qmeq+GOazMgd70TnBmdYP3L0OsOGL4AmvSPukknnwwvveSCWY0cCVlFgk1nZECnTi7AVceOMHo03HKLizH/+OM+8JWnZA4Z0JPv+s6mTq09pLx7HGuWfFVsXisw5j//Iul7z+LLrYPoct6b1E8ree9UVUMWL9/dxRkgdQcKgMeAP5jZz3yaS6oJrMT5A8sEFgBjzOwLSdOA6Wb2gqR/Ap+ZWakBCNLT023hwji6T68ufDAO1r8Cfe6BpbfAgZDIVarpJiRUAzqf7VzJN+wac5OmT4czz3RzH1dcAddcAzNnwoUX/jyw1uGHw7x5ULNmzM3yJAlffbiYFsuH0bBuDh9umUDnUbfQ7hedfri+5O130dJb6NNmHss3D6bDhLerrMsUAEmLzOxnX/wTLiaFSJpL8WJyOHCbmQ0Pzm8ILt0LZAGtzCy/aL6S8GISI/Zsgte7QX4u7jtCEVIawMnLoF7HuJq1YoVblfXii8455N69zqV8UTp0gHUx8AzjSW6yvt3CFy/fy+Cmj1JDBazPPgRDpNTIo1PTr9mc05oVtW5iyNiJ1Emt2is6ihOTqtKRbwuEBunKDNKaAjvMLL9IelgkXShpoaSFWUXHPDzR4aDWcOSLhBUSgPxdcRcSgG7d4Pnn3RzKmDHhhQRgfdhQcB5PyTTv0JJjrr6f7KNWM3/75WzN60lWXg827e3Luzvvp9HZqzlm4mVVXkhKIiUehUiaBbQKc+kmMyuL84tws1RWQnpYzGwSQRji9PT0ytElS0bajITUjrA7zFf81A7xtyeEww6DSZPg7bfD90A6JNY8TxWndZe2tL4qeu5WqhJx6ZmY2TAzOyzMT1m9KGUC7UPO2wEbgW1AmqSUIumeRNPnbqiZ+tO0mqkuvRJw992QWsS81FSX7vF4yk9VGeZaAHQNVm7VBs4CZpib8JkDjA7yTQCi6ObPU2E6j4NBk1wPBbnjoElR2UcSDcaNcz2Ujh2dz62OHd35uMphnsdT5Uj4BLykXwMPAs2BHcASMxsuqQ3whJn9Msj3S+ABoCYw2czuDtK7AC8ATYDFwHgz21dauX4C3uPxeMpPpV/NFW+8mHg8Hk/5qeqruTwej8dTifFi4vF4PJ6I8WLi8Xg8nojxYuLxeDyeiKm2E/CSsoCKOs5ohtvjUp3wda4e+DonP5HWt6OZNS+aWG3FJBIkLQy3miGZ8XWuHvg6Jz+xqq8f5vJ4PB5PxHgx8Xg8Hk/EeDGpGJMSbUAC8HWuHvg6Jz8xqa+fM/F4PB5PxPieicfj8XgixouJx+PxeCLGi0k5kTRC0gpJX0u6PtH2RBtJ7SXNkfSlpM8lXRGkN5H0X0mrgmPjRNsabSTVlLRY0r+D886SPg7q/GIQ/iBpkJQm6WVJXwXtfXiyt7Okq4K/6+WSnpdUN9naWdJkSVslLQ9JC9uucvwjeJ8tldS/ouV6MSkHkmoCDwMjgR7AGEk9EmtV1MkHrjGz7sAQ4LKgjtcDs82sKzA7OE82rgC+DDm/D7g/qHM2cH5CrIodfwfeMrNfAH1wdU/adpbUFrgcSDezw3DhLM4i+dr5KWBEkbTi2nUk0DX4uRB4tKKFejEpH4OAr81sjZnl4eKonJpgm6KKmW0ys0+DzztxL5i2uHpODbJNBU5LjIWxQVI74GTgieBcwPHAy0GWpKqzpIbA0cCTAGaWZ2Y7SPJ2xoUqPyiIzpoKbCLJ2tnM3gO+K5JcXLueCjxtjo9wkWtbV6RcLybloy2wPuQ8M0hLSiR1AvoBHwMtzWwTOMEBWiTOspjwAHAtUBCcNwV2mFl+cJ5sbd0FyAKmBEN7T0iqRxK3s5ltAP4KfIsTkRxgEcndzoUU165Re6d5MSkfCpOWlGurJdUHXgGuNLPvE21PLJE0CthqZotCk8NkTaa2TgH6A4+aWT8glyQa0gpHME9wKtAZaAPUww3zFCWZ2rk0ovZ37sWkfGQC7UPO2wEbE2RLzJBUCyckGWY2PUjeUtj9DY5bE2VfDDgCOEXSN7ihy+NxPZW0YDgEkq+tM4FMM/s4OH8ZJy7J3M7DgLVmlmVm+4HpwFCSu50LKa5do/ZO82JSPhYAXYPVH7Vxk3czEmxTVAnmCp4EvjSz/wu5NAOYEHyeALwWb9tihZndYGbtzKwTrk3fMbNxwBxgdJAt2eq8GVgvqVuQdALwBUnczrjhrSGSUoO/88I6J207h1Bcu84Afhes6hoC5BQOh5UXvwO+nEj6Je5ba01gspndnWCTooqkI4F5wDJ+nD+4ETdvMg3ogPunPMPMik7yVXkkHQv8wcxGSeqC66k0ARYD481sXyLtiyaS+uIWHNQG1gDn4r5gJm07S7odOBO3anExMBE3R5A07SzpeeBYnKv5LcCtwL8I066BqD6EW/21GzjXzBZWqFwvJh6Px+OJFD/M5fF4PJ6I8WLi8Xg8nojxYuLxeDyeiPFi4vF4PJ6I8WLi8Xg8nojxYuLxFIOkDpJ2SWqTaFvCEdh2eKR5oo2k4ZLmRXB/z8CTcZ1o2uWJLV5MPAlH0lxJ+yTtlJQjaY2kZyQNSKRdZvatmdU3s0q5Izqw7UNw+2Mk5ZeUJx4E+xbux+1tKEz7taT1krZIuqRI/n8qCHNQiJl9jtvv8ft42OyJDl5MPJWFO82sgZk1Ao4D1gEfSfp1gu3ylI+TcJsg58APYRseB36F87p9l6SWwbVjgV7Ag2GeMxn4H0n+HVVF8A3lqXSY2Tozuxl4Gngw+LaLpKaSnpa0SdJmSVMlNSm8T9I3km6WC+61S9IySb0ljQmC/+QE3nFTQu6ZEnxr3inpC0ljQ651kmSBe3ok3SZptqR75IIPbQ12VIcl5P6JklZK2iHpNUktQvKUVqfLJa0N7Nsg6Z6QaybpyGAYbiZQM6j3LkkTQvOE3HO6pM+C38VnoWIt6Zzg93S5pExJ2ZIeCwShrJwGzLIfd0M3Bfab2RIzWwesAjpKOgi383qimRWEec57QCugbznK9iQQLyaeyswLOFcXhf6jMoDGuMBk3XHuIp4pcs8E4NIg32fAq7ieTh/ct+BTgN+G5H8f98JKA+4AnlLJAc+OxrmjaIP7tn2jpCNKqcfvgvs64FzUPBtyrdg6SToUuBcYZWYNgJ6E8QUXDMONBA4Ew1r1zWxq0XzB3EkGzjtwU5ybnOclDQ7J1hFoCRwMDATOwPkrKyv9cf6uCskCDkgaIKkz0An4GrgLeM7Mvvz5IyBwZ7IqeJ6nCpBSehaPJ2FkBsemwbfv4cChZpYNIOlq4CtJrUOc000qfEFJeg4YBwwxs1wgV9Jc3EvyOQAzezKkvBck/QHn1yj0hRjKSjP7Z/D5Y0lLgHTggxLqcXvgWBFJfwRWhUzqF1snnP8oAT0lrQuCV31UQjmlcS7wipnNDM7fkPQqcB7O9xrAHuBPZnYA+FrS7KB+GWUsozHwQ8gCMzNJY3BRHWsDlwCHAMcAR0t6ACcYq4HLg4BshXyP85flqQL4nomnMtMuOG7nRzfZa0Ourw6OoS60Qz2e7sZ9W88qktYAQFINSXdIWhEM++zA9WCal2BTUY+quYXPK4FvwnxuRyl1MrM1ODG8ANgo6X1JJ5VSVkm0xzl0DGU1P/39bQ2EpJCy1C+UbKBhaIKZzTOzI81sEPAGMAlXp7OBxmZ2NO73WjSeSkN+HjHQU0nxYuKpzJwJbABW8GM0uE4h17sEx9BIceVhDM5r7Om4l1oabmgsXMCgSOgU5nMmZaiTmU03sxNxw1/TgNckpYYpI9y8Q1HW4wJDhdKFiv/+wrEYN2RXHDcDb5jZYpxwzw/S38NF9QRALsRD1+B5niqAFxNPpUNS+2Bi+xzgiiA+9UbgbeBvktLkoub9DZhZ0fgLuG+++bhx/RqSzsO94KLNLZJaysVdvw+YbWYbS6uTpG6SRgTisR8XZtYILxybcRPwRcUilKeA0+X2gdSUNBL4DTAlajV1rs5PCHdBUm9cpMM7gqTVwInBgoiRuLmUQo7GuU/3YlJF8GLiqSzcEqxY+h73LfUQYKiZvRKSZzywE/gq+NmBm9yuKFNxcwVf43pAPXCxXKLNs8Fz1+PmDcaHXCupTrVx+zU2BemXA6eb2d6iBZjZSuAR4JNg1djZYfLMxy1Q+CtuOOp/cbE7IpmHKcp/gHy5Zb8/EKwIewK4JCRWyGO4Om7Hifg9IbecBzxYzEovTyXExzPxeGKEpE64+ZD2ZpZZcu7kQdII4MZgLqQi9/fArcLrXZWDVFU3vJh4PDGiuoqJp3rih7k8Ho/HEzG+Z+LxeDyeiPE9E4/H4/FEjBcTj8fj8USMFxOPx+PxRIwXE4/H4/FEjBcTj8fj8UTM/wPFHqgDJfOdiQAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "x0 = [12, 38, 63]\n", "x1 = [12, 25, 63]\n", "yw = mwarp1d.warp_landmark(y, x0, x1)\n", "plot_warp_results(y, yw, x0, x1)" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "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.6.9" } }, "nbformat": 4, "nbformat_minor": 2 }