From 376230d6cff8997d1bb1cb9d359a572a02b74922 Mon Sep 17 00:00:00 2001 From: Luke Naylor <l.naylor@sms.ed.ac.uk> Date: Mon, 29 May 2023 00:51:56 +0100 Subject: [PATCH] Jupyter notebook plot tests --- plot_playground.ipynb | 94 +++++++++++++++++++++++++++++++++++++------ 1 file changed, 82 insertions(+), 12 deletions(-) diff --git a/plot_playground.ipynb b/plot_playground.ipynb index 36aa023..17b76f8 100644 --- a/plot_playground.ipynb +++ b/plot_playground.ipynb @@ -2,7 +2,7 @@ "cells": [ { "cell_type": "code", - "execution_count": 2, + "execution_count": 1, "id": "05f8d418", "metadata": {}, "outputs": [], @@ -13,7 +13,7 @@ }, { "cell_type": "code", - "execution_count": 10, + "execution_count": 2, "id": "b46f8a6a", "metadata": {}, "outputs": [ @@ -26,10 +26,10 @@ "$\\displaystyle \\text{Chern Character:} \\\\ \\begin{array}{l} \\mathrm{ch}_{0} = 3 \\\\ \\mathrm{ch}_{1} = 2 \\ell^{1} \\\\ \\mathrm{ch}_{2} = -2 \\ell^{2} \\end{array}$" ], "text/plain": [ - "<pseudowalls.chern_character.Chern_Char object at 0x7f0034ca7310>" + "<pseudowalls.chern_character.Chern_Char object at 0x7f1e3f2d1150>" ] }, - "execution_count": 10, + "execution_count": 2, "metadata": {}, "output_type": "execute_result" } @@ -41,7 +41,7 @@ }, { "cell_type": "code", - "execution_count": 15, + "execution_count": 3, "id": "8ed6de82", "metadata": {}, "outputs": [ @@ -57,7 +57,7 @@ "-2/3" ] }, - "execution_count": 15, + "execution_count": 3, "metadata": {}, "output_type": "execute_result" } @@ -69,7 +69,7 @@ }, { "cell_type": "code", - "execution_count": 12, + "execution_count": 4, "id": "edf2389e", "metadata": {}, "outputs": [ @@ -85,7 +85,7 @@ "16" ] }, - "execution_count": 12, + "execution_count": 4, "metadata": {}, "output_type": "execute_result" } @@ -96,7 +96,7 @@ }, { "cell_type": "code", - "execution_count": 16, + "execution_count": 5, "id": "1dcae94c", "metadata": {}, "outputs": [ @@ -109,10 +109,10 @@ "$\\displaystyle \\text{ Twisted Chern Character for $\\beta={ -\\frac{2}{3} }$ } \\\\ \\begin{array}{l} \\mathrm{ch}_{0} = 3 \\\\ \\mathrm{ch}_{1} = 4 \\ell^{1} \\\\ \\mathrm{ch}_{2} = 0 \\ell^{2} \\end{array}$" ], "text/plain": [ - "<pseudowalls.chern_character.Twisted_Chern_Char object at 0x7f002672b710>" + "<pseudowalls.chern_character.Twisted_Chern_Char object at 0x7f1e34aa1ed0>" ] }, - "execution_count": 16, + "execution_count": 5, "metadata": {}, "output_type": "execute_result" } @@ -123,9 +123,79 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 7, "id": "aa77740a", "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "<html>\\(\\displaystyle \\text{Chern Character:} \\\\ \\begin{array}{l} \\mathrm{ch}_{0} = 1 \\\\ \\mathrm{ch}_{1} = 0 \\ell^{1} \\\\ \\mathrm{ch}_{2} = 0 \\ell^{2} \\end{array}\\)</html>" + ], + "text/latex": [ + "$\\displaystyle \\text{Chern Character:} \\\\ \\begin{array}{l} \\mathrm{ch}_{0} = 1 \\\\ \\mathrm{ch}_{1} = 0 \\ell^{1} \\\\ \\mathrm{ch}_{2} = 0 \\ell^{2} \\end{array}$" + ], + "text/plain": [ + "<pseudowalls.chern_character.Chern_Char object at 0x7f1e34534cd0>" + ] + }, + "execution_count": 7, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "u = Chern_Char(1,0,0)\n", + "u" + ] + }, + { + "cell_type": "code", + "execution_count": 37, + "id": "c711b1ce", + "metadata": {}, + "outputs": [ + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnQAAAFbCAYAAABcexOxAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAAB2fklEQVR4nO3dd1hT1xsH8G/YioIDRBBR3OBEXFj3QNFqa1t3XXWPuqsI7lG1zqo4WhW1zv5aa63S1gmolYqIExURlC2CylAZkvP745gAZUjIuDfJ+3mePDcJd7zcG8PrueecV8IYYyCEEEIIIVrLQOgACCGEEEKIciihI4QQQgjRcpTQEUIIIYRoOUroCCGEEEK0HCV0hBBCCCFajhI6QgghhBAtRwkdIYQQQoiWMxI6ADGSSqWIj49HxYoVIZFIhA6HEEIIIXqIMYb09HTY2dnBwKDkNjhK6IoQHx+PmjVrCh0GIYQQQghiYmJgb29f4jqU0BWhYsWKAPgJtLCwEDgaPbJ6NbBmTd7zKVPUfsiTJ4ERI/jzI0eAPn3UfkhCiCrcvAl07gwEBAAtWggdDVEBuqSFpaWloWbNmvK8pCSU0BVBdpvVwsKCEjpN+vZbwMgIWLkSWLAAMDcHvv5arYf88ksgJATYvBmYPBm4cQNwdFTrIQkhqlChQt6Svqd1Al3S4pWm+xcNiiDiIZEAy5fzZA4Apk8HtmxR+2HXrgXatQNevQIGDways9V+SEIIIUSlKKEj4iKRAKtWAZ6e/PWMGcCmTWo9pIkJcOwYULkyEBycl08SQkTM1hZYsoQviU6gS6ocCWOMCR2EWPj4+MDHxwe5ubkIDw9Hamoq3XIVCmPAokU8uQOAdeuAuXPVesiTJ4FPPsl73q+fWg9HCCGElCgtLQ2WlpalykcooSvCh05gdnY23r17J0BkeoYxYNcu4Icf+Ovp04HRo9V2OCMjI8yfb4LNm3lr3c2bgIOD2g5HCFFGWhpw9Srg5kYdrnQEXdLCKKFTUkknMDs7G/fu3YNUKhUoOqIuBgYGqF+/Mbp2NcH160D79oC/P2BsLHRkhJBCbtwAXF35qKaWLYWOhqgAXdLCFEnoaJSrgt69ewepVIratWujXLlyQodDVOTt27d48uQJDAze4dgxE7i4AP/8w+/6ymZSIYQQQsSKEroyKleuHMqXLy90GETFUlNTUalSJrZvN8WXX5pj7Vqge3egZ0+hIyOEEEKKR6NcCcnn1KlTOHDgALKyjmHChNcAgJEjgaQkgQMjhBBCSkAJHSH5mJubw8TEBGlpaViyJAuNGwOJiXwsBnWbJERETE2BunX5kugEuqTKoYSOkHxMTU1h+v7bxMwMOHqUL//8U+3T4RFCFNG4MRARwZdEJ9AlVQ4ldISUoEmTvETO0xMIChI2HkIIIaQolNDl4+PjA2dnZ7Ru3VroUIiITJzIS4K9ewcMGgSkpAgdESEEt28D1tZ8SXQCXVLlUEKXz9SpUxEWFobg4GChQyEiIpHwuY3r1wdiYoCxY/mcx4QQAb17ByQn8yXRCXRJlUMJHSGlYGHB672amAC//w74+AgdESGEEJKHEjpCSsnFhZeUBYA5c3hpMEIIIUQMKKETifv37wsdgmCSkpJw8OBBHD58GHPmzCmxTm5wcDDevHmjwegK+vproF8/IDub96vLyBAsFEIIIURO9AldYGAg+vXrBzs7O0gkEpw4caLU2165cgVGRkZo0aKF2uJThc2bN8unytBHFy5cwKtXrzBs2DDk5OTg3Llzxa67d+/eD5Zce/78OdavX6/qMAHw/nS+vkCNGkB4ODBtmloOQwj5kAYNeH2+Bg2EjoSoCF1S5Yg+oXv9+jWaN2+Obdu2KbRdamoqRo4cie7du6spMtU4c+YMrK2tUadOHQCAv78/jh07hr1792LQoEG4cePGB/fx77//YsOGDVi6dCl69eqFgICAQus8ffoUO3fuxJ07d1T+OxRn/vz5uHTp0gdjGTJkCKa9z4wSExPRqFGjIvcXHx8vT+yB4n9va2trdOrUCT5q6uhWtSpw+DBgYADs3w/89JNaDkMIKUmFCoCbG18SnUCXVDmiT+g8PDywcuVKfPbZZwptN3HiRAwbNgxubm5qikw1tmzZgiFDhshff/7558jMzMRXX32Ftm3bYtCgQSVu/+bNG5w4cQJz5szB0qVLMX78eHh4eCAuLk6+Tk5ODry8vDB+/HjcunVLbb9LfpcuXcK+ffuQm5tb4P2SYjl9+jQGDhyI2rVrF7nPQ4cOYfjw4QA+/Hu3adMG165dQ3Jysup/OQCdOgFLlvDnU6YAjx+r5TCEkOLExgKzZ/Ml0Ql0SZUj+oSuLHx9ffH48WMskf3F/YCsrCykpaUVeCgqIQHIzOTPo6KAGzcKPhISCm9z+/Zt2Nvbw9DQUP6ev78/vvjiC/nrjA900oqIiMCaNWsQEREBAOjVqxfevn2LK1euyNdZv349PvroIyQlJcHW1lbh301RaWlpuHXrFpycnAr9rLhYgoODYWNjg4EDByIsLKzI/UZGRspbMkvze/fv3x/79u1T4W9WkLc30LEj70c3fDiQk6O2QxFC/ispic/6TYWWdQZdUuXoXEL36NEjeHp64tChQzAyMirVNqtXr4alpaX8UbNmTYWPu2tX3v8qFi0CXF0LPnbtKrzN+fPnC01i3LRpU5ibmwMA/vjjD6xdu7bE4zZt2hRXrlxB3bp1AQDR0dEAgPr16wPgyerWrVsxaNAgnD17Ft26dVP4d1PUDz/8gPHjxxd6v7hYrl69irFjx8LT0xOdOnUCK2KSt+vXr6NVq1by1x/6vQGgc+fO8PPzU+nvlp+hIXDwIGBpCfz7L7BsmdoORQghhJSodBmPlsjNzcWwYcOwbNkyNFCgV+WCBQswe/Zs+eu0tDSFk7qJE4HKlfnzFSt4s3F+RTWMxcbGokuXLoXev3btGk6dOoXWrVtj8ODBJR5XIpGgffv28tdr1qzBzJkz4eLiAgD4888/0bRpU6SmpsLJyUne/0xdTp06hd69exc5yKO4WNzc3HD7A1OD//zzz1i4cKH89Yd+bwCwsrLCYzXfC3Vw4JMODx4MfPst0LMn0LmzWg9JCCGEFKJTCV16ejquX7+O0NBQeSd7qVQKxhiMjIxw5syZIluo8hdkL6v8CZujI398SEZGRpEjNtu0aYM2bdpg586d6NixIy5evIgKpeglunfvXlSvXh3fffed/L2LFy+iVatWiIqKQo8ePUr1u8h8//33iIqKKvbnnTt3xoABA+Sv4+Pj8eLFC3z88cdFrl/WWHJycpCVlQULC4sif17U7y2j7gQW4OXA/vqLj3798ktetkaW3BNCCCGaoFMJnYWFRaFRnNu3b8eFCxfwyy+/wLE0WZYGWVlZ4eXLl/LXQUFB+OSTTxAUFARHR0d06dIFkydPxl9//VWgX11R/Pz8kJubi3Xr1iEzMxOJiYmoXbs2IiIi0KdPH4WTOQCYMWOGQuv//fffSExMxJo1awAA4eHhOHz4MN6+fQsPD48yx3Lq1Cn07du3yJ8V93vLlPa2u7K2bAEuXwYePQImTAB+/plPcUIIURMrKz4iycpK6EiIitAlVY7oE7qMjAx5x3cAiIqKws2bN1GlShU4ODhgwYIFiIuLw4EDB2BgYIAmTZoU2L5atWowMzMr9L4YODk5yft+ATz5qF+/PmrUqAGADwIwNjZG8+bNAQABAQEoX758oX53gYGBiI+PR79+/ZCYmIigoCDY2Njg9evXiIiIKPK2bXH7UsaYMWMKvN65cyeGDRuGLl264N69e8XG8iFnz54tctqa4n5vWUKXk5ODihUrlul3UVSFCnwqEzc34Jdf+FQmI0dq5NCE6CcHB6rBp2PokipH9And9evX0bVrV/lrWV+3UaNGYd++fUhISCiQFGmT3r17Y/z48ZgzZw4AoFWrVpg9eza2bNkCQ0NDXL58GadOnZJ39P/hhx+Qnp6OkydPyvcRGRmJjz/+GOnp6QX2nZqaitDQUIwfPx7JycmwtLTEnj17MGHCBBgYGBS5L1WJjY3Fli1bkJiYiA0bNuDNmzcwNzcvNpaSJCcno2rVqoXWK+n3lgkODtbIIBCZVq34wAhvbz7hcMeOpbv1TggpgzdvgAcPgEaNgPLlhY6GqABdUuVIWFFDCvVcWloaLC0tkZqaWqjf1ps3b3D//n04OTmhvAo+cRMmTMDy5ctRvXr1Uq3/008/YcSIEaXef2pqKnx9fVGuXDkMHTq0wO+j6L6UVVIsxdmyZQt69OgBZ2dnhY+3cOFC9O/fH23atPngurLr+uTJE7x69Qrp6ekYOXIkqlSpotAxc3P5oIgrV4B27YDAQMDYWOHQCSEfcuMGn0IgJARo2VLoaIgK0CUtrKR85L90btoSbbN48eJSV8F48eJFkVN6lMTS0hIzZ87ExIkTC3wYyrIvZRUXS0nu3btXpmQuLS0Nz58/L1Uyp0qGhsChQ3wqk6AgmsqEEEKIZlBCJzB7e3sMGDCgVPOl3b59u0BVCWWocl/qEhkZWWAaktJijGHTpk1YuXKl4gctRam1D6lVC/jxR/7822+BixeV3iUhhBBSIkro8vHx8YGzs7NKBwqUhqurK/r06fPB9bp06QITExOVHFOV+1KXOnXqYNKkSQpvl5KSgsmTJ8Pa2lrxg16/zpvWlDRwIDBuHMAYn8okJUXpXRJCCCHFooQun6lTpyIsLAzBwcFCh0KUYGVlhWrVqpV9B1eu8MnklLR5M9CwIRAfz6cyod6qhKiQgQFQsSJfEp1Al1Q5dNoIyU/W5y40FNi+XaldmZvzqUyMjIDjx/nEw4QQFWnRAkhL40uiE+iSKocSOkLya9EC6NCBP1+5EtiwQandtWzJS8EBwPTpQL4pFQkhhBCVoYSOkP9q2xZ4P5kz5s4FNm1SanfffAN06gS8fs370+XkqCBGQvRdWBjQuDFfEp1Al1Q5lNCJxP3794UOgeTXogXwfhJrzJ7Na3uVkaEhrxxhaQn8+y9v+COEKCkzk//lz8wUOhKiInRJlUMJnQhs3rwZpqamQocBAHj+/DnWr18vdBjiMGcOL/sAADNmAFu3lnlXDg7Ajh38+cqVwNWrKoiPEEIIeY8SOoGdOXMG1tbWqFOnjvy9jIwMfP7556UuaRYcHIyZM2di//79mDBhAqKiosq8jbW1NTp16gQfKqgHSCS8A5ynJ389fbpSt1+HDgWGDwekUl7n9fVrFcVJCCFE71FCl48Q89Bt2bKlwAS/e/fuxbp163D8+HFIpdIPbp+VlYWBAwdiwYIFGDVqFMaMGYMxY8YotU2bNm1w7do1JCcnl/0X0xUSCZ8dWNZSN3s2oEQL5rZtQI0afHDEN9+oKEZCCCF6jxK6fJSahy4hIe/Gf1QUrziQ/5GQUGiT27dvw97eHoaGhvL3vvrqKyxToF5UYGAgLCwsYGNjA4AnY0FBQXj27JlS2/Tv3x/79u0rdRw6TdZSt3gxf/3NN8DatWXaVaVKgOy07tgB/P23SiIkRP/UqQP8/jtfEp1Al1Q5lNCpyq5dQGwsf75oEa8wnP+xa1ehTc6fP690a+CTJ08KFJA3NDRExYoVce/ePaW26dy5c6nKkekNiYQXZl26lL/29OQtd2XQowfw9df8+ZgxwIsXqgmREL1SqRLQvz9fEp1Al1Q5lNCpysSJgL09f75iBRASUvAxcWKhTWJjY5WraAAgOTkZZmZmBd4zMzPDq1evlNrGysoKjx8/Vio2nbRkSd7Ect7ewPLlZdrNmjW8ikRCAjB1qgrjI0RfJCYCq1fzJdEJdEmVQwmdqtjaArIkydGRzyib/2FrW2iTjIwMlCtXTqnDWlpagv2nplRGRgasrKyU3kYikSgVm85auDCvdW7JkjIldeXLAwcO8ClNjh4Fjh1TcYyE6Lr4eMDLiy+JTqBLqhxK6ARkZWWFly9fKrWPRo0aISkpSf46KysL6enpqFWrltLbGBkZKRWbTluwIK8f3ZIlZbr92qYN//ICeCsd/a+UEEJIWVFCJyAnJ6dST00CAAEBAYUGbHTq1AlJSUmIfd9/LyAgAK1bt5YnZ2XZBgBycnJQsWLFMv1eemPePH5/AOC3X8swUGLhQj6HcUoK8NVXwH8aTgkhhJBSoYROQL1798bly5cLvHfkyBFMnjwZAODp6Ynt+QrE//DDD1gh67/1npGREfbv349Vq1Zh//79OHr0KA4dOqTUNgCfp65bt24q+T11mqdnXukHT0+FpzQxMQEOHgRMTYE//yxy7AwhhBDyQXRPTUBWVlawtrZGYmIiqlevDgAYOnQohg4dih2ysgL5HDp0CD/99FOh93v06IEePXoAAEaNGqX0NgDg5+eHwYMHK/5L6SNvbyA3l996/eYb3jFu1qxSb964MR8kMWsWL07RrRvQoIEa4yVEF1SqBHzxBQ2J1CF0SZVDLXT5CDGx8OLFi7Ft27ZSrfvixYtCgxnUsU1aWhqeP3+ONm3aKLSdXlu8mE9XA5Sp9uv06UD37sCbN8CIEUBOjhpiJESX1KkD/O9/NGmZDqFLqhxK6PJRamLhMrK3t8eAAQNKNefb7du3C1SVKA1Ft2GMYdOmTVhJFeQVt2xZ3iiHGTMUSuoMDPiEw5UqAdeu5XXNI4QUIzubz/2ZnS10JERF6JIqhxI6EXB1dUWfPn0+uF6XLl1gYmKi0L4V3SYlJQWTJ0+GtbW1Qsch4JMPr1yZV/t1xgyFar/a2wOyErorVgChoWqIkRBdcfcuULMmXxKdQJdUOZTQkQKsrKyUnuxYrxVV+3XdulJvPnQo8PnnwLt3wKhRQFaWmuIkhBCiUyihI0TV/lv7dd48PuqhlJvu2AFYWwN37vC7uIQQQsiHiD6hCwwMRL9+/WBnZweJRIITJ06UuP7x48fRs2dPWFtbw8LCAm5ubvibKqATTZPVfpVlZAsW5E1v8gHW1sDOnfz52rXAv/+qKUZCCCE6Q/QJ3evXr9G8efNSjwQNDAxEz5494efnh5CQEHTt2hX9+vVDKHVIIkJYvBhYtYo/X7So1E1un30GDBsGSKX81uvbt2qMkRBCiNYT/Tx0Hh4e8PDwKPX6mzdvLvD622+/xe+//44//vgDLi4uKo6OkFLw8uJz03l6AkuX8tY72e3YEmzdCly8CDx8yCtKbNig/lAJ0RotWgCZmYCxsdCREBWhS6oc0bfQKUsqlSI9PR1VqlQpdp2srCykpaUVeBCiUvPn5w2OWLIkr9WuBFWqAD/8wJ9v3gwEBakvPEK0joEBL7FioPN/xvQGXVLl6Pxp27BhA16/fo1BgwYVu87q1athaWkpf9SsWVODERK9MXdu3uCIhQtLVfv144+BL7/kt16/+or/75UQAiA8HOjShS+JTqBLqhzR33JVxpEjR7B06VL8/vvvJU7FsWDBAsyePVv+Oi0t7YNJ3Vvq1KRTNHY958/nZcK8vfktWCMjXu+rBN9/D5w9C9y/DyxfzmdFIUTvZWQAAQF8SXQCXVLl6GxCd+zYMYwdOxb/+9//5DVLi2NqagpTU9NS7dfIyAgGBgZ48uSJCqIkYsIYQ25urvoP5OXFJ5pbsoS32hkaAjNnFrt6lSp8KpPPPgO++47PU+fqqv4wCSGEaA+dTOiOHDmCr776CkeOHEHfvn1Vum8TExM0btwY7969Q2pqKk6dOgVzc/NSJ4REvHJzc/Hu3TvNHGzxYp7UrVgBzJrFk7qvvy529QEDgMGDgWPHgDFjgOvXAQWLhhBCCNFhok/oMjIyEBERIX8dFRWFmzdvokqVKnBwcMCCBQsQFxeHAwcOAODJ3MiRI/H999+jXbt2SExMBACUK1cOlpaWKonJxMQEJiYmyMzMxOvXr5GTk0MJnQ7J0lR5hmXL+O3Xb78Fpk/nHeVmzCh29a1bgfPn+YTDa9aUaqAsIYQQPSH6hO769evo2rWr/LWsr9uoUaOwb98+JCQkIDo6Wv7zXbt24d27d5g6dSqmTp0qf1+2viqZmprCwsICaWlpyKZqwjrFwsJC/Um6rParVMoztJkzeYKXrz9nftbWwJYtfH66lSt5q13TpuoNkRDRcnAAfvyRL4lOoEuqHAljjAkdhNikpaXB0tISqampsLCwKHHd169fa65Fh2iMqakpzM3NNXMwxvikw7KpTNau5eXCiln100+BkyeBli35VCY0ZxMhhOgmRfIR0bfQiZ25ubnm/vAT3SSr/WpkxG/Dzp/P+9d5eRW56s6dwKVLwI0bfJCEt7cAMRMitORk4MQJ/j8cKyuhoyEqQJdUOTo/Dx0hWkEi4VUkli/nr729eZJXBFtbfusV4Pnf3buaCZEQUYmOBsaP50uiE+iSKocSunx8fHzg7OyM1q1bCx0K0VeLFuVNNLd4MU/yijB8ONCvH5CTA4wezRv0CCGE6C9K6PKZOnUqwsLCEBwcLHQoRJ8tWMDvpQK8CU7WapeP7NZrpUpASAiwcaNmQySEECIulNARIkbffAOsX8+fF1P71c4O2LQpb5VHjzQYHyGEEFGhhI4QsZoz54O1X0eNAnr25DVex4/nM6AQohcqVAA6d+ZLohPokiqHpi0pgiLDhAlRu2+/zRvKun59odqvUVFAkybAmzfArl3AhAkCxEgIIUTlFMlHqIWOELHz8uJ96QBe+3Xz5gI/dnTMuyP7zTdAXJxmwyNEEFIpkJVFzdI6hC6pciihI0QbLF7MR8ACvPbr1q0Ffvz110DbtkBaGjBlCp+AmBCddvMmYGbGl0Qn0CVVDiV0hGiLZcv4CFiA137N11JnaAjs3s2rRpw8Cfz2mzAhEkIIEQYldIRoC4mE31uVJXWzZgEbNsh/3KQJLzIB8Ba71FQBYiSEECIISujyoYmFiejJkjrZ7de5cwuMfvX2BurVA+Lj+cBYQggh+oESunxoYmGiFSQSPtmwbKCEpyewciUA3v9k507+to8P8O+/AsVICCFEoyihI0RbLV4sT+SwaBEvE8YYuncHRozgAyMmTODlwQjROU2aADExfEl0Al1S5VBCR4g28/bOm3x42TJeMgK8a12VKsDt24VmOSFEN5iYAPb2fEl0Al1S5VBCR4i2mz8/r0zYihXA8uWwts57a+lS/r9eQnRKZCQwcCBfEp1Al1Q5lNARogvmzClU+3X0aKBDB15BYvp0QaMjRPVevQJ++YUviU6gS1qEZ89KvSoldIToiv/UfpV8txY7dgBGRsCJEzQ3HSGEaJVnz4CPPy716pTQEaJL5s/PGyjh6Ykmf2/AvHn85bRpvJIEIYQQkXv+HOjeHQgPL/UmlNDlQ/PQEZ3g7c07zgHA3LlYWGmbfG46Ly9BIyOEEPIhyck8mbt3D6hevdSbSRijqo//lZaWBktLS6SmpsLCwkLocAgpm8WL+SAJABemHUf3bQMgkQD//AO0aydwbIQoKzER8PUFxoxR6I8eES+6pABSUoAePXhB2+rVkXb6NCxdXUuVj1BCVwRK6IhOYIy31q1eDQAY3fou9gc3RosWQHAw71tHCCFEJGS3We/cAapVAwICkGZnV+p8hG65EqKrZGXC3t9nXRfcBZXLvcXNm7yKBCFa7dUr4ORJGhKpQ/T6kj57BnTtypO56tUBf3+gUSOFdkEJHSG6TCLhgyQWL4Y1krHm7QwAvLBEfLzAsRGijMhI4JNPaNIyHaK3lzQhAejShfeZs7PjyZyTk8K7oYSOEF0nkfAqEsuXYxx2oy2CkJ4OzJ4tdGCEEKLn4uJ4MvfgAS+TERAANGxYpl2JPqELDAxEv379YGdnB4lEghMnTnxwm4CAALi6usLMzAx16tTBTlm1ckL02aJFMPh2FXZgMgyQi2PHgLNnqAstIYQIIjaWJ3Ph4YCDA0/m6tUr8+5En9C9fv0azZs3x7Zt20q1flRUFPr06YOOHTsiNDQUXl5emD59On799Vc1R0qIFliwAC7fDcPX2AoAmDLkBbKyBI6JEEL0TWws0LkzEBEB1K7Nk7k6dZTapejHuXl4eMDDw6PU6+/cuRMODg7Y/L4iuZOTE65fv47169fj888/L3KbrKwsZOX7q5ZGs68SLcAYQ3p6OlJSUkp8vHz5Erm5uZBKpfLHu9pxMH5yBBEvjVDf5iVsG1YEYwxSqRQGBgawsLBApUqVUKlSJVhaWhb5XPba2toa5cuXF/p0EH1jZgY4O/Ml0Ql6c0nj4/kAiMhInsRdvMhb6JQk+oROUVevXoW7u3uB93r16oU9e/YgJycHxsbGhbZZvXo1li1bpqkQCSkVqVSKhIQEPHr0CBEREfJlREQEkpKSkJKSgpycHKWPE5MKxFxTbh/VqlWDo6MjHB0dUbt2bflzR0dHODg4wMTEROk4CSnA2Zl3Iic6Qy8uaUICT+ZkLXMqSuYAHUzoEhMTYWNjU+A9GxsbvHv3DsnJybC1tS20zYIFCzA7Xw/xtLQ01KxZU+2xEgIAz549w/3794tM3N6+ffvB7c3MzFC1atUiH1ZWVqhcuTKMjY1hYGAAiUQCAwMDGBgYAJBg0ZQXePC8GjriEuaNSobBoEHIzc1FamoqUlNT8erVK/ky//P872VlZSEpKQlJSUn4999/C8VnYGCAGjVqyBM8JycntGzZEi1btkTVqlXVcEYJIUSEnj3LK+fl4KDSZA7QwYQOACQSSYHXsrmT//u+jKmpKUxNTdUeFyHp6ekICQnBtWvXcO3aNQQHByM6OrrY9Q0NDeHo6Ih69eqhfv36qFevHurWrQs7OztYWVmhatWqSt3urFcPcG0pxSX2CZbv74IuzR4oPPz11atXiIqKKvLx5MkTvH37FjExMYiJiUFgYGCBbWvVqiVP7lxdXdGyZctC/yEjpEg3bwKdOgGBgUCLFkJHQ1RApy+pbNLg+/f5aNYLF3gLnQrpXEJXvXp1JCYmFngvKSkJRkZG1BpANCo7Oxt37twpkLyFhYXhv8VZJBIJ6tatK0/Y8i9r1apVZDcBVWnRApg4yQA7dgBfYytC57jAyMgImD691PuoVKkSXFxc4OLiUuhnjDEkJSXJE7zIyEjcuXMHISEhiIiIwNOnT/H06VP89ttv8m3s7OwKJHgfffQR/dslhUmlQHo6XxKdoLOXVFbOSzbP3IULQN26Kj+MziV0bm5u+OOPPwq8d+bMGbRq1UqtfxgJeffuHYKCguDn54fAwECEhIQgMzOz0HoODg5o06YN2rRpg9atW8PV1RUVK1YUIGJu5Urg558Z7qY0xQ5Mxtcz+OTDiiR1xZFIJLCxsYGNjQ3a/aeAbGpqKm7evImQkBDcuHEDN27cwIMHDxAfH4/4+HicOnVKvg9XV1e4u7vD3d0dbm5u1CePEKIdkpOBnj2B27d5BYgLF4D69dVyKNHXcs3IyEBERAQAwMXFBRs3bkTXrl1RpUoVODg4YMGCBYiLi8OBAwcA8GlLmjRpgokTJ2L8+PG4evUqJk2ahCNHjhQ7yvW/qJYrKa1nz57hr7/+gp+fH86cOYNX/6lZU7lyZbRu3bpAAlddhFWnd+4EJk8GKpu9QUSmPargJbBhg8ZnH87IyMCtW7fkCd61a9cQFhZWYJ0KFSqga9eu8gSvfv36xXanIDrsxg3A1RUICQFathQ6GqICOndJ/1ub9eJFPvJDAQrlI0zkLl68yAAUeowaNYoxxtioUaNY586dC2zj7+/PXFxcmImJCatduzbbsWOHQsdMTU1lAFhqaqqKfguiK969e8euXr3KFi1axFxdXQt9LitXrsyGDh3KfH192cOHD5lUKhU65FLJyWGsaVPGAMa+bv0PfwIwtmaN0KGxuLg4tn//fjZ8+HBmbW1d6JzXqlWLjR8/nv3vf/9jL168EDpcoikhIfwzGhIidCRERXTqkiYmMta4Mf+FbG0Zu3+/TLtRJB8RfQudEKiFjuT36tUrnD59Gn5+fvj777+RkpJS4OctW7ZEnz590KdPH7Rp0waGhoYCRaqc8+d5Nw9DQ+D2RB84b5/Gf7ByJeDtLWxw70mlUty6dQtnzpzBmTNncPnyZWRnZ8t/bmBggG7dumHIkCH47LPPULlyZQGjJWr15g0vl9SoEUDzIOoEnbmkCQlAt278l6lRg7fMlfE2qyL5CCV0RaCEjmRlZcHPzw+HDh3CqVOnCkw8bWlpiV69esHDwwO9e/cW5S3UsvrkE+DkSX6X4GznlZAsXsR/sHQpsHgxrwsrIq9fv0ZgYKA8wct/e9bY2Bi9evXCkCFD0L9/f0H7KRJC9ERsLE/mHj0CatbkyZwSAyAooSsjHx8f+Pj4IDc3F+Hh4ZTQ6RmpVIorV67g4MGD+Pnnnwv0h3N2dsYnn3wCDw8PuLm58ZGgOigyknfxyMoC/vc/4IvHawFPT/7DhQuB5ctFl9TlFxkZiWPHjuHo0aO4ffu2/H0zMzN8/PHHGDJkCPr27QsznZ+KXg9ERwNr1wLz56t0Li8iHK2/pNHReRUgatXiyZyjo1K7pIROSdRCp18iIyNx4MABHDhwAFFRUfL37ezsMGzYMHz55Zdo1qyZ3nS8X7KE5201a/I7BuV3bADmzuU/9PLit2C14FyEhYXh6NGjOHr0KB49eiR/38LCAp9//jmGDx+OLl26aO0tcr2ncz3oiVZf0pgYoEsXnsw5OvJkrlYtpXerU4MihECDInRfamoq2717N+vYsWOBDvYVK1Zko0ePZufPn2fv3r0TOkxBvHnDmIMD78u7bNn7NzdvzhsosXSpoPEpSiqVspCQEPbNN9+wmjVrFrjeNWrUYCtWrGDPnj0TOkyiKJ3qQU8Y0+JLGhvLWL16PPg6dRiLjlbZrhXJRyihKwIldLrr/v37bOLEiaxcuXLyP+oSiYT17NmTHTx4kL1+/VroEEXh2DH+3VS+PGMxMe/f3LAhL6lbuVLQ+MoqNzeXBQQEsAkTJrBKlSrJPwMmJiZs1KhR7Pr160KHSEpLa//6k+Jo5SWNj2esQQMeeO3ajD19qtLdK5KPGCjdHkiIyDHGcO7cOfTt2xdOTk7YtWsX3r59CycnJ6xZswYxMTE4c+YMhg8frlQZLV0ycCDQsSMfdSbrQofZs4E1a/jzhQuB774TLL6yMjAwQKdOnbBr1y4kJibi4MGDaNOmDbKzs7F//360atUKHTp0wM8//4ycnByhwyWEiJmaa7MqTKWppI6gFjrdkJmZyfbu3cuaNm1aoDXu008/ZQEBAVozR5xQQkIYk0j4fzyvXs33g5Ur81rqNmwQLD5VCgoKYsOHD2fGxsYFbseuXLmSJSUlCR0eKUpMDGOzZuVrQibaTqsuaVISY02a8O9Be3vGIiLUchiah05JNChCuyUnJ2PHjh3w8fHBs2fPAADm5uYYM2YMZsyYgXr16gkcofYYOxbYuxdo0wa4ehUwkLXpL10KLFvGn3//vUrKhIlBQkICdu3ahZ07d8o/O6amphg6dChmzJiBFjpXMZwQorCUFD41ye3bgK0tEBCgtnJeNChCSdRCp50iIyPZtGnTCvSPs7e3Z2vXrqUKAmWUkMBYxYr8P6E//ZTvB1IpY15eeS11338vWIzqkJmZyQ4ePMhat25dYBBF9+7d2ZUrV4QOjzDGWHo6Y//8w5dEJ2jFJX3+nLHmzfn3no0NYw8eqPVwNChCSZTQaZeHDx+yoUOHMgMDA/kf3pYtW7LDhw+z7OxsocPTeqtX8+8uBwc+AlZOKmVs/vy8pG79esFiVKegoCA2dOhQZmhoKP989e3bl4WGhgodmn7Tyh70pCSiv6TPnuXVSLSxYezePbUfkgZFlJGPjw+cnZ3RunVroUMhpRAfH4+JEyfC2dkZR44cgVQqRa9evXD+/Hlcv34dQ4cOhbGxsdBhar0ZM/icdNHRwJYt+X4gkQCrV/MBEgCfq042aEKHtG3bFocPH0ZkZCTGjRsHQ0NDnD59Gi4uLhg0aBAePHggdIiEEHVLTOSTBt+5w2+z+vvzWdjFRO3ppRaiFjpxe/nyJfP09Cxwa/Xjjz+mFhM1OnCA/6fUwoL3BS5k2bK8lrrlyzUenyaFh4ezoUOHMolEwgAwAwMDNnr0aBYVFSV0aPpF9M05RFGivaRxcYw1bMiDq1GDsfBwjR2aWuiITnr79i3WrVuHOnXqYM2aNXj79i3at2+PS5cu4Y8//qAO62o0fDjg4gKkpfEqEoUsXgysWpX3fMkSnt7poPr16+Pw4cO4desWPvnkE0ilUuzbtw8NGjTAtGnTkJCQIHSIhBBViY0FOncGHj7kU5KocQCEsiihI6KXm5sLX19fNGjQAPPmzcPLly/h7OyM33//HZcvX0aHDh2EDlHnGRgA69fz5zt38mmXCvHyypubbvlyfitWR5M6AGjatClOnDiBoKAg9OjRAzk5OfDx8UHdunUxb948vHjxQugQdZuREWBlxZdEJ4jukkZH82QuIgKoXZsnc3XrCh1V8TTQYqh16JareJw5c4Y1a9ZMfmu1Zs2azNfXV2/Lcgmtb19+1+Hzz0tYaePGvNuv3t588IQeuHDhAnNzc5N/VitVqsTWr19PA3MI0UbR0Yw5OvLvsbp1VV4BorTolivRevfu3YOHhwfc3d1x+/ZtVKpUCevWrUN4eDhGjx5NBdUFsnYtHwvx669AcHAxK82axeemA/ht2CLv0eqerl274sqVKzh16hSaNm2KV69eYe7cuWjVqhWuXr0qdHiEkNKKi+PzzEVF8RY5f39hK0CUEiV0RFTevn0LLy8vtGjRAn/99ReMjIwwY8YMREREYO7cuTAzMxM6RL3WuDEwciR/7ulZwh3V6dOBDRv486VLgW+/1UR4gpNIJOjbty9CQ0Oxe/duVK1aFbdv38ZHH32ESZMm4eXLl0KHqDvu3QPq1eNLohNEcUkTE3kyJ7vNeuECYG8vYEClRwkdEY0LFy6gWbNmWL16Nd69e4f+/fsjLCwMmzdvRtWqVYUOj7y3bBlgYsK/586dK2HF/LVfvb2Bdes0Ep8YGBoaYuzYsXjw4AFGjx4Nxhh27dqFRo0a4dChQ2A63LdQY7KygMeP+ZLoBMEvaVIST+bEUptVQZTQEcGlpKTgq6++Qvfu3REREQE7OzscP34cv//+O+qLdDSRPqtVC5gyhT9fsACQSktYef58YMUK/nzePGDTJrXHJyZWVlbw9fWFv78/GjVqhKSkJHz55Zdwd3fHo0ePhA6PECKTnAx07w7cv89b5C5e5C10WoQSunxoYmHNYozh2LFjcHJygq+vLwBg8uTJCAsLw4ABAwSOjpTEywuoUAEICQF++eUDKy9cyKcyAXir3datao9PbDp37oybN29i5cqVMDMzw7lz59C0aVMsX74c2dnZQodHiH5LSQF69ADu3uWTBl+4ANSpI3RUilP3CA1tRKNc1S8hIYF9+umn8hGBzs7O7PLly0KHRRSwdCkfANagAWM5OR9YWSplbMGCvNGvmzZpIkRRioiIYO7u7vLPfuPGjdnVq1eFDkv7iHYWWlJWglzSpCTGmjXTWG1WRdEoVyJajDEcPHgQzs7OOHHiBIyMjLBkyRKEhobio48+Ejo8ooBZs4CqVXl3kwMHPrCyRMJHvHp65m2sR33q8qtbty7++usvHDlyBNbW1rh37x7at2+PmTNnIiMjQ+jwtEe9esBff/El0Qkav6TPnvFyXrdvA9Wr89usDRtq6OCqJ2GMeuf+V1paGiwtLZGamgoLCwuhw9EZstqrp06dAgC4uLjA19cXzZs3FzgyUlYbNvASrvb2PLErV+4DGzDGR73KpjJZtYrfv9VTKSkpmD17Ng68z4jr16+Po0ePomXLlgJHRoiOS0jgAyAePADs7PhtVhEmc4rkI9RCRzTi4MGDaNy4MU6dOgVjY2OsXLkS//77LyVzWm7qVKBmTV4dx8enFBtIJHyYrCyh8/bmr/X0/5VVq1bF/v378ddff8He3h6PHj1Cu3btsGnTJkhLHG1CkJDA/3NApdZ0hsYuqayc14MH/AssIECUyZzC1H3/VxV8fHxY7dq1mampKWvZsiULDAwscf2DBw+yZs2asXLlyrHq1auz0aNHs+Tk5FIfj/rQqU5GRgYbPXq0vL9Qq1at2J07d4QOi6jQ3r28+0nlyoy9fKnAhmvW6GVFieKkpKSwAQMGyP+teHh4sGfPngkdlnhRHzqdo5FL+vQpY3Xq8APVqsVYZKQaD6Y8RfIR0Sd0R48eZcbGxuzHH39kYWFhbMaMGczc3Jw9LaYMx6VLl5iBgQH7/vvvWWRkJLt06RJr3Lgx+/TTT0t9TEroVOPu3bvMycmJAWAGBgZs6dKlLOeDveeJtsnJYczJiX8/enkpuPGGDXlJ3bx5ep/USaVStmPHDmZmZsYAMBsbG3b27FmhwxInSuh0jtovaVQUY7Vr84PUqSNYOS9F6FRC16ZNGzZp0qQC7zVq1Ih5enoWuf66detYnTp1Cry3ZcsWZm9vX+wxMjMzWWpqqvwRExNDCZ2SDh8+zMqXL88AMFtbW3bx4kWhQyJq9Ntv/DuyfHnGFG5U2rIlL6nz8tL7pI4xxu7cucMaN27MADCJRMKWLl1K9Yv/ixI6naPWS5q/Nmv9+ozFxKjhIKqnM6Ncs7OzERISAnd39wLvu7u7459//ilym/bt2yM2NhZ+fn5gjOHZs2f45Zdf0Ldv32KPs3r1alhaWsofNWvWVOnvoU+ys7Px9ddfY9iwYXjz5g26d++OmzdvokuXLkKHRtTok0+A1q2BN2/yikOU2tdfA9u28efffsv71Om5Jk2aIDg4GOPHjwdjDEuXLkXv3r2RlJQkdGiEaJ+4OD6aNSqKD6G9eFFrynkpRO3ppRLi4uIYAHblypUC769atYo1aNCg2O3+97//sQoVKjAjIyMGgPXv359lZ2cXuz610KlGTEwMa9eunbwPkLe3N7Uq6JE//+T/+TUzYywurgw72Lgxr6Vu5UqVx6etDhw4UKC1OyAgQOiQxCEykrHhw0XfB4qUnlouaXw8nywT4C100dEq3Ln66UwLnYxEIinwmjFW6D2ZsLAwTJ8+HYsXL0ZISAj++usvREVFYdKkScXu39TUFBYWFgUeRDHnz5+Hi4sLgoKCUKlSJfzxxx9YuXIlDA0NhQ6NaEivXkD79kBmJm9oU9isWcDatfz5woXAd9+pND5tNWLECAQHB8PJyQkJCQno1q0b1q5dS6NgHR2Bgwf5kugElV/SpCRezktWm/XCBT6qVVepP78su6ysLGZoaMiOHz9e4P3p06ezTp06FbnNl19+yb744osC7126dIkBYPHx8aU6Lg2KKL3c3Fy2atUqZmBgwACwFi1asMePHwsdFhHIhQv8P8ImJkr0N165Mq+lbuNGlcanzTIyMtiXX34pbwHv06cPe/HihdBhCeftW8YePeJLohNUekmfP2esaVP+PVKjBmNa+ndJ4y10oaGhSE5OLvbnWVlZZdqviYkJXF1dcfbs2QLvnz17Fu3bty9ymzdv3sDAoOCvJWslYno615W6pKamYsCAAfD29oZUKsVXX32Ff/75B3W0sQYeUYmuXfkjOxtYubKMO/H2BpYs4c/1tPZrUczNzXHgwAH8+OOPMDU1hZ+fH9q1a4eIiAihQxNGWBhQvz5fEp2gskuakgL07AncucNrs168qJ21WRWligwyKCiITZw4kQ0cOJCdPHmywM9ycnLYwIEDC/WDKy3ZtCV79uxhYWFhbObMmczc3Jw9efKEMcaYp6cnGzFihHx9X19fZmRkxLZv384eP37MLl++zFq1asXatGlT6mNSC92HhYeHs4YNGzIAzNTUlO3evVvokIhIXL7M/1NsZMTY+3+mipNK+YhXaqkrUmhoKHNwcGAAmJWVFQsKChI6JM2jUa46RyWXNH9t1mrVGLt/X2XxCUGwaUtyc3PZjz/+yPr06cNu3bpV4GcjR44s8359fHxYrVq1mImJCWvZsmWBTsGjRo1inTt3LrD+li1bmLOzMytXrhyztbVlw4cPZ7GxsaU+HiV0JQsKCmJVq1ZlAJi9vT0LDg4WOiQiMt278+/T/8w4pBiplE84LEvq1qxRWXy6ICEhgbVs2ZIBYOXKlWMnTpwQOiTNooRO5yh9SRMSGHN25jupXp2xsDCVxicERfIRldVyffToEY4fP44nT54gOTkZf/75J/r16wepVIrw8HC0bdsWO3fuVMWh1I5quRbv9OnTGDRoEN68eYNWrVrh1KlTsLGxETosIjKBgbyyjokJ8PixEjMEMMbLhC1dyl8vXw4sWqSqMLVeRkYGBg0ahD///BMGBgbYsmULpk6dKnRYmnHjBuDqCoSEAFT7VicodUnj43lt1ocPgRo1+ACIBg3UEqcmKZSPqCKD9Pb2ZoaGhkwikRR4GBgYMDMzMzZ9+nSWnp6uikNpBLXQFW3v3r3M0NCQAWC9e/fWqmtKNK9TJ/4f5a+/VsHOVq3Ka6lbtIgmH84nJyeHjRs3Tj5YYt68eSw3N1fosNSPWuh0TpkvaXQ0Y/Xq8Y0dHBiLiFBLfELQ6C3XY8eOsfbt27MTJ06w+Pj4QqWdfvjhB3b06FFlD6MR27ZtY05OTqxBgwaU0OUjlUrZqlWr5H8wRo4cWeK8foQwxti5c3nz0pVygHnJ1q3LS+o8PSmpy0cqlbIVK1bI/40OGTKEZWZmCh0WIeoXFZVXAaJ2bf5ah2g0oZswYQJ7+4ExxmPHjmUZGRnKHkpjqIUuT25uLvv666/lfyg8PT2ZlP6QklKQShlr355/z86eraKdbt6cl9TNmUNJ3X/s379fPqF6ly5d2KtXr4QOiRD1iYzkLXIAY3Xrat2kwaWh0WlLLC0tYWZmVuI6kydPxqFDh5Q9FNGwd+/eYcyYMdi6dSskEgm+//57rF69uthJnQnJTyLJ6+62axefSUBpM2YAPj78+YYNgJcXT+8IAGDkyJHw8/NDxYoV4e/vj65du+puubCHDwE3N74kOkGhSxodzfvMRUcDDRsCAQG6PWlwKSid0JmYmODevXslruPq6oqbN28qeyiiQVlZWRg4cCAOHDgAQ0ND/PTTT5g+fbrQYREt06sX0KIF8Pp1Xh6mtClT8na2Zk3egAkCAOjZsyf8/f1hbW2N0NBQdOjQAU+ePBE6LNV7/RoICuJLohNKfUnj4ngy9+QJn7ju4kU+EELPKZ3QTZo0CV988cUHE7bc3FxlD0U0JCMjAx9//DFOnDgBU1NTHD9+HMOHDxc6LKKFJBLA05M/37JFhX97p0wBNm3iz5cvV2IWY93UsmVLXL58GbVq1cKjR4/w0UcfffA/3oRohYQEnsw9fswnC75wgU8eTJRP6Ozt7TF9+nS0bdsWU6dOxe3btwutc/v2bbx7907ZQxENePnyJdzd3XHu3DmYm5vDz88P/fv3FzososU+/xyoW5ffct2zR4U7njkzr97rokVU+/U/GjRogCtXrqBx48aIj49Hx44dcePGDaHDIqTsiqrNWuY5kXSPSkp/TZ48GRs2bICvry9cXFxQo0YN9OvXD2PHjsUnn3yCTp066c/cSFrs2bNn6NKlC65evYrKlSvj/Pnz6Natm9BhES1nZAR88w1/vn49kJOjwp1/801e69z8+cDmzSrcufarUaMGAgMD0a5dO7x8+RI9evSg7i9EOyUnAz16APfv89urFy8CtWoJHZW4qHI0xoMHD9i4ceNY5cqV5XPRubi4MH9/f1UeRu30cZRrXFycfLoWGxubQpU+CFHG27eM2djwwWj79qnhAIsX541+/f57NRxAu6WmprJ27doxAKxq1ars7t27QoekvJQUxn76iS+JTij2kj5/zliLFvzft60tY+HhgsQnBEEqRfzXixcvYGxsjIoVK6pj92qlb5UiXrx4gc6dO+Pu3btwcHDAuXPnUL9+faHDIjpmzRpgwQKgWTPg5k3ev05lGOMjXtes4a/XrwfmzFHhAbRfamoq3N3dce3aNdja2uLy5cuoow8Fy4l2k91mvXsXqFYN8PcHnJyEjkpjFMlHVHLLtShVqlTRumTOx8cHzs7OaN26tdChaMzr16/x8ccf4+7du7Czs4O/vz8lc0QtJk4EzM2B27eBv/9W8c4lEuDbb4GFC/nruXP5ayJnaWmJP//8E02aNEFCQgK6d++OuLg4ocMqu+fP+Wjn58+FjoSoSKFLmpAAdOnCkzlbWz41iR4lc4pSW0KnjaZOnYqwsDAEBwcLHYpGZGdn4/PPP5f3mfv777/h6OgodFhER1WuzJM6IK8hTaUkEmDFCj7qFQC8vfmUJjRPnVyVKlVw5swZ1K1bF0+ePEHPnj2RnJwsdFhlExMDTJvGl0QnFLiksbG8IPT9+3zgQ0AA0KiR0CGKGiV0eio3NxcjR47E33//jfLly+P06dNo0qSJ0GERHTdrFmBszL+bg4LUdJBFi/IyxmXLeGJHSZ2cra0tzp07B3t7e9y/fx+9e/dGamqq0GERkichgSdzjx7xgQ+BgXy+OVIiSuj0EGMM06ZNw7Fjx2BsbIzjx4/Dzc1N6LCIHrC3B778kj9fu1aNB5o/H9i4kT9fvZqPhqWkTq527do4e/YsrK2tERISgn79+uHNmzdCh0UIN24cEBnJ55kLDATozlGpUEKnh5YvX46dO3dCIpHgp59+Qq9evYQOieiRb77hd0dPnAAePFDjgWbNArZt4883bABmz6akLp9GjRrh77//hoWFBS5duoSBAwciR6VzyhCioNhYvkxMABo04Mmcg4OwMWkRSuj0zN69e7H0fakkHx8fDB48WNiAiN5xcgJkc1Wrfdq4qVN5IVnZwTw9KanLx8XFBadPn0a5cuXg5+eHSZMmQU0TH6hexYqAuztfEu339CkqTv4S7vgbFeu8H81K5bwUorZpS7SZrk5b8vfff6Nv377Izc2Fl5cXVq1aJXRIRE9dugR06gSYmfEO0FZWaj7gjh28XBjA+9jJBk4QAMDJkycxYMAASKVSLFy4ECtWrBA6JKJPZAMgIiN5X7mAACrn9Z4opi0h4nLz5k188cUXyM3NxZdffomVVPuSCKhDB6BVKyAzE9i5UwMHnDw5rzlwxQqq/fof/fv3x873F2LlypXy56KWmwukpfEl0V4JCXyeuchI5DrWQ9rvF5FbjZK5sqCETg9ER0ejT58+yMjIQNeuXbFnzx5IVDqrKyGKkUh4FzeAd3PLytLAQWfMANat48+p9msh48ePx5IlSwDwKZxOnDghbEAfcusWYGnJl0Q7/ac2662tgbB0rkGXtIwooctHFycWfvXqFfr06YOEhAQ0btwYx48fh4mJidBhEYKBA3kXmWfPgCNHNHTQuXOp9msJlixZgvHjx0MqlWLo0KG4fv260CERXZWSUrg2K91mVQoldPno2sTCOTk5+OKLL3Dv3j3Y2dnhzz//RKVKlYQOixAAfD666dP5802bNDhWwdsbWLyYP581C/j+ew0dWPwkEgm2b98ODw8PZGZmYsCAAXhOlRiIqj1/zlvm7twBqlfnyRyVoVMaJXQ6bMaMGTh//jzMzc3h5+eHmjVrCh0SIQWMHw+UK8fLgV2+rMEDL13KC8sCwMyZdPs1HyMjIxw5cgQNGjRAbGwshg0bhlzqp0ZUJTER6NqV3yq3sQHOn6dJg1WEEjodtXfvXuzYsQMSiQSHDx9G8+bNhQ6JkEIqV86baFg2ZZxGSCTAqlV5LXXz5/PBEgQAr/t6/PhxlC9fHufOncNi2XkiRBlxcbw26717gJ0dH83q7Cx0VDqDpi0pgrZPW3Lr1i20a9cOmZmZWLFiBRbKCpYTIkK3bgEtWgBGRsCTJwJMPbVqFSD7N7JwIZ/ShAYNAQCOHj2KoUOHAgB+//139JdNICgGOTnAq1dApUr8/j0Rt+hooFs34PFjPlnwhQtA3boFVqFLWpjOTVuyfft2ODo6wszMDK6urrh06VKJ62dlZcHb2xu1atWCqakp6tati71792ooWmGlpqbiiy++QGZmJjw8PODl5SV0SISUqHlzoGNH4N074IcfBAjA2ztv9OvKlby1jv6fCwAYMmQIpr/v6Dhy5EhEREQIHFE+xsaAtTX95dcGUVF8nrnHj3kZr8DAQskcQJdUWaJP6I4dO4aZM2fC29sboaGh6NixIzw8PBAdHV3sNoMGDcL58+exZ88ePHz4EEeOHEGjRo00GLUwGGP46quvEBERAQcHB/z0008wMBD9JSYE06bx5a5dQHa2AAHMnQts2cKfr1vH+9VRUgcAWLduHdq3b4/U1FR8/vnn4qn5+vgxLzny+LHQkZCSPHrEZxF/8oT3lQsMBGrVKnJVuqRKYiLXpk0bNmnSpALvNWrUiHl6eha5/p9//sksLS1ZSkpKqY+RmZnJUlNT5Y+YmBgGgKWmpioVu6Zt3bqVAWDGxsYsKChI6HAIKbXsbMZsbRkDGDt8WMBAdu3iQQCMTZ/OmFQqYDDiERsby6pVq8YAsK+++krocLiQEH6dQkKEjoQUJyKCMTs7fp2cnBiLjy9xdbqkhaWmppY6HxF18012djZCQkLg7u5e4H13d3f8888/RW5z8uRJtGrVCt999x1q1KiBBg0aYO7cuXj79m2xx1m9ejUsLS3lD20cDXr79m3MnTsXAP8fddu2bQWOiJDSMzbmI14B4McfBQxkwgRA1j1jyxa6/fpejRo1cPToUUgkEuzduxdHNDZxINFaT5/yPnPx8UDjxrw2K80zp1aiTuiSk5ORm5sLGxubAu/b2NggMTGxyG0iIyNx+fJl3L17F7/99hs2b96MX375BVOnTi32OAsWLEBqaqr8ERMTo9LfQ93evHmDIUOGICsrC3379pX3eSFEm4wdCxgY8CmpwsMFDGTMGH7vF+C3X2mEJwCga9euWLRoEQBg4sSJ4upPR8QlNpZPTRIdDTRsyKcmqVZN6Kh0nqgTOpn/lqlijBVbukoqlUIikeDQoUNo06YN+vTpg40bN2Lfvn3FttKZmprCwsKiwEObzJ49G/fv30f16tXh6+tLZb2IVnJwADw8+PPdu4WNBRMm5PWpW7mSar++t2jRInTs2BHp6ekYMmQIsgXp8EhELSGBt8xFRfGBD+fP8/nmiNqJOqGzsrKCoaFhoda4pKSkQq12Mra2tqhRowYsLS3l7zk5OYExhtjYWLXGK4Tjx49j1/vWhAMHDsDa2lrgiAgpuwkT+NLXV0P1XUvy9dcFa7+uXy9sPCJgZGSEQ4cOoUqVKggJCcEC2eTMQqhRA9iwQYB5bkixZLVZHz3iAx8uXFDo+tAlVY6oEzoTExO4urri7NmzBd4/e/Ys2rdvX+Q2H330EeLj45GRkSF/Lzw8HAYGBrC3t1drvJoWExODcePGAQDmzZuHnj17ChwRIcrp04fPN5qcDPz+u9DRoGDt12++yWu102M1a9aEr68vAGDjxo04ffq0MIHY2ACzZ1Prj1jkr81qb8+TOQcHhXZBl1RJah+ioaSjR48yY2NjtmfPHhYWFsZmzpzJzM3N2ZMnTxhjjHl6erIRI0bI109PT2f29vbsiy++YPfu3WMBAQGsfv36bNy4caU+piKjSoSSm5vLunXrxgCwVq1asaysLKFDIkQlFi3iI9169BA6knwWLswb/bpxo9DRiMKMGTMYAGZlZcUSExM1H8CLF4z9/DNfEmE9e8ZYs2b830f16oyFh5dpN3RJC1MkHxF9QscYYz4+PqxWrVrMxMSEtWzZkgUEBMh/NmrUKNa5c+cC69+/f5/16NGDlStXjtnb27PZs2ezN2/elPp42pDQbdmyhQFg5cuXZ+Fl/MdDiBhFRvK/CxIJY0+fCh3Ne1IpYwsW5CV1q1cLHZHgMjMzWbNmzRgA1q9fPybV9BQvNMeFOMTH8ylJZMlcWFiZd0WXtDCdS+g0TewJ3cOHD1m5cuUYALZt2zahwyFE5bp25V/sK1YIHUk+Uiljy5blJXVLluj9PHW3bt1iJiYmDADbvXu3Zg9Of/2FFxPDWP36/DrY25e5ZU6GLmlhOjMPHSns3bt3GDVqFN6+fYvu3btj8uTJQodEiMqNHs2X+/aJaBo4iYRPYbJmDX+9bBng5SWiADWvWbNmWPm+j+GsWbMQFxcncEREY5484RUgHj0CatfmFSDq1xc6Kr1GCV0+Pj4+cHZ2RuvWrYUOpVjr169HUFAQLCwssHfvXirtRXTS558DFSrwEkCXLwsdzX/Mnw9s3syfr1nDe3HrcVI3e/ZstGvXDunp6ZgyZQqYHp8LvRERwZO5qCigXj0gIIDXaCWComwgn6lTpyIsLAzBwcFCh1KkO3fuYPH7SU6///57OCg4gogQbWFuDgwcyJ+/H1ApLjNmADt28OebNwNTpwJSqaAhCcXQ0BC7d++GsbExTp48iV9++UUzBy5XDnBx4UuiOQ8e8GQuJgZo1Igncyr6W0SXVDkSRv+dKiQtLQ2WlpZITU0VzSTDubm5aNeuHa5fv47+/fvjxIkTNIEw0WmXLwMdO/LkLjGRt9iJjq8vL3HBGE/qtm7lt2b10JIlS7B8+XLY2triwYMHovnuJCr06BHQuTOfPLhpU+DcOaoAoWaK5CPUQqcldu/ejevXr8PS0hI7d+6kZI7ovI8+4ndzXr8Gjh8XOppijBkD7N/PkzgfHz5XnZ7+H9nLywv16tVDQkICli9fLnQ4RNWiongFCFkyd+ECJXMiQwmdFkhJSYGXlxcAyP8HTIiuk0iAkSP58/37hY2lRCNG5NV+3bBBb2u/mpqaYuvWrQCAzZs34+7du+o9YGgoYGrKl0S9YmJ4BYjYWH6b9dw5wMpK5YehS6ocSui0wKJFi/DixQs0adIEU6ZMETocQjRmxAi+vHiR1/kWrfHj+e1WQK9rv/bu3RuffvopcnNzMW3aNPUOkGAMyM7W2xZRjUlI4MmcbADE+fNqa5mjS6ocSuhELjQ0FDt37gQAbNu2DUZGRgJHRIjm1K4NdOnCv+APHhQ6mg+YNi2v3qse137dtGkTypUrh4CAABw5ckTocIgy8tdmrV2b32a1sxM6KlIMSuhEjDGGr7/+GowxDBkyBJ07dxY6JEI0Lv9tV9H/z33OHL2v/Vq7dm14e3sDAObMmYO0tDSBIyJlkpIC9OxZsDZrzZpCR0VKQAmdiP3888+4cuUKypcvj3Xr1gkdDiGC+OILPo1BeDgg0hmFCvL2BhYu5M9nzOD96vTM3LlzUa9ePSQmJmL16tVCh0MU9ewZbxq/fRuoXp3fZqV55kSPErp8xDSxcFZWFhYsWAAAmD9/Puzt7QWOiBBhVKwI9O/Pnx89KmwspbZ8OU/sAGDuXL3rU2dqaoqNGzcC4Ldgnzx5ovqDODkBd+/yJVGduDg+Ncndu4CtLW+Za9BAI4emS6ocmoeuCGKYh27jxo2YM2cObG1t8ejRI5ibmwsSByFi8PvvwKef8u47MTGA1hRIWbUqr7XOy4sndnoy5RBjDD169MCFCxcwZMgQ6k+nDZ4+5VOTREby26sXLvCBEEQwNA+dlnvx4gVWrFgBAFixYgUlc0Tv9e4NWFoC8fHApUtCR6MAb++8W67ffsv72OnJ/6ElEgk2bNgAiUSCo0eP4urVq6o9wNOnwLhxfEmUJyvnFRkJ1KnDa7NqOJmjS6ocSuhEaNWqVXj16hWaNGmC0bIq5YToMVNT4LPP+HOta+iZPZtPOgwAmzbpVZmwFi1aYMyYMQB4zVeV3hBKSQH27OFLohxZOa/oaKBhQ57M1a6t8TDokiqHEjqRiYyMlE/OuW7dOhgaGgocESHiMHQoX/7yC5CTI2wsCpsyhf+lkkh4Ddjx4/UmqZPdZQgKCsLPP/8sdDjkv+7dyyvn1aQJr81ao4bQUZEyoIROZJYuXYqcnBz06NEDvXr1EjocQkSja1c+n2lKCp9oWOt89RXw00+AoSGwdy9P8vTg9qudnR3mzZsHAPD29kaO1mXjOiw8nM8zl5QEuLjwf1g2NkJHRcqIEjoRCQsLw8H3s6euXr2a6rUSko+RER8YAYi4tuuHDB/OkzqJhJcLmz1bL5K62bNno1q1anj8+DH27dsndDgE4H3lunXjU5S0aMGnJlFDOS+iOZTQicjSpUvBGMOnn36KVq1aCR0OIaLz+ed8eeIEkJsraChlN3Qov/0KAJs384ETOp7UVahQQV6PetmyZcjMzFR+pzY2gKcntSiVhaw2a1wc4OwMnDkDVK4sdFR0SZVE05YUQYhpS27duoUWLVpAIpHg1q1baNq0qUaOS4g2yc7mX/avXvF+2x07Ch2RErZv5wMkAD5v3aJFwsajZpmZmWjQoAFiYmKwceNGzJo1S+iQ9FNCAh8AEREB1K/P+8zZ2godFSkGTVtSRkJOLLx48WIAwODBgymZI6QYJiZ5kwxr7W1XmSlT8qY0WbxY52u/mpmZYcmSJQCAb7/9Funp6crtMD0d8PfnS1I6z58DPXrwZK52bX6bVUTJHF1S5VBCl8/UqVMRFhaGYA3XF7p27RpOnjwJAwMDLF26VKPHJkTbyKYv+fVXHbhTOXu2XtV+HTVqFOrXr4/k5GRs3rxZuZ09esRHyjx6pJLYdF5yMq/NGhbGR7GKsDYrXVLlUEInAqtWrQIAjBgxAg0bNhQ4GkLEzd0dKF+edwO6eVPoaFTA2zuvTNiMGcDatcLGo0ZGRkZYtmwZAGDLli14+/atwBHpiYQEPjXJrVu8z8KFC1SbVQdRQiewO3fu4OTJk5BIJPLarYSQ4pUrxxsaAODUKWFjUZkVK/L60Hl6AkuW6EDzY9EGDhyIWrVqITk5GQcOHBA6HN0XE8OTubAwXjvP319jtVmJZlFCJ7A1a9YA4F9y1DpHSOn068eXf/whbBwqI5HwgRGrV/PXy5cD8+bpZFJnZGSEmTNnAuA1q6V6MsGyICIj+cihR4+AWrX4SKJGjYSOiqiJViR027dvh6OjI8zMzODq6opLpSzmeOXKFRgZGaFFixbqDbCMHj9+jKNHjwIAtc4RooA+ffgyOBhITBQ2FpXy9AS+/54/X78emDZNJytKjB07FpaWlggPD8fp06fLthNjY94XzNhYtcHpigcPeDL39CmvyRoYCNStK3RUJaJLqhzRJ3THjh3DzJkz4e3tjdDQUHTs2BEeHh6Ijo4ucbvU1FSMHDkS3bt311Ckivvuu+8glUrRp08f0SadhIiRrS0gG4xe1nxAtKZPB374gbfabd8OjB2rxZPuFa1ixYqYMGECAGCDbKSvopo2BWJj+ZIUdPs2n5okPp7PMxcYCDg4CB3VB9ElVY7oE7qNGzdi7NixGDduHJycnLB582bUrFkTO3bsKHG7iRMnYtiwYXBzc9NQpIp5/vw59u/fDwDw9PQUOBpCtM/HH/Olztx2zW/8eODAAV4mbN8+YORInUvqpk+fDiMjIwQEBCAkJETocHTHzZtAly58ihIXF5pnTo+IOqHLzs5GSEgI3N3dC7zv7u6Of/75p9jtfH198fjxY/mcRx+SlZWFtLS0Ag9127t3L7KystCqVSt06NBB7ccjRNfIErrz5/mEwzrnyy+BY8d4zbPDh4EJE3Tq9qu9vT0GDx4MgM8BqrA7dwB7e74kXFgYHzH08iXQti0fzapF5bzokipH1AldcnIycnNzYfOfOiA2NjZILKbjzKNHj+Dp6YlDhw7ByMioVMdZvXo1LC0t5Y+aap6bRyqVYteuXQCAKVOmUM1WQsqgRQvA2hrIyACCgoSORk0+/5wncwYGwN69fFoTHRooMfV9pYwjR47gxYsXim2ck8NLV+XkqCEyLRQRwScNTk4GXF2Bv/8GKlUSOiqF0CVVjqgTOpn/JjyMsSKToNzcXAwbNgzLli1DAwWGZS9YsACpqanyR0xMjNIxl+TMmTOIiopCpUqV5P9DJYQoxsAgb/qSM2eEjUWtBg7kt10lEmDbNmDBAp1J6tq1a4cWLVogMzMT+/btEzoc7RUdzWuzJiQATZrwZM7SUuioiIaJOqGzsrKCoaFhoda4pKSkQq12AJCeno7r169j2rRpMDIygpGREZYvX45bt27ByMgIFy5cKPI4pqamsLCwKPBQp507dwLgs6aXL19ercciRJfJemP8/bewcajdiBGArN/w2rV51SW0nEQiwZQpUwAAO3bsoClMyiIhgSdz0dF8frmzZ4GqVYWOighA1AmdiYkJXF1dcfbs2QLvnz17Fu3bty+0voWFBe7cuYObN2/KH5MmTULDhg1x8+ZNtG3bVlOhFys2NhZ/vO/FPXHiRIGjIUS7yVroQkL4nSadNnEisHEjf754cV4dWC03bNgwWFhYICIiAufOnRM6HO1SVG3W6tWFjooIRNQJHQDMnj0bu3fvxt69e3H//n3MmjUL0dHRmDRpEgB+u3TkyJEAAAMDAzRp0qTAo1q1ajAzM0OTJk1gbm4u5K8CANi9ezekUik6d+4MJycnocMhRKvZ2fEpDhgD9CIXmDUrr3Vu7lxA2XqoImBubo7Ro0cD4HOOllr9+sDFi3ypj54/L1ib9fx5PqJAi+n7JVWW6BO6wYMHY/PmzVi+fDlatGiBwMBA+Pn5oVatWgCAhISED85JJxZSqVTeT0Q2BxMhRDk9evDlxYvCxqEx3t6Alxd/LkvwtLxPnexuxalTp/Ds2bPSbVSxIp+eo2JF9QUmVnFxfJ45WW3Wc+eAOnWEjkpp+nxJVUHCmJZ/E6hBWloaLC0tkZqaqtL+dP7+/ujatSssLCyQmJiIcuXKqWzfhOirP/4A+vfn3YcePhQ6Gg1hjCdyixfz1/Pn87JhWjxivl27dvj333+xYcMGzJ49+8MbxMXxQSLTpvEWKn0RFcX7zEVF8Ra58+d1pjarvl7SkiiSj4i+hU6XyFrnBg8eTMkcISrSsSPPY8LD+cT4ekEiARYtyutHt3Yt8PXXWj1Pney2q6+vL0rVzvDsGbBmDV/qi4cP+Qc+KoqX8bp0SWeSOUA/L6kqUUKXj4+PD5ydndFaVlNIhTIyMvDLL78A4KNbCSGqUakSnxAf4JPi65XZs4GdO3mC5+Oj1WXChgwZAjMzM9y9e5cqRxRFVs4rLi6vnFft2kJHRUSEErp8pk6dirCwMAQHB6t837/++itev36NevXqFTlClxBSdl268KXe9KPLb+LEgmXChg3TytIZlSpVwoABAwDwVjqSz7Vr/EOelJRXzsvOTuioiMhQQqchBw4cAMBb56gyBCGq1bUrX/r7CxqGcL78Evjf/wBjY+Dnn4EhQ4B374SOSmFjxowBwCtHZGVlCRyNSFy7xkf+vHwJtG+vdeW8iOZQQqcBiYmJ8H//l2b48OHCBkOIDvroI7589Ig3YuilAQOAkycBU1Pgt9+AMWO0rk9dt27dYGdnh5cvX+LMh8p/VK3KbzHr8iS6t28DvXsD6em8hU4Ly3kpQh8uqTpRQqcBx48fh1QqRZs2beDo6Ch0OITonMqVebciALh6VdhYBNW7N2+pMzICDh4Epk7VqilNDA0NMXDgQADAsWPHSl65Vi1g926+1EXh4XyeuZcvATc3Ppy7QgWho1IrXb+k6kYJnQbIvpiobish6iPrmnrlirBxCK5fP+Cnn/hAiZ07gW++0aqkbsiQIQCA33//HW/fvi1+xbdvgXv3+FLXPHnCpyZJSgJatAD8/HQ+mQN0+5JqAiV0ahYfH49Lly4BAL744guBoyFEd3XowJeXLwsbhygMGQL8+CN/vmEDsHy5sPEooG3btqhVqxYyMjLg5+dX/Ir37/NC9Pfvay44TUhI4H3mYmOBRo2AM2d0+jZrfrp6STWFEjo1+/XXX8EYg5ubGxwcHIQOhxCd1bEjX16/Tv/DB8A7I8lKgy1dmlcHVuQkEon8bsYHb7vqmuRknsw9fgw4OvIKENbWQkdFtAQldGr2v//9DwAwaNAggSMhRLc5OvKZHHJygH//FToakZgxI6/265w5wPr1wsZTSrKE7o8//sDr168FjkZDEhOBbt0K1malcglEAZTQ5aPqiYWTk5Nx5X2Hns8++0wl+ySEFE0iyetHRwldPl5evP4rwPvTLVok+j51Li4ucHR0RGZmJs6ePSt0OOr39ClvYr5zB6henbfM0QA6oiBK6PJR9cTCfn5+kEqlaN68Od1uJUQD2rblS0ro8pFIeCvd6tX89cqVvOVOxFOaSCQS9O/fHwBw8uTJ4lYCTEy0un4tgLxyXhERfHjnpUu875we0pVLKhRK6NTojz/+AAD069dP4EgI0Q9t2vDltWvCxiFKnp68PBgAbN3K+9iJePLhTz75BABw6tQp5BZVzszFBcjKyqv7po1u3uTJXEwMT+IuXwbq1RM6KsHowiUVEiV0apKdnY2///4bACV0hGiKqytgYMDLXcbFCR2NCE2ZUrBM2JAh/C+oCHXo0AGVKlXC8+fP8a8uNrn+8w+fLPj5c57BBAYC9vZCR0W0GCV0ahIQEID09HRUr14drVq1EjocQvSCuTmf9gCgVrpijRjBJx82MQF+/RXo3x8Q4cADY2Nj9OnTBwCfk66Q+/eBli21c46Lc+f4pMGpqbzMyYULNJoV2n1JxYASOjU5deoUAKBv374wMKDTTIimyPrRUUJXggEDgFOngPLl+TxnAwaIsqVOdttV1n2lgLdvgdBQ7ZujJiCAT/785g3g7q7z5bwUoa2XVCwo01CTc+fOAQA8PDwEjoQQ/eLqypc3bwoahvj17MmTOXNz4OxZYNgw0fWp69mzJwwMDHD//n3ExsYKHY7ygoOBjz8GMjP58uRJfv4JUQFK6NQgPj4eYWFhkEgk6Nq1q9DhEKJXWrTgy9BQQcPQDh99BJw4wW+/Hj8OjB8vqtGvlStXlndZ0frpS+7d47V2MzKArl35bW9TU6GjIjqEEjo1OH/+PADA1dUVVapUETgaQvRL06Z8YMSzZ3yuVvIBPXoAR4/mDZSYNUtU89R1794dAODv7y9sIMqIjOQtoi9eAK1bA7//DpiZCR0V0TGU0OWjqomFZbdbe/TooYqwCCEKKF8eaNCAP6fbrqU0YACwdy9/vmULsGyZsPHkI7vL4e/vD5Y/0XR0BH7+WfwT8MbH82QuIQFo3Bj480+gYkWhoxIlbbmkYiVhTET/FROJtLQ0WFpaIjU1FRYWFgptyxiDvb094uPjcf78eXTr1k1NURJCijN0KG90Wr2aT79GSmnrVmD6dP58wwZg9mxh4wHw+vVrVK5cGTk5OXj8+DHq1KkjdEill5ICdO7Mb7fWqcMnDbazEzoqokUUyUeohU7FwsPDER8fDzMzM7SX1SEihGiUrB8dtdAp6OuvgRUr+PM5c4BvvxX89qu5uTnavJ8xOiAgIO8Hz54BGzfypRjFxQGdOvFkzs6OT1VCyVyJxH5JxY4SOhW7fPkyAKBNmzYwoz4ShAhCNhcdzWdVBt7evN6r7PmcOYIPlJD957jABMNxcTw2Mc4gHR7OB5yEheUlc3Qf8YPEfEm1ASV0KnblyhUAwEcffSRwJIToL2dnvnz4UHQzcYifRAIsXw5s2sRfb9oEfPWVoCeyXbt2AICgoCDBYii10FCgQwfg6VOgfn3gyhXAyUnoqIge0IqEbvv27XB0dISZmRlcXV1x6dKlYtc9fvw4evbsCWtra1hYWMDNzU1egksTKKEjRHi1agHlyvG5cqOihI5GS82cCezfz0e/7t8PfP65YDO+yhK6O3fuICMjQ5AYSiUgoGA5r8uXgdq1hY6K6AnRJ3THjh3DzJkz4e3tjdDQUHTs2BEeHh6Ijo4ucv3AwED07NkTfn5+CAkJQdeuXdGvXz+EamBSqufPnyM8PBwA4ObmpvbjEUKKZmDAa50DdNtVKSNHAr/9xqfYOHmSz6OWmqrxMOzs7FCzZk1IpVJcv35d48cvlZMngV69gLQ0PhDi4kWgWjWhoyJ6RPQJ3caNGzF27FiMGzcOTk5O2Lx5M2rWrIkdO3YUuf7mzZsxb948tG7dGvXr18e3336L+vXrF106RsX++ecfAICzszPNP0eIwGS3XcPChI1D6/Xrx8tTWVjwAvKyFigNa/u+ppv8tqulJY/N0lLjsRTy00/AZ5/xJuH+/YG//hJHXFpGTJdUG4k6ocvOzkZISAjc3d0LvO/u7i5Pnj5EKpUiPT29xAQrKysLaWlpBR5lIfuiodGthAhP1m2JWuhUoFMnwN+ftzjdvAl4eADp6RoNQZbQXZMV6a1bl7eK1a2r0TgKOX4cGD0ayM0FRo0Cfv2VJg0uI7FcUm0l6oQuOTkZubm5sLGxKfC+jY0NEks5BfyGDRvw+vVrDBo0qNh1Vq9eDUtLS/mjZs2aZYr31q1bAHiFCEKIsOrV48vHj4WNQ2e4uPA+YlZWQEgI8OmnvEVKQ1q2bAkg73sWOTm8pTAnR2MxFHLxIp/0UCoFxo7lkzMbGQkXj5YTwyXVZqJO6GQkEkmB14yxQu8V5ciRI1i6dCmOHTuGaiX0ZViwYAFSU1Plj5iYmDLFKfuiadasWZm2J4Sojiyhi4gQNg6d0qgRr3RQoQJw4QIwfDhvmdKA5s2bAwAiIyP5XZQ7d3iL4Z07Gjl+ITduAJ98AmRn80obO3fyzpukzIS+pNpO1J8+KysrGBoaFmqNS0pKKtRq91/Hjh3D2LFj8fPPP3+wBJepqSksLCwKPBSVnJyM+Ph4AEDTpk0V3p4Qolqy2zbPnvF66ERFWrUCTpwATEz47cUpUzQy+XDVqlVRo0YNAHy0q6AePeIDRNLTga5dgcOHqWWOCE7UCZ2JiQlcXV1x9uzZAu+fPXu2xH5qR44cwejRo3H48GH07dtX3WECAG7fvg0AqFu3LipSnT5CBFepEiDrOku3XVWse3fg0CE+Z90PPwCLF2vksLJWuptClgCJjwfc3fOmJjlxgvrMEVEQdUIHALNnz8bu3buxd+9e3L9/H7NmzUJ0dDQmTZoEgN8uHTlypHz9I0eOYOTIkdiwYQPatWuHxMREJCYmIlXNQ+3pdish4iNrpaOETg2++ALYvp0/X7mS135Vsxbva7rJ+9Fp2vPnfGqSJ0/4Pf0//+SjfwkRAdEndIMHD8bmzZuxfPlytGjRAoGBgfDz80OtWrUAAAkJCQXmpNu1axfevXuHqVOnwtbWVv6YMWOGWuO8e/cuAEroCBETWUIXGSlsHDpr0iReVQIA5s4FPD3VWiZM9v0q+77VqMePgfbtgbt3AVtb4OxZ4ANdfwjRJK246T9lyhRMmTKlyJ/t27evwGt/f3/1B1QE2YTCjWSzmRJCBCcbsF7GcU6kNBYu5LdeFy0C1q7lhTj37OF97FRM9v368OFDoHlzPsmxubnKj1NIcDDQty9voatdm88zRxUgVE6Tl1QXib6FTlvIEroGDRoIHAkhRIYSOg2QSHhSt3cvLxN28CBPfso4n2dJ6tevDwB48eIFUl694rc7DQ1VfpwC/PwKlvO6ehVo2FC9x9RThoaauaS6ihK6fHx8fODs7IzWrVsrtN2rV6+QlJQEIO8LhxAiPEroNGjMGOCPP3jzyrlzfDLi9yP/VaV8+fLyeUIfnjvH+7M9eqTSYxSwdy+v/PDmDR8IERAAVK+uvuPpuUeP1H9JdRkldPlMnToVYWFhCA4OVmi7R+8/fba2tjTClRARkSV0sbHCxqE3PDzyKkrcugW4uam8VEfD961jD+/fB86cUU/FCsZ438CxY/k8eyNH8mSVvt/VKj1dfZdUH1BCpwKy260NqRmeEFGRJXTPnvH5X4kGtGrFb0vWrw9ERwMffQRcv66y3cu+Z8OfPlXZPgtgDJg2DViyhL/28gL27VNLn0BCVIkSOhV4/H5OhHqyqekJIaJgbc3/DjMGJCQIHY0eqVMHuHIFaNsWePmS96mLilLJrmXdWiLUdR99xQo+HYtEAvj4AKtW8eeEiBwldCogmzbFwcFB4EgIIflJJPzuHwC87+ZKNMXamt8/a96cn/w+fYAXL5TerWzKquhS1vNWyP79eS1zO3bwKhiEaAlK6FSAEjpCxEs2VdizZ8LGoZcsLIDTpwF7e+DBA17zNCtLqV3Kvmejnz8Htm3Lu6+urHPngHHj+HNPT2DiRNXsl5RazZqqvaT6hhI6FYh53/Rfkz6FhIgOtdAJrEYNPvWHhQUQGAiMHq3U5MOyhC4xKQlZ48bxlkBl3b4NfPYZ8O4dMHQov81KNM7aGpg6VTWXVB9RQqckxhi10BEiYtRCJwJNmwK//QYYGwNHj/KBBmVUtWpVlCtXDgAQu22b8rdxY2P57eD0dKBzZ8DXFzCgP41CePGCT2Oogjvzeok+tUp68eIF3rx5AwCwt7cXOBpCyH9RC51IdOvGK0gAvKLEypV8tIqCJBJJ3m3XuXN5XdWyio4Gevfm1S2cnHjSaWpa9v0RpTx5AowYodwl1WeU0OVTlomF499PnFm1alWYmZmpKzRCSBlZWfFlcrKwcRDwv9YrV/LnixYBQ4YAr18rvBvZf57jlInF359PsXLvHq/N6ucHVK6szB4JERQldPmUZWLh58+fAwCqyZoBCCGiIvsb/eqVoGEQGS8vPi2IkRHw88988uHISIV2Yf2+k9XzshyfMWDLFqBHj4LlvKg2K9FylNApSZbQWVMvTkJEqVIlvqSETiQkEmDyZODiRd7B8c4d3lJ25kypdyH7D7TCCd3bt3xQxowZvALEl1/y+fLeT4VCiDajhE5J1EJHiLhRQidSHToAISF5kw97eADffVeqfnWy/0AnWVvz2rGlER0NdOwIHDjAq79v2sSfvx9gQYRnbg60a1f6S0oKooROSUnve1pTCx0h4kQJnYjVqMEL3o8dy6cymT+/VP3q5C10bm5AaUouyvrLhYQAVavy1sCZM6kChMg0bMjvflMVzbKhhE5JdMuVEHGTJXQvXwoaBimOqSnw448F+9V17sxvjxZD3kJXmqHLJ08W7C8XEsJH3BKiYyihU9LL938lqlSpInAkhJCiVKzIl69fl2mWDKIJ+fvVWVnxpGvGjGJXl33fvgwKAm7cKH6/UVHAyJG8v9yQIcDly9RfTsRu3OAfhZIuKSkeJXRKysjIAABUlP3VIISISvnyec8zM4WLg5RChw584mGJhLfaHTpU5Gqy79v0kvaVnQ0MHgykpvKRtAcOFPwwEKJjKKFTUno6/0qpUKGCwJEQQoqSv8/7+znAiZh17w4sXsyfT5wIPHxYaBVZQpdR0n7mzQOCg4EqVXiSaGys+lgJERFK6PIpy8TCsoSOWugIEScjo7y/5ZTQaYlFi4CuXfl98oEDC/Wny5/QsaLuo//2G/D99/z5/v0AlWUkeoASunzKMrEw3XIlRPxkd9pK6GdPxMTQkN9urVaNz1M3ahSvufr+UTEtDQAgBfAmOrrAz3DjBjBmDN/P3LnAxx8L93sQokFGQgeg7czNzVGhQgW65UqIiJUvz7tSUQudFrG15Umduzvwv//xx3vlAVR4v3zz2Wcoctqydu2Ab7/VTKxEJZydgUePACqLXjaU0CkpNDRU6BAIIR9gYsKX2dnCxkEU1KMHsHEjvwWb7+JJ8IEBEU2aAMeOUb85LWNmBtSrJ3QU2otuuRJCdJ6hIV/m5gobBymDmTOB9HQgK6vg48ED3r/uwYPCPwsJoX5zWigqildji4oSOhLtpBUJ3fbt2+Ho6AgzMzO4urri0qVLJa4fEBAAV1dXmJmZoU6dOti5c6eGIiWEiBEldDro5Ut+S5ZmjNYZdEmVI/qE7tixY5g5cya8vb0RGhqKjh07wsPDA9HR0UWuHxUVhT59+qBjx44IDQ2Fl5cXpk+fjl9//VXDkRNCxIISOkKIrhN9Qrdx40aMHTsW48aNg5OTEzZv3oyaNWtix44dRa6/c+dOODg4YPPmzXBycsK4cePw1VdfYf369RqOnBAiFrKE7t07YeMghBB1EXVCl52djZCQELi7uxd4393dHf/880+R21y9erXQ+r169cL169eRk5NT5DZZWVlIS0sr8CCE6I4GDYCmTalQACFEd4k6oUtOTkZubi5sbGwKvG9jY4PExMQit0lMTCxy/Xfv3iE5ObnIbVavXg1LS0v5o2bNmqr5BQghonD8OHD7Nq8ARXSErS2wZAlfEp1Al1Q5ok7oZCQSSYHXjLFC731o/aLel1mwYAFSU1Plj5iYGCUjJoQQola2tsDSpfTXX4fQJVWOqOehs7KygqGhYaHWuKSkpEKtcDLVq1cvcn0jIyNUrVq1yG1MTU1hamqqmqAJIYQQQjRM1C10JiYmcHV1xdmzZwu8f/bsWbRv377Ibdzc3Aqtf+bMGbRq1QrGNMkkIYQQQnSQqFvoAGD27NkYMWIEWrVqBTc3N/zwww+Ijo7GpEmTAPDbpXFxcThw4AAAYNKkSdi2bRtmz56N8ePH4+rVq9izZw+OHDlS6mPKbtHS4AhCCCGECEWWh8jykhIxLeDj48Nq1arFTExMWMuWLVlAQID8Z6NGjWKdO3cusL6/vz9zcXFhJiYmrHbt2mzHjh0KHS8mJoYBoAc96EEPetCDHvQQ/BETE/PB3EXCWGnSPv0ilUoRHx+PihUrljj4AuDZc82aNRETEwMLCwuljtu6dWsEBwcLvg8xxaKq86tr50VV+6Hzq779iO27QVX7EUssYju/YjkvqtoPnV/17qO055cxhvT0dNjZ2cHAoORecqK/5SoEAwMD2NvbK7SNhYWF0h96Q0NDUexDbLEAyp9fXTwvdH7FHwsgnu8GVe1HTLEA4jm/YjsvdH7FHwtQuvNraWlZqn2JelCEvpk6daoo9qGq/agqFlXQxfNC51d9+1DlflRBTL+TmGJRFV08L3R+1bcfMZ3b/OiWq5LS0tJgaWmJ1NRUlWXsJA+dX/Wi86s+dG7Vi86vetH5VS91nF9qoVOSqakplixZQvPYqQmdX/Wi86s+dG7Vi86vetH5VS91nF9qoSOEEEII0XLUQkcIIYQQouUooSOEEEII0XKU0BFCCCGEaDlK6AghhBBCtBwldIQQQgghWo4SOgU9efIEY8eOhaOjI8qVK4e6detiyZIlyM7OLnE7xhiWLl0KOzs7lCtXDl26dMG9e/c0FLX2WLVqFdq3b4/y5cujUqVKpdpm9OjRkEgkBR7t2rVTb6Baqiznlz67pffy5UuMGDEClpaWsLS0xIgRI/Dq1asSt6HPb/G2b98OR0dHmJmZwdXVFZcuXSpx/YCAALi6usLMzAx16tTBzp07NRSpdlLk/Pr7+xf6nEokEjx48ECDEWuHwMBA9OvXD3Z2dpBIJDhx4sQHt1HFZ5cSOgU9ePAAUqkUu3btwr1797Bp0ybs3LkTXl5eJW733XffYePGjdi2bRuCg4NRvXp19OzZE+np6RqKXDtkZ2dj4MCBmDx5skLb9e7dGwkJCfKHn5+fmiLUbmU5v/TZLb1hw4bh5s2b+Ouvv/DXX3/h5s2bGDFixAe3o89vYceOHcPMmTPh7e2N0NBQdOzYER4eHoiOji5y/aioKPTp0wcdO3ZEaGgovLy8MH36dPz6668ajlw7KHp+ZR4+fFjgs1q/fn0NRaw9Xr9+jebNm2Pbtm2lWl9ln11GlPbdd98xR0fHYn8ulUpZ9erV2Zo1a+TvZWZmMktLS7Zz505NhKh1fH19maWlZanWHTVqFPvkk0/UGo+uKe35pc9u6YWFhTEALCgoSP7e1atXGQD24MGDYrejz2/R2rRpwyZNmlTgvUaNGjFPT88i1583bx5r1KhRgfcmTpzI2rVrp7YYtZmi5/fixYsMAHv58qUGotMdANhvv/1W4jqq+uxSC50KpKamokqVKsX+PCoqComJiXB3d5e/Z2pqis6dO+Off/7RRIg6z9/fH9WqVUODBg0wfvx4JCUlCR2STqDPbuldvXoVlpaWaNu2rfy9du3awdLS8oPnij6/BWVnZyMkJKTA5w4A3N3diz2XV69eLbR+r169cP36deTk5KgtVm1UlvMr4+LiAltbW3Tv3h0XL15UZ5h6Q1WfXUrolPT48WNs3boVkyZNKnadxMREAICNjU2B921sbOQ/I2Xn4eGBQ4cO4cKFC9iwYQOCg4PRrVs3ZGVlCR2a1qPPbuklJiaiWrVqhd6vVq1aieeKPr+FJScnIzc3V6HPXWJiYpHrv3v3DsnJyWqLVRuV5fza2trihx9+wK+//orjx4+jYcOG6N69OwIDAzURsk5T1WeXErr3li5dWmSHz/yP69evF9gmPj4evXv3xsCBAzFu3LgPHkMikRR4zRgr9J4uKsu5VcTgwYPRt29fNGnSBP369cOff/6J8PBwnD59WoW/hXip+/wC+vvZBRQ7v0Wdkw+dK33//JZE0c9dUesX9T7hFDm/DRs2xPjx49GyZUu4ublh+/bt6Nu3L9avX6+JUHWeKj67RiqNSItNmzYNQ4YMKXGd2rVry5/Hx8eja9eucHNzww8//FDidtWrVwfAs3BbW1v5+0lJSYWycl2k6LlVlq2tLWrVqoVHjx6pbJ9ips7zq++fXaD05/f27dt49uxZoZ89f/5coXOlb5/folhZWcHQ0LBQa1FJn7vq1asXub6RkRGqVq2qtli1UVnOb1HatWuHgwcPqjo8vaOqzy4ldO9ZWVnBysqqVOvGxcWha9eucHV1ha+vLwwMSm7odHR0RPXq1XH27Fm4uLgA4H0YAgICsHbtWqVjFztFzq0qpKSkICYmpkACosvUeX71/bMLlP78urm5ITU1FdeuXUObNm0AAP/++y9SU1PRvn37Uh9P3z6/RTExMYGrqyvOnj2LAQMGyN8/e/YsPvnkkyK3cXNzwx9//FHgvTNnzqBVq1YwNjZWa7zapizntyihoaF6/TlVFZV9dhUaQkFYXFwcq1evHuvWrRuLjY1lCQkJ8kd+DRs2ZMePH5e/XrNmDbO0tGTHjx9nd+7cYUOHDmW2trYsLS1N07+CqD19+pSFhoayZcuWsQoVKrDQ0FAWGhrK0tPT5evkP7fp6elszpw57J9//mFRUVHs4sWLzM3NjdWoUYPObREUPb+M0WdXEb1792bNmjVjV69eZVevXmVNmzZlH3/8cYF16PNbOkePHmXGxsZsz549LCwsjM2cOZOZm5uzJ0+eMMYY8/T0ZCNGjJCvHxkZycqXL89mzZrFwsLC2J49e5ixsTH75ZdfhPoVRE3R87tp0yb222+/sfDwcHb37l3m6enJALBff/1VqF9BtNLT0+XfrQDYxo0bWWhoKHv69CljTH2fXUroFOTr68sAFPnIDwDz9fWVv5ZKpWzJkiWsevXqzNTUlHXq1InduXNHw9GL36hRo4o8txcvXpSvk//cvnnzhrm7uzNra2tmbGzMHBwc2KhRo1h0dLQwv4DIKXp+GaPPriJSUlLY8OHDWcWKFVnFihXZ8OHDC03zQJ/f0vPx8WG1atViJiYmrGXLliwgIED+s1GjRrHOnTsXWN/f35+5uLgwExMTVrt2bbZjxw4NR6xdFDm/a9euZXXr1mVmZmascuXKrEOHDuz06dMCRC1+sile/vsYNWoUY0x9n10JY+973hFCCCGEEK1Eo1wJIYQQQrQcJXSEEEIIIVqOEjpCCCGEEC1HCR0hhBBCiJajhI4QQgghRMtRQkcIIYQQouUooSOEEEII0XKU0BFCCCGEaDlK6AghhBBCtJyR0AEQQoguS09Px48//oiwsDBkZ2fD2dkZc+bMoYLxhBCVotJfhBCiJtevX8fWrVsxd+5cNG3aFIwxDBw4EObm5ti/f7/Q4RFCdAjdciWEEDW4c+cOdu/ejT179qBp06YAAIlEgu7du+PgwYNIT08XOEJCiC6hW66EEKJiWVlZWLFiBQ4cOAAjo4Jfs7GxsZBKpcjOzhYoOkKILqJbroQQomKbNm2Cq6srOnXqVOhnzZo1g5GREW7cuCFAZIQQXUW3XAkhRMWCg4OLTOb8/PwQFhaGHTt2CBAVIUSX0S1XQghRoVu3bqFly5YAgOjoaKxatQpSqRRZWVm4ePEifv/9d7Rt21bgKAkhuoZa6AghRIX8/f3RrVs3AMDr168h69ViZGSE58+fU985QohaUB86QghRoTFjxmD37t0wNDQs9LMhQ4YgKCgIT5480XxghBCdRi10hBCiQlKptMhkDgDq1KmDp0+fIjIyUsNREUJ0HSV0hBCiImFhYWjUqFGxP09JSQEAZGRkaCokQoieoISOEEJUJDAwEB06dCj25zdv3oSRkRFq166tuaAIIXqBEjpCCFGRwMBA1K1bt8ifJScnIyQkBN26dYOFhYWGIyOE6DpK6AghREUuXbqEly9fFvmz3bt3Izc3F97e3hqOihCiDyihI4QQFYiMjESjRo1w/PjxQj+LjY3FmjVrMGfOnCInHCaEEGVRQkcIISoQGBiIefPm4cWLF/D19ZXPPxcWFobevXtjwoQJWLduncBREkJ0Fc1DRwghKjBx4kRs3LgR5cuXx+HDh3H69GkYGhqiQoUKGDduHFxdXYUOkRCiwyihI4QQFRgxYgR++uknocMghOgpuuVKCCFKio+PR40aNYQOgxCixyihI4QQJQUEBKBz585Ch0EI0WOU0BFCiJKeP3+Ojz76SOgwCCF6jPrQEUIIIYRoOWqhI4QQQgjRcpTQEUIIIYRoOUroCCGEEEK0HCV0hBBCCCFajhI6QgghhBAtRwkdIYQQQoiWo4SOEEIIIUTLUUJHCCGEEKLlKKEjhBBCCNFylNARQgghhGi5/wONQuO/vuhJwgAAAABJRU5ErkJggg==\n", + "text/plain": [ + "Graphics object consisting of 5 graphics primitives" + ] + }, + "execution_count": 37, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "alpha = stability.Tilt().alpha\n", + "beta = stability.Tilt().beta\n", + "\n", + "mu = stability.Mumford().slope\n", + "\n", + "pr = 5\n", + "\n", + "coords_range = (beta, -5, 5), (alpha, 0, 5)\n", + "\n", + "p = (\n", + " implicit_plot(stability.Tilt().degree(u), *coords_range , rgbcolor = \"red\")\n", + " + implicit_plot(stability.Tilt().degree(v), *coords_range )\n", + " + line([(mu(v),0),(mu(v),5)], linestyle = \"dashed\", legend_label = r\"$(3,2\\ell,-4\\ell^2/2)$\")\n", + " + line([(mu(u),0),(mu(u),5)], rgbcolor = \"red\", linestyle = \"dashed\", legend_label = r\"$(1,0,0)$\")\n", + " + implicit_plot(stability.Tilt().wall_eqn(u,v)/alpha, *coords_range , rgbcolor = \"black\")\n", + ")\n", + "p.xmax(1)\n", + "p.xmin(-2)\n", + "p.ymax(1.5)\n", + "p.axes_labels([r\"$\\beta$\", r\"$\\alpha$\"])\n", + "p" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "2887e82d", + "metadata": {}, "outputs": [], "source": [] } -- GitLab