{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Scripting sequential manual warps\n", "\n", "The **mwarp1d** GUI supports sequential manual warps for arbitrarily fine control over local warping effects. This control is replicated in the **SequentialManualWarp** class, as demonstrated below." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZMAAAEJCAYAAABR4cpEAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOydeXhcdbn4P28m+74nzb523/cVKTuKoFwBQRRQQQURvfq7FwVEUbgqKlyvgIJsIoIIQlEQKLRA27SFdN+bpVmbZl+afZnv748zE9I0yySZmXNOcj7Pc56ZOXPO+b6TzMw77y5KKSwsLCwsLCaCj94CWFhYWFiYH0uZWFhYWFhMGEuZWFhYWFhMGEuZWFhYWFhMGEuZWFhYWFhMGF+9BdCL2NhYlZGRobcYFhYWFqZi165ddUqpuMH7p6wyycjIID8/X28xLCwsLEyFiJQOtd9yc1lYWFhYTBhLmVhYWFhYTBhLmVhYWFhYTJgpGzOxsLCYevT09FBRUUFnZ6feohiewMBAUlJS8PPzc+l4S5lYWFhMGSoqKggLCyMjIwMR0Vscw6KUor6+noqKCjIzM106xzBuLhF5SkRqROTgMM+LiPxORApFZL+ILB7w3A0iUuDYbvCe1BYWFmais7OTmJgYS5GMgogQExMzJgvOMMoEeAa4ZITnLwVyHdstwGMAIhIN3AusAJYD94pIlEcltbCwMC2WInGNsf6dDOPmUkp9KCIZIxxyBfBnpfXM3yEikSIyDTgX2KiUagAQkY1oSukFz0psDGpaOnlldyUdXT3EtheQ3bSDWXE2ooMDtAMCIyA4BkLiIHEuhCXqK7CFqWjr6uXlXRXUt3ZpO0RYmxPL8sxofQWzMByGUSYukAyUD3hc4dg33P6zEJFb0Kwa0tLSPCOll6hv7eKPHxbz+vb9fEm9yed9dpDpcwoAe4lgF/BhiFk10VmQvhoyPwVZ6yH0rEJWi8lAbzdU7YPUZeM6vaO7j7/sKOWxD4poaOvG+SNVKfjdewWsy43lPy+czqI0ywkwHioqKrjttts4fPgwdrudyy67jAcffBB/f/8zjjt58iTf+c53ePnll0e83qc//Wn++te/EhkZOWZZfvKTnxAaGsoPfvCDMZ87EDMpk6FsLjXC/rN3KvU48DjA0qVLTTsVbEtBLd9+bidX2t9mk/8/CFLtSOY5MPtOGtIu4g/5Lfx5ewk9fXYeuDSNa2YHQ2s1VO6C0jw48i/Y8xftYonzYfYVsOh6y2qZLJRuh399F2qPwo1vQsaaMZ1eVNvKl/+0k5PNnWcpjYFK5vOP5nHj6gzu/exsy3U0BpRSXHnllXzrW99iw4YN9PX1ccstt3DXXXfx4IMP9h/X29tLUlLSqIoE4M033/SkyC5hJmVSAaQOeJwCnHTsP3fQ/ve9JpWXySus44FnN/C6/8Ok28sh41y4+H8gYTYA0cCPPp3M19dlcucrB/jvN8ro85/HdStWaxbJ6tvB3qf9ai3aBIXvwqafweYHYPolsOQGyLkAfGy6vk6LcdDbDf/+L9j1NIQlafuq9o5JmZyoa+Pax3dgV4oXb1nJyqyYM54P8rdx8zlZXLcijQffPsYzeSUopfjJ5XMsheIimzZtIjAwkJtuugkAm83GQw89RGZmJpmZmWzevJnOzk7a2tp46qmnuOyyyzh48CDt7e3ceOONHD16lFmzZlFSUsIjjzzC0qVL+9tDtba2cumll7J27Vry8vJITk5mw4YNBAUF8cQTT/D444/T3d1NTk4Ozz33HMHBwW57XWZSJq8D3xaRF9GC7c1KqSoReRt4YEDQ/SLgh3oJ6Ul2FNfz62df4gW//yE80B8ufwFmXApDfIjjwwJ57PrFfPO5Xfzo1QP4+ghXL3PoYh8bJC/WtnN+APVFsPtZ2PtXOPYGhKdolsqSGyA8ycuv0mLcHHxFUyQrvgnn3QO/WwjVh10+vbReUyS9dsULN69kRmLYsMeGBPhy72dn42cTnthyAh8f4ceXmctC+ek/D3H4ZItbrzk7KZx7PztnxGMOHTrEkiVLztgXHh5OWloavb29bN++nf379xMdHU1JSUn/MY8++ihRUVHs37+fgwcPsnDhwiGvX1BQwAsvvMATTzzB1VdfzSuvvML111/PlVdeyc033wzA3XffzZNPPsntt98+sRc8AMMoExF5Ac3CiBWRCrQMLT8ApdQfgDeBTwOFQDtwk+O5BhH5GfCx41L3OYPxk4kjVS088syfec72KwJDo/G5YQPE5ox4ToCvjceuX8Itz+3iv/+xn/AgPy6ZO4QrKyYbLrwP1t8Nx97UFMsHv4Qtv4H518Ca70DcDA+9Mgu3UbIVgqI0S9XHB+JnQ80hl05tbOvmuid20tnbN6oicSIi/OjTs+i1K57eVkJkkD93XJA70Vcx6VFKDal0nfsvvPBCoqPPTnDYunUrd9xxBwBz585l/vz5Q14/MzOzX9EsWbKkXyEdPHiQu+++m6amJlpbW7n44ovd9Io0DKNMlFLXjvK8Am4b5rmngKc8IZcR6Omz88e/vsQTPg/gG5WG7cYNEJHi0rmBfjYe//ISrv7jdn706gGWZkQRGxow9MG+/jDnc9rWcAJ2PAq7n4O9f4GZl8G5d0LiPPe9MAv3UroV0lZrigQgYQ7kP6W5NUdxW/7kn4eobunklW+tZta0cJeXFNEskoa2bn63qYDzZsYzLyViIq/Ca4xmQXiKOXPm8Morr5yxr6WlhfLycmw2GyEhIUOep30Fjk5AwCefb5vNRkdHBwA33ngjr732GgsWLOCZZ57h/fffH98LGAYj1ZlYDMOT7+7je82/RAXH4vu1t1xWJE4C/Wz8+qoFtHb2cu8G136pEp0Jn34QvncQzvkvOLEF/rAW/nY9nDowjldh4VGaK6Gx5Mz4SMIc6O3UfhiMwNuHTrFh70luPy+XBaljzwYSEe67fC6xof784O/76O61j/kaU4nzzz+f9vZ2/vznPwPQ19fH97//fW688cYRYxhr167lpZdeAuDw4cMcODC2z+Hp06eZNm0aPT09PP/88+N/AcNgKRODc6SqhYRt95DqU0fQF58edyrv9IQw7rgglzcOVPHG/irXTwyJhfPugu/uh0/dCcUfDlAqQzYrsNCD0jztNn2AMonXkjJGcnU1tXdz16sHmT0tnFvXZ497+YhgP/7nynkcqz7N7zcVjPs6UwER4dVXX+Xvf/87ubm5TJ8+ncDAQB544IERz7v11lupra1l/vz5/PKXv2T+/PlERLhuBf7sZz9jxYoVXHjhhcycOXOiL+NslFJTcluyZIkyOt29fepXv/qpUveGq/a3fz7h6/X09qnLfrdFLb7vHVV3unN8F2lvVGrTA0o9kKLUveFKvXSDUvXFE5bNYoK8fof2P+nr/WRfV5tS90Zo/69huOOF3Sr7h2+og5VNbhHje3/bo7J++IY6UOGe67mbw4cP6y3CuOnt7VUdHR1KKaUKCwtVenq66urq8uiaQ/29gHw1xHeqZZkYmA3v7+CbrY/SGLOYoPP/e8LX87X58OurFtDS2cOv3zk+vosERcL6H2qWyjn/D46/DY8sh3fugY6mCctoMU5Kt0HayjNjI/7BWnJF9dAW5EcnGnht70luXZ/DnCT3xDnuvWwOMSH+3LPhoMs+fgvXaG9vZ+3atSxYsIDPf/7zPPbYY2cVOeqJpUwMSmdPHyFb/wd/sRN5/dNgc0+uxIzEMK5dnsbf88spb2gf/4WCouC8u+H2XTDvKsj7P/i/JVoxpN3ymXuV1hqoO67VEQ0mfjbUDJ0e/NuNx4gLC+DWc8fv3hpMRLAf371gOnvKmnj/eK3brmsBYWFh5Ofns2/fPvbv38+ll16qt0hnYCkTg/LvdzdysX0LdXNvQqIy3HrtW8/NwcdH+P2mwolfLDwJPvco3PK+9it4w23w1EVwcs/Er23hGv3xkrVnP5cwRwvAd7edsTuvqI4dxQ3cdm42gX7uLVD9wpIUUqKCeGjjccs6mUJYysSAtHf3EvfRL+nwCSH5M3e6/fqJEYF8aUUaL++uoLS+bfQTXCFpIdz0FnzuMS2r6PH18K/vQfukK/kxHqXbwC9Y+x8MJn42oLTWKg6UUjy08TiJ4YF8cbn7e9T5+/rwnfNy2V/RzHtHatx+fQtjYikTA/Luv//BWrWbhsW3ae4kD/Ctc7Pxswm/e88N1okTHx9YeJ3m+lrxTdj1rOb62vWM5fryJKV5kLocbENMxEtw1FIMqITfWljHxyWN3HZejtutEiefX5xMekwwv7WskymDpUwMRltnD+l7HqTBFkPqxd/12DrxYYFcvyKdV/dUUFzb6t6LB0bApb+Ab26B+Fnwzzvgmc9AnZUy6nbaG6D60NAuLoCoTM1qqdbSg5VS/HbjcZIjg7h66djqlcaCn02zTg5XtfD2oWqPrWNhHCxlYjC2vfVXFnCc0yt/oGXjeJBvnpuNv68Pf/yg2DMLJMyBG9+Ay3+vBYEfWw0fPgh9vZ5ZbypSvhNQQwffQbMW42b215psL65nT1kTt63PIcDXs808r1iYRFZsCI++70brdxLwve99j4cffrj/8cUXX8zXv/71/sff//73+e1vf+s1eUJDQ91yHUuZGAi7XRF94EnqfWJIP+9mj68XGxrA5xcls2FfJc3tPZ5ZRAQWfxm+/THM/Axs+rlmpTSVeWa9qUbVfkCGjpc4SZjT7+Z6bnspUcF+XLl4yJE/bsXX5sONazLYX9HM3nIrbdzJ6tWrycvTkibsdjt1dXUcOvRJYWleXh5r1oze6Vkphd1A7mOXlYmIZInIj0TkEcfjGSKiT3ObScru3TtY2reP6ulfGtr/7QGuX5lOZ4+dv+8qH/3giRAaD1c9A1c+oblcHlsLB0af02AxCtUHtYFn/kP3cwI0ZdJeR83JMt45XM3Vy1I9FisZzOcXJRPib+O57aVeWc8MrFmzpl+ZHDp0iLlz5xIWFkZjYyNdXV0cOXKEWbNmcf7557N48WLmzZvHhg0bACgpKWHWrFnceuutLF68mPLyckJDQ/n+97/P4sWLOf/886mt1VKyi4qKuOSSS1iyZAnr1q3j6FEtCePEiROsWrWKZcuWcc8997jtdblUvCAiFwL/ADajdfa9DYgF7kabzW7hBk5/+Chd+JF96be9tuacpAiWpEfxlx2lfHVNJj4+Hm4hPv9qSFkG/7gFXvkaFL4Hn/4VBIzepdZiCKoP9c+yGRZHW5UPtn2IXUVx/Yp0LwimERbox+cXJ/NSfgV3fWYW0SHGKbLj33e6v89c4jwtXjgCSUlJ+Pr6UlZWRl5eHqtWraKyspLt27cTERHB/PnzCQ4O5tVXXyU8PJy6ujpWrlzJ5ZdfDsCxY8d4+umnefTRRwFoa2tj8eLF/OY3v+G+++7jpz/9Kb///e+55ZZb+MMf/kBubi47d+7k1ltvZdOmTdxxxx1861vf4itf+QqPPPKI2166q5bJL4CrlFKXA32OfbuBxW6TZIpzsrqa5c1vczz2IgIiEry69ldWpVNS386WwjrvLBidCTf9W2sguf9F+MM6qNjlnbUnE93t0FAMCXNHPi5KUx5Hjh5i/Yx4UqM9G4sbzFdWZdDda+elfA9bvybCaZ04lcmqVav6H69evRqlFD/60Y+YP38+F1xwAZWVlVRXa4kM6enprFy5sv9aPj4+XHPNNQBcf/31bN26ldbWVvLy8rjqqqtYuHAh3/jGN6iq0nrybdu2jWuv1Zq0f/nLX3bba3K1rDpbKfWW474CUEp1iIh3fDFTgCP/fozzpYv4C9w3rMZVLpmbSGyoP89tL+VT0700E97mqzWQzF6vWSlPXQTrfwRrvmtNeXSV2iOA+iT9dzgcUxdDO2v48krvWSVOpieEsSIzmud3lnLzuixsnrZ+XWUUC8KTOOMmBw4cYO7cuaSmpvKb3/yG8PBwvvrVr/L8889TW1vLrl278PPzIyMjg87OToBhW9Q7ERHsdjuRkZHs3bt32GPcjauWSbmInPHzR0QWACVul2gK0tXTQ27JCxQGzCZh5iqvrx/ga+OaZalsOlpNReMEWqyMh/TV8M2tMOuz8N598OcroOWkd2UwK4503/7uwMPh60+TTxS5gU3e+7EwiC+vSqe8oYMPjltFjKBZJv/617+Ijo7GZrMRHR1NU1MT27dvZ9WqVTQ3NxMfH4+fnx+bN2+mtHT4mJPdbu+fE//Xv/6VtWvXEh4eTmZmJn//+98BLVi/b9++/rVffPFFALe2ondVmfwO+IeIXA/YROQ/gL8AD7lLEBG5RESOiUihiJxV9i0iD4nIXsd2XESaBjzXN+C5190lk7fYvekV0jhFzxLPZ3ANx3UOP/oLH+mQZRUUCV94Gq54BCp3w2NroPzj0c+b6lQf1mpIojJHPKyg+jSlvVEsjGjzfExsGC6ek0h8WAB/2WFl8QHMmzevPxYycF9ERASxsbF86UtfIj8/n6VLl/L888+P2DI+JCSkfxTwpk2b+PGPfwxoiuLJJ59kwYIFzJkzpz+I/7//+7888sgjLFu2jObmZre9JnG1OlVEbgFuBzLRLJKHlVJ/cosQIjbgOHAhUIE2gvdapdSQHepE5HZgkVLqq47HrUqpMSVLL126VOXn509McDex81dXMKN9F+F3FeHjN8wURC9w49MfUVDdypb/Wq/blw51hfD8F6C1Gq5+DnIv0EcOM/DMZdDTATe/N+Jhv3zrKAvzvs35cS343v6Rl4QbWo7HPyxm54/OH37ap4dxZkpNJkJDQ2ltdXPhsYOh/l4isksptXTwsS6nBiulHldKzVNKhSql5rpLkThYDhQqpYqVUt3Ai8AVIxx/LfCCG9fXjfr6Wha0baMo4WJdFQnA5xYmU9nUQX5po35CxObA196BmBx44RrY/5J+shgZpVzK5LLbFa/vPYlPRDK+rWMYiuYBPrcwmT67GttwNgvT4JIyEZGk4TY3yZEMDEz1qHDsG0qWdDTraNOA3YEiki8iO0Tkc8MtIiK3OI7Ld+Zi682R9/5CoPQQt+YreovCRXMSCPa38eqeSn0FCY3XKufTVsE/boYdj+krjxFprYaOhlEzuT4uaaCyqYOktGzoaoHOFi8JeDYzEsOYNS1c//fXJMNTVslYcdUyqUD7sh9qcwdD+VSG8799EXhZKdU3YF+aw+y6DnhYRIYc0OCwrpYqpZbGxekTiBxMeME/qPBJIm3eOXqLQrC/LxfNTuCN/Sfp6u0b/QRPEhgOX3pZC8y/dSe89zPt17iFhnPg1SiZXK/trSTY30ZOjsPn3qLvF/nnFiaxt7yJE3Vu6lY9DqzGk64x1r+Tq8okE8gasK0D3gZuHNNqw1MBpA54nAIMl9LzRQa5uJRSJx23xcD7wCI3yeVRyk4cZ273AU6lX661HTEAn1uUTEtnL+8fM4Dl5hcIVz0Li78CW34N//ou2HVWckbB2QV4hEyurt4+3thfxcVzEgmIdny8dFYmly9MQgQ27NVHjsDAQOrr6y2FMgpKKerr6wkMDHT5HJfqTJRSg/PSSkXkBjRX03OuizgsHwO5IpIJVKIpjOsGHyQiM4AoYPuAfVFAu1KqS0RigTXAr9wgk8cpff9Z0kSRfu5NeovSz9qcWGJD/XltTyUXz0nUWxyt5uSzv4PgWNj6W+ho1Fqy+OobX9Kd6kNa/Uhw9LCHbD5aS0tnL1csTIIIrUaBZn2VybSIIFZmxvDankruOD/XI/UOI5GSkkJFRQVGcXMbmcDAQFJSXO8sPZFZsO2AWyqglFK9IvJtNGvHBjyllDokIvehDa93pvteC7yozvxZMQv4o4jY0SytXwyXBWYklN1OStnrHPOfzYz04dP+vI2vzYfL5ifx14/KaO7oISLIAHWpInDBvRAcA+/cBZ3NcM3zEOCebqempPrQqC6uDXsriQ31Z21OLFrjCtHdMgGtX9d/vbKffRXNLEyN9Orafn5+ZGaOnEptMT5cDcBfN2i7GfgnsNNdgiil3lRKTVdKZSul7nfs+/EARYJS6idKqTsHnZfnyDJb4Lh90l0yeZJj+/LIVGW0Tr9Sb1HO4nOLkunutfPWQYNl3az+tjbJ8cQW+PPlugaTdaWvB+qOjahMmjt6eO9IDZfNT8LX5qM1Dg1NMIQyuWReIv6+PrxmBeInFa7GTO4ftH0HKMV9MZMpR/3OF+lRNqafr38W12AWpESQERPMv4yYwrnwOrjmOTi5FzbcOjWD8vWF0Nc9ojLZeLia7j675uJyEpGsu5sLIDzQj/NnxvPGgSrs9in4/5ukuKRMlFKZg7Z5SqmblFL6vzNNiLLbSat+l2NBCwiL8m5TR1cQES6em8j2onrPzTmZCDM/Axf+FI78E/J+p7c03seFNipvHzrFtIjAM91I4UmGsExA6wdXe7qLPeU61jRZuBVrOJYOFB3OJ1VV0ZH9Gb1FGZZL5iTSa1dsOmbQkaurvg2zr4B3fwLFH+gtjXepPQZig9jcIZ9u7+7lw+O1XDwn8cwAd3iKZpkYwJpbPzMeP5tYI30nEcMqExHZKCLvjLZ5U9jJQvXOl7ErIWvd1XqLMiwLUiJJCA/grYOn9BZlaES0Xl4xOfDyV6dWc8i6YxCVMWxG2wfHaunqtXPRnEFWb0Qy9LRpCQw6Ex7ox+rsWN46eMpK050kjGSZbAW2ubBZjJGEyo0c959FTGKa3qIMi4+PcNHsRD44XktHt0FrOwLC4Jq/aP2pXv7q1JktX3sc4mYM+/Tbh04RFezH8oxBacPhjviJQVxdF89JpKyhnaOnTustioUbGDY1WCn1U28KMlUoLz5Cjr2Yndn/qbcoo3LJ3ESe21HKhwW1xqg5GYq4GfDZ/4V/fB02/UyLpUxm+nqhoQimXzTk0929dt47WsMlcxK1LK6BhDtqBporR5+B4gUunJ3AXa8d4O1Dp5g1LVxvcSwmyJhiJiISJCKpIpLm3Dwl2GSlfJvWuDB1tXFdXE6WZ0YTEeTH20Z1dTmZfxUsuRG2PQzH39ZbGs/SVKplcsUObZlsL67ndGfv0Mo/wtHuziCWSVxYAEvSoozrSrUYE67WmWSJSB7QitZ+/sSAzWIMRJW9RbEtk6RM47fB9rP5cP6seN49Uk1Pn11vcUbmkl9Awjx49RuTO35Sd1y7jZ0+5NNvHzpFsL+NtbmxZz8ZmgjiYxhlApr1e/TUacrqvTyUzcLtuGqZ/B6tqeMC4DQwH3gN+JqH5JqU1FaVMaP7CDXJQ7sojMglcxJp6exlZ3GD3qKMjF8QXPUMdLfBpvv1lsZz1B7TbofI5LLbFRsPV7N+RjyBfkOMPrb5agrFALUmTpwW1NuHLOvE7LiqTFYAX1dKHQRQSh0CvgH8P08JNhkp2vISPqJIXPkFvUVxmXOmxxHkZzPHhz02B5bfAnuf/6QWY7JRV6BVsged3YZkT3kjtae7zs7iGkhEsqEsk9ToYGZNCzfH+8tiRFxVJnagw3G/VUQigQbAipmMgcATG6mSeNJnnjWkzLAE+tlYlxvLpqM15kjhXPd9rX39uz/RWxLPUHdsWBfXe0dqsPkI586IH/58AxUuOrlwdgK7yxppbOvWWxSLCeCqMjmE1o0XtH5cD6HNhbdiJi7S2dHGjPbdlMeuQ3zMVSu6fmY8lU0dFNQYYwjPiARHw9r/hIJ3tB5ekwmltJjJMMpk87FalqZHjdyc00CFi07Wz4jDruDDAquTr5lx9VvtO2iWCGiurWRgKZqry8IFju18iyDpJnD2pXqLMmbOnaENEtt8tEZnSVxkxTcgPBk23gN2gycOjIXWGq3gcIgak1PNnRypamH9zBGsEtDcXL0dWit/gzA/JZLoEH/zvL8shsTV3lz7lVIHHPeLlVIXKaVWKqW2j3auhUbHoTfpVH5MX2E+ZTItIoiZiWFsPmaSD7tfEJx3N5zcA8fe1Fsa91E3fPD9fcf/Zv1ILi4wXOEigM1H+NT0OD44Xkuf1fjRtLiaGvyuo/W862O3LPpRdjupdVs4HryIwGBzzuA4b2Y8+SWNtHQasPHjUMy7WnPpfPS43pK4j/604LMtk01Ha0iODGJ6wijvr4GFiwZi/cx4Gtt72FfRpLcoFuPEVTfXNuDnQJWIPCYiy9wtiIhcIiLHRKRQRO4c4vkbRaRWRPY6tq8PeO4GESlwbDe4W7aJUl50gGRVTUfmBXqLMm7Wz4yn167YWlCntyiuYfOFpTfBiQ+09iOTgdrj4B/6iXXhoKu3j22FdZw7I270yYWRjvG9zeUeEnJ8nJMbi4/A+5ary7S46ua6VymVBfwHEAJsFpGDIvI9dwghIjbgEeBSYDZwrYgM1V/7b0qphY7tT45zo4F70dKXlwP3Okb5GoaTH20AIHXZFTpLMn4WpUYSHuhrLr/24hvA5g8f/0lvSdxD3XHNxTVIYeSXNNLW3Te6iwsgJF77mzSVeUjI8REZ7M/itCg2mcWVanEWY0orUkptUkp9BW1cbynwazfJsRwodMRjuoEXAVe/eS8GNiqlGpRSjcBG4BI3yeUWQso2UeqTSlKmccbzjhVfmw/nTI/j/eO15hloFBoHcz4P+16ArknQTHCYTK7NR2vwt/mwOidm9Gv4+EBEquGUCWjW78HKFmpaOvUWxWIcjLU31wIReRg4AswD/sdNciSjVdg7qXDsG8x/iMh+EXlZRFLHeK4utJ5uYkbnfqriz9FblAmzfkY8tae7OHTSRONyl90MXS2w/296SzIxuk5rQfMhlMmmYzWsyIom2H/Yvq1nEplqODcXfJI1+P5xK0XYjLgagL9dRHYD24EE4HogXSl1t5vkGMrRO/jn7z+BDKXUfOBd4NkxnKsdKHKLiOSLSH5trXfesMe3/wt/6SNs3qe9sp4n+dSMOEQwT1YXQMpSmLYAPvqToWorxkxdgXY7SJmU1rdRXNvGeaOlBA8kMg2ajKdMZk8LJyE8oD8zzcJcuGqZfAV4AkhSSl2rlHpHubccugJIHfA4BTijW59Sql4p1eV4+ASwxNVzB1zjcaXUUqXU0ri4OLcIPho9R9+hTQUyfZl5g+9OYkMDmJ8cYa4Pu4hmndQegTITZ7I7lcmgGpMPHL/iR6x6H0xEGrTVaEGo2U4AACAASURBVHNgDISIcO70eLYcr6PX6I1FLc7C1QD8MqXUY0opT+XtfQzkikimiPgDXwReH3iAiEwb8PByNFcbwNvARSIS5Qi8X+TYZwhSG3dwPGQxfv6TI6t6XW4c+yqazZMiDDD3Si3ofOzfeksyfuqOa6N6ozLP2L2loI6UqCAyYoJdv1akowtSc4UbBXQP66bHcrqrl30V+k+DtBgbhujroZTqBb6NpgSOAC8ppQ6JyH0icrnjsO+IyCER2YdWkX+j49wG4GdoCulj4D7HPt2pLD5EkqqmO8388RIna3Nj6bMrthfV6y2K6/iHQMpyKH5fb0nGT30hRKWDr3//rp4+OzuK6lmXGzt6SvBAnOnBTaVuFnLirMmORQTzpKBb9GMIZQKglHpTKTVdKZWtlLrfse/HSqnXHfd/qJSao5RaoJRar5Q6OuDcp5RSOY7tab1ew2DK87Xq66Ql5o+XOFmcFkWwv818H/asc+HUfmgzkRIcSH2RNu9+APvKmzjd1cu63DG6bJ2WiQHjJlEh/sxLjmCL1afLdBhGmUxG/Es/4BRxpGTP01sUt+Hv68PKrBi2FppQmQCUfKinFOPDbtcsk5gz26hsKahDBFZnu5ASPJCwaeDja8j0YIC1ObHsKW/itJlcqRaWMvEUvT3d5LTtoixquem6BI/G2pxYTtS1Ud5goul4SYsgINycrq7TJ7XmjDHZZ+zeWljH/OQIIoP9hzlxGHxsWiNMA6YHwyeu1B1GH8hmcQYuf8uJyFoReVxE/ul4vEREJk8wwM0U7ttCOO3Ycs/XWxS3s84xEtZU1onNFzLWmlOZODO5Bri5Wjp72FveNPR4XleITDOsZbIkPYogPxtbLVeXqXC1zuQ6YAPQCTgViALu85Bcpqdx/zvYlZC9/DN6i+J2cuJDSQgPMGfcpLFE28xEfaF2O6Bb8PaievrsirU540xxN2itCUCAr40VWdFsMdv7a4rjqmVyF3CxUuo7aFMXAQ4Cczwi1SQgsmoLRb7ZRMYm6i2K2xER1uXGsa2ozlwtwzM/pd0Wf6CvHGOlvgj8QrRYh4OtBXUE+9tYnH72+F6XiEiF01XQa8zphmtzYimua6OyyVi1MBbD46oySVJK5TvuO789egGb+0UyP60tjeR0H6UuYc3oB5uUdbmxNLX3cOikieoB4mZAaKL5XF31BVq8ZED679bCOlZkRhPgO86PYGQaoKDFeLUmAOdM1ywuy9VlHlxVJkUisnrQvtXAMTfLMyko/Ogt/KSPsDkX6i2Kx1iTo/nqTeWKENFcXSc+NNcExvrCM+IlFY3tnKhrY+1YU4IH0l9rYsy4Sa7DlWqq99cUx1Vl8nNgg4jcDfiJyPeBF7BiJkPSdexdOpQ/uUsmX/DdSWxoALOmhZswbvIpaK+DmkN6S+IavV3aF/6AeInzb75uvMF3MHStCWiu1DU5sWwrrDNPl+opjqvtVF4DrkObGVIKnAd8VSll4v4UniOh/iMKguYREDiGFhcmZG1ODLvKGuns6dNbFNdxxk1OmKTepOEEKPsZlkleUT1xYQHkxk9gamd4MoiPYS0T0KrhG9t7OHpqEowPmAK4nBqslNqolPqsUmquUuozSqmNnhTMrNSdKifDXkZb0uSNlzhZlR1Dd6+d3aWNeoviOhHJEJlunqaPzkwuR42JUoq8onpWZcWMrYXKYGx+EJZk2FoT0N5fAHlFJrN+pyjDKhMRSXJl86awZqBk11sAxMw1f5fg0ViWEY3NR8gzU58ugLRVULbDHC3p68+sMSmqbaWutWvsVe9DYeBaE4CkyCAyY0PM1QduCjOSZVKBNnRquM35vMUA7EUf0EIwWfNW6S2KxwkL9GN+SoT5fjmmrYS2Wmgo1luS0akv1EbtBkYA9Cvu1dkTiJc4iUw1bMzEyarsGHaeaLBa0puAkZRJJpA1wuZ83mIASY35FAUvxNdvjC0uTMrq7Bj2VzTT2tWrtyiuk+5ITDSDq6uu8Izge15hPcmRQaRGB0382pFp2vTGPuP+71Znx9Da1ctBM033nKIMq0yUUqWubN4U1uicKisgRVXRlTI4i3rysjo7ll674uMSE/VRip0OQdHmUCb1hf3xErtdsb24ntXZE4yXOIlIBdWn9f4yKCuzrLiJWRhLb65rROQtETnouP2iJwUzI+W7tZlc8fMv0lkS77EkPQp/m4+5/NoimqurbIfekoxMR6OWxuyIlxyuaqG5o4fVOW6Il8CA9GDjxk1iQwOYkRBmrvfXFMXV3lz/DfwObfjUQ47bhx37LZyc+JBGwsmYtVRvSbxGoJ/W0sN0vxzTVmq/+lsNXGFdX6TdOlrP7yjWvlBXZbkhXgKmUCagxU0+Lmmgq9dEKehTEFctk9uAS5VS9yilnlRK3QNcijYd0S2IyCUickxECkXkziGe/08ROSwi+0XkPRFJH/Bcn4jsdWyvDz7XGyi7nbTmfE6ELsLHNrW6zKzOjuXQyRaa2o3Z52lI0hwJEkZ2dfWnBWuWSV5RPVlxISRGuGkEdGSaNgrYqbQMyursGDp77Owt89TUcAt34KoyCQb2D9p3AHBDFBBExAY8gqagZgPXisjsQYftAZYqpeYDLwO/GvBch1JqoWO7HB2oKD5EAvX0pK3TY3ldWZ0dg1KYa/7EtAXgG2hsV1d9oVZYGJVBT5+dnY54iduw+UFUxidKy6CsyIrBRzBfCvoUw1Vl8mfge4P2fRd41k1yLAcKlVLFSqlu4EXgioEHKKU2K6Wc05h2ACluWtstnNyjxUuSFk2deImT+SmRBPvb2G4mV5dvACQvMb5lEqnNfT9Q2Uxbd597UoIHEpNjeMskIsiPuckRVtzE4IxUtLhRRN4RkXeA+cD9IlIiIh+ISAlwP7DATXIkc2bNSoVj33B8DRjYyiVQRPJFZIeIfG64k0TkFsdx+bW17vWV+5ZtpYboSTWi11X8fX1YmhFtvl+Oaaugah90t+ktydDUF/Vncjm/SFdkRrt3jZgcaCgyfOPLVVkx7ClvpKPbipsYlZEsk63ANse2BU15PA1sctw+4NjvDobKcxyyPFlErgeWAg8O2J2mlFqK1j/sYRHJHupcpdTjSqmlSqmlcXET6Lg6+Lp2O+mn91AWvnjSjeh1lZVZ0RTUaNXZpiFtlZYaW5E/+rHeRimtqDJaeyvvKK5nRkIYMaEB7l0nJht62rXZJgZmZVYMPX2K3WUmat0zxfAd7gml1E+9KEcFkDrgcQpwVvK7iFyANqjrU0qp/m8tpdRJx22xiLwPLAK8ZruXFewnnSaK0yZ/P67hcNYDfHSigU/PmzbK0QYhdZl2W/6R1k3YSLTWQHcrxGTT02cnv6SRq5d6wLPrbCBZX6j1LTMoSzOi8BFNqTrHH1gYi7HUmdhEZKaIrBORc5ybm+T4GMgVkUwR8Qe+CJyRlSUii4A/ApcrpWoG7I8SkQDH/VhgDXDYTXK5xKl9Ws/LaQsn7/yS0ZiXHEGwv60/fdUUBEZoX6ZVe/WW5GwaHL+ForPZX9FMR09fv8J2KwOViYEJC/RjXnKEud5fU4xhLZOBiMhi4B+AYzwb4rjtAybcN0Qp1Ssi3wbeRpve+JRS6pCI3AfkK6VeR3NrhQJ/d1T/ljkyt2YBfxQRO5py/IVSyqvKxLc8T4uXZE3dKcZ+Ni1uYroP+7SFxszo6q8xyWbHPu1vutzd8RLQRgH7BRs+CA+a9fv0thI6uvsI8p9a6fdmwFXL5GHgVSACaAHC0ayEG90liFLqTaXUdKVUtlLqfse+HzsUCUqpC5RSCYNTgJVSeUqpeUqpBY7bJ90lk0ty2+2kn949peMlTlZmRXO8upV6M8VNkhZqo2vbDJaJVl8IPn4QkcrOEw2eiZcA+PhocRmDWyagKZPuPjt7rLiJIXH1228ecKdS6jQgSqlW4L+wJi1SXrifWJroS506/biGY0XmJ3ET0zDNkZBoNFdXQ5FWX4IP+SUNrMjygFXiJCb7E7eagRkYN7EwHq4qk54B95tFJN6xL9H9IpmLqn3vApC0cPLPLxmN+SkRBPmZLG7iVCYnDaZM6oshJpsDlc20d3soXuIkJgcaS6CvZ9RD9SQsUKs3MVVx7BTCVWWyC3BGl98HnkMrLBxcFT/lsJXlUUvUlKwvGYwWN4ky14c9MAKis4xlmdjt/WnBTsXskXiJk5hssPcavkcXaK6uveVN5hoVPUVwVZl8HdjnuP+fQBnQBdzkCaHMgjNeUhpmxUucrMyK4Vj1aXPFTaYthJP7Rj/OW5yugt4OLfhe3MD0hFBiPREvcWKSjC7Q4nLdfXar3sSAuPQNqJSqVEqVO+7XK6VuVkp9USl1zLPiGZuKogPE0UjfFK4vGczAehPTkLQQmsug3SAyO77Ue6OyyC9p8KyLC0ylTJZmRDviJgb5X1n0M1I7lRUD7q8ebvOOmMbk5F4tXjJtwfk6S2IcnHGTnWZSJv1xkz36yuHEEQw/2h1He3dff2KDxwiOhqAoUyiT8EA/5iRZ9SZGZKQ6k3eBMMf9rcMco9DqQqYktvI86ogkNWe+3qIYBmfcxFRN+QZmdOUY4IdBfRH4BrKlWnNteTSTy0lMjimUCWiurmfzSuns6SPQb8p+/RiOkcb2hg247zPMNmX/k8puJ61lN6Vhi6x4ySBWZEZzrPo0jW0mmW8SFKW1Yq8ySNykoRiiMtlZ0khOvIfjJU5M0D3YyYpMrd5kb7k138RIjPotKCK+IvKGiLhpIs/k4GTJEeJpoNeqLzmLFc64iZnmwk9baJz04Poi7NFZ5Jc0ur9L8HDEZENLpXE7KA9gWWY0IrDTipsYilGViVKqF1gC9HpeHPPgjJckzjOAW8RgzE+JIMDXx1wf9qSF0FSqfxDe3geNJ6gLSKW1q7dfMXscZxC+odg7602AiCA/ZiWGs/OEiVypUwBX/TPP4cYRvZOC0m00Ek7ajEV6S2I4AnxtLEqLNNeHvT9uorOrq7kc+ro52q2NSPCeZeJQJnUF3llvgizPjGZ3WSPdvcaewzKVcFWZLAZ+5ZjP/q5zaJZjcNaUJKV5NydCFljxkmFYkRnD4aoWWjqNXVXdz7SF2q3eysQRt9jRFEVGTDAJ4V7yLkdnAwJ1x72z3gRZmRVNZ4+dA5VW3MQouNQ1GPjQsVkAVaXHmEYtpSlTumZzRFZkRaPeg/ySBs6bmaC3OKMTHA1hSVBzRF85HG6md06FsGKul1xcAP7BWieA6kPeW3MCLHekS+8obmBJupesN4sRcUmZeHlQluGp2Pse04C4uVa8ZDgWp0XhZxN2FptEmQDEz4RanZVJfRF232AKW0O51RspwQNJmA01Xp3eMG6iQ/yZnhDKzhMN3LZeb2kswDjDsUyFKtlKCyFkzl6mtyiGJdDPxoKUSHaYqXgxbibUHtd3HnpDEY2BKYB4th/XUMTP0Syj7nbvrjtOlmdGs6ukgd4+K25iBFxSJo7hWEVoEwzfd2yb0QobpxxJTbspDp6Pj23Kltm4xIqsaA5WNtPWZZJEwLiZWk+splL9ZKgvokQlkhwZREpUsHfXTpgNyg61R7277jhZkRlDW3cfh0626C2KBQYajiUil4jIMUeQ/84hng8Qkb85nt8pIhkDnvuhY/8xEbnYXTINRe3JElJUFZ1JKz25zKRgRWYMfXbFrlKTNOWLn6Xd6vVl2teLaiplX3uMd6reBxPvmBRqEleX829kqqzBSYwhhmOJiA14BLgUmA1cKyKzBx32NaBRKZUDPAT80nHubLSZ8XOAS4BHHdfzCKV7tHnvMXPO89QSk4Yl6VHYfMQ8H/bY6dqtXsqkuQyx93KkO46Vnu7HNRTRmeAbBNXmUCbxYYFkxYaYq55pEmOU4VjLgUKlVLFSqhttVsoVg465AnjWcf9l4HzRhsFfAbyolOpSSp0ACh3X8wh9xVtpVUFkzrUsk9EICfBlXnKEeT7sQZGOjC6dlEm9lslVYk/UxzLxsWlJCDXmyOgCzTr5qKSBPrvSWxRTULDnQ3Y/+Fkqi92faGKU4VjJQPmAxxWOfUMe46jKbwZiXDwXABG5RUTyRSS/trZ23MIeD1uOr5//uM+fSqzIjGZfhYmGGemZ0eXoFtwamkZatJfjJU7i55jGMgEtCH+6s5ejp6y4iSvU7X+LxW0fEhwW6fZrj3c4VinuHY4lQ+wb/FNjuGNcOVfbqdTjSqmlSqmlcXFxYxRRY8Xtz7L4B6+P69ypyIqsaHr6lHmGGemY0aXqC2kjkNzMbDSjWwcSZkNbDbSO/8eWN3G25zeN9aszIVU7KfFJIypumtuvPd7hWLe4eThWBZA64HEKcHK4Y0TEFy0ZoMHFcy10wjnMyDQfdh0zujpOFXDCnsiKbB3iJU7iHaFKk7i6kiKDSI0OMk9cTkd6e7rJ7jhIddRij1zf1dTgQhH5kYgM6T5yAx8DuSKSKSL+aAH1wT//XwducNz/ArBJKaUc+7/oyPbKBHKBjzwkp8UYCQ/0Y3aSiZryxc3UbnUIwvfVaWnBHh+GNRIJc7VbE7m6VmTG8NGJBuxW3GREig9sJ0Q68c1a55Hru+rmuh+4CDghIv8WkatExM9dQjhiIN8G3gaOAC8ppQ6JyH0icrnjsCeBGBEpRHO13ek49xDwEloNzFvAbUopkzjopwbLM2LYU9ZEV68J/i1xM7RbbyuTvh6C2yup8U0iOy7Eu2sPJDQOQuJMY5mAFpdrbO+hsLZVb1EMTcPhTQCkL7nII9d31c31tFLqXGAWkA88CFSJyP+6SxCl1JtKqelKqWyl1P2OfT9WSr3uuN+plLpKKZWjlFqulCoecO79jvNmKKX+7S6ZLNzDiqxounrt7K9o1luU0dEro6upDBt9+Mbn6BcvcRI/23SWCcBOa5TviARW7qBckohNTPPI9cfU8lYpVaSUugdYCezEaktv4QLLMxzFZWb5sMfN8HpGV22p9uUdmza4vEoHEuZoDS/tJrAkgdToIKZFBJqrdY+X6evtJatjP1VRSzy2xlh7c10hIq8BJUAk8A1PCWYxeYgK8WdmYhg7zfJhj5/l9YyuyqKDAOTOXuC1NYclfraWhNBYorckLiEirMiMZmdxA1oY1WIwJYc/Ipx2fDLWeGwNVwPwv0XLkPo9WkxjvlJqjVLqTx6TzGJSsTwzml2ljfSYoSmfDhldrVXHOU0Q2WkZXltzWBIc1pFJ2tGDNiq6rrWL4jrjjx3Wg9qDWrwkZdGFoxw5fly1TFLR+nClK6V+qJQyxwQdC8OwIjOG9u4+DlaaIG6iQ0aXX/MJGvxT8LEZYNha3CwQHzh1QG9JXMY5kfIjs1i/Xsa/cgcnJYHE1ByPreFqAP4qpdS/lVIm+FlpYUSc7dRN4epyZnR5aVDWqeZOEnsq6Y3M9Mp6o+IfrCmUk7v1lsRlMmNDiAsLME9czosou53Mtn1URnimvsSJAX4GWUwF4sICyI4LMceHPSgSQhP625t4mo+KqkiRWsKSZ3plPZdIXgSVu8EkMQgRbf7LzhNW3GQwpcd2E0ULKt1z8RKwlImFF1meGUN+SaM5mvJFZULDCa8sVXD0MDZRxKbN8sp6LpG0GDoa9J3tMkZWZkZT1dxJeUOH3qIYiur97wGQvOACj65jKRMLr7EyK5rTXb0cNsMwo2jvKZO6cs2d5hPrOX/2mEl2uEQqzePqWpnlmAtvlm4LXsK3PI9TxJKUMcOj67iazTXkeF4R8UxdvsWkpP/DbgZXV1QmnD4JPZ79lVvT0klQS4n2IDrLo2uNifg5YPM3VdwkJz6UmBB/c7y/vISy28lo3UN5xGLEx7O2g6tX/9cw+ze4SxCLyU9CuDbMyBQf9mhHMLzRs26eHScayJBT9PmHQ7COPbkG4+sPifOgco/ekriMiLAyK8aqNxlA2fG9xNAM6Ws9vparyuSs/g4iEgZY2V0WY2JFltaUz/BxE6eV0OhZV9eO4npybNX4xGSD3m1UBpO0GKr2mqYSHrTWPZVNHVQ0WnETgFP73wUgaaFn4yUwijIRkQIROQ4EicjxgRvaQKrNHpfQYlJhmrhJlMMy8XDcZEdxPbl+NUhMtkfXGRfJi6G7FeoK9JbEZZyu1O1msH69gG/ZNqqJISnD88kdvqM8/3M0q+QxtM7BTuzAKWCTh+SymKQ4P+w7T9QzLyVCZ2lGIDgaAsKhoXj0Y8dJzelOKmqbiAmsASMqkyRnEH6XNoHSBOTGhxId4s/O4gauXpo6+gmTGGW3k966l5LwZSR4OF4CoygTpdSzACJyVCm1w+PSWEx6EsIDyXTETb6+zkAB58GIQFSGR91cO4sbSJUafLBDtAGVSWwu+IdqQfhFX9JbGpfQ4ibR5ojLeZiygv2k00SxF+Il4HoF/A4RCRSReSKyeuDmaQEtJh8rs7TiMuPHTTybHryjuJ5Z/o7xuEa0THxsMG2hqdKDQWvdU9nUQXlDu96i6MqpfRsBSPJwfYkTV1ODLweq0ObAbx2wbfGcaBaTlZVZMZzu7OVIlQniJk1lHgtA7yiuZ210k/bASGnBA0leBNUHobdbb0lcxlQp6B7Et3wbNUSTnOWdsQauOtJ+A/wUCFVK+QzYbBMVQESiRWSjI9i/UUSihjhmoYhsF5FDIrJfRK4Z8NwzInJCRPY6toUTlcnCsziHGRn+wx6dBfYeaK5w+6VrTndSVNvG/KB6CIrWYjRGJHkJ9HVrCsUkOOMmO4pN0AfOQyi7nfTTeygL93x9iRNXV0lQSj2slPKE3Xgn8J5SKhd4z/F4MO3AV5RSc4BLgIdFJHLA8/9PKbXQse31gIwWbiQxwhk3MfiHvb/WxP2uLmd321R10pguLifOILyJihd9fBzzTaZwJXx54X5iacLupXgJuK5M3hGRlR6S4QrgWcf9Z4HPDT5AKXVcKVXguH8SqAHiPCSPhRfQ4ib1xo6b9KcHuz+ja0dxPaEBvoS2lRoz+O4kMg1C4qFsp96SjImVWTFUNE7duEnVPq2+ZJqX4iXgujIpAV4Xkf8TkR8N3NwgQ4JSqgrAcRs/0sEishzwBwa2dL3f4f56SEQCRjj3FhHJF5H82tpaN4huMV6ccRND15uEJ2ktRTwQhM8rqmdtejDSYnDLRATSV0PpNtN0EAZYlT216018S7dQTQwpWXO8tqarymQJcAiYC1w4YHNJ7YnIuyJycIjtirEIKyLTgOeAmwbMVvkhMBNYBkQD/z3c+Uqpx5VSS5VSS+PiLMNGT1Y5gqR5RXU6SzICPjaITHe7m+tUcyfFtW1ckOj41WzU4LuTjLXQUqklI5iE3PhQYkP92V409ZSJva+PzNbdlEUs9Vq8BEYvWgRAKbV+IosopYZVOiJSLSLTlFJVDmVRM8xx4cAbwN0Da16cVg3QJSJPAz+YiKwW3iE+PJCc+FDyiur5xqcM/Ms8OgsaStx6ye3FmgJdHt6o7TCyZQKaZQJQmgdR6frK4iLOPl15RXUopRCjtarxIKVH88mkhcIM7/bhdVltiYjNUVtyjeNxsIgEuUGG14EbHPdvYIjmkSLiD7wK/Fkp9fdBz01z3ApavMU8aSdTnNXZMXxc0mDsufDRmZpl4kYXT15hPRFBfqQ4fwcZOWYC2tTFwEgo3aq3JGNidXYs1S1Tby58taO+JHXxRV5d19U6k2y0L+k3gScduy8CnnCDDL8ALhSRAjTX2S8cay4VkT85jrkaOAe4cYgU4OdF5ABwAIhFawFjYQJWZ2tz4fdXNOktyvBEZWr9qdrc547LK6pnVVYMPg1FEBIHgeFuu7ZH8PFxxE3y9JZkTKzOdrpSp5arK6BiGxWSyLR0z84vGYyrlsn/AS+ixSR6HPveByZsRyml6pVS5yulch23DY79+Uqprzvu/0Up5Tcg/bc/BVgpdZ5Sap5Saq5S6nqlVOtEZbLwDisyYxDRfqkblmj3ZnSVN7RT2dTB6pwY7ZpGt0qcpK/R5G2pGv1Yg5AeE0xSRCA7ppAy6evtJbt9Lyejlnl9bVeVyXLgfkfQWwEopZqAyBHPsrAYgagQf2Ylhhv7l2OUe2tNnAkHq7JioL7I+PESJ864SZl5rBMRYVV2LNuL67EbOQXdjRQfyCOcdnyyP+X1tV1VJi0MUhwikgRUu10iiynF6uwYdpU10tlj0JkZUemAuC09OK+ontjQAHIigdZTxs/kcpI4X2v6aEJXV0NbN8eqT+stileoP6DFSzKWXOL1tV1VJv8AnhKRFAARiQEeRnN9WViMm9U5MXT32tld1qi3KEPjG6DVmzRNfOKiUoq8onpWZ8cgTuVkFsvE5gupK6Bkm96SjIlVUyxuEnwyjxKfNGITvd9+31Vlcg/QCpShWSg1QBfwgIfkspgiLMuIxuYjxq4HiEhxS3+uotpWak93aYHhBkfNrVliJgAZa6D2CLQZ+H81iKTIIDJjQ9hu5HomN9Hd1UlOxwGqY5brsr6rLeg7lFLXoVWnLwcSlVJfVkp1elQ6i0lPWKAf81MijP3LMTxZK9qbIM7XuDo7VouXgHncXKAF4QHKtusrxxhZla3Nhe81cgq6Gyja+wHB0oV/7oTKAsfNmMojlVJ1jiwrqxeJhdtYnR3DvvImWrt69RZlaCKSoblywrUm2wrrSI4MIjU6SMuMCk2EgFA3CekFkhaBbyCUmK3eJIbTXb0cqGzWWxSP0nToXexKyFpyoS7rD6tMHG3di0fbvCmsxeRkTU4svXZl3BTOiFTo65pQrUlvn13rx5UTq1Vj1xVokwzNhG8ApK2C4s16SzImVmfHIgJbCya3qyuyaiuFfrlExCTosv5IlsndaLGSe4CngSDgr2iz4F8AAoGnPC2gxeRnSXoUQX42thYa9MMenqzdtow/brK/spnTnb2szY3VdtQXmif4PpDs86D2KLSc1FsSl4kO8WdOUjhbjPr+cgOnmxvI7T5KfYJ+w2+HVSZKqeedG1pl+meVUncrpZ5UW708ggAAIABJREFUSt0FXIZWBW9hMSECfG0sz4zmwwKDek8jHMqkefxxk60FdYhoVhjtDdDRADE5bhLQi2Q7/PHF7+sqxlhZmxPHnrJG47pSJ0jRR2/hK3bCZ+v3lexqzGQhMHjo1H7HfguLCbMuN5bi2jZONnXoLcrZRDjSLCeQ0bWloJY5SeFEh/hrVglAjMncXADxc7QWMEWb9JZkTKzLjaWnT7Fzkrak7zq2kTYVQM4SfYLv4LoyOQZ8b9C+7wLH3SuOxVRlXa42EsCQfu3gGC3wPE43V2tXL3vKmvpf4yfKxISWiY8PZJ2rWSZ282RHLUmPItDPhy1GfH+5gaT6HRQGLyAgMFg3GVxVJrcB/y0iJSLygYiUoI3Xvc1jkllMKaYnhBIfFmBMv7aIVrg4TjfXjqJ6eu2KdTkD4iU+vqZp534W2edBWy3UHNJbEpcJ9LOxPDPGuHG5CXCqrIBUdZKO1HN0lcPVOpOPgCzgLrSZIncB2Uopc83ytDAsIsLanFi2FdYZs4/SBAoXtxbWEejnw5KMKG1HXQFEZYDNz33yeZOsc7XbInNlda3LiaWwppWqZgO6UidA2cdvAJC46FJd5XC5zkQp1eIIyP/KcTu5k7YtvM7a3Fga2ro5XGXAUb7hKeMuXNxSUMvyzBgCfG3ajvoic7q4nIQnQdxM08VNnJl0hnSlTgDbic3UEE36jMW6yuHqPBMRkWtF5Dci8vjAzdMCWkwd1jrcQIb0a0ckw+kq6BtbNtDJpg6Kats4x5kSbLdrrVTMrEwAstZrlfA95mmCMTMxjNjQAGO+v8aJva+PrNZdlEYs9+qI3qFwdfXHgN8DqYDfoG1CiEi0iGwUkQLHbdQwx/UNGIz1+oD9mSKy03H+3xxTGS1MSHx4IDMTw9hixBThiBRQdk2hjAHnr+D++pKWCujtNL8yyT5Pex0maq0iIqzLNbArdRwUHcgjitNItn5ZXE5cVSZXAcuVUlcrpW4auLlBhjuB95RSucB7jsdD0TFgMNblA/b/EnjIcX4j8DU3yGShE2tzYskvaaSj22At6cNTtNsxurq2FNYRFxbAjIQwbYeZM7kGkrEGfPygYKPekoyJtTmx1BvVlToO6vf9G4CM5Z/RWRLXlUk7WsdgT3AF8Kzj/rNoc9xdwjH3/Tzg5fGcb2E8PjUjju4+O9uLDeaK6C9cdD0I32dXfHi8lnNy47QWKvBJg0ezKxP/EJhxKex9Hrrb9ZbGZdZN1yzED44b0PodB+GVH1Boy9al5fxgXFUmvwJ+LP2fCLeSoJSqAnDcxg9zXKCI5IvIDhFxKowYoEkp5XRkVwDJHpDRwkssz4wmyM/GpqM1eotyJuFjVyZ7yhpp7uhh/cy4T3bWFWhDpsIS3SygDqz8FnQ2wX7zjDWKDwtkXnIEm432/hoHzQ21TO86TG2i96cqDoWvi8d9B0gHbheRM/4LSqnpo50sIu8CQ3167nJxfYA0pdRJEckCNonIAbQJkIMZ1hkqIrcAtwCkpaWNYWkLbxHga2NNTiybj9ailMIzv1/GQWA4BESMyc21+VgNNh/5pFgRPunJZZTXNRHSVsG0BbDjD7DkJtO8pvUz4vj95kKa2ruJDDZviLVwxwaWiJ2ohfq7uMB1ZfLziSyilLpguOdEpFpEpimlqkRkGtrgraGucdJxWywi7wOLgFeASBHxdVgnKcCwHeiUUo8DjwMsXbp0ckTgJiHnzYzn3SPVFNa0kuuMNRgBZyt6F9l8tJYl6VFEBA3IU6kvhJSlHhBOB0Rg5a3w6je0NOGc8/WWyCXWz4znd5sK+bCgjssXJOktzrixH3uHRsLIXXSu3qIArhctPjvc5gYZXgducNy/Adgw+AARiRKRAMf9WGANcFgppYDNwBdGOt/CXJw7Q/slv/mYwVwRESnQXO7SoaeaOzlc1cL6GQO8tj2d0FRmzp5cwzHn8xCaADse01sSl5mfEkl0iD/vm9jVZe/rI7t5O0XhK7D5umoTeBaXE5NFZK2jtuSfjsdLRMQd9fu/AC4UkQK07sS/cFx/qYj8yXHMLCBfRPahKY//396Zx0dVnY3/+yQhhJCQBBK2sIQ9LMoOgsgqm1rAolKtaF1q61u1trWt9dVfq69t7e5GVVSUUqwrqNW6gVBXZDEgOwl7IEAIhCUBsj2/P84NDGEmmTAzuZPJ+X4+87kz555773Ny4T73nGd7RFU3OPt+CfxURHIwNpTngyCTxUXaJjchs3UiSzaFmZG0FhUX/7vFPKjOspcc3g5o/Te+exLTGAbdAjkfGXtQPSA6ShjVPY2lW/LrrYvw1jWf0pyj0C18Erf7G7R4HWYGcRKoVCAKPBSoAKpaoKrjVLWbsz3ktK9U1Vud71+o6gWq2tfZPu9x/DZVHaKqXVX1alU9FahMFvcZ3aMlK3Yc4ujJUrdFOUNSOhQXQGnN6Tg+3nSAtklxZ1yCwcMtuB7WMamOQTdDdCx88bjbkvjN6B5pHCoqYU1uoduinBcFWe9QrkKXYVPdFuU0/s5M/heYoKp3AZWpQtcBvUMilaXBM6ZHGmUVyufhFK18OhV99bOTkrIKPss+yOjMlmc7EFS+uUfSzAQgIc0Y4LPmw4FNbkvjFyO7pRElsGRzmM1+/aR53n/ZEptJSmr4eAX6q0zaqupK53vlvLAMiA6+SBYLDOiYQmJcTHjZTfysuLhyxyGKSsrPtpeAiTFp2tJ4hkUao35pYk8W/dptSfwipWks/TuksDSc/n35ScH+3XQv28LhtqPdFuUs/FUmW0Wkaj3I4Zg6JxZL0GkUHcXI7mks2WxchMMCPysuLtl8gNjoKC7u2uLsHQc3Q1qPEAnnMk1bwCU/gy3vw7b/ui2NX4zNbMk3uUfIP1a/Vsa3f2l8jFoO+JbLkpyNv8rkYeAtEbkfaCQiP8PUgQ/YZmKx+GJsj5bkHzvFN7lhkqDaj8BFVWXRxgMM7dyc+NgYzx2QvwVSawzLqr8M/aFZCvzw/npROKty5vjxpv0uS1I7orPf5wDN6dxnmNuinIW/rsFvAtcBQ4GdmBQmN6vqeyGUzdLAGdezJdFRwgfr97ktiiGmsVmmqmaZK/vAcbYfLGJi7ypr2cf2wakjJnV7pNIoDsb9GvZ9Uy+i4nu2SaRdShM+WF9/lMnJ4mP0OL6cbS1GERXtbpbgqtSmnslHqvotVe2jqperav3K8GapdyTHxzK0U/PwUSZQY5GsD9btQwQm9Gp19o6DzopwWgTPTAD6TIf0QfD+vXBou9vSVIuIMLF3az7LPsjxU7UrLeAWmz9/i3g5RXzfK90W5RxqE2fSQUTuFZEnnW1G6MSyWAwTe7dma34ROQeOuS2KoQZl8v76ffRvn0zLZnFn78ivVCYRPDMBUyP+Ksdz/7Ubw77eycTerSkpr6g3ubpK171NIQlkDp3ktijn4G+cySSMsf1yIMnZbnTaLZaQMaG3ecMPm6WIpPZGmXhxCth9qJj1e4+eu8QFRpnEJZlo8UgnJQOufAby1sD7v3RbmmoZ2DGF1ITY8Jr9+qCs5BTdj3zGlqQRxDZu7LY45+DvzORPwC2qeomqzlTVS4Cbgb+ETjSLBdokNaFv++Tw+c+e1A5Ki+HE4XN2fbjBKDyfyiS1R71JhhgwPSbDiJ/Aqhdh9b/clsYn0VHC+F6tWLLpACdLw6yGThU2f/UezSgiuveUmju7gL/KJAOoalF7BbCpdy0hZ2LvVnyTe4Q9hTVHnoecJKdIlpccXR+s30dm60QyUpuee1wkuwX7Ysz9kHEJvHM35H3jtjQ+mdC7NUUl5XyxNYwCZL1QvGYhRdqYnhfXb2WyFBhdpW0UUD8cyi31mknOm/6H4TA7Oa1MzrabHDx+ihU7DjHB26yk+BAU5Tc8ZRIdA1e9AE2awyvXm79DGDK8SwsSG8fwwbowWUr1QkV5OZ0OLmVTwlDim4ZRJm0P/FUmOcBCEZkrIg+KyFzgTWCLiNxX+QmdmJaGTOe0BLq1TAiPpa5kZzJeRZks2rAfVTOLOoeGYnz3RkIaXPMPOLoXFtwWlvEnjWOiGZPZko827qesPPzkA8jJWkIqhZT1uMJtUXzirzLpB3yNWdYa6Wy/xtQUGe98fNYssVgCZWLv1izffoiC4y5HK8e3gJi4c5a53l+/j3YpTejVxkuqlHwnX1UkByxWR/vBMPkPJrPw4t94dV5wm4m9W3OoqIQVO861hYUDhSvfoESj6XHJdLdF8YlfifBVdUyoBbFYquPyC9vw5JIc/rM2j5nDMtwTROQc9+CC46f4LPsgt4zo5L0y5MEt0Cj+TKLIhsigm2HfWvj8MeMuPOkR40YcJozJTCM+Npq31+xlWJcWNR9Qh1SUldF533usix/MgJRUt8XxSW3iTKJFZLiIzHB+x4tIk9CJZrGcoWebZvRolcjCLP8rHYaMKsrk3bV5lFUo0/qne++fvwlSu4XVw7POEYHL/wrD7oDlz8CC70NZidtSnSY+NoYJvVrx7jd7OVUWXl5dW756l1QOU9r7GrdFqRZ/40y6YFLO/4czxacmAM+GSC6L5Rym9U/n612F7CoodleQKsrkzaw9ZLZOpKe3JS4wObkaor2kKlFRMOFhuPRBWPc6/GsGnDrutlSnmdY/naMny1gaZmnpi1e+xDFtQp8xM9wWpVr8fVV6AuMa3ByorFa0FLgkUAFEpLmIfCQi2c42xUufMSKy2uNzUkSmOfteFJHtHvv6BSqTJTyZ0s/U635ztcuzk6T2JtdWWQk7C4r4elchU/v5mJWcPGpyeTVUe0lVRGDE3TDlSdi2FP4xBYoK3JYKgBFdU0lNiOXNcJj9OpQUHyPz8BLWJY+hadMEt8WpFn+VyRDgt6pagVPPRFULgeQgyHAvsFhVuwGLnd9noapLVLWfqvbDJJksBj706PLzyv2qujoIMlnCkPTkJgzt1Jw3s/a4m5Y+qR2gcGwvb2btBWCqo+jOobIglp2ZnM2AmTBjPuxfDy9MgsJz43bqmpjoKK64sC2LNx3gyInwqPC56b8vE88p4gZd57YoNeKvMjlKFcUhIm2BYDhmTwXmOt/nAtNq6H8V8J6qurzWYXGDaf3T2XawiLV7XExL78SaaOEu3lq9h6GdmtM22Yf58HSCxwYWY+IPmZfB9Qvg2H6YMzEsqjRO659OSVkF76/Lc1sUAKLXvkoeqVwwbLLbotSIv8pkATBHRNoBiEgL4FHOjYo/H1qpah6As21ZQ//vYGqpePJbEflGRP4mIj6T1ojIbSKyUkRW5ueH17qoxT8u69OG2Ogodw3xjlfW7u3ZbDtYxJW+DO9gjO9RjSClUx0JV8/IuBhuehcqyswMZfdyV8Xp2y6JTqlNT8843eRYwR4yi1aQ0+oyYmL8crx1FX+VyQPAcWAXZoZyADgF/M6fg0VkkYis8/KZWhthRaQNcAHwgUfzr4BMYDDGpuMzs5yqzlbVQao6KC0trTaXtoQJSfGNGJOZxr/X5LkXYNbMLGnl5GwkNjqKyRe08d03f7Op+R4d/g8D12h9AdzyITRJgblTYMsHNR8TIkSEqf3asmx7AXlH3E3fk7N4LtGipA6f6aoc/uJvcawTqnodkIaxn7R2Ej76lV9aVS916qBU/bwF7HeURKWyqC4X9DXAQlU9vaCpqnlqOAW84MhniWCu7N+Og8dPscQtr5tGTdCmaRTmbWdsZkuSmjTy3XfvamjTt+5kq6+kZMDNH5p6Ly/NgI9+7Zrr8LR+6ajC6yt9lxoIOaq02PIym6O6kHnhYPfkqAW1iTPpDnQHjqpqMP8Xvw3c6Hy/EXirmr7XUmWJy0MRCcbesi6IslnCkHE9W9KqWWP+uWynazIUNmpJi7IDfPeianKdHs2D4/ugrXUw9IuENLjpPRhwA3z+KMyZAAVb61yMjNSmjOiaykvLd7k2+9256n06lO1kX/eZ3gNhw5AalYmITBWRXcBG4HNMHZNdla65QeARYLyIZGPSsjziXHeQiDznIUcG0J5zk0vOF5G1wFogFVOv3hLBNIqO4tohHfjvlnx2HCxyRYZNJ5LoGHOIi7tUE5Gc5zgWtrHKxG9im8KUx+GaeaZS49OXwNfz6jwFy8xhHck7cpLFLhXNKvpkFoc0kb6X3eLK9c+HapWJiAwAXsO44U4EejnbD4FXRWRgoAKoaoGqjlPVbs72kNO+UlVv9ei3Q1XTHfdkz+PHquoFzrLZ9aoaPlFQlpBx7ZAOxEQJ87+q+9nJuj1H2FCURHpUAVHVvTTuXQ0SZWwCltrRawrc/gWkD4C374DXvue1hkyoGJfZkrZJccz7su7/fR3bt5UeRz4jK20qyc18BMKGITXNTO4G/qCqt6rqIlXd5GxvBX7v7LdY6pxWzeKY2Ls1r67M5URJ3aa/mPflTvKj0mhUfqL6B1zeahOs2Di8g83ClqR0uOEtuPQ3sOkdeGoE7Pi8Ti4dEx3FdUM78FnOQbbm1+376Y73HkMR0i+9o06vGyg1KZPhwGwf+5519lssrjBzWEeOnCjl39/UnRvnkeJS3lqzh3adujkN1Rhp92bZJa5AiYo2FRtv+QhiYmHuFfDxw1Ae+qDCGYM70Cha6tQ2pyVFdNz5BssbDyMzs2edXTcY1KRMWqiq19BUVc0Fwiu9pqVBMbRTc7q3SmDelzvrLCL+tVW7OVlawbD+jpLwpUyO5sHx/dC2f53IFfGkD4AffAr9roNP/gRzJsGBjSG9ZFpiYyb3acPrq3IpLikL6bUq2br4RZpxnNJBt9XJ9YJJTcqkpv31w83AEpGICDMv6sjaPUdYtTP06+ll5RXMW7aTgR1T6NLVeWv0pUz2Zpmt9eQKHo0TYOosU73x0DZjnP/4tyalfYi4YVhHjp0s441VdeAmXF5KwqpZbCaDoaPDtwiWL2pSFo09KylW/QCxdSGkxeKL6QPbkZoQy6OLskN+rYVZe9hZUMwPRnaGpqkQ3dhrLXjA2Eus8T009Pk23LEC+kyHT/4IT18MWz8OyaUGdkxhYMcU/r50KydLQ2uby13yLK3L9rC5553Exda/INealMkyzlRS9PZZFlLpLJYaiI+N4YejuvBZzkG+2ha67LOl5RU8/nE2F6QnMb5XK5P9Nrm9eUP2xt7VkNrDuLpagk/TVPj2Mya3l1bAvCvhlZlBTxgpIvx0fHfyjpzklRUhTEZZUkzTL/9MFj0YM+WG0F0nhFSrTFR1tKqOqe5TV4JaLL64/qKOpCU25m+LtoTsGm+symX3oRP8dHz3M0Fk7QbDri/PrWuuapa57BJX6Ok6Dm7/EsbeD9kfwawh8MmfoSx45Z2Hd2nBkE7NmbUkJ2Szkz0fPkpKeQHb+/6cxCb1c8GnAZd+s0QKcY2i+dHoLizbdogvth4M+vlLyip44uMc+rVPZnQPj5xunUZBcQHsr5J04VgeFB2wxve6olEcjPw53LEcul4KH/8f/P0io1yCQOXs5MCxU8z/aldQznkWJw6TvGoWnzKACZddGfzz1xFWmVgigu8M6UDrZnH87aMtQffsenXlbvYUVpmVAHQeZbbblp59QKXx3boF1y3JHWDGPLP0JdEw/yr45/SgeH1d1LkFw7u04KmlOUH37Nr37u9oUlHEvsG/IKFx/bOVVGKViSUiiGsUzY/GdmXFjsN8uCEYZXYMR0+W8sTH2QzqmMIl3aqkTmnW1thFtlfJ8GMj392l6zgTPT/hYdi9Ap4aDv++G44GFo/0k/HdOXi8hNmf+LCTnQcVO76k5brneCdqNJePHx+087qBVSaWiGHGoPZktk7k/jfXUVgcnIyzD7+zgfxjp3jgil7eE+51HgU7vzh7jX73MlNZMTY+KDJYzoOYWBh+J9yVBYO/D1nz4LF+8J9fmBig82BwRnO+1bcts5bksGnf0cBlPFFI0cs3sbsilfIJvye+HnpweWKViSViiI2J4s9X9+VwUQkP/XtDwOdbuvkAr67M5QejutC3vY8K1Z1HQ2kx5K4wv/ethe2fQO/6u/YdUTRtAZf9Ee5cBRdeAyueg8f6wnu/PC+l8uCU3iQ1acQ9r62hNJCMwqoULbiTuBMHmJf+ANMuqv9lna0ysUQUfdKT+J8xXVmQtYdFASx3HT1Zyq8WrKVbywTuvrSb744ZI8yS1jZnqevTv0JsIgz5/nlf2xICUjJg6pOOUrkalj8Lj/czSqW6lDhVaN40loen9WHdnqMBLXdVfD2PptlvM0uu4dbvXFNv0sxXh1UmlojjjjFdyWydyH0L13K4qPbLXarKw+9sYP/Rk/zp6r40jon23TkuCdIHGiP8wRxYvxAG32KqBlrCj+adTBT9navggquMUnmsLyy83W9D/aQ+bbjiwjY8umgLG/POY7lrzSvwzt18Xt6b9MvvpXVSXO3PEYZYZWKJOCqXuwqLS7lhznKOnKhdUsBHF2Xz6spcfjiqC/18LW950mkU7FllXFJjGsOwH52n5JY6o1Kp3JUFg2+FDW8ad+IXr4C1r9eYouXBKb1Jjo/lxjnL2V6bmjrLn4WFt7GsvAcvdfo9Vw3qGOBAwgerTCwRSZ/0JJ6eOYBN+45yw5zlHD3pn0J5YnE2jy3O5uqB7bhnQg//LtZ5NGi5eSD1nwkJLc9bbksdk9IRJv8BfrIexj4AhTvhjVvgr5nw/n2Qv9nrYS0SGjP/1qGUVSjXzl7GzoIaFErpSVj8EPznHj6qGMiTbX7PH797cUQsb1XiujIRkatFZL2IVIjIoGr6TRKRzSKSIyL3erR3EpGvRCRbRF4RkfoZPmoJOmMzWzHrugGs33OE781ZTsFx31HR5RXK44uz+ctHW/h2/3QemX4hUdVWvvKg/RCIaQJRMXDxXUGS3lKnxDeHkffAXWtg5kLoNBKWP2Mi6p+fCF/OMstgHjFM3VslMv/WoZwsK+fa2cvIOeCl7okqbPw3+veh8OlfeL18FLNb/ZrZN4+gaT2OKfGG1FXqbp8CiPQEKoBngHtUdaWXPtHAFkw+sFxgBXCtqm4QkVeBBar6sog8DaxR1adquu6gQYN05cpzLmWJQN5fl8cdL2URGxPFTRdn8P1LOpMcb945KiqUd9fm8eiiLWzNL2Jqv7b89Zp+RPurSCr54H+hUROT1sMSGRzPhzUvweqXIH+TaUtsY5wuOg6HDsOheWfW7T/B9c9/xbGTZUwfkM6dY7vRPuYIbHwb1r4GuSvYHdOBXxV/l+PplzDvliEkxjVyd2wBICKrVPWcF3/XlUklIrIU38pkGPAbVZ3o/P6Vs+sRIB9oraplVftVh1UmDYut+cd5dFE273yzl/hG0bRNbgLA8VNl5B05SfdWCfzk0u5M7N3a/xmJpeFQuBu2LTHZiXd+YWrVVNK4GeVxKRSciuLIiVKiKSdD9hGFsl06MKdkLIuaXMbt43owY3D76h066gG+lEl9mWelA54pO3OBoZjiXIWqWubRnu7rJCJyG3AbQIcOHUIjqSUs6ZKWwBPX9ueOMV158Yvtp43yIsKEXq244sK2tZ+NWBoOye1hwA3mo2qyRe9aZqLqiw8SXVxAy/ISmpVWsP3gcd6VSWQljmRfbEcGdmzOkqEdiGtUv5VITdSJMhGRRUBrL7v+V1Xf8ucUXtq0mnavqOpsnDLEgwYNCo8pmaVO6dE6kd9/+0K3xbDUZ0SgRRfzqUIc0NP5fKuu5XKZOlEmqnppgKfIBdp7/G4H7AUOAskiEuPMTirbLRaLxVKHuO7N5ScrgG6O51Ys8B3gbTUGnyXAVU6/GwF/ZjoWi8ViCSKuKxMRuVJEcoFhwLsi8oHT3lZE/gPgzDruAD4ANgKvqup65xS/BH4qIjkYG8rzdT0Gi8ViaeiEjTdXXWO9uSwWi6X2+PLmcn1mYrFYLJb6j1UmFovFYgkYq0wsFovFEjBWmVgsFoslYBqsAV5E8oGd53l4KibGpSFhx9wwsGOOfAIdb0dVTava2GCVSSCIyEpv3gyRjB1zw8COOfIJ1XjtMpfFYrFYAsYqE4vFYrEEjFUm58dstwVwATvmhoEdc+QTkvFam4nFYrFYAsbOTCwWi8USMFaZWCwWiyVgrDKpJSIySUQ2i0iOiNzrtjzBRkTai8gSEdkoIutF5MdOe3MR+UhEsp1tituyBhsRiRaRLBF5x/ndSUS+csb8ilP+IGIQkWQReV1ENjn3e1ik32cR+Ynz73qdiPxLROIi7T6LyBwROSAi6zzavN5XMTzuPM++EZEB53tdq0xqgYhEA7OAyUAv4FoR6eWuVEGnDPiZqvYELgJ+5IzxXmCxqnYDFju/I40fY0ocVPIH4G/OmA8Dt7giVeh4DHhfVTOBvpixR+x9FpF04C5gkKr2AaIxtZEi7T6/CEyq0ubrvk4Gujmf24CnzveiVpnUjiFAjqpuU9US4GVgqssyBRVVzVPVr53vxzAPmHTMOOc63eYC09yRMDSISDvgcuA557cAY4HXnS4RNWYRaQaMxKn/o6olqlpIhN9nTHXZJiISA8QDeUTYfVbVT4BDVZp93depwD/UsAxTubbN+VzXKpPakQ7s9vid67RFJCKSAfQHvgJaqWoeGIUDtHRPspDwKPALoML53QIodAqzQeTd685APvCCs7T3nIg0JYLvs6ruAf4M7MIokSPAKiL7Plfi674G7ZlmlUntEC9tEelbLSIJwBvA3ap61G15QomIXAEcUNVVns1eukbSvY4BBgBPqWp/oIgIWtLyhmMnmAp0AtoCTTHLPFWJpPtcE0H7d26VSe3IBdp7/G4H7HVJlpAhIo0wimS+qi5wmvdXTn+d7QG35AsBFwNTRGQHZulyLGamkuwsh0Dk3etcIFdVv3J+v45RLpF8ny8FtqtqvqqWAguA4UT2fa7E130N2jPNKpPasQLo5nh/xGKMd2+7LFNQcWwFzwMbVfWvHrveBm50vt8IvFXXsoUKVf2VqrZT1QzMPf1YVb8LLAGucrpF2pgec9yLAAAGfUlEQVT3AbtFpIfTNA7YQATfZ8zy1kUiEu/8O68cc8TeZw983de3gRscr66LgCOVy2G1xUbA1xIRuQzz1hoNzFHV37osUlARkRHAp8BaztgP7sPYTV4FOmD+U16tqlWNfPUeERkN3KOqV4hIZ8xMpTmQBVyvqqfclC+YiEg/jMNBLLANuAnzghmx91lEHgRmYLwWs4BbMTaCiLnPIvIvYDQm1fx+4NfAm3i5r45SfRLj/VUM3KSqK8/rulaZWCwWiyVQ7DKXxWKxWALGKhOLxWKxBIxVJhaLxWIJGKtMLBaLxRIwVplYLBaLJWCsMrFYfCAiHUTkuIi0dVsWbziyDQu0T7ARkYki8mkAx/d2Mhk3DqZcltBilYnFdURkqYicEpFjInJERLaJyDwRGeimXKq6S1UTVDUsI6Id2b4EEx8jImXV9akLnLiFv2FiGyrbrhSR3SKyX0Rur9L/aXHKHFSiqusx8R531IXMluBglYklXPg/VU1U1SRgDLATWCYiV7osl6V2TMAEQS6B02UbngW+hcm6/bCItHL2jQYuAJ7wcp45wJ0iYp9R9QR7oyxhh6ruVNX7gX8ATzhvu4hICxH5h4jkicg+EZkrIs0rjxORHSJyv5jiXsdFZK2IXCgi1zrFf4442XFjPI55wXlrPiYiG0TkOo99GSKiTnp6ROQ3IrJYRH4npvjQASei2isex98qIltEpFBE3hKRlh59ahrTXSKy3ZFvj4j8zmOfisgIZxnuPSDaGfdxEbnRs4/HMdNFZI3zt1jjqaxF5HvO3+kuEckVkcMi8oyjEPxlGrBIz0RDtwBKVXW1qu4EsoGOItIEE3l9q6pWeDnPJ0BroF8trm1xEatMLOHMy5hUF5X5o+YDKZjCZD0x6SLmVTnmRuB/nH5rgIWYmU5fzFvwFOAaj/6fYR5YycBDwItSfcGzkZh0FG0xb9v3icjFNYzjBue4DpgUNf/02OdzTCLSHXgEuEJVE4HeeMkF5yzDTQbKnWWtBFWdW7WfYzuZj8kO3AKTJudfIjLUo1tHoBXQBRgMXI3JV+YvAzD5rirJB8pFZKCIdAIygBzgYeAlVd147inASWeS7ZzPUg+IqbmLxeIauc62hfP2PRHorqqHAUTkp8AmEWnjkZxuduUDSkReAr4LXKSqRUCRiCzFPCRfAlDV5z2u97KI3IPJa+T5QPRki6o+7Xz/SkRWA4OAz6sZx4NOYkVE5OdAtodR3+eYMPmjBOgtIjud4lXLqrlOTdwEvKGq7zm/3xWRhcDNmNxrACeA/6eq5UCOiCx2xjffz2ukAKdLFqiqisi1mKqOscDtQFdgFDBSRB7FKIytwF1OQbZKjmLyZVnqAXZmYgln2jnbAs6kyd7usX+rs/VMoe2Z8bQY87aeX6UtEUBEokTkIRHZ7Cz7FGJmMGnVyFQ1o2pR5fmqYYeX7+2oYUyqug2jDL8P7BWRz0RkQg3Xqo72mISOnmzl7L/fAUeRVOLP+Dw5DDTzbFDVT1V1hKoOAd4FZmPGNBNIUdWRmL9r1XoqzTi3YqAlTLHKxBLOzAD2AJs5Uw0uw2N/Z2frWSmuNlyLyRo7HfNQS8YsjXkrGBQIGV6+5+LHmFR1gaqOxyx/vQq8JSLxXq7hze5Qld2YwlCedOb8/37eyMIs2fnifuBdVc3CKO4vnPZPMFU9ARBT4qGbcz5LPcAqE0vYISLtHcP294AfO/Wp9wIfAn8RkWQxVfP+Arx3vvUXMG++ZZh1/SgRuRnzgAs2D4hIKzF11/8ALFbVvTWNSUR6iMgkR3mUYsrMKt4Vxz6MAb6qsvDkRWC6mDiQaBGZDHwbeCFoIzWpzsd52yEiF2IqHT7kNG0FxjsOEZMxtpRKRmLSp1tlUk+wysQSLjzgeCwdxbyldgWGq+obHn2uB44Bm5xPIca4fb7MxdgKcjAzoF6YWi7B5p/OeXdj7AbXe+yrbkyxmHiNPKf9LmC6qp6segFV3QL8HVjueI3N9NLnC4yDwp8xy1F/xNTuCMQOU5UPgDIxbr+ncTzCngNu96gV8gxmjAUYJf47j0NuBp7w4ellCUNsPROLJUSISAbGHtJeVXOr7x05iMgk4D7HFnI+x/fCeOFdWJ+LVDU0rDKxWEJEQ1UmloaJXeayWCwWS8DYmYnFYrFYAsbOTCwWi8USMFaZWCwWiyVgrDKxWCwWS8BYZWKxWCyWgLHKxGKxWCwB8/8BAjL3BA4hTQEAAAAASUVORK5CYII=\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 first warp:\n", "Q = 101 #domain size\n", "warp0 = mwarp1d.ManualWarp1D(Q) #constrained Gaussian kernel warp object\n", "warp0.set_center(0.10) #relative warp center (0 to 1)\n", "warp0.set_amp(0.3) #relative warp amplitude (-1 to 1)\n", "warp0.set_head(0.0) #relative warp head (0 to 1)\n", "warp0.set_tail(0.0) #relative warp tail (0 to 1)\n", "\n", "#define second warp:\n", "warp1 = mwarp1d.ManualWarp1D(Q)\n", "warp1.set_center(0.90)\n", "warp1.set_amp(-0.3)\n", "warp1.set_head(0.0)\n", "warp1.set_tail(0.0)\n", "\n", "#create and apply sequential warps\n", "seq = mwarp1d.SequentialManualWarp()\n", "seq.append( warp0 )\n", "seq.append( warp1 )\n", "y = np.sin( np.linspace(0, 4*np.pi, Q) ) #an arbitary 1D observation\n", "yw = seq.apply_warp_sequence(y) #sequentially warped 1D observation\n", "\n", "#plot:\n", "plt.figure()\n", "ax = plt.axes()\n", "ax.plot(y, label='Original')\n", "ax.plot(yw, label='Warped')\n", "ax.legend()\n", "ax.set_xlabel('Domain position (%)', size=13)\n", "ax.set_ylabel('Dependent variable value', size=13)\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Individual warps can be accessed and applied individually like this:" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZMAAAEJCAYAAABR4cpEAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOydd3gc5dW377Nqq14sNyxbcpG7JBfJNm4SrvROqAkh4SVfIAmEhN4SAkleUjAkJG8IAUJJIBBagjHYuGODLRv3XrHcbVm9a8/3x+waWd6VVtKWWWnu65prdqceaaU98zy/U0RVsbCwsLCw6Ai2YBtgYWFhYRH6WM7EwsLCwqLDWM7EwsLCwqLDWM7EwsLCwqLDWM7EwsLCwqLDhAfbgGCRmpqqGRkZwTbDwsLCIqRYs2bNCVXt3nx7l3UmGRkZFBYWBtsMCwsLi5BCRPa7225Nc1lYWFhYdBjLmVhYWFhYdBjLmVhYWFhYdJguq5lYWFh4T319PUVFRdTU1ATbFIsAYbfbSUtLIyIiwqvjLWdiYWHRKkVFRcTHx5ORkYGIBNscCz+jqpw8eZKioiL69+/v1TmmmeYSkRdF5JiIbPKwX0TkWRHZJSIbRGRMk303i8hO53Jz4Ky2sOga1NTU0K1bN8uRdBFEhG7durVpJGoaZwK8DJzfwv4LgEznchvwZwARSQEeA8YD44DHRCTZr5ZaWHRBLEfStWjr522aaS5VXSoiGS0cchnwiho18z8XkSQR6Q0UAPNVtRhAROZjOKV/+tdic3D4MLzyClRUKMfq69lTV8WgCyroPqgegKTwcFIjIugREUFOXBy9o6KCbLFFKKGqfFxcTEpDA0dqa3EA8WFhxIWFWc7F4gxM40y8oA9woMn7Iuc2T9vPQkRuwxjV0K9fP/9YGSCOH4ennoLnnlOqqwUEIAI0iQXPJML0o3DzPuhz5jB1UHQ0UxMTmZaczMzkZHpERgbDfAs/U18P27fDyJHtv8aa8nJ+tHMnK8rK+CgxEVtd3el9cWFh9I6MJCGATqWoqIg77riDLVu24HA4uPjii/nNb35DZCt/w64E5dTUVK/u8+yzz/LnP/+ZMWPGcO2117Jlyxbuv/9+r87dt28fK1as4IYbbvDq+M6Emaa5WsPdX6y2sP3sjarPq2ququZ2735WNYCQYf58GDBA+d3vlYb8Y8hrnzPjyw38pegwW4rq+OndQvTyXoTdPJ7fbZ7CjnHjWDpqFL8dOJDhMTG8e+IEN23dSs8VKxhTWMgv9+/ncG1tsH8sCx9QUQFz5sDAgZCVBStXtv0a9Q4H39+xg7w1a9hVXc0LQ4aQFhXF6Lg4RsfF0S8qijqHg53V1eytqSEQDfZUlSuvvJLLL7+cnTt3smPHDioqKnjooYd8fq8//elPzJ07l9dff51LL73UrSNpaGhwe+6+ffv4xz/+4XObmtLY2OjX67cbVTXNAmQAmzzs+wtwfZP324HewPXAXzwd52kZO3ashiKffqoaZXdo1KBK5e9f6Ix163RjeflZxx0+rHrRRaqg+vzzZ+5rcDh0dWmpPrlvn05eu1ZZtEjDFi3SyzZs0P+eOKENDkeAfhoLX7Jtm2pKivGZjx9vrJ95pm3XcDgc+u2tW5VFi/SunTu1pL5eVVW3bNlyxnGNDocerKnR1WVluquqSh1+/ptZsGCBTpky5YxtpaWlmpKSopWVlfrSSy/pFVdcobNnz9ZBgwbpPffcc/q49PR0PX78uD788MM6Z86c09sffPBBfabZL+h73/ueRkRE6MiRI/X3v/+9vvTSS3rHHXeoqurNN9+sP/7xj7WgoEDvvvtuXbx4sebk5GhOTo6OGjVKy8rKdPz48ZqQkKA5OTn6+9///oxrf//739f3339fVVUvv/xyveWWW1RV9YUXXtCHHnpIVVUvu+wyHTNmjA4fPlz/8pe/nD43NjZWH3nkER03bpwuW7ZM09PT9d5779W8vDzNy8vTnTt3dvRX7Jbmn7uqKlCobr5TQ2ma6wPgByLyBobYXqqqh0XkY+CXTUT3WcADwTLSnyxZAhderDScU0Xy05v414SBXOIhwqZXL/j3v+HKK+G22yAsDL7zHWNfmAi5CQnkJiTwYHo6O6uqeOHwYV4+coT3T56kb1QU3+nVi/855xz6WBpLyPCPf0BJCXz2GZx7LnTvDhs3tu0aj+zdy8tHjvBYejo/8xASetfOnayrqACgTpVah4MIEey29k90jIqLY05mpsf9mzdvZuzYsWdsS0hIoF+/fuzatQuAdevW8eWXXxIVFcWQIUP44Q9/SN++fU8f/93vfpcrr7ySO++8E4fDwRtvvMGqVavOuOb//d//MW/ePBYtWkRqaiovv/zyGft37NjBggULCAsL45JLLuG5555j0qRJVFRUYLfb+fWvf81vf/tb/vvf/571M0ydOpVly5Zx6aWXcvDgQQ4fPgzA8uXLue666wB48cUXSUlJobq6mry8PK666iq6detGZWUlI0eO5PHHHz/j51+1ahWvvPIKd911l9t7BhLTTHOJyD+BlcAQESkSke+KyP8Tkf/nPGQusAfYBfwVuB1ADeH9F8Bq5/K4c1unYsMGuOAipb5nFb2e3cJn07K4NDW1xfnqqCjDocyeDbfeCu++6/64zJgY/nfgQA6cey5vDR/OsJgYHt+/n/6ff853tm1ja2Wln34qC18yfz7k5cHEiSBiTHO1xZn8+eBBnvzqK/6nd28e87KidqQIUTYb9arUOBztM9wLVNXt33rT7dOnTycxMRG73c7w4cPZv//MeoQZGRl069aNL7/8kk8++YTRo0fTrVu3NtlxzTXXEBYWBsCkSZO4++67efbZZykpKSE8vOVn8ylTprBs2TK2bNnC8OHD6dmzJ4cPH2blypVMnDgRMPSanJwcJkyYwIEDB9i5cycAYWFhXHXVVWdc7/rrrz+9Xtme+UwfY5qRiape38p+Be7wsO9F4EV/2GUG6uvhmpsaqY5qZOBz21k0LYu+drtX59rthhOZOtUYoUyaBD16uD820mbj6h49uLpHD/ZUV/N0URF/O3yYl44c4YrUVB7LyCAnLs6HP5mFrygthVWr4IEmY/KsLHjxRXA4oLVBw5rycn6wcycXd+vGnzIzW3xIcTeCOFBTw9H6egZHR5PQypdqexgxYgT//ve/z9hWVlbGgQMHGDhwIGvWrCGqySg6LCzMra5x66238vLLL3PkyBG+4xqqt4HY2NjTr++//34uuugi5s6dy4QJE1iwYEGL5/bp04dTp04xb948pk6dSnFxMf/617+Ii4sjPj6exYsXs2DBAlauXElMTAwFBQWn8zzsdvtpJ+ai6Wdkhsg604xMLDzz8182smNjGKn37GHFjJFeOxIX0dHw8stQVgY/+IF35wyIjuYPmZnsnzCBR9LTWXjqFKMKC7lq0ybWO6c4LMzDokXQ2AgzZ369LSsLKith376Wz3WocvuOHXSPiOC1YcMIb8d0VZ+oKOwi7KupodEPgvz06dOpqqrilVdeAQwR+ic/+Qnf/va3iYmJ8fo6V1xxBfPmzWP16tXMnj27Qzbt3r2brKws7rvvPnJzc9m2bRvx8fGUl5d7POfcc89lzpw5TJ06lSlTpvDb3/6WKVOmAFBaWkpycjIxMTFs27aNzz//vMX7v/nmm6fX5557bod+Fl9gOROTs2ED/OoJG0w/ynvf793uUN4RI+BnP4O33jIWb+keGcnj/fuzb8IEHktP51PLqZiS+fMhNhYmTPh6W3a2sW5tquuFw4dZVV7ObwcOJLGdowqbCBl2O3WqFPkhMlBEePfdd3nrrbfIzMxk8ODB2O12fvnLX7bpOpGRkZx33nl84xvfOOtJv63MmTOHkSNHkpOTQ3R0NBdccAHZ2dmEh4eTk5PD008/fdY5U6ZMoaGhgUGDBjFmzBiKi4tPO5Pzzz+fhoYGsrOzeeSRR5jQ9MN0Q21tLePHj+eZZ55xe6+A406V7wpLKERz1dWppmfVKcm1eu/a/R2+Xn296tixqt27qx471r5rnKqr08f27NGEpUuVRYv0mk2bdHdVVYdts+gYmZlG9F5TysuNiK5f/MLzecdrazVl2TLNX7u2xYgsd1E97viqulpXl5VpqTMKzGw0NjZqTk6O7tixI9imdAhXhJq/aUs0lzUyMTG/+mM9+zdGMPTBIp7MSevw9cLDjemukhJ4+OH2XSMpIoKfOUcqD6en8+HJkwxbtYp7d++mpL6+wzZatJ39+2HnzjOnuADi4mDAAGN064kH9u6lrLGR5wYP9sm8e5+oKKJE+Kq21hWqbxq2bNnCoEGDmD59OpktRI5ZtA/LmZiU6mr49a/AllPC3Nt7t2se2x0jRxpC/Isvtj6X3hLJERH8on9/dowfzw09e/LbAwcYvGoVLx0+jMNkXyKdnfnzjfWMGWfvaymia2NFBS8cPsxdaWmMaCIsdwSbCH2ioqhxOCj2kNgXLIYPH86ePXv43e9+F2xTOsy+ffu8zugPFJYzMSmPPltL9fEIrruviv4x0T699gMPGHknTzzR8Wv1iYripaFDKRw7lszoaL6zfTuTvvySNS2IkBa+Zf58OOccGD787H1ZWcaoxV3x198cOECszcYDPi4tlBweTrTNxmETjk4s/IflTExIZSU8+5SN8DEl/PFq35d96dMHvvc9Y8pr927fXHNMfDzLR4/m70OHsre6mrw1a/h/27dz0pr68isOB3z6qTEqcTdLlZVlRHlt3Xrm9gM1Nfzz2DH+55xzSPGy+ZG3iAjnREZSo8pJk41OLPyH5UxMyD2/r6auOIJbH6wl2cf/6C7uvx8iIuAXv/DdNUWEb/Xqxfbx47kzLY0XDh9myBdf8NdDh6ypLz+xbh2cPHm2XuIiK8tYN5/qmlNUhKry47SOa3HuSAoPJ8Y5OrE++66B5UxMRnm58tffhxM17hS/u9x/c6K9e8Ptt8Orr8KOHb69dmJ4OE8PGsS63FxGxsZy244dFKxbxzYrk97nLFtmrM87z/3+zEyjEkJTZ3Kqvp7nDx/muh496NfGnCVvcY1OalUptkanXQLLmZiMnzxbRUNJBD9+uIGYDsbBt8a99xpfNE895Z/rj4yLY9GoUbw4ZAibKivJKSzkiX37aPBj2Y2uxpo1xoNBH7dNF4wIvmHDznQm/3foEBWNjdzj5zYMic7RyZG6Op9oJ0VFRVx22WVkZmYycOBA7rzzTuqalMX3REZGBidOnOjw/QPF4sWLufjii4NtRpuxnImJcDjg9efDiRxZxuMXta1mUHvo2RNuuskoEHjqlH/uISLc0rs3W8eN44rUVB7Zt4/8devYV13tnxt2MdasgWb1D8+iaURXrcPBswcPMjs52e+lcUSEHk7tpLyDZdNVA1eCvjOjqjj89DDntTMRkQEi8qCIPOd8P0RERvjFqi7KSx9WU/VVFJfeWkeEj0KBW+P2240w5GbFUX1Oz8hI3hgxgn8MG3Z6lPKPo0f9e9NOTmUlbNvWujPJzoZDh6C4GN47cYIjdXXc3aSarj9JCQ8nHDjWwamuhQsXYrfbueWWWwCj9tbTTz/Niy++SFVVFS+//DJXXnkl559/PpmZmdx7771nXeORRx7hmWeeOf3+oYce4tlnnz3jmMrKSi666CJycnIYOXLk6ZIla9asIT8/n7FjxzJ79uzTFX937drFjBkzyMnJYcyYMezevRtV5Z577mHkyJFkZWWdvsbixYspKCjg6quvZujQodx4442nR2zz5s1j6NChTJ48mXfeecft7+DCCy9kgzNpaPTo0acrCD/yyCO88MILVFRUMH36dMaMGUNWVhbvv/8+YIQRDxs2jNtvv50xY8Zw4MAB4uLi+MlPfsKYMWOYPn06x48fb98H0xR3mYzNF2AmUI5RBr7MuW0S8JE355txMWMGfP/p5UpSrR4oqw3ofSdOVB00SLWxMTD321tVpRPXrFEWLdJvbtli2mxps7N8uZHh/sEHLR83b55x3OLFquevX699V6zQxjb2H2maCX3nnar5+d4vE6Y06pjJ9To13+HxmDvvbPn+zzzzjN51111nbR81apSuX79eX3rpJe3fv7+WlJRodXW19uvXT7/66itV/TpbfO/evTp69GhVNTLhBwwYoCdOnDjjem+//bbeeuutp9+XlJRoXV2dnnvuuXrMWTbijTfeON2LZNy4cfrOO++oqmp1dbVWVlbq22+/rTNmzNCGhgY9cuSI9u3bVw8dOqSLFi3ShIQEPXDggDY2NuqECRN02bJlWl1drWlpabpjxw51OBx6zTXX6EXNyxmo6q9+9Sv94x//qKWlpZqbm6uzZs1SVdWCggLdtm2b1tfXa2lpqaqqHj9+XAcOHKgOh0P37t2rIqIrV648fS1AX3vtNVVV/fnPf366Z0tz/JEB/2vgGlW9FHCNV9cCYzruziwANu1pYO+iWLKuKSMtPrCtdO+4A3btglaKnvqMjOholowaxc8yMnj96FFGFxbyeWlpYG7eiSgsNNatjUwGDDDW63bV80lxMTf36oUtgFVmI2zGveoc7ddNVANTgj4rK4sFCxZw3333sWzZMhITE9m+fTubNm1i5syZjBo1iieeeIKioiLKy8s5ePAgV1xxBWBU9o2JiWH58uVcf/31hIWF0bNnT/Lz81m9ejUA48aNIy0tDZvNxqhRo9i3bx/btm2jf//+ZDqrNd90001ufwdTpkxh6dKlLF++nIsuuoiKigqqqqrYt28fQ4YMQVV58MEHyc7OZsaMGRw8eJCjztF/enr6GbW+bDYb1157LQA33XQTy5cvb8/HcgbeVnUbqKrznK8VQFWrRcQ/catdkLufrgLi+eWd/omuaYmrroIf/xieew5mzQrMPcNtNh7LyGBGcjI3btnC5C+/5Of9+3N/v36EmaCcdijgEt/POafl41zRvx9tq8QxEL7Vs2eH7jtnTlvPEHZW1VHlcJAVG9suRxaoEvSDBw9mzZo1zJ07lwceeIBZs2ZxxRVXMGLEiLN6hpSVlbm1VVsINvBkozelbPLy8igsLGTAgAHMnDmTEydO8Ne//vV007DXX3+d48ePs2bNGiIiIsjIyDhdwj62lQoHviil4+3I5ICIjGx28xxgX4ctsKC6Rln4ejTJU0q5eETg+4VERRnNs/77X6POUyCZlJjIutxcvtGjBw/v3cu0dev4yl26tsVZeCO+g9GCIDVV+WJ3HZMSEshsQ8l2X9EjMpJ6VUramcQYqBL0hw4dIiYmhptuuomf/vSnrF27liFDhnD8+PHTzqS+vp7NmzeTkJBAWloa7733HmBU8a2qqmLq1Km8+eabNDY2cvz4cZYuXcq4ceM82jR06FD27t3LbmcG8T//+U+3x0VGRtK3b1/+9a9/MWHCBLcl7Hv06EFERASLFi06a2TWFIfDwdtvvw3AP/7xDyZPnuzFb69lvHUmzwLviMhNQJiIXAW8Bvis7rGInC8i20Vkl4jc72b/0yKyzrnsEJGSJvsam+z7wFc2BYqfv1JB46kIvvf94CV3fe97xvr55wN/76SICF4fNoy/Dx3K2ooKRhUW8pk17dUi3orvLrqd46DkUBg39+rlX8M8kBAWRpQIx9spxAeqBP3GjRsZN24co0aN4sknn+Thhx8mMjKSt99+m/vuu4+cnBxGjRrFihUrAHj11Vd59tlnyc7OZuLEiRw5coQrrriC7OxscnJymDZtGk899RS9Wvi92+12nn/+eS666CImT55Menq6x2OnTJlCz549iYmJYcqUKRQVFZ12JjfeeCOFhYXk5uby+uuvM3ToUI/XiY2NPd0KeeHChTz66KPe/go9405IcbcAtwEbgQpgE3Crt+d6ce0wYDcwAIgE1gPDWzj+h8CLTd5XtPWeZhLge08pU1vPaq2uD5AC7oELLlDt1y9wQrw7dlZW6uDPP1f7kiX6nwCU2A5VvBXfXWScV64ysFxL2hns4G0J+pY4WFOjq8vKtCaIf2CdpQR9R4mNjfXqOL+UoFfV51U1S1XjVHWkqr7QcVd2mnHALlXdo6p1wBvAZS0cfz3gfiwYYuw61MDhFXGMuqwKe3hw035uvBG++go++yx4NgyKiWH56NGMjI3l8k2beNkZgmlxJmvWGGtvRia1DgeHk8qIOB7d7uZXvsBVAyxYGfFWCXr/4tVfloh4lPhU9ZAP7OgDHGjyvggY78GWdKA/sLDJZruIFAINwK9V9T0P596GMcKin5+zf73lsRcroTGRn34nsBFc7rj8cqNb32uvgXPkHBS6R0ayMCeHKzdv5pbt2zlWX8+9Jvm8zEJhIfTq1br4DvBxcTG1qTVQFkZlpfEZBwO7zUaszUZxQwO9mwjRgcJVgt4CKvzQJdXbR+EijC97d4svcBdK4ElAuA54W1WbptT2U9Vc4AZgjogMdHeic3SVq6q53bv7vhpve/jPm+FEZlZy3bgg/Yc3ITbWcChvvQV+6LzaJuLDw/kwK4tru3fnvj17uGf3bqtgYBPWrIHcXO+Ofef4cWJ6GqOBAx34j1Uf/P67RURQ7XBQ1cGMeAv/09bP21tn0h9Dz3AtU4CPgW+36W6eKQKapuSmAZ5GPNfRbIrLNTpS1T3AYmC0j+zyK0s31VK+KZaCq2t9EprnC266ySit8tFHwbYEIm02Xh8+nDvOOYffHjjALdu2UW/V9WqT+N7gcPCfkyeZlGmEnBcVte+edrudkydPdtihJDun2azij+ZGVTl58iT2NhQC9WqaS1Wbx5jtF5GbMaaaXvXeRI+sBjJFpD9wEMNh3ND8IBEZAiQDK5tsSwaqVLVWRFIxMvP9VLrQtzz+txqQSB79jm+bX3WEGTOgRw94/XVjlBJswkT4Q2YmPSIjeWzfPoobGnhz+HC/F8E0M+vWGXXcvHEmS0tLKW5o4PJhccyn/SOTtLQ0ioqKfFJ2o7Kujq2qlEVFuZ2SsDAHdrudtDa0KOiIGlcFeI5hawOq2iAiP8AY7YRhRGptFpHHMSIHXOG+1wNv6JmPR8OAv4iIA2Ok9WtV3eILu/yJw6EsezuKxNxyJg1KCLY5pwkPh2uvNUKES0shMTHYFhlhoY9mZNA9IoI7du5k1vr1/Ccry2+9XszO2rXGeowX9SfePXECu83GdSOSuIP2O5OIiAj69+/fvpObse7oUW7YupXFo0aRn5Tkk2taBB9vBfjmo4RYjC/2L3xliKrOBeY22/Zos/c/c3PeCiDLV3YEitcXV1FXFMs37jJf9dybboI//AH+/W9wkyQcNL7fpw+pERHcuHUrU9etY/GoUXTrgg5l/XpITW1dfFdV3jtxgtnJyaTEhNGjR8c0E19xaWoqsTYbrx09ajmTToS3msmTzZYfAfvxnWbS5Xj2tXqIcPD4zYHPeG+NvDwYNAjeeCPYlpzNNT16MDcrix1VVdy4ZQuNXVCUX78ecnLct+ltSmF5OUW1tVyeajRZ69u3/ZqJL4kNC+Py1FTeOX7c6m3TifDKmahq/2ZLlqreoqoH/W1gZ8ThUNbNs5M6voL+qeZ7shaBK6+ERYv81+ekI8xISeEPmZl8fOoUP9u3L9jmBJSGBti0yXAmrfHuiROEAZc0cSZmGJkAXNm9O8UNDSy1Kh10GqzmWEHgnc+raThsZ/al5g2PvOIK44vrww+DbYl7/qd3b77TqxdP7N/PByHURa+j7NwJNTVGj5LWeO/ECaYmJZ2eCkxLM48zmZ2Sgt1m490u9Nl1djxqJiIyH8+5HqdR1QDVme08/PGNOrBFc991gS+45y3jxhlz8u++a2goZkNEeC4zkw2VlXxz61YKx44NSgHDQOPsjdTqyGR7VRVbq6r4fhNhpW9fKCszloQgx3zEhoUxOzmZ906c4NlBg0wTGm/RfloamSwHPvNisWgjX8yNJGFUBVl9A58F7C02mxEa/NFHUFUVbGvcYw8L498jRhAuwtWbN3eJRLj167/u694Sc0+eBOCSJv06XM0VzaCbAFzRvTtFtbUUlpcH2xQLH+BxZKKqPw+kIV2FhRtrqNkdw6xHTShGNOPKK+FPf4JPPjFHzok7+tntvD5sGBdu3MgdO3fy4pAhnfopd/16w5G0Vo3kk1OnGBIdTUb01zlMLmdy4AAMH+5HI73kkm7dCAPeOXGCvGAPlSw6TJs0ExGJFpG+ItLPtfjLsM7K7/9p9Or4yfXmSVT0xNSpkJxsTHWZmfO7deOR9HRePnKEF48cCbY5fsUVydUSNY2NLCkpYXZKyhnbXflnZtFNUiIiyE9K4l1f9B+3CDpeORMRGSAiKzDKz+8D9jZZLNrA0v+GEz20kqlDA99Rsa1ERMAll8AHH4DZq188mpHBzORk7tixgy2VlcE2xy+cPAkHD7Yuvi8rLaXa4WBWM2fSp48RqWeWaS4worq2V1eztZN+Zl0Jb0cmf8Qo6pgDlAPZwHvAd/1kV6dk3d5ayjfGce6FdcE2xWuuvBJKSmDJkmBb0jJhIrw2bBjRYWH81NmxrrPhrfj+yalTRIpQ0CwhMCLCqDRslpEJcDoHxorqCn28dSbjMZphbQJQ1c3A94B7/GVYZ+SpN4xs9x9cZ17hvTmzZkFMjPmnusBoDftwejofFRczv7g42Ob4nPXrjXVrzuTj4mImJyYS66Z+mZnCgwH6REUxLj7eciadAG+diQNw1f2oEJEkoBiwNJM28Ok8G+F9arhsrPn1EhfR0YZD+fBDCIVk8x/06UOG3c5Pd+/udNnx69dDz57G4olDtbVsrKw8Sy9xYabERReXpqZSWF7O0brQGbFbnI23zmQzRjVeMOpxPY3RF97STLzkVEUjx76IZdi0amy20Io2uvBC2L8ftpi+fCZE2Wz8esAANlRW8monE+O9Ed/nO0sWzEpOdrvf5UzM5GcvdDq+eZ1wNNmV8NaZ/AhjJALG1FYfIBdjqsvCC56bWwW1YVxzceiVTr/gAmM9d27Lx5mFb3Tvzvj4eB7au7fT5J7U18Pmza2L7x8XF9MzIoLsOPc13/r2NfqhmKmKyai4OHpHRp7OjbEITbytzbVBVTc6X+9R1VmqOkFVV7Z2roXBW/9phKhG7rgw+B0V20pamvElFirORET47cCBHKqr42VhoXcAACAASURBVK+dpIf89u1QV9fyyMShyvxTp5iVkoLNQ66N2cKDwfi8LkhJ4ePiYqvwYwjjbWjwAhG5QUTMH89qQhwOZcsiOz3GV5ISF3ojE4CLLoLly831RNsSk5OSmJiQwB+KijqFduJNJNeXFRWcqK/3qJfAmYmLZuLCbt0obWxkRVlZsE2xaCfeTnN9BjwBHBaRP4tInq8NEZHzRWS7iOwSkfvd7P+2iBwXkXXO5dYm+24WkZ3O5WZf29ZRPtlQQ8NBO+edH7pPXRdeaBR+nD8/2JZ4z11paeyuqeHDTjB9smGDEdo7dKjnYxY69ZLpLfQI6ecMmfnqK19a13FmJCcTLmJNdYUw3k5zPaaqA4CrMBpjLRKRTSLyY18YISJhwHPABcBw4HoRcVfw4U1VHeVcXnCemwI8hhG+PA54zNnK1zT8+R0j6/2OK0J3YDdhAiQlhc5UF8AVqan0jYpijpmy9NrJhg1GGZWWeoEtKSlhSHQ0vVqotdK7t3ENs1XuTwwPZ3JiInMtET5kaVM5FVVdqKrfwmjXux/4rY/sGAfscuoxdcAbwGVenjsbmK+qxap6CpgPnO8ju3zC8k/CiexfxZQQyHr3RHg4zJ5tFH4MlWntcJuNH/bpw6KSEjZUVATbnA6xYUPL4nujKstKS89KVGxOWJgxOtm/38cG+oCLUlLYWFnJgZqaYJti0Q7aWpsrR0TmAFsxWuX+ykd29MHIsHdR5NzWnKtEZIOIvC0ifdt4blA4fKqB4jWxZE8P/Rj6Cy+EI0dg3bpgW+I9t/buTYzNxjMhPDopLm69jMq6igrKGhu9aoObnm6+kQkYugnAR9boJCTxVoD/oYisBVYCPYGbgHRVfdhHdrgLPWmumv4HyFDVbGAB8Pc2nGscKHKbiBSKSOHxABWXe/Y/VdBg44ZLQ1N4b8r55xu1nczaMMsdyRER3NyrF68fPcqxEE2K27jRWLfkTJaUlAB45UwyMsw5MhkWE0N6VJSlm4Qo3o5MvgX8FThHVa9X1U9UfRoiUwT0bfI+DTjU9ABVPamqtc63fwXGentuk2s8r6q5qprbvXt3nxjeGu9/6ICYBm6bFXohwc3p0QNyc42prlDih336UKvKm8eOBduUduGK5MrK8nzM4pISMqOjOae12vQYI5PDh42OjWZCRLigWzc+LSmhPlTmUi1O460An6eqf1bVEj/ZsRrIFJH+IhIJXAd80PQAEend5O2lGFNtAB8Ds0Qk2Sm8z3JuMwU7l9vpNa6K2KjO0SF51ixYtSp0QoQBhsXGMtBu5xMzNrT3gg0boFs3Qzx3h0sv8WZUAsbIBMwXHgxG5n5FYyOfWyHCIYcpvuFUtQH4AYYT2Ar8S1U3i8jjInKp87AfichmEVmPkZH/bee5xcAvMBzSauBx57ags3BzDQ2H7EyZ3jmysAFmzoTGRli0KNiWtI2ZKSksDtEnXpf47qnn18aKCkoaGshPTPTqeunpxtqMusl5SUnYIGQdf1fGFM4EQFXnqupgVR2oqk86tz2qqh84Xz+gqiNUNUdVz1PVbU3OfVFVBzmXl4L1MzTnrx8Y8wjfvSR0o7iac+65EBsbWvkmADND9InX4YBNm1rWSxa3QS+Br0cmZtRNkiIiGJeQ0CmrPnd2TONMOiNLPxXCetYyM6vzOJPISCgoCD1nMs35xDs/xJ549+yBqqpWxPfSUgbY7fS1e/d31qePESJsxpEJGI5/dXk5p8zekc3iDCxn4idq6hwc+iKWgZNDr0pwa8ycCTt3mvfLyB2nn3hDzJm0Jr47VFlaUuL1qASMnKG0NHOOTMDQTRzAohJ/SbQW/sBrZyIik0XkeRH5j/P9WBGZ6j/TQpvXllZBRTgXzu5cjgQMZwKhNzqZmZzMqrIySkLoiXfDBkMrGTHC/f5NlZUUNzS0yZmAMdVl1oeB8QkJxIeFhZzj7+p4m2dyA/A+UAO4HIgCj/vJrpDnH3PrQZTbL4kJtik+Z9gwOOec0HQmDmBhCD3xbtwImZlGt0t3LHOG1U31Unx3kZ5u3pFJhM1GQVKSpZuEGN6OTB4CZqvqjzC6LgJsAjw8L1msWRRBzNAqMnu1UEwpRBExRieffmpEdoUKExISiAuxJ97WyqisKC2ld2QkGV7qJS4yMoyserPmcc5MTmZ3TQ17qqtbP9jCFHjrTM5R1ULna1eyYgMQ+mndfuDQqQbKNsaQU2DS/1QfMGuWUebjyy+DbYn3hNoTb0UF7N7dcrLiyrIyJiYkIJ7ihj2Qnm5Eipm1ysxMZ6fIUHL8XR1vncluEZnYbNtEYLuP7ekU/HluFTTauPr88GCb4jdmzDDWoTbVNcv5xLs3BJ54N2822ut6Gpkcrq1lb00NE9s4xQXmDg8GGBITQ9+oqJBx/BbeO5MngPdF5GEgQkR+AvwTSzNxywcfOyCqkVtnhn4JFU/06GE0ago1ZxJKT7yt1eRa6cyZmZiQ0OZrmzlxEYzSKjOTk1lYUtIpmpt1Bbwtp/IecANGz5D9wDTgO6oaYlWaAsP2zyLpNqaKhOjOHXk9fTqsWGG+Gk8tMSQmhp4RESwNARF+wwaIi/t6FNGcFaWlRIowOj6+zdfu29fQvsw6MgGYnpzMqYYG1od4+4Cugtffdqo6X1UvUdWRqnqRqobYM2lg2HSgjto9MeTlNwTbFL8zbRrU1sLKlcG2xHtEhPykJJaUluLbWqW+Z/16Qy+xefgvXVlWRm58PFGeDmiByEgjIs+sIxMwSqvA1x0kLcyNx79CETnHmyWQxoYCz8+tAuC6WZ0viqs5U6YYmdQLFwbbkrYxNSmJotpa9pt4SKXaciRXrcNBYXl5u/QSF2YtRe+id1QUQ2NiQiqUuyvT0iNNEUbTKU+La79FE+Z/CsQ1cO2kzpdf0pyEBMjLC0Fn4vwCXmLi0sdFRVBSYuhS7lhbXk6dKue2Qy9xYebERRfTkpJYGqIFOrsaLTmT/sCAFhbXfosm7F5pp1duFfbIzq2XuJg2zShJX14ebEu8Z0RsLCnh4abWTdavN9aeRiYrnOJ7R5xJerrhtBpMPCM7LTmZSucozMLcePzGU9X93iyBNNbsrNxRQ32RnQlTQyiTr4NMm2Z8GS1fHmxLvMcmwpTERJaaeGTSWk2ulaWl9Lfb6e1FMyxPZGQYn90ht63kzIGrrL411WV+2lKb61oRmScim5zr6/xpWCjy4jxjDv7GC9r/Dx5qTJxoiLkhN9WVlMSu6moO1da2fnAQ2LAB+vc3phKbo6p85kxW7Aiu8GAz6yapkZHkxMZaInwI4G1trvuAZzGaTz3tXM9xbrdwsnAhSFI9l+dFB9uUgBEdbfQ4CTln4nziNetU1/r1nqe49tfUcKSujnM7IL7D1yHHe/d26DJ+Z1pyMp+VllITSrV7uiDejkzuAC5Q1UdU9W+q+ghwAUZ3RJ8gIueLyHYR2SUi97vZf7eIbBGRDSLyqYikN9nXKCLrnMsHzc8NBA6Hsn9lNH3GVxEe1vkqBbfEtGlGWZVQSlYeFRdHXFiYKae6qqthxw7P4rsv9BIwnElYmNFOwMxMS0qiVvV0kqaFOfHWmcQAG5pt2wj45BFcRMKA5zAc1HDgehEZ3uywL4FcVc0G3gaearKvWlVHOZdLCQKfbqqh8VgUUwq6XtTJtGlGKOvSpcG2xHvCbTYmJSSYcmSyZYtRN8vTyGRVWRnRNhvZsR2rsBAZaUyl7djRocv4nSnOxmaWbmJuvHUmrwA/brbtLuDvPrJjHLBLVfeoah3wBnBZ0wNUdZGqVjnffg6k+ejePuHlj4y595sv6DxdFb1l3DijRHqoTXXlJyWxuaqKEyYrneuK5PI0MlldXs6YuDjC25Gs2JzBg83vTBLDw8mNj2eRpZuYmpaSFueLyCci8gmQDTwpIvtEZImI7AOeBDz8ubeZPpyZs1Lk3OaJ7wJNS7nYRaRQRD4Xkcs9nSQitzmPKzx+/HjHLG7GssWCLbVztej1lshII4Ex1JzJVGeG9TKTTXVt2GA45wFuAu/rHQ7WVlSQ18EpLheDBxvTXCYvBsC05GS+KC+nwsxxzF2clsraNg/2XOZHO9yJDG7/vEXkJiAXyG+yuZ+qHhKRAcBCEdmoqrvPuqDq88DzALm5uT7793E4lIOroul7bjU2W9eJ5GpKQQE88AAcPw7duwfbGu/IjY8nUoQVZWVcYSKjN2zwXEZlc2UlNQ4H49pRj8sdgwdDZSUcPmyUVzErBUlJ/Pqrr1hRVsaslJRgm2PhBo/ORFV/HkA7ioC+Td6nAWdFv4vIDIxGXfmqejqmU1UPOdd7RGQxMBo4y5n4i3nrqnEUxzA1v6r1gzspBQXGeulSuOqqoJriNVE2Gzlxcaw2kbCrakxzXX21+/2rnMl7eT50JmBMdZnZmUxKSCAMWFJSYjkTk9KWPJMwERkqIlNEZKpr8ZEdq4FMEekvIpHAdcAZUVkiMhr4C3Cpqh5rsj1ZRKKcr1OBScAWH9nlFa99Yvi1b87umqMSgLFjITYWFi8OtiVtIy8+njUVFaYpc37okBEV50l8X11eTnJ4OAOjfRN+nplprM2um8Q5dZPFlghvWrzNMxmD8aS/BVjsXBYBC3xhhKo2YIQZfwxsBf6lqptF5HERcUVn/QaIA95qFgI8DCgUkfVOm36tqgF1Jp8ttWFLrWX6yK6nl7iIiIBJk0LTmVQ0NrK9yhyjSlfme0uRXOPi49vcWdETaWlgt5vfmYAx1bWqvJxKK9/ElHg7MpkDvAskAmVAAsYo4du+MkRV56rqYFUdqKpPOrc9qqofOF/PUNWezUOAVXWFqmapao5z/Tdf2eQNDodStCqavuNrsNm6Vn5JcwoKYNMmQzcJFVxC9mqT1H5qqSZXZWMjmysrfSa+g6HLZGaGjjNpUGWlyQImLAy8dSZZwP2qWg6IqlYA92J1WuTjDTU4TkYyeao5pkmCSVPdJFQYGhNDXFiYaXST9euNMifuktu/LC+nEXwmvrsIhfBggEmJiYSBNdVlUrx1JvVNXpeKSA/ntl6+Nym0eO1jSy9xkZtrhLQuWRJsS7wnTISxcXGmGZmsWwejR7vft9rH4ruLwYNh925zVw8GiA8PZ6ylm5gWb53JGmCm8/Vi4FWMxMLmWfFdjs+WCrZudV0yv6Q5IaubJCSwrqKCuiD3zKiqarmMyqrycvpGRdGrA5WC3TF4sOFIzN7bBL7WTaos3cR0eOtMbgWcs7ncDXwF1AK3+MOoUMHhUA58EU3f8dVdXi9xUVAAGzfCiRPBtsR78uLjqVNlY2VlUO3YtMkoozJqlPv9q8vKfD4qgTPDg81OQVIS9VadLlPilTNR1YOqesD5+qSq/o+qXqeq2/1rnrmZv9Gpl+RbeomLUNRNXF/QwdZN1q0z1u6cSXF9PbtrahjnQ/HdRSg5k0mJidiwdBMz0lI5lfFNXk/0tATGTHPyqlMvuXGWpZe4yM01ytKH0lRXht1Ot/DwoOsm69YZwnt6+tn7/KWXAHTrBsnJoeFMEpy6yRLLmZiOlsqpLABcf7me+ugpEOZTi0KI5U69ZHa2pZe4iIwMPd1ERMhLSDCFM8nJAXcpJK62tWP94ExEQieiC4yprjlFRVQ3NhId1mW/fkxHS21745u8tnlYuuwn6XAoB1bZ6ZNn6SXNyc83dJNQ6m+SFx/P5srKoCXEORxGwqInvWRteTmZ0dEkhrf0/Nd+QsmZ5Dt1k88t3cRUtKqZiEi4iHwoItbjdxMWb63FcTyKSVZ+yVnkO0twLvNnaVAfkxcfjwMjlyMY7N5tFFz06EwqKhgTF+e3+w8eDAcOGBFlZmdyYiIC1lSXyWjVmThLnYwFTB6FHlhe/djo9379zMggW2I+xo2DqKjQyjc5LcIHyZm4xHd3YcHF9fXsq6lhjB+muFy4RPhdu/x2C5+RGB7OqLg4U3bJ7Mp4Gxr8Kj5s0dsZWLrU6Pd+8Ziu0+/dW6KiYMKE0HImvaKi6BUZyfqKiqDcf/16CA+H4c37iwJfOm3y98gEYHuIxGfmJyWxsqyM2iDnBll8jbfOZAzwlLM/+wJX0yxn46wuyf4v7PTOrbL0Eg/k5xtP26H08JgdGxu0XJN162DYMKPoYnPWOkdLo/04MsnMNIT4rVv9dgufkp+YSI3DEfRwbouv8daZLMXorPgqRpOsz5osXY7PttfQeMTOhCmWXuKJ/HxDVP4shP5CsmJj2VJVRUMQnnbXrWtZL0mPiqJbRITf7h8ba3R23LjRb7fwKZOdxcuWhNLTSifHq9CQADfKMj2vflID2LnW0ks8MmGCUV5lyRK48MJgW+MdWXFx1Dgc7KquZmhsbMDue+IEHDzouYzK2vJyv+olLrKzQ8eZpEZGMjI2lqUlJTzkLjHHIuCYpTlWSLFoCUhcA1eOs/QST8TEQF5eaGXCZzkdSKCnulxl592NTMoaGthRXe1XvcRFVpbRD7662u+38gn5iYl8VlpKvaWbmAJTNMcKNfZ+YafH2CrCwyy9pCXy86Gw0Ah5DQWGx8RgI/DOpKVILldAQCBGJllZxtTkloC2lms/+UlJVDocrA1S0ITFmZimOZaInC8i250i//1u9keJyJvO/V+ISEaTfQ84t28Xkdm+sskda/fWUl9kJ2+K9TTUGvn5RjXaFSuCbYl32MPCGBwTE3Bn8uWX0KcPpKaevW9tACK5XGRlGetQmeqa4tRNllr5JqbAFM2xRCQMeA64ABgOXC8izYMkvwucUtVBwNPA/zrPHY7RM34EcD7wJ+f1/MLLnxhzANdM908mcmdi4kQICwutEOGs2Fg2BPhJd+1aGDvWw77ycnpHRvq87Lw7Bg0yoslCxZn0iopiSHS0lbxoEszSHGscsEtV96hqHUavlMuaHXMZ8Hfn67eB6WI0wr4MeENVa1V1L7DLeT2/8OliIKaBb0yM8dctOg3x8TBmTOg5kz01NVQEqFNUZSVs22b8ntzh78z3poSFwYgRoeNMwJjqWlZaSqNakZXe8MriCtJmlbJ4S43Pr22W5lh9gANN3hc5t7k9xpmVXwp08/JcAETkNhEpFJHC4x1oVJ42pRJ7pNexC12a/HxYtSp0RF2XCL85QHVF1q8HVffOpKqxkS2VlQHRS1xkZRk1wkKF/KQkyhobg5ZsGmq8Preeg/MT6Z3k+++v9jbH2o9vm2O5U7KbP2p4Osabc42Nqs+raq6q5nbv3r2NJhps/mcSB+a5adBt4Zb8fKirgy++CLYl3pHtHAVsDNCX09q1xtqdM9lYWYmDwOglLrKy4OhR6MCzVkCZaukmbWLNZ2FEDahiyDm+T2tob3Os23zcHKsI6NvkfRpwyNMxIhKOEQxQ7OW5FkFi8mQjszpUproy7HZibbaAifBr10KPHnDOOW72OTPfAz0ygdCZ6kqz2xlgt1vJi15QU+fg5JcxDBpf55frexsavEtEHhQRt9NHPmA1kCki/UUkEkNQ/6DZMR8ANztfXw0sVFV1br/OGe3VH8gEVvnJTos2kpRkhLyGSr6JTYSRsbFsCKAzGTPGfQ+TLysqSAkPp28AxHcX2dnGOlScCRhTXUtLSnBYukmLvPlZFVSHM73APykN3k5zPQnMAvaKyEcico2I+Ky2g1MD+QHwMbAV+JeqbhaRx0XkUudhfwO6icgujKm2+53nbgb+hZEDMw+4Q1WD05TCwi35+bBypTHdFQpkxcWxsaIC9fOXU00NbN7sWXxfV1HBqLg4xJ2n8RM9e0L37qHnTIobGtgcKglNQeKtBUYc1Xdn+6ebiLfTXC+pagEwDCgEfgMcFpFnfGWIqs5V1cGqOlBVn3Rue1RVP3C+rlHVa1R1kKqOU9U9Tc590nneEFX9yFc2WfiG/HxDgF+9OtiWeEdWbCwnGxo44mfvt2mTkYfjzpk0OBxsrKxkVAD1EhchJ8K7dBNrqqtFVi8PI6JfNdnp/hnptknSV9XdqvoIMAH4AqssvYUXTJlirENFN8kOUFmVlsT3ndXV1Dgc5ATJmWzebGTDhwLpdjt9o6KsfJMWqGtQjq2NYeD4Wr/do621uS4TkfeAfUAS8D1/GWbReUhNNfIXQkU3yXJ+gfs7eXHtWkNTysg4e986572DMTLJzjY6Lu7Z0/qxZkBEyE9KYklJid+nJkOVd7+ogopwCvxYTdFbAf73GBFSf8TQNLJVdZKqvuA/0yw6E/n5Rjn6AOUCdohuERH0jIhgm59zTVoS39dXVBApwtCYwCfHhlpEFxhTXcfq69keCn2Hg8CbTr3k237SS8D7kUlfjDpc6ar6gKru8JtFFp2S/HyoqPh6asfsDI6JYYcfMy3r6w1doiXxfURsLJG2wCfHjhgBNtvXBShDgfykJMDSTTzxxXIh/JwaxmcG2Zk4he+PVDVEZlEtzMZU5/A6VKa6BkdHs8OPT7lbt0JtbcvOJBh6CRjtA4YPD52ACYBB0dH0joy0dBM3OBzK4cIYMsb7voRKU6yaIBYBoVcvo894qIjwg2NiOFpfT6mf5uVaEt+P1NZytL4+KHqJi7w8w5mEigQhIkxNTLR0Ezf8Z001WhJBvp+7T1nOxCJg5OfDsmXQGAJZQIOjjcZnO/00Olm71miVm5l59r71ziiyYDuTEydg//6gmdBm8pOSOFhXx54a/z6Bhxr/nG+EuN80y39TXGA5E4sAkp8PpaVfdxY0M4Odwre/dJPVq41RiTtJxBXJlRPA1sHNycsz1qE01VXg1E2sqa4zWbFUCOtZy9Sh/q2k4G00l9sBkohM8a05Fp2Z/HxjHQpTXQOjoxHwi25SV2c0xBrnoVHCuooK0qOiSIrwWZGJNpOdDZGRoeVMhsbE0CMigsWWMzmNw6EUrYomfXw1Npt/Kyl4OzL5r4ft7/vKEIvOT1qa0YBp8eJgW9I6UTYbGXa7X0YmmzYZ4rvr6b85651lVIJJZKRRU21VCFW5ExEKkpJYbOkmp5n7ZTV6KpKp+f6/l7fO5CyXJiLxgBXdZdEmCgqMiK5Q0U38MTJxfUG7G5lUNTayvaoq6M4EDPvWrAmNz8pFflISB2pr2WvpJgC8/omhl3zTj/klLlp0JiKyU0R2ANEisqPpgtGQapHfLbToVBQUQElJaNR+cuWa+Popd/Vq6NbNfeb7JmcPEzM4k7w8Izdou68aTQQASzc5k8+WCGE9aikY5v/K0601Mn8CY1TyZ4zKwS4cwBFgoZ/ssuikNNVNRo8Ori2tMTg6mvLGRo7W1fm0B/uqVcZTv7vM99Piu0mcCRjOb/jw4NriLcNiYuju1E1u6d072OYEFYdDObg6mvRJ1dhsQXYmqvp3ABHZpqqf+90ai05PWhoMHGjoJnfdFWxrWqZpRJevnElFBWzZAlde6X7/+ooKEsLCyLD7f1qiNYYMgbg4w5ncfHPrx5sBV50ul24SyPL9ZmPeumocxTFMzQ9MiRlvM+A/FxG7iGSJyMSmi78NtOh8uHQTs1eldeWa+FI3WbvW+Lk9RXJtqKwkO8A9TDwRFgZjx4ZWRBcYU11f1dayr4vrJq99YlQI/uaswDRX8zY0+FLgMEYf+OVNlmX+M82is1JQAKdOmb+QYF+7nSgRdvowosslvruL5FJVNlRUBDW/pDl5eUaNrlBpbAaWbuJi+RIbtu61nDciMKNcb6O5fgf8HIhTVVuTJayjBohIiojMd4r980Uk2c0xo0RkpYhsFpENInJtk30vi8heEVnnXEZ11CYL/+LSTcweIhwmwsDoaJ+GB69eDenpRt/35uyvqaGssZFsE+glLvLyDEdidsfflOExMaR28XwTh0M5uCqafgHIL3HhrTPpqapzVNUfk2/3A5+qaibwqfN9c6qAb6nqCOB8YI6IJDXZf4+qjnIuIVTrtGvSty8MGGB+ZwLOiC4fTnOtWuU5v8TVdz7bZCMTCL18k/zExC7tTD7eUIOjODD5JS68dSafiMgEP9lwGfB35+u/A5c3P0BVd6jqTufrQ8AxoLuf7LEIAAUFRkRXKOgmu6qrafRBePDx47Bvn2e9ZH1FBQKMNJEzycgw+sKvWBFsS9pGQVIS+2tr2efHNgJm5rWPDb3kpgDpJeC9M9kHfCAifxCRB5suPrChp6oeBnCu3UwAfI2IjAMigd1NNj/pnP56WkQ8/vZE5DYRKRSRwuPHj/vAdIv2ct55hm5i9nyTwTEx1KnylQ/EXJeQ3dLIZGB0NHHhrUXsBw4Ro31AqLQOcHFesjFbvqiLjk6WLTLyS6aPDFxUoLfOZCywGRgJzGyyzPDmZBFZICKb3CyXtcVYEekNvArc0qS3ygPAUCAPSAHu83S+qj6vqrmqmtu9uzWwCSbnnWesF5o8U+l0RJcPnnBXrza+nMeOdb9/fUWFqaa4XEydCl99ZYyqQoXhzjpdC7ugM2loVIq+iCb93MDpJdB60iIAqnpeR26iqh6djogcFZHeqnrY6SyOeTguAfgQeLhpzotrVAPUishLwE87YqtFYOjTx+hvsnAh3H13sK3xzOlck6oqZqekdOhaX3wBw4ZBfPzZ+yobG9lVXc1NPXt26B7+oGljM3dZ+2ZERJiWnMzCU6e6XL7J+4XVaEkM5xUEtoWx1yXoRSTMmVtyrfN9jIhE+8CGDwBXStTNuCkeKSKRwLvAK6r6VrN9vZ1rwdBbNvnAJosAMG2aoZvU1wfbEs/0iIggISyswyMThwNWroSJHjKzNldWophLfHcxciQkJ4feVNe0pCQO1dX5tf2yGXl9nhHHfcsFgU189TbPZCDGl/Rc4G/OzbOAv/rAhl8DM0VkJ8bU2a+d98wVkRecx3wDmAp8200I8OsishHYCKRilICxCAGmTTMywtesCbYlnhER+tvtHU6A277dqEnmyZmsd5ZRMVNYsAubDaZMCUFn4tRNFp46FWRLAsvKxTbC+9QwaYgJnQnwB+ANDE3C9Ry5GOhwPxNVPamq01U107kudm4vuwZjTAAAIABJREFUVNVbna9fU9WIJuG/p0OAVXWaqmap6khVvUlVKzpqk0VgKCgw1otMXi403W5nfwediSsa6txz3e/fUFlJvEnKqLhj6lTYuRMOH279WLMwwG6nX1RUl9JN6hqUI4UxDJoY+Ox/b53JOOBJp+itAKpaAiS1eJaFRQt07w5ZWeYX4TOcI5OOVA9esQJSUgydyB0bKirIio3FZtK5/aa6Sajg0k0WnTqFo4v0N3ljeSVUhDNzeuDv7a0zKaOZ4xCRc4CjPrfIoksxbRosX240izIr6XY75Y2NlDQ0tPsaK1caoxJ3bXpV1YjkMuEUl4vRo42ij6HkTADOS0riZEMDG50JoZ2dNz4xJo5uuzAm4Pf21pm8A7woImkAItINmIMx9WVh0W6mTYOaGvjcxDWpXVNP7dVNioth61bPesmB2lpKGxtNVZOrOeHhMGlSaDoTgEVdRDdZtSSMqAFVjOwbGfB7e+tMHgEqgK8wRijHgFrgl36yy6KLMHWq8bRu5qkulzNpr27icpQe9RITi+9NmTrVaDl84kSwLfGevnY7mdHRXUI3qahxcHJtLEMmBacqp7cl6KtV9QaM7PRxQC9V/aaqdu0azxYdJikJxowxtzNJ7+DIZMUKo5y7x57vzimYLBOPTOBr3WT58uDa0VamJSWxpKSEBrPX7ukgryyqhJowLpweHN3N6zwTAFU94YyysmqRWPiMadOMp/cKk8bhpYSHExcW1m5nsnIl5OQYmoM7NlRU0N9uJ95EZVTckZcHdntoFOhsyrTkZMoaG1lj1j8wH/HWJw0gyq3nB14vgRacibOs+57WlkAaa9E5mTEDGhrMOx8vIqRHRbG/HVECDQ1G5rsnvQScDbFMPioBiIoy8k3mzw+2JW1jWlISAswvLg62KX5lzaIIYodVMbBnRFDu39LI5GEMreQR4CUgGvgHRi/4fwJ24EV/G2jR+Zk82fiiMvOXVEY7Exc3boTKSs96SU1jIzuqqkyvl7iYOdNoO3zwYLAt8Z7UyEhGx8UxvxOL8EXFDZRvjCUnP3hdzDw6E1V93bVgZKZfoqoPq+rfVPUh4GKMLHgLiw4RHW088X7ySbAt8Ux7ExdXrjTWnkYmW6qqcGDOMirumDnTWC9YEFw72sqslBRWlJVR3oHwbjPz57lV4BCuuSB4U6XeaiajgOZNpzY4t1tYdJhZs8z9xJtht3OqoYGyNn4ZrVgBvXsb3RXd4YrkygqRkUl2ttEl0syjSHfMTE6mQbXTtvL9zzwHRDfynRnBeyjx1plsB37cbNtdwA7fmmPRVTH7E297woNVjUKWkycbpefdsaGyErvNxqBoX9RM9T82m6FxzZ9v/sZmTZmUmEi0zdZpp7q2LYuie24lCdFtiqnyKd7e+Q7gPhHZJyJLRGQfRnvdO/xmmUWXIjvbKK9i1ife9oQH790LRUVf97x3x8bKSkbGxhJm0jIq7pg5E44dC62+8FE2G1MTEzulM1m5o4b6r6KZdF5jUO3wNs9kFTAAeAijp8hDwEBV/cKPtll0IVxPvAsWGE/0ZqM9I5MlS4y1q6ClOzaYtCFWS7hGkWZ1/J6YmZLC1qoqinzQNdNM/OU/xs/z7YsD16LXHV6PiVS1zCnIP+Vcl/rTMIuux8yZcPSoOZ94u0dEEG2ztWlksmQJpKbC8OHu9x+tq+NYfX3IRHK56NPH+JlCzZnMcpak72yjk4XzBVtqLZeMDe5UqVfSv7Px1HVALnBGnzhVvc0Pdll0QZo+8WZnB9eW5ogI6W0MD1682Mga96iXuMT3EBuZgPFZ/eUvRl01k1bNP4uRsbH0ioxk/qlT3NK7d7DN8QkNjUrR5zEMOK8Smy00RiZ/Bv4I9AUimi0dQkRSRGS+iOx0rpM9HNfYpDHWB0229xeRL5znv+nsymgRgqSlGW1tzfrE25bExf37jaU1vQRC15nU1IRWaRURYUZyMgs6UUn6Nz6rQksjmDUz+Jqbt87kGmCcqn5DVW9puvjAhvuBT1U1E/jU+d4d1U0aY13aZPv/Ak87zz8FfNcHNlkEiZkzjUx4M05rtyVx0Vu9pHdkJN0jQ+/5Jz8fIiJg7txgW9I2ZiYnc7y+/nRny1Dntf8aSYq3XxL8aEBvnUkVRsVgf3AZ8Hfn679j9HH3Cuf02zTg7facb2E+Zs+G6mpzllZJt9s5UV9PZWPrUTOLFxvNsEaO9HxMqJRRcUdcHFxyCbz8MlRVBdsa73HpJh91ktIqKxdEED20Migl55vjrTN5CnjU+eXta3qq6mEA57qHh+PsIlIoIp+LiMthdANKVNWVSVYE9PGDjRYBoqDAmIM34xNvWyK6liwxsvrdNcMCaHA42FJZGTLJiu646y44dQpefTXYlnhPr6gocuPjmXvyZLBN6TB7j9VTtj6WMdODV0KlKd46kx8B9wGnRGRH08Wbk0VkgYhscrNc1gZb+6lqLnADMEdEBgLunJvHyVARuc3pkAqPH7cKH5uRmBijivCHHwbbkrPx1pkUFcGePS3rJTurq6lVDdmRCRjJmGPGwDPPmDOc2xMXpqSwsqyM4vr6YJvSIZ55zyih8q3LglPYsTneFnJ5oiM3UdUZnvaJyFER6a2qh0WkN0bjLXfXOORc7xGRxcBo4N9AkoiEO0cnacChFux4HngeIDc3N4T+/LsWF15ojEx27oTMzGBb8zXeJi56pZc4xfdQCwtuigjceSfcfLMRNDErRCr1XditG4/v///tnXd4VGXWwH8nFUhISIEo3SgqKCKKooiuoiBgQRAFRRdF7H5Y97O7roUPuytWRARXsayKYltQFtu6FlSCICpBBEIvCaRA6vn+eO/IGKYlM8ncGd7f88xzZ249b+5kzn1PXcncrVsZnZcXbXEazbvvK5JZzQUnuOOBJNSkxRn+XhGQYTYw1nk/Fni7/g4ikiUiqc77XOAY4EdVVWA+MDLQ8ZbYYsgQs3SbqWuvlBRSRIIqk48/hszMwOHNi8rKSBLhwFbR6T0RKUaNgrw8MzuJFfq0bk1ucjLvx7DfpKZW+fWTVnQ5tpyUpOhHckEDkhZFpL+ITBGRd5zPh4vIcRGQYRIwUESWYaoTT3LO30dEpjr7dAcWiEgBRnlMUtUfnW03AteJSCHGh/JcBGSyRJH8fDjwQPcpkwQROgepHqwK8+aZ/JLERP/nKigr44CWLUn151SJEVJT4YorzL36+edoSxMaiSIMzs7mg61bqY0l+5wXL31ajpakMHSoOxQJhKhMRORczAxiJ+BRIArcFa4AqrpFVU9U1W7OcquzfoGqjnfef6GqPVW1l7N8zuv4X1X1SFXdT1XPUtWGdzCyuI6hQ80TvmMNcg0dU1MpCpBrUlhoanKdfHLg8ywqL6dXDJu4vLnsMkhJgYceirYkoTM0O5vN1dUsKC2NtiiN4vm3ayBBufqM6IcEewj1sehWYJCqTgA8tUIXAwc1iVSWPZ6hQ6Gqyn294TulprI6gDLx9GQJpEy2VlezurIybpRJu3Zw6aUwbZppIxALnJydTQLEbFTXgg+TSe9Zzv57Rz8k2EOoyqS9qi5w3nvmhTVAgIm8xdJ4+vc3uQxuM3V1Sk1lTWWlX/PInDnGTLfffv7P4SmjEi/KBOCOOyAtDW68MdqShEZ2cjJHZ2TEpN9k8eoqyn9Mo+9J7opGC1WZLBeR+r3i+mH6nFgsESc11VQRfv99d4WddkpNpRZTpLE+VVUwf37wqKbfI7liOCy4Prm5cOut8O677ptN+mNoTg4LSkt93ks389gskyU67gz3zEogdGVyD/C2iNwGJIvI9Zg+8GH7TCwWfwwdCqtWuauKcMdUU0zPl6nriy+grCy4v6SgrIy2ycnsFYNlVAIxYQJ07gw33BAbjbNOyckB4N0YM3W9/04CCbmVjOrnrkjAUEOD38IkC/YFVmJKmIxT1Q+aUDbLHs7pp5tchlmzoi3JLjo5uSarfUR0zZkDSUkm6TIQBWVl9EpPp2kKSkSPFi1g4kT4/nt48cVoSxOcQ9LS6NqiBbNiKIF5S2ktaz5Po/ugHSQmuOv705B+Jh+q6mmqerCqnqKqLq3taokX8vKgXz+XKZMAM5O5c+HooyEjw//xNXV1LC4vp1ccmbi8Oecc6NvXJDOuWBFtaQIjIgzPzeXD4mK219QEP8AFPDyrHHYmcv5I94WUNyTPpLOI3CQijzvLrk0nlsViGD4cCgpMeRI3kJWURMuEhN3CgzduhO++C+4v+b2MShw5371JSICZM837kSPdWf3ZmxG5uVSpxkzhx1ffUCSjmquGuu/7E2qeyWCMs/0UINNZLnXWWyxNxvDhZumW2YmI+AwP9vRgCcVfAvEVyVWf/Hx44QWjXK+5JtrSBObozEzaJSfHhKlrR2Udv/47jX0GlJOWGrszkweAi1T1WFU9X1WPBcYBMZSmZIlF8vOhVy/3KBPwnWsyZw7k5JjCh4EoKC8nWYTuMV5GJRinnQY33WS6Mb7wQrSl8U+iCMNyc3lv61Z2htBaIJpMfq8cLUvi7BHu8pV4CFWZdAVeqbfuVaBzRKWxWHwwfLiJlNqwIdqSGDq1aPEHB3xNDXzwgTFxBSqhAmZm0r1VK1JivIxKKNx9tyl2edllxlTpVkbk5lJWW8u8kpJoixKQf7xeCy1ruX64O2e1oX6jPwaOr7fuT8AnkRTGYvHF8OEm1+Rtl5Tw7JiayrqqKmqc+Nf//Ac2b95lkgvEorKyuPWX1CcpCV55xTQJGz4c3OqWGJCVRUZiIm+62NRVU6v8OLcV7fuXkZvuzlzxUJVJITBLRGaIyN9EZAbwFvCLiNzieTWdmJY9mZ49jbnLLaauTqmp1AHrnGS3WbNMkqWn2rE/tlRXs6aqKm4juXyRlwevv256vJx3njvzT1ISEjglJ4fZW7b8/oDgNqb/u4K6LSmcPtxFGbz1CFWZHAp8hzFrHecsv8P0FBnovPz2LLFYwkHEPNnOmwfbtkVbmj+GB6saZTJwoCn/Eog9wfnui6OOgsceM6bAm292V0UDDyNyc9lcXc3nbviC+WDKK9WQVMeNZ7n3QSSk5liqekJTC2KxBOKss0xV2jffhAsvjK4s3lnwLReaLP077gh+XDzW5AqVSy+FhQvh/vthxw549FH/LY2jwZCcHNISEnh540aOd/rEu4XqGuW72a1oe0wpXXMzoy2OXxqSZ5IoIv1EZJTzuZWIuKf+sSWuOfJI2Hdfd2RWe2YmRZWVzJplfhRPPz34cQXl5eQlJ9MuzsqohIIIPPkkXHstTJ5sTF5uKomVlpjIGbm5vLZpE5UuM3VNfqec2s0pjBrjLrnqE2qeyb6YkvPvs6v51CDg2SaSy2L5AyIwZowppLhmTXRlyUxKIj0xkdU7dzJrlqlw3LZt8OO+LS3lsNatm15Al5KQYGaXkybByy+b8GFnsuYKzsvLo6Smhg9cVqvrmRm1kFbDX89x93cn1JnJZExocDbgqXv8MXBsuAKISLaIfCgiy5zlbnNMETlBRBZ6vXaKyBnOtukissJr26HhymRxJ2PGGHv7K/WD1JsZEaFjaipLlymLF4cWxbWjtpYfy8s5fA9WJmAeCm68EZ57Dj76CE480UTCuYGTsrJol5zMi26JQQeKS+v4ZW46+w4uJTc9JK9E1AhVmRwJ3KuqdTj9TFS1BGgTARluAuapajdgnvP5D6jqfFU9VFUPxRSZrADmeu3yF892VV0YAZksLmT//eGII9xj6vpxjkk8POOM4PsXlJVRCxy+B/pLfDFunPF/FRTAsccav1O0SUpIYHS7dry7ZQsl1e7oFXL3zHLYkcglf3aRg8kPoUq4nXqKQ0TaA5FQ4cOAGc77GUCwf82RwAeqWhGBa1tijDFjjCN3yZLoytEpNZUN8zLp3Ru6dg2+v6c97J4+M/Fm2DBTHHPdOjjmGHd0aRyTl0elKm+4ZLo080VIyNvJNUPd/70JVZm8CUwTkY4AIpIDPMruWfGNIU9V1wE4y3ZB9h+N6aXizb0iskhEHhGRVH8HisglIrJARBZscnGCksU/o0ebLPOXXoquHK3WplG9pDUjR4XmFP22rIx2ycm/R4JZDMcdB598YqoIHHss/Pe/0ZXniNat6dayJS+5wNS1bE0NG75Ip/cZFaQkxc/M5HagDFiFmaFsBCqBiaEcLCIfichiH69hDRFWRPYGegJzvFbfDBwIHIHx6fhtHKqqU1S1j6r2aRuKx9TiOvLyTAfGmTOjmwD369ttIEE5caT/fvDefFtayuGtW8ddD5NI0KuXqSKQnW18KO+9Fz1ZRIQxeXl8XFJCUZRLHt8xrRzqhBvGxUb0X6jNsXao6rlAW4z/ZC+n4GNIf21VPcnpg1L/9TawwVESHmWxMcCpzgZmqervBk1VXaeGSuB5Rz5LHHPeebBypXmijQZ1dfDNG62gz1aqcoIrkwrrfA9Kfj58/jl0776rSGS0QofHtGuHAtPXr4+OAJhAk9kzUkg9sIxRR7g3UdGbhuSZ7A/sD2xX1UjaiGYDY533Y4FAFZjOoZ6Jy0sRCcbfsjiCsllcyJlnQlaWyVuIBvPnw6Y1iXDy+t36mvjCOt9DIy8PPv0Uxo+H++4zfpRly5pfjv1ateKkrCyeWbcuauVVpr6/g4rlLRk2vipmZrNBlYmIDBORVcBS4D+YPiarPKG5EWASMFBElmHKskxyrttHRKZ6ydEV6MTuxSVfEpEfgB+AXEy/eksc07KliQaaNQvWrm3+60+fDpmZCv23+Oy4WJ9vHed7HzszCUpaGkyZAm+8AcuXQ+/eMG1a85dgubJ9e4oqK6PWH37iIzWQWcUjl8TOdyagMhGRw4B/YsJwTwZ6OMu5wGsicni4AqjqFlU9UVW7OcutzvoFqjrea7/fVLWDE57sffwAVe3pmM3OU1UXpUFZmorLL4faWni2mdNmt283P3SjRwsZLSU0ZeI43ztY53vIjBgBixaZygcXXQSjRkFxcfNd/9ScHDqlpvJEFJ5WFhXW8Nv8dA4ZtZ32rZOb/fqNJdjM5BrgPlUdr6ofqepPznI88H/Odoul2dl3Xxg82DRfas6UgNdfN7Wlxo6FLi1asGLHjqDHWOd74+jY0XSwnDTJzEJ79TJmsOYgKSGBS9u356PiYn6uaN4shOseMt+pide0aNbrhkswZdIPmOJn27POdoslKlx5pclRaM4+J9Onm+TJo46Cg9LSWFxeHnD/itpalljne6NJTDQZ8198Ycr8n3AC3H578zxAjN97b5JFeKoZ6/eUlyvzZ7Yk6/gSTukeWz62YMokR1VX+9qgqkVATuRFslhCY8gQ6NKl+RzxCxbAZ5/BxRebsiA909JYWVnJ9poav8cUlJVRh/WXhMsRR8D335sZ4T33mJyUpk5czUtJYWTbtkxfv57yZmrpe8fUCuq2J3HZVe4u6uiLYMok2HY7b7dEjcRE4zuZPx9++KHpr/fAA5CRAZdcYj73dJpcBZqdeJzvNpIrfNLTjTP+1VehsNA45++4A5oyHeTKDh3YVlvbLGHC1dXwzCNJJO5Xxu2nRaJSVfMSTFmkendSrP8CYiObxhK3XHwxtG4Nd93VtNf59VfjL7n8cqNQYJcy+SGQMrHO94hz9tmwdKmphnD33XDooca30hT0y8igX0YGk1atYmcTz07ufGon5StTOft/y2mZ5M7WvIEIpky+ZFcnRV+vL5tUOoslCNnZcM015oe+oKDprvPww6an+YQJu9Z1adGC1omJ/OCnjrqq8nFJCUdnZFjne4Rp2xZeeAHmzDGlWAYNgpEjI18wUkS4a599KKqsZOq6dZE9uRcVFfDQPYkkHrydpy6ITe9BQGWiqser6gmBXs0lqMXij2uvhcxMuPPOpjn/pk3GvHL++dC+/a71IsLBaWl+ZybLd+zgt507GZid3TSCWRg0CBYvNn6U99+HAw+Ee++FECK2Q2ZAmzYcl5nJxFWr2NFEs5O/PLCTyk3JjLujgsxkd5ea94f7q4dZLEHIyoLrroO33oJvv438+Z94woQDX3/97tt6OspEfWTVfegkRgx0WRvYeKNFC7j1VmP6GjoUbrsNDj7Y9JyPBJ7ZybqqKp5pgryT4mKY8lASyUdv5aHhuRE/f3NhlYklLrj6aqNUIj072bwZHnvMtOXt3n337T3T0iiuqWGtj0JSHxYX0zk1lW4tbXfr5qBLF2PunDPHBGcMHWoi/iIR9fWnNm0Y0KYN/7dqVcQju666s5KaskSu+mslrZNic1YCVplY4oTMTLjhBnj3XRO+GyluvhlKS2Gin/rYPZ0orfp+k1pV/l1czMCsLOsvaWYGDTLZ8w8+aEraH3IIXHZZ+O2e/9a1Kxurq3lwtc9siUbx2efKzMdTaDFkI3efFKz7hruxysQSN0yYYBpVXXSRMUuFy5dfwtSpxsF/0EG+9/FEdC2q5zdZUFrKttpa6y+JEikpxixZWGiSW597zlRNmDCh8fXc+rdpw+h27bh35Uq/QRcNoaQEzjinFvbaycOPKGmJsRfB5Y1VJpa4IT3d/GgsW2bs5uFQWwtXXAEdOphcBn9kJyfTPiVlNyf8h1u3IsCJbWIvXyCeyM01ZsqffzatC5580pS7v/rqximVyfvtR1ZSEhf89BPVYVQUVoXzLq5h67oE+t1XxGXd8hp9LrdglYklrhgwwJg0HnnENFxqLE89ZTKuH37Y5LEE4pD09N2eVD8sLqZ3ejq5KTYVyw3k55tZ5i+/mNbPTzxhZipXXw0NsVrlpqTw1P77811ZGQ+EYe56bpry3utJpF60ktfO6hwXplCrTCxxx/33Q+fOcOGFJn6/oRQUwC23mI6OZ50VfP+eaWksraj4/Um1rKaG/27fbqO4XEh+vpm9/vILnHuuUSr5+XDBBaE76ke0bcuotm2587ffWNQIc9eLL8KllwK9i3ny9hZxk9BqlYkl7mjd2uSFFBaaUuYNKbexcqWJAMrMhOefNzW4gtEzLY0qVZY5jppPtm2jWtX6S1yMR6kUFhpz5j//acKJBwyAl18O/p2Z3K0bOcnJDF60iGUNeGJ54gmTr1R3SAmDn1zLhR32CnMk7sEqE0tcMmCAMWvMmQPDh4emULZsMWXtd+yAf/3LlEAPBe+yKnWqvLlpEy0SEjjGU3fF4lq6doW//91kzt97L6xYYWYsHTqY3KWlS30f1zYlhY969aJalRMWLmR5kIiPnTtNLsxVVwH9NtPvyZX8s88BcWHe8hB1ZSIiZ4nIEhGpE5E+AfYbLCI/i0ihiNzktX4fEflKRJaJyKsiYo3UFsB0Y5w61SiGESMgkEXip5/MjGTFCpg923/0li+6p6WRCDyzdi09vv6aaevXMyI3lxYxHp2zJ5GTY0yby5fD3LnmYeTxx6FHD+jf3/jgliz5Y8fHg9LSmNerFzvq6jhh4UJ+8lEJQdX0YjnoIGXiRJDB6+n78Cr+1edg0mM4p8QnqhrVF9AdOAD4GOjjZ59EYDmQjykuWQD0cLa9Box23j8NXB7KdQ8//HC17Bk8+6wqqObkqN5zj2pJya5tJSWq112nmpSkmpGhOmtW465x8NdfK/Pn62HffKOvbNig1bW1kRHeEjU2bFC9/37VHj3M9wdU27dXPfdc1aefVl2yRLWyUvW77ds1+7PPNHH+fB23dKmuqKjQNWtUH3tM9aijzHGt8iuUBxZq3wULdFt1dbSHFhbAAvXxmyra3M2V/SAiHwM3qOoCH9uOBu5U1ZOdzzc7myYBm4C9VLWm/n6B6NOnjy5YsNulLHHKV1+Z+k3vvmt8KtnZUF4O27aZQoHjx5vt7RqZN1ZQVsa2mhqOzcyMK9OFxbBqlalMPHeu6fboXZE+IwOycpTSxCqKa2rQWoGilqBCSn4FVacWsdeZm7lt3y6M33tvUhOibhAKCxH5VlV3syLFyjyrA+Adh1cE9MU05ypR1Rqv9R38nURELgEuAejcuXPTSGpxJX37wjvvmHDfp582Nuy0NKNYRo2Cww4L7/y9bL+SuKZzZ5MMe9FFZo5SWGhCz4uKTCHQzZuFqqpUdtQlsaxiB0nDNtFx4DYy9q3imIwMLm3fl5ZxbvZsFmUiIh8BvsIWblXVUJqu+nrU0wDrfaKqU3DaEPfp08cdUzJLs9K7t+kbb7E0FhHo1s28dicRSHdesV0epaE0izJR1ZPCPEUR0Mnrc0dgLbAZaCMiSc7sxLPeYrFYLM1IrBjvvgG6OZFbKcBoYLbjDJoPjHT2GwuEMtOxWCwWSwSJujIRkeEiUgQcDbwnInOc9e1F5H0AZ9ZxFTAHWAq8pqqefNUbgetEpBDjQ3muucdgsVgsezquieZqbmw0l8VisTQcf9FcUZ+ZWCwWiyX2scrEYrFYLGFjlYnFYrFYwsYqE4vFYrGEzR7rgBeRTcDKRh6ei8lx2ZOwY94zsGOOf8IdbxdVbVt/5R6rTMJBRBb4imaIZ+yY9wzsmOOfphqvNXNZLBaLJWysMrFYLBZL2Fhl0jimRFuAKGDHvGdgxxz/NMl4rc/EYrFYLGFjZyYWi8ViCRurTCwWi8USNlaZNBARGSwiP4tIoYjcFG15Io2IdBKR+SKyVESWiMjVzvpsEflQRJY5y6xoyxppRCRRRL4XkXedz/uIyFfOmF912h/EDSLSRkReF5GfnPt9dLzfZxG51vleLxaRl0WkRbzdZxGZJiIbRWSx1zqf91UMjzm/Z4tEpNE9R60yaQAikgg8AQwBegDniEiP6EoVcWqA61W1O3AUcKUzxpuAearaDZjnfI43rsa0OPBwH/CIM+Zi4KKoSNV0/B34l6oeCPTCjD1u77OIdAAmAH1U9WBMW8TRxN99ng4MrrfO330dAnRzXpcATzX2olaZNIwjgUJV/VVVq4BXgGFRlimiqOo6Vf3OeV+K+YHpgBnnDGe3GcAZ0ZGwaRCRjsApwFTnswADgNedXeJqzCKSARyeW5E/AAAH60lEQVSH0/9HVatUtYQ4v8+Y7rItRSQJaAWsI87us6p+Cmytt9rffR0GvKCGLzGda/duzHWtMmkYHYDVXp+LnHVxiYh0BXoDXwF5qroOjMIh/hpcPwr8L1DnfM4BSpzGbBB/9zof2AQ875j2popIGnF8n1V1DfAgsAqjRLYB3xLf99mDv/sasd80q0wahvhYF5ex1SKSDrwBXKOq26MtT1MiIqcCG1X1W+/VPnaNp3udBBwGPKWqvYFy4sik5QvHTzAM2AdoD6RhzDz1iaf7HIyIfc+tMmkYRUAnr88dgbVRkqXJEJFkjCJ5SVXfdFZv8Ex/neXGaMnXBBwDnC4iv2FMlwMwM5U2jjkE4u9eFwFFqvqV8/l1jHKJ5/t8ErBCVTepajXwJtCP+L7PHvzd14j9plll0jC+Abo50R8pGOfd7CjLFFEcX8FzwFJVfdhr02xgrPN+LPB2c8vWVKjqzaraUVW7Yu7pv1V1DDAfGOnsFm9jXg+sFpEDnFUnAj8Sx/cZY946SkRaOd9zz5jj9j574e++zgb+7ER1HQVs85jDGorNgG8gIjIU89SaCExT1XujLFJEEZH+wGfAD+zyH9yC8Zu8BnTG/FOepar1nXwxj4gcD9ygqqeKSD5mppINfA+cp6qV0ZQvkojIoZiAgxTgV+BCzANm3N5nEfkbMAoTtfg9MB7jI4ib+ywiLwPHY0rNbwD+CryFj/vqKNXHMdFfFcCFqrqgUde1ysRisVgs4WLNXBaLxWIJG6tMLBaLxRI2VplYLBaLJWysMrFYLBZL2FhlYrFYLJawscrEYvGDiHQWkTIRaR9tWXzhyHZ0uPtEGhE5WUQ+C+P4g5xKxqmRlMvStFhlYok6IvKxiFSKSKmIbBORX0XkHyJyeDTlUtVVqpquqq7MiHZk+y+Y/BgRqQm0T3Pg5C08gslt8KwbLiKrRWSDiFxeb/+nxWlz4EFVl2DyPa5qDpktkcEqE4tbuFtVW6tqJnACsBL4UkSGR1kuS8MYhEmCnA+/t214FjgNU3X7HhHJc7YdD/QEJvs4zzTgf0TE/kbFCPZGWVyHqq5U1duAF4DJztMuIpIjIi+IyDoRWS8iM0Qk23OciPwmIreJae5VJiI/iMghInKO0/xnm1MdN8nrmOedp+ZSEflRRM712tZVRNQpT4+I3Cki80RkopjmQxudjGqfeB0/XkR+EZESEXlbRNp57RNsTBNEZIUj3xoRmei1TUWkv2OG+wBIdMZdJiJjvffxOuZMESlw/hYF3spaRC5w/k4TRKRIRIpF5BlHIYTKGcBHuisbOgeoVtWFqroSWAZ0EZGWmMzr8apa5+M8nwJ7AYc24NqWKGKVicXNvIIpdeGpH/USkIVpTNYdUy7iH/WOGQtc4exXAMzCzHR6YZ6CTwfO9tr/c8wPVhvgLmC6BG54dhymHEV7zNP2LSJyTJBx/Nk5rjOmRM2LXtv8jklE9gcmAaeqamvgIHzUgnPMcEOAWsesla6qM+rv5/hOXsJUB87BlMl5WUT6eu3WBcgD9gWOAM7C1CsLlcMw9a48bAJqReRwEdkH6AoUAvcAM1V16e6nAKecyTLnfJYYICn4LhZL1ChyljnO0/fJwP6qWgwgItcBP4nI3l7F6aZ4fqBEZCYwBjhKVcuBchH5GPMjORNAVZ/zut4rInIDpq6R9w+iN7+o6tPO+69EZCHQB/hPgHH8zSmsiIj8BVjm5dT3OyZM/SgBDhKRlU7zqi8DXCcYFwJvqOoHzuf3RGQWMA5Tew1gB3CHqtYChSIyzxnfSyFeIwv4vWWBqqqInIPp6pgCXA7sB/wJOE5EHsUojOXABKchm4ftmHpZlhjAzkwsbqajs9zCrjLZK7y2L3eW3iW0vSueVmCe1jfVW9caQEQSROQuEfnZMfuUYGYwbQPIVL+iarnnfAH4zcf7jgQZk6r+ilGGFwNrReRzERkU5FqB6IQp6OjNcv7499voKBIPoYzPm2Igw3uFqn6mqv1V9UjgPWAKZkznA1mqehzm71q/n0oGu3cMtLgUq0wsbmYUsAb4mV3d4Lp6bc93lt6d4hrCOZiqsWdiftTaYExjvhoGhUNXH++LCGFMqvqmqg7EmL9eA94WkVY+ruHL71Cf1ZjGUN7k0/i/ny++x5js/HEb8J6qfo9R3F846z/FdPUEQEyLh27O+SwxgFUmFtchIp0cx/YFwNVOf+q1wFzgIRFpI6Zr3kPAB43tv4B58q3B2PUTRGQc5gcu0twuInli+q7fB8xT1bXBxiQiB4jIYEd5VGPazCq+Fcd6jAO+vrLwZjpwppg8kEQRGQKMAJ6P2EhNqfMTfW0QkUMwnQ7vclYtBwY6ARFDML4UD8dhyqdbZRIjWGVicQu3OxFL2zFPqfsB/VT1Da99zgNKgZ+cVwnGud1YZmB8BYWYGVAPTC+XSPOic97VGL/BeV7bAo0pBZOvsc5ZPwE4U1V31r+Aqv4CPAl87USNne9jny8wAQoPYsxR92N6d4Tjh6nPHKBGTNjv7zgRYVOBy716hTyDGeMWjBKf6HXIOGCyn0gviwux/UwsliZCRLpi/CGdVLUo8N7xg4gMBm5xfCGNOb4HJgrvkFhuUrWnYZWJxdJE7KnKxLJnYs1cFovFYgkbOzOxWCwWS9jYmYnFYrFYwsYqE4vFYrGEjVUmFovFYgkbq0wsFovFEjZWmVgsFoslbP4fC3BksFFjATwAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "w0 = seq.warps[0]\n", "w1 = seq.warps[1]\n", "\n", "yw0 = w0.apply_warp(y)\n", "yw1 = w1.apply_warp(y)\n", "yw01 = w1.apply_warp(yw0)\n", "\n", "#plot:\n", "plt.figure()\n", "ax = plt.axes()\n", "ax.plot(yw0, 'c', label='Only first warp')\n", "ax.plot(yw1, 'b', label='Only second warp')\n", "# ax.plot(yw01, 'orange', label='Sequentially warped')\n", "ax.legend()\n", "ax.set_xlabel('Domain position (%)', size=13)\n", "ax.set_ylabel('Dependent variable value', size=13)\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The main value of the **SequentialManualWarp** class is to apply manual warps (that were constructed in the GUI) to other data. For example, if the original data represent the time series of a three-component force vector, and if the warps were created using only the X component time series, then the **SequentialManualWarp** class can be used to apply those warps to the Y and Z components. This is demonstrated in the \"Applying GUI warps\" example." ] } ], "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 }