Skip to content
Snippets Groups Projects
Figures-notebook.ipynb 119 KiB
Newer Older
{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
Luke Naylor's avatar
Luke Naylor committed
   "id": "51c1d049",
   "metadata": {},
   "outputs": [],
   "source": [
    "from pseudowalls import *"
   ]
  },
  {
   "cell_type": "code",
Luke Naylor's avatar
Luke Naylor committed
   "execution_count": 2,
   "id": "cbd6285e",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/latex": [
       "Chern Character:\\begin{array}{l} \\mathrm{ch}_{0} = 1 \\\\ \\mathrm{ch}_{1} = 0 \\ell^{1} \\end{array}"
      ],
      "text/plain": [
       "<pseudowalls.chern_character.Chern_Char object at 0x7ff55f6a1c30>"
Luke Naylor's avatar
Luke Naylor committed
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "O = Chern_Char(1,0)\n",
    "O"
   ]
  },
  {
   "cell_type": "code",
Luke Naylor's avatar
Luke Naylor committed
   "execution_count": 3,
   "id": "32366efd",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/latex": [
       "Chern Character:\\begin{array}{l} \\mathrm{ch}_{0} = 1 \\\\ \\mathrm{ch}_{1} = 1 \\ell^{1} \\end{array}"
      ],
      "text/plain": [
       "<pseudowalls.chern_character.Chern_Char object at 0x7ff55f6a1840>"
Luke Naylor's avatar
Luke Naylor committed
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "O1 = exponential_chern(1,1)\n",
    "O1"
   ]
  },
  {
   "cell_type": "code",
Luke Naylor's avatar
Luke Naylor committed
   "execution_count": 4,
   "id": "54385028",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/latex": [
       "Chern Character:\\begin{array}{l} \\mathrm{ch}_{0} = 1 \\\\ \\mathrm{ch}_{1} = -1 \\ell^{1} \\end{array}"
      ],
      "text/plain": [
       "<pseudowalls.chern_character.Chern_Char object at 0x7ff55f6a2680>"
Luke Naylor's avatar
Luke Naylor committed
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "O1inv = exponential_chern(-1,1)\n",
    "O1inv"
   ]
  },
  {
   "cell_type": "code",
Luke Naylor's avatar
Luke Naylor committed
   "execution_count": 5,
   "id": "f0768fcd",
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "text/latex": [
       "Chern Character:\\begin{array}{l} \\mathrm{ch}_{0} = 1 \\\\ \\mathrm{ch}_{1} = 2 \\ell^{1} \\end{array}"
      ],
      "text/plain": [
       "<pseudowalls.chern_character.Chern_Char object at 0x7ff55f6a2a10>"
Luke Naylor's avatar
Luke Naylor committed
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "O2 = O1 * O1\n",
    "O2"
   ]
  },
  {
   "cell_type": "code",
Luke Naylor's avatar
Luke Naylor committed
   "execution_count": 6,
   "id": "d5c5f903",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/latex": [
       "Chern Character:\\begin{array}{l} \\mathrm{ch}_{0} = 0 \\\\ \\mathrm{ch}_{1} = 1 \\ell^{1} \\end{array}"
      ],
      "text/plain": [
       "<pseudowalls.chern_character.Chern_Char object at 0x7ff55f6a2d10>"
Luke Naylor's avatar
Luke Naylor committed
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "O_x = Chern_Char(0,1)\n",
    "O_x"
   ]
  },
  {
   "cell_type": "code",
Luke Naylor's avatar
Luke Naylor committed
   "execution_count": 7,
   "id": "f5e93bd3",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnYAAAD4CAYAAABhcat4AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAA9hAAAPYQGoP6dpAABOTUlEQVR4nO3dd3hb5fXA8e+VrvbwiLM3CWHPhBlWoEAptKWbQkPhB6W0zDLLptAWCimEVSC00Ja2UEah7EJDEqA0lLBSZkMgw0mc2PHQHle6vz9ey7Zs2bFlyxo+n+fxY1tXunpt2dLRe973HM00TRMhhBBCCFH2LMUegBBCCCGEGBoS2AkhykI8HufPf/4z3//+95k9ezZjxozB7XZjs9nQdT3rY//99ycUChV7yEIIMew0ScUKIUrde++9x/z58/nKV77CXnvtxaRJk3A6nVgsFubPn89OO+3EFVdcAYDL5WLGjBlFHrEQQhSHBHZCiJIWCoWYNWsWJ598MjfeeGPWsXvvvZezzjqLtWvXMnHixCKNUAghSocEdkKIkvb+++/zn//8h1NPPRVN0zouN02T7bffnlmzZvHcc88VcYRCCFE69GIPQAgh+rLrrruy66679rh88eLFrF69mgsuuKAIoxJCiNIkmyeEEGXp73//OwBHHnlkkUcihBClQ1KxQoiyNHPmTGKxGPX19cUeihBClAyZsRNClJ36+npWr17NoYceWuyhCCFESZHATghRdpYuXQrAwQcfXNyBCCFEiZHATghRdpYsWQLAUUcdVeSRCCFEaZE1dkKIsmKaJpMmTWLMmDG88847xR6OEEKUFJmxE0KUleXLl7Nx40Z+/OMfF3soQghRciSwE0KUlbvvvptp06Zx8sknZ12+detWrr32WuLxeJFGJoQQxSeBnRCibCxfvpwnn3ySxx9/HIfD0XF5OBzmC1/4Aps3b866XAghRhpZYyeEKFlvvfUWP/7xj3n77bcxTZNUKgWA1WrtuI6macycOZMLL7yQ0047LavtmBBCjDQS2AkhStbmzZtpbGzs9bimaYwePZoxY8YM46iEEKJ0SWAnhBBCCFEhZI2dEEIIIUSFkMBOCCGEEKJCSGAnhBBCCFEhJLATQgghhKgQEtgJIYQQQlQICeyEEEIIISqEBHZCCCGEEBVCAjshhBBCiAohgZ0QQgghRIWQwE4IIYQQokJIYCeEEEIIUSEksBNCCCGEqBAS2AkhypJpmgQCAUzTLPZQhBCiZOj9v+pbhRuFEEIMUDgcYt68w1i2bCler7fYwxFCiHazi3rvAwjshBAit7fffps//vGPfPTRxzQ1NbFgwQLmzTus1+uvWLGCH/7wzB6XP/bYY0yfPq1QwxRCiIongZ0QYtCi0SizZs3iK1/5ChdffEm/b/e3vz2Ox+Pp+L6mpqYQwxNCiBFDAjshxKDNnTuXuXPnDvh2tbW1+Hy+fl03kUiQSCQ6vg+HwwO+PyGEqHQS2AkhiubEE08kkUgwffp0TjvtdPbZZ06v133ggQdYtOi+YRydEEKUHwnshBDDrq6ujiuvvIIdd9yJZDLBs88+x49+dCaLFi1i7733znmbU089lZNOOqnj+3A4zJe+dOxwDVkIIcqCBHZCiGE3bdo0pk2b1vH97rvvzubNDTz44IO9BnZ2ux273T5MIxRCiPIkdeyEECVht912Y9269cUehhBClDUJ7IQQJeGTTz6hrq6u2MMQQoiyJqlYIcSgRSIR1q/vnG3buHEDn3zyCX5/FePHj+OOO+6ksXEL1113HQB/+ctfGD9+AjNmbEcymeS5555n8eKXufnmm4r1IwghREWQwE4IMWgffvhhVsHhW265FYDjjjuOn/3sWpqammhoaOg4nkwmWbhwIY2NW3A4HGy33Qxuu20hBx100LCPXQghKolm9rvRorQUE0KUjlAoxKGHSksxIUSpKW5LMVljJ4QQQghRISSwE0IIIYSoEBLYCSGEEEJUCNk8IYQoK4888giPPPII6XS62EMRQoiSI5snhBBlSTZPCCFKk2yeEEIIIYQQQ0ACOyGEEEKICiGBnRBCCCFEhZDATgghhBCiQkhgJ4QQQghRISSwE0IIIYSoEBLYCSGEEEJUCAnshBBCCCEqhAR2QgghhBAVQgI7IYQQQogKIb1ihRBlRXrFCiFE76RXrBCiLEmvWCFEaZJesUIIIYQQYghIYCeEEEIIUSEksBNCCCGEqBAS2AkhhBBCVAgJ7IQQQgghKoQEdkIIIYQQFUICOyGEEEKICiGBnRBCCCFEhZDATgghhBCiQkhLMSFEWZGWYkII0TtpKSaEKEvSUkwIUZqkpZgQQgghhBgCEtgJIYQQQlQICeyEEEIIISqEBHZCiLJkteocfvg8rFZrsYcihBAlQzZPCCGEEEIMGdk8IYQQQgghhoAEdkIIIYQQFUICOyGEEEKICiGBnRBCCCFEhZDATgghhBCiQkhgJ4QQQghRIfRiD0AIIXL68EP429/gJz8BhwNeegmWLoWmrTB+HHznBNht12KPUgghSorM2A1GWxukUtu+XigEyWThxyMGLhyGRGLb14vGIBot/HgEpNNw332wYAGccgps2QKnnw4rV0LdaNiyGZ5/AU4/Dd54o9ijFUKIkiKBXVfr6+Hdd8E04Z134I9/hKefgVis53U3NcBf/wqWfvwKTRNuvx3i8SEfsuimsRH+8x/19YcfwoN/gieeUMF1d62tcP8DoPXjMbRocNddEMxxHjF0DAOuvRbWrYN7F6nH7cwz4cKL4NJL4eKL4LHHYL99wUjBLbeo/y8hhBCABHadnngCHn4Ixo6D666DFStg9Wq47mfqhaXrzFzSgDvvhJNPBk1Tl7W0qBf+mxfABRfAs892vuD4fPDVr6rbiMJ5+WW4+26YMgVuuAFeeQXWrYVf/Qq+//3sGTfThFtvhZPng63LioRoTM0GLViQfW6HA+bPh18vkECiUEwTrv+5+l+69lowkvDTn8J3vpOdcnU44MorVbD96WrYsKFoQxZCiFIja+wAfvc7NfN21llw3nlwxRUwbZo6NnMmLLwN/vUvOOQQddljj8Jhh4HTqb5PpVRQd8EF4HZDQwOceKJ6gfre9zrP879V8L//waxZw/0TVr4nnoTPVsNll8GFF6rU3e67q2N77QVXXQ3PPw9f/7q67J//hB13gqqqznMsXgzvvw+ffZZ7JnbsWJW2Xb4cDjig4D/SiPO3v6k1dE/8DaxW9UYoHFb/S91NmAA77QQffAj19TBp0rAPV/RixQrYbTd4800183rQQerNVsayZbDjjur/qb/SafUmOvNG+tHH4KtfAbt9aMcuRAWQGbuXXlIpu1NOUamfPfboDOoAvvkt0HX46CP1fToNzzyrAruMjz9WT1Yff6y+HzcOjjwSHnww+76+eHTPy8TgrVgBzzwN558PD/8Vams7gzqAL34Raqo7H0NQ1zvu2OzzHHGECuxnzuz9vr74RXkMC6GtTS1XOOlE9fg1NalA++qrO99AdTd5svqcTg/fOEXvUqnON8mLFoHDqd7cfv/76vEF+Ohjtd64v0Hd1q0qm3LqqepcGUcdCfffP/Q/gxAVYGQHdskk3Hwz/PjHav1bLArfPyX7Oi4n1NV1vnh89BFMnJCdvqurg/Hjs9dq+f3qyaxr2u7ww+H1f6tUrhgapgk33qjS5RYLbKiHs8/Ovo7FAhMndT6GDQ1gplWKfKAOOEDNumZeqAQAb7/9Nueffz5HH/1FZs+ew5IlS7d5m8bGps5vFi+GUBi+/BX1fU0NPPMMHHxw7yfIPJ61tfkPXAydO+5Qj9c/XoRjjoF95qgsyISJ6s0XwNNPqefB/vL74dhje25Sq6pSKfitW4du/EJUiJEd2P373+oJYsYMNStw2WXg8/a8XjSiUj8AK//bc0Zn7Fi10WKvPTsve/ddmDu3M3UA6r6mT4O1a4b25xjJPvxQvZOfPVv9rn/6UxVodxcJw8SJ6uv/5ngM+8tmg9lzsmf/BNFolFmzZnHppZf06/obNmygpaW584IXXoBJE1UZE1CpWH0bK0W2bFFvsLbbLs9RiyGz6lMIBlW2Y889s/+/pk9TSxgiEbU5pj8bzjJstt7TrYcfDq+8OohBC1GZRvYau8WLt/2isHEjtAVg3/3U9/X1MHNG37dZvhzWroU//KHnsZoaWL8+/8BCZHv5ZZg2ve8Xi0hEPR77ZR7DDTB6dP73Wdv+GO6/f/7nqDBz585l7ty5/b7+448/zqmnnqq+aWmFt9+Gffft/x3G4yqoP/DAni/8sZja+JSh67DDDmrWJ7NcAmDUKLVsQgze88+p/y+7HY75Yvaxxia1rnj16s43V0Nh553hgQeG7nxCVIiRHdht2qQCrb784x8w77DOmYRwSO3K683mzXDPPWq93vjxPY/7fCrlJIbGpk3bvs7ixbDLLuqFANRj2HXTxED5fGpRv8jbypUrsWcCss0NYKJmZ/rr1VchkYRvfKPnsVgMPvtc/R8aSbjkEhXYGQb87Qn49+vwrW/DQXMlsBsqq1apTUrdtbZCOqXeQL/2Gng8nceSSbXWlRy7zEePVutZ++L1qvMLIbKM7MDO44FPPlHrtLqmTDMCAfj73+E3v+m8zG5XM0C5BAJqAfjNN6snpqYmNSvQ9dyx2LZTTKL/PB5Y/alat2jL8XuNx9Vmh+uu63wc+noM+yMWU2t/RN62bt2KmVl/minevXlz/25smvDQwzBnjpqx6666Gr58nKqBd/ttsPfsztu5nPCXv6jriKEzapQqIN19TeTvf68Ca01Ty126lhyy2WD+9/K/z3i87zfZQoxQI3uN3YEHqqLEr7zS81g8Dr/4paqXlVlfBzBmjArguovFVLX8iy/uTPM98WTPgLEtAGPHDNmPMOIdeKAqGvzU33seS6XgppvgRz9S5RUyensM+6stoM4hBiWR6fiRmdn+9FOVttuWF/4Baz6Ha67O/YYs47hjwWJVfxtJA+6+B047TYK6Qvjq8fDQQ/DqayqAbm2FP/0JjvgCbL+9us7UqWoZxECl0+r5tfubsQ0bs8uoCCGAkT5j98Vj4Pd/UMHbxRerxbgWi1rv89I/VS207buthdthR3jyiezLDAMuv1w9+fziF+qJzTDA2u3XmzRg1f86n+jE4M2dqx6jm25Wwfixx4LdoerRPfsMfO1rqoRNVzvsqILuXGIxtdtuS6NaDO71ZgcPpgkffgD/d2rBfqSRYNSoUTQ0NOD3+9Vmlz12h/dWwh9+Dxdd1PsNP/4Y7rgdFi7MfsOVi8+n1ns99hgEgnDCCdteeiHys9eecMON6rnxhedVuaFvfANcrs7rjBkD9evVGy6rtX/nbW5Wu6OPPhpefFGt49tlF3Xsny9tO10rxAg0sgM7nxd+ex/89ndqq/6dd6pFvsccA1ddmTtlOntvdSwY7CyXsb4+d8mFPfbM/v7N/8DkKfmV2RC52e0qVX7//Wpn8v33q13ORx6pdjnnStXM2l6t6+pe2PZf/4IlS1Qwt9OOcNttqtDqV7/aeZ1Vq9TaLpkpGJTdd9+dFStWMCtTrPv441Vg9/DDMH26KiTdPaB+6SX1GN92e883XL351rdUED99Wuc6WVEYBx+kPvpy9NGq9ElmI9O21NaqDj/dpVLQtFUVqRZCZNFMs7/9kd4q7EjKyS9+Abvumv2C3x9XX62e0I49dtvXFYV1770qcD/ttIHd7vbbVSov14vNCBaJRFi/fj0AJ554Ehdc8BPmzJmD31/F+PHjuOOOO2ls3MJ1110HqHInN910E7fddps6gWnClVepsiegZvAOPkT9rhsa1ML7vfdS9Qq7LsDfloceVjN2kyaqQF0UVzoNf/gjnPjdwa2Pe+YZ9aZr6tShG5sQQ2Z2Ue99ZK+xy9cPfpB7XV5fIhFVOuXoowszJjEwJ56kZg4G0rXAMNQC8W99q3DjKlMffvghJ554EieeeBIAt9xyKyeeeBL33HMPAE1NTTQ0NHRcf+LEiZx88sk0NjaqCzQNrr8Obluoyshs2KgCshdeAKcDFt6qWsUNJKh76ilVVugHp6vZ2Pr6ofpxRb4sFvjeSWpjy2DMmSNBnRC9kBm7fL38slpE/7Xjt31d04QFv4Zvf0s9Ga1dK09KpeDtt1Wgdsop/bv+PffA3IOyG9KL0vTSP8Fug0MPVbtuv/QlNVN+/vnFHpkQouLJjF15OvxwmDC+Z6ubXBo2w9e/poK5QEA1mQf1rrW/cbUYenvvrdZB9qf0SSCgNmpIUFf6liyBeEwFdaDKanzta2oGLxYr7tiEEKLAZMZuuGV+3ZoGF1ygeh2ecYYq29FX6QYhRN+efhr++U+Vdr3+550dED5fA/ctUj1Mv/AFtZNy3mFFHKgQorIVd8ZOArtiWr4cFi1SuwF33UUCPCEGo6Ghcwbd6VRFc0F1CenaocDlyr2LXQghhoQEdiObacIbb6gA7+OP4dlnpdZWqeqtQ4kQQpH/ESEodmAna+yKTdPULsDf/Q7++lcV1EUicO65KqUka/BKQzCkulhEZY2WEDkZBtx6q2qlKIQoGgnsSoWmweTJ6uvmZrWx4tzz1I5NCfCKr349PPU0XHiBBHdCdGcYqk7nX/+qingLIYpGArtSNGmSmsG7607Veufc8+CXvyz2qEa2nXZSDeXfWwkv/qPYoxGitLzxBix+GW64AQ44oNijEWJEkzV2pS6zBs/lUj1PP/pIzejJJoviqK+HiRPV717WE4mRruv/QPcWfUKMWLLGTvQlswYv08j++eclRVtMkyapx+Qf/1DrIKUumhipDAOuuAIeekh9L0GdECVBArty85OfZKdoTzkFPv+82KMaeerq4O13VC1CCe6KItncDC+9RFIW6w8/w4CrroLFi2HMmGKPRgjRhQR25abrLtq77gSvr7M8yvp6mcEbLrNnd665k+CuKJKbNvHRTy8juXlzsYcysmSCupdfVmt/jzii2CMSQnQhgV25ygR4d90J1dVqF+33vicp2uHUNbi75ZZij2bESU+dyvfaP4thdP/9EtQJUcJk80Sl6Fro+L2VsOuu8MMz1A41WeBfWO++CxMnwei6Yo9kRAmFQhx66GEsW7YUr9db7OGMHMEQfPC+emMphMhBNk+IodA9RWvR4A9/7Ny9KTN4hbPnniqoa2pS5R6kzt2wsKxaxb/bP4sCMwxYsADWrwefV4I6IUqYBHaVJhPg3X+/eiIGWP6GpGiHw4YN8MyzUsR4uJgm9vbPooAyxYcffVQ2aglRBiSwq1Sapt5ZA3g92btoJcArjD326FxzJ8FdwTzyyCN885vf5Kqrri72UCpfJqhbvFjNRh9ySLFHJITYBlljN1J0X4N35ZXwteOLParK9NZbKojea0+44w5Z41ggkbffxv2DM4jctwj33nsXeziV6dprVe3MG26Aww8v9miEKBPFXWOnF/XexfDJpGj3208FeLvtpi5/7DEYP146WQylzG7ZLVvkdyrK21FHqVk6CeqEKBuSih1pMgGex6Nm8ZYskRRtIcyeDccco36fTz8jde4KID11Kt9Cyp0MOcOAJ5+EdFq94ZOgToiyIoHdSKZpcOed3TpZnCprw4bSxo1w441SxLgQHA4+a/8shkhmTd0NN8DHHxd7NEKIPEhgN9J1L5MyZw64nGqm6c0VMoM3WBMnSoeKAkmtqedKIPTOp7Stb6NtfRuB+gDBhhCRpgjR1hiJcAIjbpBOpYs93NLXdaPEL38JO+9c7BEJIfIga+yEkgnwMvWpVqyAM38Eu+4CZ5wha/AGI7Pm7tzzVHB3yy3gdBZ7VGUvsbGJrwGNkVDHZelUmnTcIG6kMdPd3pRoGlbdgtVhRXfq6E4dm8uGRZf3tz2COukoIUTZksBO5DZnjprBW7RIBSS77gJnnw377FPskZWnTHD30EPFHknF8Yx2455c1eNyM22SMlKYRpfPyRRGzCDaHO0I/Cy6pSPQ0506NrcNi3UEBns2uwR1QlQACexEbt130S5aBJ9+qgK7eBzsdpnBG6jZs9UHwOdrYNw4lfYWBaFZNHS7DnawYetx3EgYGLH2j2jXYE/D5tFx+BzYvXasNuvwD364GIb6W9x+Jvzs2mKPRggxBCSwE33rGuCl29cp3XorfPSRpGjzlUyq2c+pU+DXt0hwVyS6XVeBn7/zMiNukAwniQfjhBrCQAjdqWP32rF7HdhcFfSUmUm//ut1eOZp8PmKPSIhxBAYgfkGkRdNUztnQZU/6LqL9vXXZZPFQNhscN3PpEPFIKX91TwAmDW1Q3ZO3aHjqnVRPbWaUbNq8U3wYbVbiTZHaV3TwtZVWwltCZFKpobsPoui65q6a66WoE6ICiKBnRi4ffft3EVr0eCcc2HdumKPqrxk1txJcJe39KhR3AmYo+sKcn6L1YKzyol/op9Rs0ZRNaUKu9dBrCVG86fNtK1vIxFKFOS+C6p7mzCpUydERamgvIIYVl1TtKtWwdSpkEqpF4wvfUlStP2RCe4uvwLWr4NZs4o9orLwyCOP8Mgjj6DHE6pxTyQCXm9B71PTNOweO3aPHe9YD7G2OLGWKG3r27DarDhrnDirnOWxw7axEd59V4I6ISqU9IoVQ6exCS69RM1CSZmU/ovFVPmTpAEpQ0qh9FPji68x+rLzi9orNhlJEm2JEg+omTtnlQNXnUut3Ss1hqH+xlzOzr85IUQBFLdXbBm8vRRlY3RdZ4o2swbvp5cVe1SlL/MCe/11UsS4zNjcNpWq3b4Wz2g3iXCCltWtBBtCpI0SKoqcSb9eeIFaDytBnRAVSwI7MbS6drL4zV1w3LHq8g0bpBfttnz1q9KhokxZdAvuOje1M1SAF2+LsfXTZsKN4eJ3vei6pu4b35AZdCEqnAR2ojA0Ta2/O/hg9f3zz7fvoj1FArzedN1QIcFdWdIsWkeA56pxEmmK0ry6WdXIK8bfvHSUEGLEkcBODI/TTstO0Z5yijQZz6VrcPf0M8UeTWmzWtkMoJfeejaLbsE71kvtzBocPgehzSGaP20mHowP70BefU2COiFGGNk8IYaXaapOFr/9LVxzDUyerHbp1dVJiqirtWthyhT1OzFN+d3ksGnVJo474cssW7YUb4F3xQ6WETcIbwmTCCVwVjnxjPUUtm1Z17+ZNWtg2rTC3ZcQopvibp4ovbe6orJl1uDtv7/6PmnAqafCqFGyi7arqVPV56VL4W9/g1/dNOI6VKSNNKlkmrSRIm2YmOk0Ztps/4BwUxiA0JYQhAFNlSWx6FYsetfPFrQi/03pDp2qyVVEW2OEN4dIhBN4x3lx+BxDf2eGAddcq/o9f+14CeqEGGEksBPFpVvhyitVL9pzz5MyKd35fPDW22o3Y4W2H0slU509W2MG6WSaVDLV3re1k2bRsj70NWt4DjA/WU1i1k4q4DNNzFT3JISG1W5Bd+hYHVZ0Z/vnIpQkcVU7sXtshBpCBOoDQz57ZyYNjEsvR//XMjh8HvIfJMTII6lYURoyKdpFi6C1FR59tLOF2Uj31lsq6N1j94oI7oyE6seaCCVUINdeFkSzairosulY7RYsNgtWmxWrzZqz8G9vdexM0yRtpEkbZvtsXxojliKVMDDiRkfgp1k17G47No8Nm9uG7hjeQC8ze6dZtCGZvTOTBsbFl6G//grBcy6DeYfjm+gr+mylECOPpGKFyO5k0dKigrpPP4Wf/xx+8IORPYOX2VBx7nlw801ql2MZMU2TRChBIpQkGU6091nV0F06jioHNpetPaAbmkBe07T2gBByPcWlkilS8RTJaJJEOEm8IQyYaFYNh9eB3WfH7rUXPCDqPnvnqnXhGePJ635N0yR+y504Xn8F45rrcBxyGIH6IGwM4psgwZ0QI4kEdqK0aBrUtjd1NwywWCRFCyq4u/NOGDe22CPpt2QkSawtTjwYw0yZWG1WbB47Hq+aIcs3/ZhpKTYlHOGWPG6fmQW0e+14RoOZNklGkiQiCRLBBLG2GJpFw+ErfJBntVnV2rvmKKHNIYy4gX+if0C/G9M0CW4Mkph3HNY5e2M74hAAfBNNghuCaBYN33hfQcYvhCg9kooVpS2Tor33Xlj5XzjnbFUqZSRra4NF96nfRYl1EEglU8RaY8Tb4qSSKSy6BUeVA4ffic05tO8jC9VSLBkzSATjxANxUonOn8FZ7SzourxEKEFgQwCLbsE/yd+v1LCZNEjc+GtC876Md48ZPdK50dYYoU1BPKM9uOvchRq6ECKLpGKF6F3XFO0bb8D06eryf/4TXK6ROYO3bh08+SR8/hnccktJBHfJmEGsOUqsLQ6a6pnq9Xuxe+x5nc9IqE0Uaq2c2h1rGmnS6TSYKt5vNb1cAVyT8hFf04rFqqmZLoumvtatWO0WrHbrgNK8NqeOzanjGe1RP1drjFhLjOjWKHavHVeNC7s3v5+rL3avnerp1QTWB2hd04p/or/P+8msqbO//gr+fffB5tu5x3Vc1U7SyRThxjBWh7Uwu3CFECVFZuxEebrkElj88shN0XbdUFHE4C4RShDZGiEZSWLRLbhqXTirnf1OJZqmiRFr3xUbV7tiU3Eja0esZlElSyy6RZ23vazJ5rWb+c7ZJ/DsI8/gcXlIp9KkUyakTfV1116tmobusGK1W9WaPpeO7tT7nWI10yaxtjixlihG3EB32vCMdhckwEun0iq1GkrgGePBParnTFvXjRLG1T/DdtwX+zxn2/o2EuEk1dOqh3zmVAjRXXFn7CSwE+Wp6y7a91aqAO/WWzvX540EmeBuzz3hjtvVesRhkowmCW8Jk4wk0Z06rloXDr+jX4FSKpnq2EiRCCc6gjirXZUi0R0q6MrsitUsuc+5efn7PHfWKXz7scfwTJ/W47hpmqQSqayPTBAJqoCvza1jd9ux++z93hWbCCUIN0YwYsmCBXimaRJuDBPdGsVZ7cQ7ztvxuzVNk+TFl2N75eV+BXWgAtPWNa2YaZPqadU5dxkLIYaKpGKFGLiuKdr//AdefBFqatSxDz6AnXeu/Bm8zG7ZdeuHLagz4gaRxgjxYByr3Yp/kr9f6T0jbhBri5EIJkglUgDoThuuWhd2j13NnvUSwPXGEmjlVCDS0gw5AjtN01SQ2C1gU7OEquRKMpok3BRRqUq7FYffgd3n6HNWy+5VGyoyAV7b+jZ0lw3vOO+QzYZpmoZ3jBerXSe0KYiZNvFNUBsgghuDmHvsj3bYof0K6kDNevon+2n5vIXAhgBVU6pkp6wQFUpm7ERl+fRT+M4JIzNF+49/wCGHFqTOnZlun0FqjmHRNdyjPTir+p6hSxtp4sE4sdY4Riypdpr6Hdg8Nuwe+6CL8g7V5gkzbZIIqx2x8WAcM6128GY2TGxrfV48GCe8JUwqkcJZ48Iz2j2k7cJigRjBDSFsTg37a0sIzzm03wF1d8lIkta1bTirnfjGl3YbNiHKl8zYCTF0ZsyAu+7M7mRx5o/ggP2LPbLCamiA666HPf4+5EWME6EEoYYQqWQaz2g3rlpXn7NryahBtDlCPJAAwO614Zvow+F1DHhWbjhkSps4fA68ppdkOEk8ECeyNUqkKYrDb8dV68LmsuW8vcPnwO61E22OEm6MEA/E8Izx4qoemsfA6XdCKoVx8eW4Vi7Heu+OOHyj8zqXza1mFkMNQXSnFVeNa0jGKIQoHbLQQlSWTIr2d79TAZ7VCm8sV8eShlqbV4nGjVNp2fdWqvZj0digT5k20gQ2BGhb34bFZqFmRjXuOnevwdlDDz7MscccxyHzDuaM887g082rGLV9LVWTq3D6nVm3W7FiBbNnz+nx8fnnawY97sHQNA27145vgo9R29fiGePGiBq0rmml5fMWYoHcv1dN03CPclM7owa7x05oU5DWda3txZgHx0waWH92LZ6Vy2n5wcXER0+h34mWHFw1TpzVLkINYZLR5KDHJ4QoLTJjJypT1zV4hqEu+9ODsHRp5aZou3aoGGRv2UQ4odZymSbe8b4+Z58SoQTPPvksC+9cyLmnn8Oc/efw1AtPc/GVF/Poo48xfvy4Xm/7t789jsfj6fi+JrNOsh/SXi9PAEf5/f2+zUBYrBbco9y4R7mJB+NEW6IENwSJbo3irnPnTIVabVb8E/3E/XFCDSFaPm/BN96Xd5mRrN2v11yHZ4g6SnjHeTBiBsGNQWqm15TkTKoQIj8yYycqm6aBrT2Ftttuagbv3PNUkeN//avyZvAywZ3VCmZ629fvJrMbs21dG1a7lZrpNb0GdUbcoG19G23r23jsmcf58rHHcdIPTmKH3Xbg4osvYuzYcTz22GN93l9tbS11dXUdH9Yc/YHTqTRG3CDUEqJpYxON65toXN9Ik2Hj50DEoTo3RFtjJEIJkl36zw4Vh89B9ZRqqqdWg6YRqA/QsqaVZCT3jJfD56Bmeg02l41AfYBg+waIgTBNk+CGNtXj9prrsB17NA6fA99EL/FAnFBDOO+fR9M0fBO8pBJpwo35n0cIUXpkxk6MHHPmqBRtpkzKuefBb38Le+1Z7JENrdmzYe+9VVBbXw91df2qc5dKpghsCGJEk7jr3CrtmmNGKJ1KE26MEGtRGymcY1ys+ux/nP7D07Kut//++7Ny5Xt93ueJJ55IIpFg+vTp/N+p/8eeu+6patpFVV27VDLdEXz/6fE/8acn/txxWwewHRDb1AZ+LUeQrmG1WbA6rB39aHWnPqhSHza3jZpp1cSDcSJNEVrXtuLwO/CO9fY4r0W3qHZhLTHVLixq4Jvk61f3CjNpEHrnY+K+cThuvAGbv/Pxc/qdmCkGvU5Od+h4RrsJN4bVppZe1hAKIcqLBHZiZOmaol25EnbfXV1+8wI48IDKSdFqmkpBn302TJiwzSLGyWiSQH0AgOqp1djcuV/kY20qSDHTdGykaNraRCqVprZbDcFRo2rZunVrzvPU1dVxxeWXM3PaTCKBCC+8+AI/PuvH3HTFTey20+7oTh2bx47LbkXTNay6lTPPP5PTz/tBx8PT8vK/mfyzy4joAdw7zsTMFCZOpkkZKUzDJJVUdesiWyMdM2YW3YLdo2rX2T32vNKQmc0W0dYY4S0hmj9r7nXDhKvGic2tZ3eU6KMjRyb96nlnBfaHH8fhz31OI24QalBlWvLt8OEa5SIeTEhKVogKIoGdGJk0DfbYQ30dDsNHH8LDD1dWmRRdh6uuUjOTF1zQa3AXD8YJbAiiO3SqJvtzzmiljTShhhDxYByHz4FnrKdHGZDus3umafa4LG2kiYcS1FhrOHiXQ8A0sVZZ2eW8XWkONPPU0qeYd/y8nDOFNrKDzagze6ZKs2hYLap9WPfrgmpTZkSNjhp2sbYYoGH32jpq0w2k9Rioll0Or53Q5hChTUHibTG84709ZuV0h071tGoCGwO0rQvgHefJOdPWvaOEY3zvBbe9Yz2kEgaBDQGqp1Xn1cc2k5Jt+ayVcGMY71gpgSJEuZM1dkJ4PJ27aC0WFQidfU5lrL/LrLl7b6UK7mLZuzojWyME6gM4vHaqp1blDOoSoQQtn7eQCCfwTfThn+TPCoCqq6uxWi09Zueam1s6ZvGS0STBjUG2rmomtClIOmXiqXNTs10NtTNr8Y3zssfee7C+vr5ghXN1u46zyol3rJea7WqomVGDd6wHM20SagjT/GkLgfoAiXBiQOe16Bb8E/1UTa4ilUjR8llre9DY83pVk6tw1jgJNYQINoSyjg+0TZimafgn+LFYLATrg6RT+a0rzKRko83RXtcMCiHKhwR2QkBnivb++1WA98Wj1WVtbeW/ySIT3H3yP1i9uuPiYEOI8JYwrlEu/JP8PdJwpmkSbAjRtr4Nq8NKzXY1qqZaNzabjZ122only9/IuvyNN5az86xdaPm8hdb2jQaeMW5qZ9ZSM02VTunaFeKTTz6hrq5uiH/43ul21Qqtemo1o2bVqp2icYO2dW00r24m2hwdULBk99qpnVGLw2cnuDFIcFOoR1kSTdPwjfPiHecl1hJVmypME9M0Cb27CsvKd/rdJgzag8rJflLJFMGNwQH9/F25RrnQnba8NnkIIUqLpGKF6CoT4GW8+BLceGP5p2hnz4ZnngaXC1IpgutbiEVNvOO8OVOCmRp2yYiBd6wXV23fC/RPOukkrrrqanbeeSd22203Hn3oUTZtauDI/b+ARbfw5+f/THNrM9dddx0Af/nLXxg/fgIzZmxHMpnkueeeZ/Hil7n55psG9GMlYEgeD4vVgqvGhavGRTKSJNocJbQ5TLgxjKvWhXtU7/X7utIsGv6JfqLuKKGGMEYsiW9izw0TrhpV5Dm4MUQ6ngQzRcI7GvtDj+EYWz2gsesOHd8EH4H6AOHGMJ7Rnm3fqPu4u6Zkm8J4x0hKVohyJYGdEH355jdg8qTsThYXXaRKp5QblwrOktf9Aseadeg3LsgZ1CVjBsH6AOl0muqpVb1upOjqqKOOorW1lUWL7mNrUxNTJ03jxmtvZMf9dkS367S0tdLQ0NB5H8kkCxcupLFxCw6Hg+22m8Ftty3koIMO6vePY0yfwQHAsu237/dt+sPmtmFz20glUx3dJ2KtMdx1bpzVzn6lil017TNgGwK0ft6Kf1LPDRPOKiekVUcJLRHHcftCHHl2q3D4HHhGewg3hrF77P16zLrTHTruOheRpiiuGteA1xsKIUqD9IoVoj9Ms7NMyjnnwF57QWsrVFWV3Qxe8vU3sF50Aakdd0G/8zY0d2dwFw/GCW4MqkK7k/39fnFPhBKEt4Qx4gZ2rx3PGE9WmrUQNq3axHEnfJlly5bi9RZuhimVTBFpjBBri2G1WXGPcedMSeeSTrXPfIYNfBO9Wbfruqau5YyLsR59pEqJ5/n3ZJomrWvbMI00Ndvlt8PVTJts/XQrDq8D3wRfXuMQQhS3V6yssROiP7q2KttrLxXonXd+WRY6th24H6mbf4314w8wzj4PMxIFINoaI1AfwOa2UT2tul9BXTqVJrgpSNv6NrBoVE+rpmpyVcGDOgBr/Xr+BFjWri3s/dis+Cb4qJ5eg9VhJbhB/bz9aRdmsaoNEw6/XXWtaFG/6x4dJb7zZRKh5NAUHTbShLbkdx7NouGp8xBri5GMGXmPRQhRPBLYCTEQXWdTfnRm2XaysM3dvyO4i93/Z7WmbFMQZ7Uz50aKXBKhBC2ftRAPxPGO81EzrXpYi9xqyQQ7AcTjw3J/NqdO1eQqfBN9JKNJ1Ts2x+7X7jRNrbtz1rgINYQIbQkRfeyZjqAu01HCO95LrDU6qE4Qul3HO8ZDrCX/Ha7OGidWm5WIdKQQoixJKlaIweiaol27rnODQplIfvAJm0NukhGDuh3q+pV+S6fShLeEibXGsHls+Mb7irIeq/HF1xh92flE7luEe++9h/W+u9b1s3vt+Mb7+tXRItwUpvGDRqy6hbHeMPY9dso6HmmKqHpyvWxq6Y+hSMnG2mIENwb7LFYthOiNpGKFKF9dU7R//pMK6pqa1A7aMpjBM8ZPwWrT8Xz2Ps7rL8eM9j0DZSQM2ta2EWtTs3TVU/qXsq00Ft2Cf5I/a/ZuW6lLM2lgu+F6fP99HatDx5g4tcd13HXu9pm9cN4zbkORknVWOdGdet63F0IUjwR2QgwFTYNx49TXgYBq51XiKdpYQLUH8030UbPjaKzvvY1x1rm9BneJcILWNa2YaZOa6dW4avLbwTlYjzzyCN/85je54dZfFuX+u3L6ndRMr8GiW2hd00o8mDstnFlTZ3t1Cf5ptVRNrSLcGO5Yc9eVd6xHtSDbEOjXOr5chiIl6xntwYgme/2ZhBClSVKxQhRC1xTteyth/vfg/POLPaoOyUiS1rVtOPx2/BP96rJ/Lcd68YVqt+xdt6O5OgO3aIuqy2Zz6/gn5m47Ntwa3lnFrad/l6ueeRrv+PFFHYuZNgluDBIPxnHXubNqyfXWUSLYECLWEsU3wadKn3SRSqZoXdOKxWalempVXjtlTdOkdU0bYFIzvSavn6t1bStmyqR6enXBOoIIUXkkFStE5emaor3rTjjuOHX5myuKPoNnJAza6tvQXXrWmrquGyqSV17b0YEg2BAi1BDCWeOkakrutmPFYHq9/BPAV/yyHJpFwz/Jj7vOTaQpQmBDoKOjROKmhTnbhPnGeXFWOQluDPWYVbParPgn+TGiBqHNee5w1TS8Yz0YMaNfmzxy8YxR3Thk1k6I8iEFioUopO6dLF54Hp78e9E6WaSNNIF1AVWGI0fNNNvc/Un++lZCcTva+jasuqV912v+i/kLRWtt4SRAa2mBAtaxGwjPaE97SZQQ6XQbmqZhHP5l/PvOwXbkYT2u7x3vxUikaKtvo2Z6TdZ6RZvLhnech1BDCLvHhsPnGPB4bG4bdq+dSGMEh98x4Fk3m8uG7rIRbY71u3afEKK4SuOttxAjxZVXqhk8i6V9Dd6psHHjsN19cFOQdDqNf3Lv6VTbAfvi3WdnQv9rIP3rW/BW6yUX1AFYW5q5ANAaG4s9lCxOvxPfOBfccSdt73yOZ5epOYM6ULNqmQA7UB/s0VvWVePC7rUTagiRNvrft7Yr92gPqWSKWGt+s3auWidGNEkyKnXthCgHEtgJMZwyM3j3368CvNF1MGqUOrZ6dUFTtJGtERKhBL7xPXuXdhcPxHEGGvAvfwn9iks7ihiLbTOTBtafXUvV0qfxNK8n3hbvEbB1pXbYVmHEjJwFir3jvJimSWhzKK/x2Jw6zionkaZIR3p9IBw+BxbdQizHRg8hROmRwE6IYsgEeAsWgMMBjY1w4kkF20WbjBqEt0Rw1bq2mdILbQkRa43hPWL/nB0qRO+6d5Twfv0o1VFiU99Bmc2l4x3nIdYaJRbInlmz2qx4x3qJB+J5r3Vzj3aTNkwiWyMDvq2mabhqXcTa4nnPGgohho8EdkKUgro6uG1hdieL118fklOnU2mCGwLoTiueMZ4+rxttjRHdGsUzxoOr2pm1ocI4+zzMpKTjemOm0yR/emWPjhK+CV5ibTEiTX0HVX2lXZ1Vzo5j+cy6WW1WnDVOos3RvIKzzK7dfNO5QojhI4GdEKWg+y5aqxWeeUYdS6cHNYMX3hImnUrjm+jrc/F8MqpmlpxVTtyj3B2XZ4K72O770rYpv8CiEEyXm2WA6ek7WB2WsZgmwU0h4jvt1RHUZTirnLhGuQg3hrc54+YdpzaB5Eq7esd58551A/DUuTFNctbO2xaLbsHht0tgJ0QZkMBOiFLSNcC76mp12Qsv5J2iTUaSxFpjeMZ4+lxXl0qmCNQH0J063vE9d5ja5u6P4wcnk4wYhB9/bpsdKoZDatx4LgDMiROLOg4zaRB99GnigTj27xyfFdRleMd4sXvtBDcGMeK9z3r2lXa12qy4RjmJbI3mVbjYoltwVjmJtcb6XPPXG2e1k1QyRSKcGPBthRDDRwI7IUqRpkGmQPD4Cdkp2n4GeKZpEmoIoTttOKt7L1VhmiaB+iAA/km9z+rZPXaqPEncC2/os0PFsDEMqts/F0tmTZ1rwS+oSm7tc/2ib4LqqRtYHyCd6j0d2jXt2j1t6h7lRrOoWdh8OGucpI10Xmv17B47VrtVZu2EKHES2AlR6vbaMztFe+55sGTJNm8WbY5ixA2847x9pmAjTRGMmIF/UtU2+77aJ48nddMCteauSMFdpqXYzRf8kMWA5bPPhn0M0LOjhH33Hfu8vsVqwTfJRzqV3mZg5h3nxUybhLuty7NYLXhGqxm9ZHTgrcJsTh2bW9Wly4ezykk8kCiZdLwQoicJ7IQoB11TtPfcDQcfoi6///6cM3ipZIpwYwRnjQubq/cUbDJqEGmK4q7r+3pdZW2oKEJw9+1vf5vHHnuMy35y+bDeb1e9tQnbFt2u4xnjJdYaIxHqPaVptVlxj3ITa4lhJLJnJJ3VDnSHTrgx/1k7I5okGRv4TKfdbwdMSccKUcIksBOinGga7LMP2HRIGvDvf+dM0UYaI2gW8Ix293oq0zQJbQqiO62463q/Xi6Z4M5MJAmsaRpRMzimaRKsbyUdDA0oqMtw1TixeWyqWHQfKVlXrQuLrhHZkj1rp2karlEukuH8igZ31qUbeECu23WsdiuJoAR2QpQqCeyEKFc2HRYtyk7Rnn46RiROrC2Gp86Dxdr7v3ikKYIRT+Ed3/du2V7vfu7+cM89JDQnwXc/HRF17sykQeitT4jHgFtvHXBQl+Eb71Op1j5SsppFwz3aQzzYM+3q8Duw2qxE86xL56x2EmuL9RlY9sbusxMP9V10WQhRPBLYCVHOupdJOfRQ4mEDjCT6O2/2uskiGeuSgnXm3zLa7nVQNdGL62eXVnwR40z61X3pOfhH6TgG0TvVarN2pmT7SGs6q1TaNbQ5OwDMzNrFg/Eeqdr+cFY7wVQdRgbK4XVgpkwMaTEmREmSwE6ISpAJ8E4+GYffgev9N7FdcgHJ784n+c+lPQK8SGMYq80y4BRsLnafEy65ZNg7VBhTpnEIkJ45s+D31XVNXfqCi3CMqR70OV01TnSnre9ZO03DPdrd3qs1e9bOWeXEoluIbh3479tqs2L32oi1DTyw0106mlUjHsqvC4YQorAksBOiwugOHe8JXyG5YCFYLNguvYjkd+eT+NcbgKptlwglcI9255WCzaVHh4rhCO6sVsIAlsI+jXVvE5arTl2+PGPcGDGjRxuxrhy+TNo1+3eqWTIp1Xh+KVWvAyOaHHAnCk3TcHgdss5OiBIlgZ0QlUjTsM07CNuf/9AR4EU/WU/L5620frQB3WHtaBM1VDqCu9X/I/TKWwXfUGHdtIE7AUt9fcHuwzRNQm99hPWdFUMe1IGqDWfz2Ig0Rvpcs5ZJu3YvTKxSqibxwMCDLLvXBkC8j925vd7WZyeVSOWVBhZCFFb+i2uEEKWvPcDjsLm4Qgla17fhuu0mnLE2kmf8ANsRh6o07hCxzd2fxKNPEGtOkVrXStU4N5qz96K9g6HFYhwARCL5tdjaFjNpENwYIF41HvvDj+MYX1uQ+/GM9tK6poVYaxxXTe5g21nlJLwlTLQ5indsZ2cQq82KzWMj1hrr9ba9sdqs6E4biWAcVx8FrHOxe+yARjKURK+VlxEhSonM2AkxEmgadp8Dm8tG6stfBZu1I0Wbaw3eYNjH1FA12Y/9t7/B+PE5xe9QkYdM+tV+8y/wT/IVLKgDsLl0HH4HkaZwr7N2mkXDWdPeDqzbTKiz2okRS/bZqqw3Dp+dRDg54NlVzaJh99ry6mAhhCgsCeyEGCGS0SRGzMB9zDxsf2pP0Vot6Ff8lNZ3PxvSF2m7x47tmCOK2qEiX13X1FmPOarPNmFDxVXrJm2k+yxa7KpxYabNHo+Tw+dAs2p5tfqyee1g5ldw2OayYeRR5FgIUVgS2AkxQkSbo+27Ie2da/D+9AeM3/8Js2YUgdWNJH949pDN4BWqQ0WmpdgNt/5ySM7XVSE3SvTF5tLRXTaizb1vOrHarOguW48SJZqm4fDlt5nB5tSx2qwkQgNvT6Y7dcy0KevshCgxEtgJMQKkjTTxQAJn93VYmoZtp5nUTKumymNAKDCkKdquwV3snt8PyYaKTEuxS66/hRsBc8yYQZ8T1EaJ2MNPDHtQl+GqcZKM9J1SdfodJEI9d7LavXZSyVRe6Vib104ij9Ilenv9Q6lnJ0RpkcBOiBEgk6brayesfftpWSla26UXEb/0qkGnaG1z9yd1z32Ejjqetvq2Idsta1ZV8ShgVlcP/lymSXBjkNCcw0je89thD+pAdZPQrBrRPlp92X2qV2v3x8TuUbOw+cza2T020kZ6wDNvFt2CRbdIOlaIEiOBnRAjQKw1hrNK9QjtU5cUbXLBQhLzjiJQHyCwZMWgZvBsu+9M1bRazHfexTj7/CGpc6eFghwDEAgM6jxm0iB5yRWYS5bin+THvvdugx5bPjKtvuJtsV43UWR2wfZIx1o07B5bXqVLbC5V9iSfAE136hLYCVFiJLATosIlYwapZAqHfwCbANoDPN8xh1A1uQrbq4sHnaK1e+x465xY310xJEWMrY1b+DlgaWjI+xyZNXW2ZYtx1ziGZaNEXxx+J2baJBnufc2bw+cgGTF6FCXOt+DwYGbeZAOFEKVHAjshKlwiGEezaNg8trxub/facV11SVaKNvnd+ST++/GAz1WUDhW9KNZGib7YnDoW3UK8j5SqehxNkpFkt8vVmrfurcf6Q3fqea2Vkw0UQpQeCeyEqHCJUBK7xz649mHdUrS43QQiOi2ft5L435oBzeB1De6Sl11V8A4VuZimSeLGBSUV1GXYfY4+NzPodhX8dQ/sOi7PN7DLMxUL+aVxhRCFISXDhahgqWQKI5bEVesbmhNmOlnMOwh/KEF4fRP6maeQHD8JTj+9350sbHP3J7ngFsJhK1p9G1WTqtAsQ9cBoy+ZjRLGYV/Gsu++2I4+fFjut78cXjuxlijJqIHNlfsp2u6x9wjsIP+ZN5vLRiQdwUgY6Pb+vyx0pHGjBvgHfLdCiAKQGTshKpiqT6a1t4AaWnavnZodx5O69udZKdrk4mX9msGzHbgfngN3J9kSIvarhQOuc2c6HKwEcLn6f5ukQXzBbSQ2NuHZc2bJBXWgUq2aRet71s6t1rZ1X09nc9tIRo0++87mPN8gZt50p04qkdr2FYUQ2xQIBFi0aBHf+MY32G677fB4POi6jq7rOBwOqqqqmDFjBsceeyz33nsv8XjP5wkJ7ERl27ABnnmm8/t0Gt55Bx5+GB78E7z2GhiVm0Yyokl0p3Xbu2Hz1T1Fa7Fg/u5+Wta0Eg/Gthng2T12qlPNOJ5+bMBFjFMTJnEqkJ48uV/Xz6ypczz6EFXBDUXfKNEbTdPaA7TeU6o2d+71dDaXDUwTIzawQGswGygsumXAGzaEENlM0+T222/npJNOAuDCCy9k06ZN2O123nnnHV599VWeeeYZbr31VvbZZx9efPFFzjzzTA477DBSqez/d0nFisqUNGDRvfDss3DNNeqyJUvhtoWwvj77ujO2g9tvh3HjhnuUBZeMJrG5h362rodMivawuSQamiGuEXnxNSyP/n6bKVrb3ruTvPnXWC++EOOsc9Hvuh3NNbCm9NvSY6PE4QcP6fmHmu7Uibb0vrFEt+ugaT1myqwOKwCpRO9p3L7uMxUf+MybBHZCDN4ZZ5zBpk2bePrpp9E0ja1btxKLxdhtt93YbbfsEkz/93//x7JlyzjiiCNYvnw57777LrNnz+44LjN2ovK0tsKPzoSPP4YHH4TZs+GGG+Du38A118Lzz8Fv7oLZe6vrr/4MLrxQzeZVEDNtkkqkOtJsw0LTsI8fRc20arxTRmWnaPsok9Kj/Vii99mqTEuxm889jbcAy//+1+eQzHQa45LLS3KjRG90p46ZMkklew+0dEfPQMxiVTNv+aRG8w3QMrcbaPpXCKE8/fTT/PWvf+X666/v2OS2Zs0aALbffvuctzn00EPZZZddAHC73VnHJLATlaWxCU47HSZOhIULobYWrr0WXn0V7r0X9toTxoyB/faDu34D+++vbvfxJ/Cv14s48KGXSasNa2DXhW2fvXp0sojc9+deO1lkgrv4DrvT1hDudbdspqXYZT+5fJtjME2T4KYQsZm7lE1QB52PWV/pWN1hzdlCzGKzDntgB8isnRB5mjx5Mh999BF77bVXx2Vr164FYIcdduj1dqFQiBkzZrDjjjtmXS6pWFE5AgE1UzduLFx9NVit8PTT8PwLKgVbU5N9fZsOl1wM3/gGmMDr/4KDDyrGyAtCBXYautNavEF0SdEml/6LxJhpJOsDuF97Efv0iT1StLa5+2PuuTfR9QHCf38Jz9EHo7n7vzmiKzNpEH3mn8R3nIP/5BOwleiaulysNmvnmrdedptaHVZigZ5Bsu6w5rlWztox8zaQ0jidgZ2JNb9SiUKMaHvuuWePy9577z0A9t1335y3Wb58OZ999hlPPvlkj/9XmbETlcE0VTC3dSv8/OcqqNvUADffDHNmw0G9BGxTp8IO7e921q0bvvEOAyNuoDutg6tfN1TaA7zqXSZRNcmPffmrvaZo7R47Vb4U7gXX5V3EOLOmznXDNVTRWrIbJfqyrTVvukMHs2e61mq3YuS1Vk79neTTuULdTnbGCjFU3nzzTTRNY7/99utxbO3atcyfP59bb72Vr371qz2OS2AnKsOLL8Krr8Hpp3fOzP32PojH4cwf9X3bWe1rGELhwo5xmKWNdOF2ww6C3efA9tu7e3SyiG9p7bzOxLGkblqQV4eK7hsl7DvNHPofYhhkZtB6o/WSArXoFjDNHi3H+nN/6nwDWyvXdcZOCDE03nzzTXbYYQdq2l/P1q1bx9KlS7n00ks56qijuOeeezj//PNz3rb0nvWFyMf994PTCV3fvdTUqM0Te+3Z9229XvXZ5y3Y8IpBBXZFTMP2pVuZFGPn3QhsTdLyeQuJV5eDaW6z/ZgxcRJfBdLTpnVcVoptwvJl0bU+AzuLNfcMm8XaHmgNOLDLnG9gM2+apqFZNZmxE2KIrFq1iqamJo444ghaW1s5/PDDmTp1KocffjiTJ09m5cqVHHHEEb3eXgI7Uf5WfQqfroYDD+gM0gDOPhtmzdr27YNB9XnixMKMr0hSyVTHi3XJag/wXFdfStXkKvSPVmI//+yOFK3twP1I3fxrCEcIrN6SvaHC7qAewK7KuZimSXDdVsytzWUf1EH2mrfcx3PPlGntgZ05wDhrMDNvVt0qM3ZCDJHnnnsOgGOPPRZN05g3bx6apmGaJg6HA4ej76UlEtiJ8vfCC+rzHnvkd/tNm9Tn9q3jlcA0TcyUWbozdjnYvXZ8xxzSI0WLkcL83e9I2DwE/vt5RxFj6+bNXA9omzZhJg1C760mnrBi3nVX2Qd1sO01b5mZMjPVPRWrbtf98m0Z1MybRRtQv2AhRO+ef/55XC4X8+bNo6qqiquuuorTTjsNgF/96lekt1GaSwI7Uf4W/1N9Hj9+4LdNJOC//wWrBQ4u7aK1A5EJBkp+xq677p0srBZib/6XcGMUpwvc1/20o0OFFgnxJYCWFrVR4tJz8I924PAPbXHjYunPmjeL1UJqiFKxQN4bbTQNqWMnxBCIRqMsW7aMo48+Gqez87nsqquuwmq1snr1ap588sk+zyGBnShv6TRs3Ki+7laksV/e+A/EE3DIIT3LoZSxTMoy8yJfdroEeI4zTgENrIvuhkQC6wcrMX58jtoYA1jvuQ/99Vcwzz0fR52vuOMeQlr7ZKvZx7tzzWLpMVOmWdpn7HqpA9jnfbanewZ+Q1TJICHEoDz11FPEYjFOPvnkrMunTJnC17/+dQBuvPHGHrdbsWIFr776KiCBnSh34QhkZiZaWwd++8ceU59PP33IhlQKOl6by2zCrgdNw+53UjOtGtvRR4Dfj8Uw0P/7Dr47bwHA8f57FbGmrruO2bNtBUw5j+f5wOcZoOUdEAohsvzxj39k1KhRHHvssT2OXXTRRYDaMfvnP/+54/KHHnqIr33ta0yfPh3oZ4Fi0zQJh0NDMWYhhoQ3s0mia0/RVavgmGP6f5KVK+G11+C7J0C3yt39YRoG0U8+IT11KgCWVat6zJ6kp04FhwOtoQEtEMi+fU0t5ug6iESw1HfrX6vrpLfbTp33s8/AyC44m540CdxutMYmtJbm7PP6/Ri1o4m0tmD844Ps13gNklPUefVN69G6te4y6sZgerxYAq1Yu5/X7cYYPQ5SBrb6njX/kpOngcWCvnkjWiyWdSxVW0fa50cLB9GbGrPP63BgjFMbV2xrP+txXmPCJEybnXQYzK98F/tOn1C75HlcTZsB+O/Uw1j/igGvPNvjtmXPVOVhek2pm2B1WnH4stPPZjqNs8qJzTOwPsFmKo2jyoFr08CKQptpE7vXjjVUPms6hSiUdDqAz+cb8NKGzZs38+KLL3LWWWdht/f8391333350pe+xHPPPccPf/hDNm7cSHNzMy+88AJLlixh0qRJ6opmP7S1tZmo93HyIR9F/3A4HNl/oHvvbZpgmrNmmWYq1Z8/adOMRk1zjz1Mc//9TTMS6d9tumtpMd/qMq4YqHF0+di5/dh9OY79sv3YoTmOre9y3vU5jh/afuyXOY7d135s5xzHYl3O+1aO499sP/aTHMf+3n6sLscxE0xf+/EXchz7cfuxk3Ice73LmHKdd0b7sQdzHHsVzHFF+juUD/mQD/no7aOtrW3ALylPPPGEef7555vBYLDX66xfv96cPXu26XA4zNmzZ5t33HGHGYvFsq6jmea2589N02THHXfkzTff3NZV87LPPvsU5NyFOm8gEGDy5MmsX78ev7+Xfj95KtSYy/XcvZ03693QL38JV1yhvn7gATjllL5Pmk7Dt7+tZvgWL4a6up7XSSTgttsgFIKLLoKGBpW2jcXg5JNhxgzSsRiR//6XdHsvP8t775GZsZs/fz4PPvigOuZyoa1bh9bcbQZszBjMCRMgGMSyenX2/dtspNt36Vo++ACSyezzzpgBPh/axo1oW7Zkn7e2FnPKFIhGsXzySfZ5NY10++5hyyefQLSzNtz8+fP5wyuvQE0N2pYtaJm1i5nz+v2Y220HyaQaU/df6267gdWqfpZMCZn28/5xyRLMujpobsbSvcOHx0O6vdG15d13e553xx3B6ST5v/9hNHbO9jU1NfGFk0/m5Q8/ZGIBStWU2/NSOZ67kM+lII/hcJxbHsPc8pmxGyr9SsVqmobVai3IgwYU7NyFHDOA3+8f8vOX4++5kOfu13lPOAF+9jMVjJ17LsyY0fsO1y1bYP58sFjglVegqir39e64QwV/118PX/86fOELcPHFUF8PBx4I77+PpbYW7z77dN6my31+4vHgPeSQzmO77tr7+P3+vmvoHXBA7+f1+3tPI/v9MHZs7+ftOvb2c/vb08r4/TCzj44NXcfQXZdG1pnz+trTyvj90KWg8IDOO2dO1rfJ+noaUE+g5fT8UY7/h4U+dyGeS0Eew+E8tzyGpaPfmyfOOuusgg2iUOcu5JgLpRx/z4U8d7/Ou912cMMN6utgEI48En76U3jrLUil1GXvvKMu22cfOOooePbZ3oM6UDttp05VgaDHA5dcArreGZQ89dTgx52Hin0MS1A5/j7K9dyFIo/h8J27UMrxMSy2fqViRbZAIEBVVRVtbW0VG/GXpbvvhp//vLP8SVdVVSpFe8klMGFC3+cxTbWxYtddYdQoWLiwM72bSKhA71e/ggsuGOIfQAxEfX19RwqoY9GwKCvyXFr+5DEsPf1KxYpsDoeDa665ZpttPcQw+9GP4LTT4PHHVWAWj6tU5L77qnSms5+FazVNdbF4+21oa4PDDus89vbbaodqf1qViYLK/P/J/2H5kufS8iePYemRGTshenPLLXD77bBmTedlV18Nd92lZgXliayoZKZACCF6kgLFQvRm6dLsBfvhMNx/v1rPJ0GdEEKIEiSpWCFySaXUrtnJk6GpCXw++MEP4DvfUZ+FEEKIEiSBnRC5vPuuWl+3dCncc49aV3faaXDEEcUemRBCCNErCeyEyGXpUlXuZM891YcQQghRBmSN3SCsWbOG0047jenTp+NyuZgxYwbXXHMNiUSi2EMTg/XyyzBvXrFHIXK466672HnnndmnW4FlUT5eeeUVvvzlLzNhwgQ0TePJJ58s9pDEAN1www3ss88++Hw+xowZw/HHH88n3TvdiKKQwG4QPv74Y9LpNPfeey8ffPABt956K/fccw+XX355sYcm8tXSotbRvfCCWmMnj2XJOeuss/jwww8L2g5IFFY4HGaPPfbgzjvvLPZQRJ6WLVvGWWedxfLly3nppZcwDIOjjjqKcDhc7KGNeFLuZIjdfPPN3H333Xz22WfFHorIh2l29GUFVOsxXVYslCIpd1IZNE3jiSee4Pjjjy/2UMQgNDY2MmbMGJYtW8YhfbUGFAUnr1hDrK2tjdra2mIPQ+RL08BuL/YohBCirLS1tQHI618JkFTsEFq9ejV33HEHZ555ZrGHIvL0/vvvF3sIQghRVkzT5IILLuCggw5i1113LfZwRjwJ7HK49tpr0TStz48VK1Zk3Wbjxo188Ytf5Fvf+hann356kUYuMro/hjvvvDPz58/ntttu4x//+AevvvoqwWAw6zaPP/44Bx10UJFGLIQQ5enss89m5cqVPPTQQ8UeikBSsTmdffbZnHDCCX1eZ9q0aR1fb9y4kXnz5nHAAQewaNGiAo9O9Ed/HkObzZb1ffdATwghRN/OOeccnnrqKV555RUmTZpU7OEIJLDLqa6ujrq6un5dd8OGDcybN4/Zs2fzwAMPYLHIJGgpGMhjKIQQYmBM0+Scc87hiSeeYOnSpUyfPr3YQxLtJLAbhI0bN3LYYYcxZcoUFixYQGNjY8excePGFXFkYiDWrVtHc3Mz69atI5VK8e677wIwc+ZMvF5vcQcnRAUKhUJ8+umnHd9//vnnvPvuu9TW1jJlypQijkz011lnncVf/vIX/v73v+Pz+WhoaACgqqoKl8tV5NGNbFLuZBB+//vfc+qpp+Y8Jr/W8nHKKafwhz/8ocflS5Ys4bDDDhv+AYl+kXIn5Wvp0qXMy1EA/Pvf/z6///3vh39AYsA0Tct5+QMPPMApp5wyvIMRWSSwE0KUJQnshBCiJ1kQJoQoK9JSTAgheiczdkKIsiQzdkII0ZPM2AkhhBBCVAgJ7IQQQgghKoQEdkIIIYQQFUICOyGEEEKICiGBnRBCCCFEhZDATgghhBCiQkhgJ4QQQghRISSwE0IIIYSoEBLYCSGEEEJUCOk8IYQoS6ZpEgwG8fl8vTYkF0KIkUYCOyGEEEKICiGpWCGEEEKICiGBnRBCCCFEhZDATgghhBCiQkhgJ4QQQghRISSwE0IIIYSoEBLYCSGEEEJUCAnshBBCCCEqxP8DL7aVSS/3Cy8AAAAASUVORK5CYII=\n",
      "text/plain": [
       "Graphics object consisting of 21 graphics primitives"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
Luke Naylor's avatar
Luke Naylor committed
    "def plot_central_charge(chern, name, argument = None, radius = None):\n",
    "    Z = stability.Mumford().central_charge(chern)\n",
    "    x = Z.real()\n",
    "    y = Z.imag()\n",
Luke Naylor's avatar
Luke Naylor committed
    "    \n",
    "    if not argument:\n",
    "        argument = arctan(y/x) if x != 0 else pi / 2\n",
    "        if argument <= 0:\n",
    "            argument += pi\n",
    "            \n",
    "    if not radius:\n",
    "        radius = (pi-argument/2)/pi\n",
    "    \n",
    "    return point(\n",
    "        Z,\n",
    "        marker = \"o\",\n",
    "        size = 600,\n",
    "        rgbcolor = \"white\",\n",
    "        #markeredgecolor = \"purple\",\n",
    "        zorder = 100\n",
    "    ) + point(\n",
    "        Z,\n",
    "        marker = name,\n",
    "        size = 500,\n",
    "        rgbcolor = \"red\",\n",
    "        zorder = 101\n",
    "    ) + line(\n",
    "        (0, Z),\n",
    "        rgbcolor = \"red\",\n",
    "        linestyle = \"dashed\",\n",
    "        zorder = 99\n",
    "    ) + disk(\n",
    "        (0,0),\n",
Luke Naylor's avatar
Luke Naylor committed
    "        float(radius),\n",
    "        (0, float(argument)),\n",
    "        alpha=.2,\n",
    "        fill=False,\n",
    "        thickness=1,\n",
    "        rgbcolor=\"purple\"\n",
    "    )\n",
    "\n",
    "p = sum(\n",
    "    plot_central_charge(chern, name)\n",
    "    for chern, name in [\n",
    "        (O1, r\"$\\mathcal{O}(1)$\"),\n",
    "        (O, r\"$\\mathcal{O}_X$\"),\n",
    "        (O1inv, r\"$\\mathcal{O}(-1)$\"),\n",
    "        (O2, r\"$\\mathcal{O}(2)$\"),\n",
    "        (O_x, r\"$\\mathcal{O}_p$\")\n",
    "    ]\n",
    ")\n",
    "\n",
    "xmax = (2.5)\n",
    "xmin = (-2.5)\n",
    "ymin = (-0.25)\n",
    "ymax = (1.5)\n",
    "aspect_ratio = (1)\n",
    "\n",
    "p += polygon(\n",
    "    [\n",
    "        (xmax + 1,0),\n",
    "        (xmin - 1,0),\n",
    "        (xmin - 1,ymax + 1),\n",
    "        (xmax + 1,ymax + 1)\n",
    "    ],\n",
    "    rgbcolor = \"yellow\",\n",
    "    alpha = 0.2,\n",
    "    zorder = 102\n",
    ")\n",
    "\n",
    "p.xmax(xmax)\n",
    "p.xmin(xmin)\n",
    "p.ymin(ymin)\n",
    "p.ymax(ymax)\n",
    "p.set_aspect_ratio(aspect_ratio)\n",
    "\n",
    "p.axes_labels([r\"$\\mathcal{R}$\",r\"$\\mathcal{I}$\"])\n",
    "p.show()"
   ]
  },
  {
   "cell_type": "code",
Luke Naylor's avatar
Luke Naylor committed
   "execution_count": 8,
   "id": "50d8265f",
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAb4AAAHWCAYAAAAIKuzeAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAo8ElEQVR4nO3de3Cc5WHv8d/77l3S7kqyLpbxlYsvsgkQY7DpxKfg1sUkDpMMHdK0DnTAlCRMbuW045NpQjLTcmaaUGaKoIGcJgWcCadxoGG4JMwE45wiAyVQwNjGOLYl25Jl3XZ13V3tvuePteUY27Jsr/bZ1fP9zOysd/Xuu7+V2PzyvLfH8TzPEwAAlnBNBwAAoJgoPgCAVSg+AIBVKD4AgFUoPgCAVSg+AIBVKD4AgFUoPgCAVSg+wADP85RMJsX1I4Dio/iAQvI86aGHpOXLpffekw4ckP7wD6UvflH6m7+RPvYxyXHkuK5izzwjx3FMJwas43DJMqBAPE/667+Wdu+WfvYz6dAh6frrpUcfldatyy+TzUp33CH9279JDQ1SW5sUCpnNDViG4gMK5W//Vnr5ZWnbNsnvl667Trr2Wumf//nk5RIJad68/P0rr0irV5vJC1iKTZ1AITz7rPRP/yRt3iyFw9I//mN+U+e3vnXqsvG4tGZN/t979xY3JwCKD9i2bZvWr1+vWbNmyXEcPfPMM2d9zXvvvXfiQTotfelL0l/+pXTZZdLwsPTYY9IPfiDV159+BZdckr/PZi/8AwA4JxQfrDc0NKQrrrhCDz300KSW37dvn37xi1+ceOJXv5IOHswXnyRFItIHH0gbNpx5JccLb+bM80wN4Hz5TQcATFu3bp3WHT/4ZBL+5V/+5eTlf/ITqbJSuvrq/GPHye/jm0h7e/7+4x8/x7QALhTFB5yj1tZWfeUrX8k/GBuT/uM/8gernK3sjvM86T//M1+Us2ad/LOxMemNN05+7tprJdfNPz82ln9uxgxp4cIL+yCApdjUCZyjzs5OhcPh/IPe3vw+vUBg8itobZUOH86f2/dR6XT+FIcvflH64z/OH/ziOPmy/MUvpJtvzh85mkgU5sMAFmLEB5yH7PF9dJlM/v7gwXw5TeaE9AcflC6//PT7ACsqpFtvlYJB6bOfla66Kr/OsTEplZJ++1tp9uyCfQ7ARoz4gHM0c+ZMJZPJ/IP6eqmqKj/y++1vz/7il1+Wnn9eevzxiUeJ69fnC+6RR6RcTrrvvvyRo5QecMEoPuAcrVq1Sq2trfkHwaD0yU/m//33fz/xC/fuzR/5+e//Ll155cTL+v3SX/1VviC/9S3pz/5Mmj//QqMDEMUHaHBwUG+//bbefvttSfnTFd5++221tbVJkjZt2qQvfOEL48vffffd2rZt24kVbNyYv3/6aenv/i4/Qvt9npcf5f3pn+aPAJ3sEaR33imNjOQ3fy5der4fD8BHcMkyWG/r1q26/vrrT3n+tttu049//GPdfvvt2r9/v7Zu3Tr+s9dff13XXHPNiYW/+U3pH/4h/+/mZumWW/JHbB4+LP3yl9IVV+R/PmPG5IM98UR+f2A6Lb3zzuT2HwI4K4oPKJTWVumBB06cjjB/vnTDDdJtt+VPdzgXTz+dHynW1+ev5bltm/SJTxQ8MmAjig8oNS+9lB8p3nZbvvyuuCI/ivzpT00nA6YF9vEBpWTrVmnHjnzpSfnNm1/+srRli9TRYTQaMF1QfEAp+OlP8+ftrVmTP1/vuLY26c038+fxbdgg/fjHxiIC0wWbOoFS0NV1ovCCQamxMf/v0VHp6NETy4VC+QlsAZw3ig8AYBU2dQIArELxAQCsQvEBAKxC8QEArELxAQCsQvEBBmR6e3XkJz9RprvbdBTAOhQfYEDqvffU+Od/rtT775uOAliH4gMAWIXiAwBYheIDAFiF4gMM8ObM0Z2SRmsalBpIKT2Yzt+G0sqMZJRNZ5Uby4krCgKF5zcdALBJS0uLWlpalEll9KGk/9nrKJ3umfA1juPI8Tlyfa4c98S/Xb8rX9AnN+DKFzj53mG2duCMuEg1YMDh1nf0reuu0D/u/EDRBQukY99Cz/PkZT15OU+5bE5e7vSPc9mccpmcspmscplTR4au/0QJ+sN++UP+/H3YTzHCeoz4AAN8hw/ph5IGuzrkX3zZBa8vN3aiBH//PpvOKpVMaTg1PF6OjuPIH/bLF/KNl6E/7FcgEpDjUoiY/ig+YBpw/flNn4qc/uee5ymbziqbympsdCx/S41ppHdE2XR2fDl/2K9gZVCBisD4jTLEdEPxARZwHCe/uTPkVygWOulnXs7T2OiYMsOZ8dtI78j4CPGkMqw8VoZsKkUZo/gAyzmuMz66O87zPI2NnL4MHddRsCqoUDSkYDRIEaLsUHyAAV5FpbZKurqqynSU03Kc05dhZjij9ED+1IuBjgF5hyhClB+O6gQMOLjroOYsmaNEIqFYLGY6znn5/SI8fi6il8sXYSgWUjgeVigeki/gMx0VOAkjPsCEXE7BY/flynEcBSuDClYGVTWz6qQiHE2Mqv9AvyQpWBlUuDpfgoFIYOKVAkVA8QEGBHbuUErS4DvvSKtXm45TEB8twtxYTqOJUY32j2qgY0DJQ8n8wTXxkMLVYQWrgmwShREUH4Ap4fpdVcyoUMWMCnk5T6mBlFKJlEb7RjXUNSTX5ypcE1akNqJQNHT2FQIFQvEBmHKO6ygcDyscDys+N67McEaj/aMa6R3RcPewfEGfIrURRWojbA7FlKP4ABTd8SNGo7OiSg+mxwtwsHNQgUhgvAR9QQ6MQeFRfACMClYFFawKKjYnplQypZGekfF9gqFoSJEZEUVqIlxBBgVD8QEGZC5bpNmS3m9uNh2lZDjOic2huWwuvym0Z0T9+/uVPJjM7y+sr5A/xP9s4cLwXxBgQjCoQ8fucSrXd+LAmLHUmIaPDuc3hR4ZVDgeVkV9hUKxEEeF4rwwES1ggK/9gP6vJGffPtNRSp4/5FdsdkyNH2tU9fxqZTNZ9X7Yq673ujTYOajcWPmeCwkzGPEBBrjJpP5U0mAiYTpK2XBcZ3wUmB5Ka/josAY6BjRweECR2ogqGys5IhSTQvEBKDvHT5SPzY5puGdYQ11DGu4ZVjgeVlVTlYKVbELGmVF8AMqW63dV1VilyoZKjfSOaLBzUN27uhWKhlQ1s+qUKZgAieIDMA04zonNoKP9oxrsHFTPnh4FKgKqmlmlcHWYA2EwjoNbAAOyDY3aJMlrajIdZdoJV4dVt7hOMxbOkOtz1fe7Ph3dcVTD3cNiMhpITEsEGDEdpiUqF+mhtAY7BzXaPyp/2K/orKgiNRHTsWAQIz6giFpaWtTc3KxP37ha6yWpv99woukvWBlU7SW1ql9SL1/Qp77f9al7V7dSAynT0WAIIz7AgCNbXlDjLTdp8JVXVDVNpiUqF6mBlAYODSg9lFYoFlLsothJM81j+uPgFgBWCUVDCi0OabR/VMlDSR3deVSRmoiiF0W5HJol+CsDsNLxWeFHekY0cHhAR3ccVUV9haKzonJ97AWazig+ANZyHEcVdRWK1EY01DWkwc5BjfSOKHZRTBV1FabjYYrwf2sAA7xQSDskKRw2HQXKXw6tamaV6pfWKxQLqf9Av7p3dSsznDEdDVOAg1sAAzidobSlB9NKtCWUGcmosr4yv/nTzzhhuuAvCQAfEawKqm5JneJz4hrpHVHXji4Ndw+bjoUCofgAAwI7dyghyX3nHdNRcAaO46iyoVL1S+sVjofZ/DmNUHyACbmcYsfuUdp8AZ+q51erblGdvJyn7l3dGjg8wOXPyhjFBwCTcHzzZ9XMqvwsEDsZ/ZUrig8AJslxHEVnRVW3uE6SGP2VKYoPAM5RoCLA6K+MUXyAAWMXX6qPS8otXGg6Cs4To7/yRfEBBniRiN6SpAquDlLuThn97erWWGrMdCxMgOIDDPAdPqSHJDnt7aajoAB+f/TnZT117+zWSO+I6Vg4A4oPMMDt69WXJTk9PaajoICOj/5C8ZD69vWp/0C/vBybPksNF6kGgAJyfa5qFtQoFAvlL3s2lFH1gmoFIsz5VyoY8QHAFKiYUaH6JfWS8ge+cMmz0kHxAcAU8Yf9qltcp0htRP0H+tW3r0+5LFfrMY3iAwzIzajTA5K8+nrTUTDFHNdR9bxq1SyoUSqRyh/1OcpRnyZRfIAB2ZlN+mtJ3kUXmY6CIonURlS35MQ5f6OJUcOJ7EXxAQY4Q0NaKUmDg6ajoIj8ofymz2BVUL0f9mrwCH9/Eyg+oIhaWlrU3NysL938R2qV5H74oelIKDLX56rmkhpVzaxS8mBSffv6OOWhyJiBHTDgyJYX1HjLTRp85RVVrV5tOg4MGekdUf+BfgUiAdVcUiNfwGc6khUY8QGAIZHaiOoW1Smbzqp7Z7fSQ2nTkaxA8QGAQcev9uIL+tSzu4dLnRUBxQeY4Pfr6LF7wBfwacbCGYrURtS3r4+DXqYY3zrAgMyiJZojKbFsmekoKBGO66h6frXcgKvkwaRymZxis2OmY01LFB8AlJDYRTH5Aj4l2hPKjeUUnxeX4zimY00rbOoEDPDv2a09ktydO01HQQmqbKhUzYIajfSOqG8vpzsUGsUHGOCk07pUklIp01FQoiK1EdVeWqvUQEo9H/QoN8Y1PguF4gOAEhWKhVS3qE5jqTF17+5WNp01HWlaoPgAoIQFKgKqW1wneeIC1wVC8QFAifOH/JqxaIZcv6ueD3oovwtE8QEGjM2brz+RlLv4YtNRUCaOn+vn+l117+5WZiRjOlLZovgAA7yqqH4lSTHO08LkuX5XMxbOyF/l5YMeyu88UXyAAW7XEX1bktPZaToKyozrdzXjMsrvQlB8gAG+o126TxQfzg8jvwtD8QFAGXJ9x0Z+AR8HvJwjig8AytT4yO94+aUov8mg+ACgjB0vP8fnqOeDHmUznOR+NhQfYEAuHteTkrzqatNRMA0cP+BFntS7p1e5LJc3mwjFBxiQnT1XGyR58+ebjoJpwhf0qfayWmXTWS5sfRYUH2BCalSXSNLoqOkkmEYCkYBqL61Veiit/v398jzK73QoPsCAwId79KEkd9cu01EwzQSrgvkpjfpGlGxPmo5Tkig+oIhaWlrU3Nysu756l+komMbC1WFVz6vW0NEhDXQMmI5TchyPsTBQdEe2vKDGW27S4CuvqGr1atNxME0Ndg4qeSip6nnVqqirMB2nZPhNBwAATI2qmVXKZrLqP9AvN+AqHA+bjlQS2NQJANNYfE5c4XhY/fv6ubrLMRQfYEBm6eVyJOWuvNJ0FFigekG13ICr3r2c4ydRfAAw7bk+V7WX1iqXyanvd33Wn+ZA8QEG+Pft1auS3D17TEeBJfwhv2ouqVF6IK2BQ3Yf6UnxAQY4w8NaJUlDQ6ajwCKhaEix2TENHhnUcM+w6TjGcFQnAFiksqFSmZGMEgcS8of9ClYGTUcqOkZ8AGCZ+Ny4ApUB9e3tUzZt32wOFB8AWMZxHNVcXCM5svJgF4oPMGDsotn6C0m5uXNNR4GlfAGfai6uUWY4o4HDdh3sQvEBBnjVNdosSbW1pqPAYsHKoKKzohrsHFQqmTIdp2goPsAAt7dHX5LkdHebjgLLVc2sUigWUt++Pmtmb6f4AAN8HYfVIsk5eNB0FEA1C2rkOI7699kxhx/FBwCWc/2uqhdUKzWQ0mDnoOk4U47iAwAoFA0p2hTVwOEBpQfTpuNMKYoPACBJqmqqUrAqqL7f9Sk3Nn0vZk3xAZIefvhhLViwQOFwWMuXL9dvfvObMy67detWOY5zym3Xrl2Tfj+vskq/lKRo9MLDAwVy/Pw+z/PUf6DfdJwpQ/HBek899ZS+9rWv6Zvf/KbeeustfeITn9C6devU1tY24et2796tjo6O8dtll1026fccm79AN0rKXXLJBaYHCssX8Ck+N67R/lGN9I6YjjMlKD5Y74EHHtAdd9yhO++8U0uWLNGDDz6oOXPm6JFHHpnwdQ0NDZo5c+b4zefzTf5Ns1lFj90DpSZSE1GkJqJEe2JanuJA8cFq6XRab775ptauXXvS82vXrtWrr7464WuvuuoqNTU1ac2aNXr55ZfP6X0Du95XUpL77rvnGhkoivjcuCQp0ZYwnKTwKD5Yrbu7W9lsVo2NjSc939jYqM7OztO+pqmpSY8++qi2bNmin//851q0aJHWrFmjbdu2nfF9UqmUksnk+G14xN4pYVAeXL87bTd5Mi0RoPxO/d/ned4pzx23aNEiLVq0aPzxqlWr1N7eru9973tavXr1aV9z//336zvf+c7446sk/fbCYwNTKlIT0WjNqBLtCQWjQfkC57A5v4Qx4oPV6urq5PP5ThnddXV1nTIKnMjKlSu1Z4LZ1Ddt2qREIjF+2/LEz847M1BM03GTJ8UHqwWDQS1fvlwvvfTSSc+/9NJLuu666ya9nrfeektNTU1n/HkoFFIsFhu/VUQqzjszUEzTcZMnmzphvW984xvasGGDrr76aq1atUqPPvqo2tradPfdd0vKj9YOHTqkxx9/XJL04IMPav78+Vq6dKnS6bSefPJJbdmyRVu2bJn0e2YWLla9pL1Ll07FRwIK6vc3eYZiIbn+8h4zUXyw3q233qqenh5997vfVUdHh5YtW6bnn39e8+bNkyR1dHScdE5fOp3Wvffeq0OHDikSiWjp0qV67rnndNNNN03+TQMBdR+7B8pBfG5cXe91KXkoqep51abjXBDHs+FS3ECJ6fjV/9Mbf/IJXf/WW4peeaXpOMCkDB0dUqItobrFdQpWBk3HOW/lPV4FypQ7MKBPS3KSSdNRgEmrqKtQoCKgRFuirKcvovgAAJPiOI7ic+PKDGc03F2+56JSfACASQtWBlUxo0IDhwbKdgYHig8AcE5is2OSpOSh8txUT/EBBmQbZ+obkrxZs0xHAc6Z63cVnRXVcPew0kPlN2ktxQcYkKur1z9J8hoaTEcBzktFfYUCkYCS7eU36qP4AAOcRL9ukaS+PtNRgPNy/ECX9FBawz3ldaALxQcY4D/Yrn+X5B44YDoKcN6CVUFFaiIaODwgL1c+pzdQfACA8xadFVUukyur0xsoPgDAefOH/YrMiGigY0C5bHmc3kDxAQAuSHRWVF7O09CRIdNRJoXiAwzwwuH8RLSRiOkowAXzBXyqbKjU4JFBZTNZ03HOiuIDDBi75DItl5T7vZncgXJW1Vglx3E02DloOspZUXwAgAvm+l1VzazS8NFhjaXGTMeZEMUHGBB4/z2NSnL/+79NRwEKprKhUq7f1cDhAdNRJkTxASZ4nkLH7oHpwnEdVTVVaaR3RGOjpTvqo/iAImppaVFzc7Pu+updpqMAU6KirkK+gK+k9/UxAztgwJEtL6jxlps0+Morqlq92nQcoKAGjwxq4NCAGpY1yBf0mY5zCkZ8AICCqqirkOM6GjxSmqM+ig8wIHPJpVoqTmfA9OT6XFU2VGq4e7gkJ6ul+AATwhG9L3ECO6atyoZKSdJQV+ldzYXiAwzwHTqoxyQ5bW2mowBTwvW7qqir0FDXUMldw5PiAwxw+/t0pySnt9d0FGDKVDVWyct5JTdzA8UHAJgSvqBPkRkRDR0ZUimdQEDxAQCmTFVjlbKZrEZ6R0xHGUfxAQCmjD/sVygW0vDR0tncSfEBBmTr6nS/JK+hwXQUYMpV1lcqPZRWZjhjOookig8wItfYpP8lyZs1y3QUYMqF4iH5gj4NHS2NUxsoPsAAZ2hQ/0OSBkr7KvZAITiOo8r6So30jpTEqQ0UH2CAf/8+bZXk7t1rOgpQFJEZEclTSZzaQPEBAKacL+BTuCZcEge5UHwAgKKorK/UWGpMqWTKaA6KDwBQFMGqoAKRgPGDXCg+wAAv4NdBSQoETEcBiqqivkKj/aPKZrLGMlB8gAFjC5dojqTc0qWmowBFFamNyHEco1dyofgAAEXj+lyFq8MUH2Ab/wc71S7J3bHDdBSg6CK1EWWGM8qMmLmSC8UHGOBkxjRbkjKlcQknoJhC8ZBcn2ts1EfxAQCKynEcRWojFB9gg5aWFjU3N+uur95lOgpgVKQ2omw6q/Rguujv7XilNDsgYIkjW15Q4y03afCVV1S1erXpOIARR949olAspOp51UV9X0Z8gAFj8xfoDyXlLrnEdBTAmEhtRKN9o0WfnZ3iAwzwKqv0iiRFo6ajAMZUzKhQLptTKlHcS5hRfIAB7pEO/YMk5/Bh01EAY/xhvwKRgEb7R4v6vhQfYICvu1ubJDldXaajAEaFq8MaTRR3cyfFBwAwJhQPKTeWU2aoeOe0UnwAAGOClUH5Ar6ibu6k+AAARoXiIY0mKD5gWstV1+iHkrzaWtNRAOPC1WGNjY5pbHSsKO9H8QEGZC+arY2SvLlzTUcBjAtFQ3Jcp2ijPooPMGF0RM2SNGJuahagVDiuo1AsVLT9fBQfYEBg74faIcndvdt0FKAkhONhpQfTyo3lpvy9KD4AgHGheEiSirK5k+IDABjnC/gUiASUHpj62RooPgBASQhGg0WZpojiA0xwHKWO3QPIC0VDGkuNKZvOTun7UHyAAZnmZQpLyl1xhekoQMkIVgUlSamBqZ2tgeIDAJQE1+8WZT8fxQcY4N+7R2+K0xmAjwpGg4z4gOnIGR3VxyVOYAc+IhQNKZvOTul+PooPAFAygtGp389H8QEASobrcxWomNr9fBQfUEQtLS1qbm7WXV+9y3QUoGQFq6b2fD7HK+Z87wAkSYdee09fW3m5/s/+/YrNm2c6DlBSRnpH1LevTzOvnCnXV/jxGSM+wAAvXq2fSVJNjekoQMkJVAQkSZnhzJSsn+IDDHC7j+rrkpyuLtNRgJLjD/vluI4yQxQfMG34jnTqAUnO4cOmowAlKVARYMQHALAHxQcAsEqwMqix1Jhy2cJPTEvxAQBKzlQe4ELxAQbkolH9QpIXi5mOApSkqTzAheIDDMjOna+bJXkXX2w6ClCypmo/H8UHmJDJqO7YPYDTC1QElBmh+IBpIfDBLh2V5O7YYToKULL8Yb+yqawKfYExig8AUJL8Ib88zyv4FEUUHwCgJPnDfknS2OhYQddL8QEASpIv6JPjOsqmGPEBACzhD/kZ8QHTQWZxs2KScpdfbjoKUNL8YYoPmB58Pg0cuwdwZr6QT2Mpig8oe/79+/SiJHfvXtNRgJLmD/uVTWfl5Qp3SgPFBxjgDA3qTyRpYMB0FKCkjR/ZWcBRH8UHAChZvmB+d0Ahj+yk+AAAJcsX8MlxHGUzFB8AwBKu31UuU7h5+fwFWxOAs2ppaVFLS4uiw6O6WtL/nj3bdCSg5LkBt6AjPscr9NU/AZzVwV0HNWfJHCUSCcWYkw+YUO/eXsmTai+tLcj62NQJGOD09+nPJam313QUoOT5Aj728QHlzn/ooJ6U5La1mY4ClDw3UNh9fBQfAKCkHR/xFWrPHMUHAChpbiBfVbmxwoz6KD4AQEnzBfInsRdqcyfFBxjgVVSoVZIqK01HAUqe689XVaEOcKH4AAPGFlyi6yTlLrvMdBSg5Dk+R5IKdqFqig8AUNIc91jxZSk+oGwFdrwrT5L79tumowAlz3EcOa7DiA8AYA/HdZTLcnALAMASrs9lxAcAsIfjOuzjAwDYw/Gxjw8oa5lLL9OlknKLF5uOApQF1+eyjw8oa6Gw9kpSOGw6CVAWOKoTKHO+g216QpKzf7/pKEBZoPiAMucmEvoLSU5/v+koQHlwJBVo2nSKDwBgFYoPAFDyHMcp2Hx8/sks5HmeBgYGCvKGgI1SqZRSqdT440RPtxol9R3p03Bbl7lgQJkY7hlWbiynUDI04XLRaFSO40y4jONNokKTyaTi8fi5pQRwRjMl/ZWkH0jqNJwFmE4SiYRisdiEy0yq+CYz4luxYoXeeOONc0s4Saz7hGQyqTlz5qi9vf2sf9zzUW6/j6led6HW+9ERX0dHh6655hq99857uuiiiy54/R917apr9VrrawVf71SuuxwzJ5NJzVswTwf2HZiS72M5/k6mOvPZvo+TGfFNalOn4zhn/aP6fL4p+cOz7tOLxWJTsu5y/X1M1bqnMrMkxWviqq6tLvh6fT7flKx3KtddjpmPT5BaXVs9Zf/9ldvvZKozF+L3XLCDW7785S8XalWs26By/X1M1brL8W8o8bsu5rqnUjn+Tsoh86Q2daJ0HN/fOpnt2ChdBw8eHN9kPXv2bNNxcJ74PpYnTmcoM6FQSN/+9rcVCk18ZBNK2/G/H3/H8sb3sTwx4gMMYKQAmMOIDwBgFYoPAGAVig8AYBWKDwBgFYqvDPT19WnDhg2Kx+OKx+PasGGD+s8ync3tt98ux3FOuq1cubI4gYFp5OGHH9aCBQsUDoe1fPly/eY3vznjslu3bj3le+c4jnbt2lXExDibSV25BWZ9/vOf18GDB/Xiiy9Kku666y5t2LBBzz777ISvu/HGG/WjH/1o/HEwGJzSnMB089RTT+lrX/uaHn74Yf3BH/yBfvCDH2jdunV6//33NXfu3DO+bvfu3ScdrVtfX1+MuJgkiq/E7dy5Uy+++KK2b9+ua6+9VpL02GOPadWqVdq9e7cWLVp0xteGQiHNnDmzWFGBaeeBBx7QHXfcoTvvvFOS9OCDD+qXv/ylHnnkEd1///1nfF1DQ4Oqq6uLlBLnik2dJa61tVXxeHy89CRp5cqVisfjevXVVyd87datW9XQ0KCFCxdq48aN6upi+htgstLptN58802tXbv2pOfXrl171u/eVVddpaamJq1Zs0Yvv/zyVMbEeaD4SlxnZ6caGhpOeb6hoUGdnWee0GbdunXavHmzfv3rX+v73/++3njjDd1www0nzRAA4My6u7uVzWbV2Nh40vONjY1n/O41NTXp0Ucf1ZYtW/Tzn/9cixYt0po1a7Rt27ZiRMYksanTkPvuu0/f+c53Jlzm+PQbp5tiw/O8CafeuPXWW8f/vWzZMl199dWaN2+ennvuOX32s589z9SAfT76PZvou7do0aKTdj+sWrVK7e3t+t73vqfVq1dPaU5MHsVnyD333KPPfe5zEy4zf/58vfPOOzpy5MgpPzt69Ogp/090Ik1NTZo3b5727NlzzlkBG9XV1cnn850yuuvq6jqn797KlSv15JNPFjoeLgDFZ0hdXZ3q6urOutyqVauUSCT0+uuv65prrpEkvfbaa0okErruuusm/X49PT1qb29XU1PTeWcGbBIMBrV8+XK99NJL+sxnPjP+/EsvvaSbb7550ut56623+N6VGPbxlbglS5boxhtv1MaNG7V9+3Zt375dGzdu1Kc+9amTNqksXrxYTz/9tCRpcHBQ9957r1pbW7V//35t3bpV69evV11d3UlfYAAT+8Y3vqEf/vCH+td//Vft3LlTX//619XW1qa7775bkrRp0yZ94QtfGF/+wQcf1DPPPKM9e/Zox44d2rRpk7Zs2aJ77rnH1EfAaTDiKwObN2/WV77ylfGjyz796U/roYceOmmZ3bt3K5FISMrPUvzuu+/q8ccfV39/v5qamnT99dfrqaeeUjQaLXp+oFzdeuut6unp0Xe/+111dHRo2bJlev755zVv3jxJUkdHh9ra2saXT6fTuvfee3Xo0CFFIhEtXbpUzz33nG666SZTHwGnwbREgAFMSwSYw6ZOAIBVKD6giFpaWtTc3KwVK1aYjgJYi02dgAFs6gTMYcQHALAKxQcAsArFBwCwCsUHALAKxQcAsArFBwCwCsUHALAKxQcAsArFBwCwCsUHALAKxQcAsArFBwCwCsUHALAKxQcAsArFBwCwCsUHALAKxQcAsArFBwCwCsUHALAKxQcAsArFBwCwCsUHFFFLS4uam5u1YsUK01EAazme53mmQwC2SSaTisfjSiQSisVipuMAVmHEBwCwCsUHALAKxQcAsArFBwCwCsUHALAKxQcAsArFBwCwCsUHALAKxQcAsArFBwCwCsUHALAKxQcAsArFBwCwCsUHALAKxQcAsArFBwCwCsUHALAKxQcAsArFBwCwCsUHALAKxQcAsArFBxRRS0uLmpubtWLFCtNRAGs5nud5pkMAtkkmk4rH40okEorFYqbjAFZhxAcAsArFBwCwCsUHALAKxQcAsArFBwCwCsUHALAKxQcAsArFBwCwCsUHALAKxQcAsArFBwCwCsUHALAKxQcAsArFBwCwCsUHALAKxQcAsArFBwCwCsUHALAKxQcAsArFBwCwCsUHFFFLS4uam5u1YsUK01EAazme53mmQwC2SSaTisfjSiQSisVipuMAVmHEBwCwCsUHALAKxQcAsArFBwCwCsUHALAKxQcAsArFBwCwCsUHALAKxQcAsArFBwCwCsUHALAKxQcAsArFBwCwCsUHALAKxQcAsArFBwCwCsUHALAKxQer9fX1acOGDYrH44rH49qwYYP6+/snfM3tt98ux3FOuq1cubI4gQFcML/pAIBJn//853Xw4EG9+OKLkqS77rpLGzZs0LPPPjvh62688Ub96Ec/Gn8cDAanNCeAwqH4YK2dO3fqxRdf1Pbt23XttddKkh577DGtWrVKu3fv1qJFi8742lAopJkzZxYrKoACYlMnrNXa2qp4PD5eepK0cuVKxeNxvfrqqxO+duvWrWpoaNDChQu1ceNGdXV1Tbh8KpVSMpk86QbADIoP1urs7FRDQ8Mpzzc0NKizs/OMr1u3bp02b96sX//61/r+97+vN954QzfccINSqdQZX3P//feP70eMx+OaM2dOQT4DgHNH8WHaue+++045+OSjt//6r/+SJDmOc8rrPc877fPH3XrrrfrkJz+pZcuWaf369XrhhRf0wQcf6LnnnjvjazZt2qREIjF+a29vv/APCuC8sI8P084999yjz33ucxMuM3/+fL3zzjs6cuTIKT87evSoGhsbJ/1+TU1Nmjdvnvbs2XPGZUKhkEKh0KTXCWDqUHyYdurq6lRXV3fW5VatWqVEIqHXX39d11xzjSTptddeUyKR0HXXXTfp9+vp6VF7e7uamprOOzOA4mFTJ6y1ZMkS3Xjjjdq4caO2b9+u7du3a+PGjfrUpz510hGdixcv1tNPPy1JGhwc1L333qvW1lbt379fW7du1fr161VXV6fPfOYzpj4KgHNA8cFqmzdv1uWXX661a9dq7dq1+tjHPqYnnnjipGV2796tRCIhSfL5fHr33Xd18803a+HChbrtttu0cOFCtba2KhqNmvgIAM6R43meZzoEYJtkMql4PK5EIqFYLGY6DmAVRnwAAKtQfAAAq1B8AACrUHwAAKtQfAAAq1B8AACrUHwAAKtQfAAAq1B8AACrUHwAAKtQfAAAq1B8AACrUHwAAKtQfAAAq1B8AACrUHwAAKtQfEARtbS0qLm5WStWrDAdBbAWM7ADBjADO2AOIz4AgFUoPgCAVSg+AIBVKD4AgFUoPgCAVSg+AIBVKD4AgFUoPgCAVSg+AIBVKD4AgFUoPgCAVSg+AIBVKD4AgFUoPgCAVSg+AIBVKD4AgFUoPgCAVSg+AIBVKD4AgFUoPgCAVSg+AIBVKD6giFpaWtTc3KwVK1aYjgJYy/E8zzMdArBNMplUPB5XIpFQLBYzHQewCiM+AIBVKD4AgFUoPgCAVSg+AIBVKD4AgFUoPgCAVSg+AIBVKD4AgFUoPgCAVSg+AIBVKD4AgFUoPgCAVSg+AIBVKD4AgFUoPgCAVSg+AIBVKD4AgFUoPgCAVSg+AIBVKD4AgFUoPgCAVSg+oIhaWlrU3NysFStWmI4CWMvxPM8zHQKwTTKZVDweVyKRUCwWMx0HsAojPgCAVSg+AIBVKD4AgFUoPgCAVSg+AIBVKD4AgFUoPgCAVSg+AIBVKD4AgFUoPgCAVSg+AIBVKD4AgFUoPgCAVSg+AIBVKD4AgFUoPgCAVSg+AIBVKD4AgFUoPgCAVSg+AIBVKD6giFpaWtTc3KwVK1aYjgJYy/E8zzMdArBNMplUPB5XIpFQLBYzHQewCiM+AIBVKD4AgFUoPgCAVSg+AIBVKD4AgFUoPgCAVSg+AIBVKD4AgFU4gR0wwPM8DQwMKBqNynEc03EAq1B8AACrsKkTAGAVig8AYBWKDwBgFYoPAGAVig8AYBWKDwBgFYoPAGCV/w8M6IFocQvblwAAAABJRU5ErkJggg==\n",
Luke Naylor's avatar
Luke Naylor committed
      "text/plain": [
       "Graphics object consisting of 4 graphics primitives"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "plot_central_charge(O, r\"$\\mathcal{O}_X$\", radius = 0.9)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "1cfa37b3",
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAbwAAAHUCAYAAABBF53oAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA8p0lEQVR4nO3dfZAcZ50n+G++VGbW+0tXv0nWm18kuW3Y9RkZyRMwO/hGawOG8cRw9nIrRjNYcxB2sODxRJyPuMFwO+c/Doxvd9osxgtext45YsfAHWMjzneggwjb2DD2ApYtybKk1ku3urveX7OyMvP+yOqWW1K3uqWszKrK7yei41G3qrJ/anX3t54nnxfBtm0bREREA070uwAiIiIvMPCIiCgQGHhERBQIDDwiIgoEBh4REQUCA4+IiAKBgUdERIHAwCPymG3bKJfL4BJYIm8x8Ig8VqlUkEwmUalU/C6FKFAYeEREFAgMPCIiCgQGHhERBQIDj4iIAoGBR0REgcDAIyKiQGDgERFRIMh+F0DUjyzTgm3ZsE0btmVf8P7FPmaZFmADlZqz/q54oggrYUEQBQiCAAiAIAgQxGX+LAmQQhLEkAhRFp3nENGqMfCIzmOZFsyW6bzp5uKf23obZsuE1bYueQ1BFCBKohNmkhNc3/77b+M/Pv0fYVomAMBsmTDqBmA7u6/AhhOQC3+2naBczkL4XaplMBI5BJv7G1HA2JaNdrN9QZAtBJxlngs0QRAgKZLzpnbakLQkyM4PtsUe2zLK5TKSySRKpRISicSl6+0EoGVasAwLpmFe0Jot589W27pgyzJRFiGrMuSwDFmTEQqHIIdlSCHp8r+IRH2IPTwaaGbLhNEwYNQNtBttp9Xbi38viOcCTYkqkNIXhpvfFoY7JdGpJ4TQio+32p0e6rsDUXd6k418Y7HXKMriYvi9OwhFibf2aTAx8GhgmIYJo2agVWvBqBkwGsbi8KMoiQhFQlCTKmKRGGRNhqzKEOXB++Uuys49vosFo23bTq+20YbRcF4E6GUd9bn6Ys9QUqQlQajEFMgqf1VQ/+N3MfUl27IXg22hNQ3n3pgUkhCKhhAdiSIUCSEUDkFS/O+p9QJBEJzhTVWGltIWP25bNtp6e0kQNvINmK1zX1MlrkCNq1DiDEDqT/yupb5hNAzoZR16WUer0oJt24s9t/BQGEpUQSga6olhyH4jiAJCYefFQRjhxY9bpoVWtYVWpYVWtYVivgiAAUj9id+l1LOstrUYcHpZh2mYEEQBalxF4qqE84tWkzkLsYtESYSW1KAlnd4gA5D6Gb8rqWfYtg2j9q5eXK0FAE6vIxOGmlChxBRnbRr5Yi0BqCZUaCkNakLl/xn1BAYe+cpsmWiWmovDlJZpQZRFqAkVqeEU1ITKIcoetlIA6mUd9VwdgihAS2nOW1Jj+JFvGHjkOattoVFooJFvoFVtQRAEZ5LJaBRqQkUoEuIwZZ86PwDbzTYahQaaxSYK7xSc8Etq0NIMP/IeF56TJ2zLRrPURCPXgF7WAQBqQnWGKpNqoNZ+rXXh+aBo6200C000Cg0YdcO5H5tQEU4H73uA/MEeHnWNbdtoVVto5JxX+JZpQYkqSFyVQDgTHsg1cLQ8WZURG4shNhZbDL9msYnCsQLDjzzBwCPXLezo0cg3YBomZFVGdDSKcCbM2XsEYGn4mS3TGfYsdMJPEKClNUSyEahx1e9SaYDwtw+5wmyZaOQbqOfqaDfbEGUR4UwY4YyzPo5oOZIiITYaQ2z0XPjV5+vIHc5BVmVEhiOIDEU4IkBXjPfw6IroFR212RqaxebibLyFJQSceHJxQb2Ht1Z6RUd9vo5moQkA0FIaIsPs9dHlYw+P1sy2bTQLTVTPVmHUDciajNSmFLS0xnsv5Bo1rkKNq7A2WKjn6ud6fZqM6EgUkaEIZ3nSmrCHR6tmtS3U5+uozdZgGibUhIrYaAxqgq+414I9vMv37hEFURIRyUYQGY7w3jCtCr9L6JLazTZqszXUc3XABsJDYWdj5vDKx9QQuW2h12e2TNTmaqjP11E9W4WW1BAdifLFF62IgUfL0is6amdraJaaEGURsbEYosNRTh4g30mKhMT6BOLjcTTyDdRma8gdySEUCSG+Lr648J3o3Rh4tIRt24u/QIy6gVA4hNSmFMKZMO+XUM8RRMEZ1sxGoJd1VKYryL+dhxJVEBuPMfhoCQYeAegEXa6BynQFZsu5Pzd03RCHiKhvqAkVakJ1gu/MueCLr4vz+5gAMPAIQLPYRPl0Ge1mG+F0GLFrY7w/R33r/ODLHclBiSmIjzP4go6BF2B6RUfldAWtWgtqQkV6SxqhCIOOBsNC8DVLTVSnq+eCb12ca/kCioEXQEbDQOV0Bc1SE6FIiEOXHpmcnMTk5CRM0/S7lEBZOL2hWWo6Pb7DDL6g4jq8AGnrbVTOVNDINyCrMuLr4winw36XFThch+evheAz6gbUuIr4ujiUGLe/CwL28ALANExUZ6qoz9UhyiKSG5OIZCPc+osCabHHV2yiMl3B/KF5hDNhJK5K8LDhAcfAG2CWaaF2tobq2SoEQUB8XRzRkSiXFxABi6ew13N1lE+VMffGHOLr4ogM88XgoGLgDahGvoHSyRJs00Z0JIrYWIwLxokuIjIUgZbSUDldQelkCfX5OpIbkxzmHEAMvAHT1tsoTZWgl3WE051hGoXDNEQrEaVzQ/2lqRKHOQcUJ60MCNu2UTtbQ2W6snifjrtM9CZOWul9C8OcsMFhzgHCHt4AaNVaKJ0owWgYiI3GEF8X5306oivAYc7BxMDrY5ZpoXKmgtpsDaFICMPXD3PhOJFLOMw5eDik2aeapSZKUyVYbevc7EsOufQFDmn2p/OHOaMjUb9LojViD6/PmIaJ8skyGoUG1ISK5MYkD78k8sD5w5zNUhOpzSn29voIf1P2kfp85xWmAKS3pBHOcJcUIi8tDHNqKQ3F40XMHZxDalMKWooTxPoBA68PmIaJ4rEi9IqOSDaCxPoE19QR+UhNqBieGEbxRBH5o3lEshEkNyQ5WazHMfB6XLPURPF4EYIgYGjrEDe7JeoRoiwic00G9fk6SidLaFVaSG1JQYlyJmevYuD1KNu2UTldQfVsFVpSQ2pzir06oh4UyUagxBUUjxWRO5RDbDyG2FiMk8h6EAOvB5ktE4V3CjDqBhJXJRAbjfldEhGtQFZlDG0bQnW6isqZCvSSjtSWFCeU9Rh2GXpMs9jE3ME5mIaJoW1DDDuiPrGwQXt2exZW28LcwTnUc3W/y6J34cuPHmFbNsqny6jN1hBOh5HclIQo8fUIUb9Rogqy12dRPllG8XgRzWITqU28JdELGHg9oK23UXingHajjeTGJKLDXNBK1M9ESURqcwpqUkVpqoS5g3NIX5PmhBaf8SWHzxr5BuYOzsG2bGS3Zxl2RAMknA5jeGIYkiIhdziHRr7hd0mBxh6eT2zLXtyUNpwJI7mRQ5hEg0gKSRjaOoTSVAmFYwUYDcPZ4J2zOD3HwPOB2TKRfzuPtt5GanMKkaGI3yURURcJooDU5hTksIzyqTLajTZSW1J8kesxfrU9ZtQNzL81D8u0kN2eZdgFyOTkJCYmJrBjxw6/SyGfxEZjGLpuCK1qC/NvzaOtt/0uKVB4WoKHmsUmCscKCIVDSF+T5qazAcXTEqjdbCN/NA/LsJC+Js0dlDzCHp5HarM15I/moSU1DG0dYtgRBZisychuzyIUDSF/JI/abM3vkgKB9/C6zLZtlE856+tiozEkruIreiJyli5krs0sHjdkNAwkNyY5maWLGHhdZFs2Cu8UoJd1rq8jogsIgoDEVQnIYRmlEyW0m21krslwkXqX8KvaJaZhYv7QPPSKjsy1GYYdES0rMhTB0LYhmLqJuTfnYDQMv0saSAy8LjAanZmYhjMTU03whjQRrUyJKshuz0KUReQO5dCqtfwuaeAw8Fyml3XkDuUgyiKy12cRCof8LomI+oSkOIvU5bCM3OEc9Irud0kDhYHnovp8Hfm381BiCrLbspyJSURrJkoihq4bghJTkH87j2ap6XdJA4OB55LqTBXFE0VEshGkr0lDEDnTiogujyAKyFybgZpQUThaQKPAPTjdwMBzQeVMBeXTZcTXxTmtmIhcIQgC0lenoaU1FN4p8Gw9F3BZwhUqny6jOlNFYn0CsTEe1kpE7hEEZw9OQRRQPF6Ebdmc8X0FGHhXoHyqjOrZKpIbkoiO8JuQiNwnCAJSm5zQK02VYJs2X1xfJgbeZSpNlVCbq3FBORF5IrnBOUKsfLoM27IRXxf3u6S+w8C7DAthl9qUQiTL0w6IyBvxdXEIorAYetyqcG0YeGtUOtkJO55jR0Q+iI3FIEjO8KZlWpwotwYMvDVY2AQ6tYlhR0T+iQ5HFyeyAEBqU8rXevoFA2+Vyqc7E1Q2JjmMSUS+W3jRXTxehBSSeE9vFRh4q1A5U3GWHlyV4AQVIuoZkaEILMNC+XQZYkjk76dLYOBdQnWmisp0xVlnN8qpwETUW2JjMZiGidJUCaIsIpwO+11Sz2LgraCeqzs7qIzHue6FiHpWckMSVttC8VgRoixCjfOElovh1mLL0Cs6SidKiGQjHBsnop6X2pyCEldQOFqAUed5ehfDwLuIdrONwtEClLiC5Mak3+UQEV2SIAjIXJOBrMnIHcmhrbf9LqnnMPDOYxomckdykBQJ6avTXN9CrpmcnMTExAR27Njhdyk0oBZOWRBlEfkjeZiG6XdJPUWwbdv2u4heYVs2codzMFsmstuzkBSeZ0fuK5fLSCaTKJVKSCS4Uwa5z2yZmH9rHmJIxNDWIYgS+zYAe3hLFI4VYDQMZK7NMOyIqG9JioTMdRmYuonC0QLYr3Ew8DrKp8poFptIb0kjFAn5XQ4R0RUJhUPIXJtBq9pC8ViRoQcGHgCgNldbPOZHS2l+l0NE5AolpiB9dRqNQgOV0xW/y/Fd4AOvWWqiNFVCdCTKM+2IaOBoKQ3JDUlUz1bRKDT8LsdXgQ48o2Gg8E4BWlLjMRtENLCiI1GEM2EUjxfRbgZ3uUJgA880TOSP5CFrMpcfENHAS21KQVZl5I/mYZmW3+X4IpCBZ9s28m/nAQHIXJuBIDLsiGiwCaKA9NVpWIaF0omS3+X4IpCBVz5VRrvRRuaaDKQQlx8QUTDImozU5hQahQaqZ6t+l+O5wAVes9hEbbaGxFUJLj8gosDRUhpiYzGUT5WhV3S/y/FUoALPbJkoHi9CS2mckUlEgRVfF4caV1F4pxCo7ccCE3i2baPwTgGCJCC1OeV3OUREvhEE536eIAqB2oklMIFXOVOBUTeQvjrNfeWIKPBEWUT66jSMuoHyqbLf5XgiEL/59bKO6kwV8fVxKFHF73KIiHqCElWQ2JBAbbaGRn7wF6UP/InnpmGicMxZXB4b5anlRETvFh2OwqgZKJ4oQg7LCIUHdzLfQPfwbNtG8VgRgsD7dkREy0luTELWZOd+njW49/MGOvCqM1XoFR2pLSmI8kD/U4mILpsgCkhvScM0TJRPD+79vIFNAb2io3Kmsjj9loiIlidrMhLrnft5enkw1+cNZOBZbQvFY0WocRWxMd63o0t7/PHHsWXLFmiahptvvhm/+MUvln3sgQMHIAjCBW9vvfWWhxUTuS86EoUaV1E8URzI/TYHMvCKx53DDlNbUtwUmi7pe9/7Hj7/+c/ji1/8Il577TV84AMfwB133IGpqakVn3fo0CFMT08vvl133XUeVUzUPanNKdimjfLJwRvaHLjAaxQaaJaaSG1KcZ9MWpVHH30Un/70p3Hvvffi+uuvx2OPPYYNGzbgG9/4xorPGxkZwdjY2OKbJPH7jfqfpEhIbEignqujWWz6XY6rBirwLNNCaaqEcDrMk8tpVVqtFn79619j9+7dSz6+e/duvPjiiys+96abbsL4+Dhuu+02/OxnP1v2cbquo1wuL3kj6mWRoQi0lOYMbbYHZ2hzoAKvfKoM2EBiAw9zpdWZn5+HaZoYHR1d8vHR0VHMzMxc9Dnj4+N44okn8Oyzz+L73/8+tm3bhttuuw0///nPL/r4Rx55BMlkcvFtw4YNrv87iNyW2pQCAJRODs5RQgOz8Fyv6KjP15HcmORQJq3Z+fd6bdte9v7vtm3bsG3btsX3d+3ahZMnT+KrX/0qPvjBD17w+IceeggPPPDA4vvlcpmhRz1PlEUkNyRROFYYmFGzgejh2ZaN0lQJSkxBdJinINDqZbNZSJJ0QW9udnb2gl7fSnbu3IkjR45c9O9UVUUikVjyRtQPwpkwtKSG0lRpIGZtDkTgVWeqMHUTyY1Jv0uhPqMoCm6++Wa88MILSz7+wgsv4NZbb131dV577TWMj4+7XR6R75Ibk7AteyA2mO77IU2jYaA6U0VsLDbQe8BR9zzwwAPYs2cP3ve+92HXrl144oknMDU1hc985jMAnCHJ06dP47vf/S4A4LHHHsPmzZtxww03oNVq4emnn8azzz6LZ5991s9/BlFXSIqE+Pq4MyEwE+7rjTz6PvBKUyVIqsQF5nTZ7r77buRyOXzlK1/B9PQ0brzxRjz//PPYtGkTAGB6enrJmrxWq4UHH3wQp0+fRjgcxg033IDnnnsOH/7wh/36JxB1VXQ4imahidKJEoYnhiGI/bm+WbD7+OS/2lwNpakShrYO9fWrDgqWcrmMZDKJUqnE+3nUN9p6G3NvzCE2FkN8Xdzvci5L397DMw0TldMVRLIRhh0RUZfJqozoaBTVs1WYLdPvci5L3wZe+WQZEIDEer5CJiLyQmwsBkEUUDlT8buUy9KXgdcsNtEoNJDckOSxP0REHhElEfF1cdRzdRh1w+9y1qzv0sK2bZROlqAmVIQzYb/LISIKlEg2glA41Jc7sPRd4NXn6jBbJpIbuOaOiMhrgiAgcVUCrWoLjULD73LWpK8CzzItVKadiSqy1vcrKoiI+pKaUKElNVROV2Bb/TPRv68Cr3a2Btu0ER/vzymxRESDInFVAmbLRG225ncpq9Y3gWcaJqpnq4iORCEp3ByaiMhPsiYjMhxBdabaN0cI9U3gVWeqEASBO6oQEfWIhdG2flmm0BeB19bbqM/VERuLcRkCEVGPEGURsfEYanM1GI3eX6bQF+lROVOBKIuIjvDoHyKiXhIdiUJW5b44TaHnA89oGGjkG4iNx/p2w1IiokG1sExBL+tolpp+l7Oing+8yukKZFVGJBvxuxQiIroILaVBjas9fy+vpwNPrzivGOLr4xAE9u6IiHpVbDwGo25AL+t+l7Ksng68yukKQpEQwmluIUZE1MvUuAolqqAy3bu9vJ4NvGaxiVatxdMQiIj6RGwshla1hVa15XcpF9WzgVc5U4GaUKEmeNYdDYbJyUlMTExgx44dfpdC1BVaSoOsyajOVP0u5aJ68sRzvawjdySH7LYslJjidzlEruKJ5zTI6rk6iseLGJ4YRigc8rucJXqyh1c9W0UoEmLYERH1mXAmDEmRerKX13OBZzScWT6xUW4hRkTUbwRBQGw0hka+gbbe9rucJXou8GqzNUghCVpa87sUIiK6DJFsBKIsona2t05S6KnAs9oWGrkGoiNRrrsjIupTgiggOhJFfb4O0zD9LmdRTwVeba4GCOCuKkREfS46EoUgCj11Xl7PBJ5t2ajN1hAZivBEBCKiPidKIiLDEdTn6rDM3jgvr2eSpZFvwGpbPBGBiGhAREeisC0b9bm636UA6KHAq83WoCWdRYtERNT/pJCESDaC6tkqbMv/Jd89EXh6WYfRMBAdZe+OiGiQREeisNoWmkX/jw7qicCrzdYQioSgxrmNGBHRIJE1GUpMQT3n/7Cm74HXbrbRLDV5746IaEBFshHoZR1my98lCr4HXvVsFVJIQjjDI4CIiAZROB2GIAqoz/vby/M18CzTQiPfQGQ4woXmREQDShAFhDNh1HN1+Hlega+B1yw2YVs2IkNcaE5ENMgi2QjMlolWxb+z8nwNvEa+ATWuQlIkP8sgIqIuU6IKZE32dVjTt8AzDRN6Wee9OyKigIhkI2gWm7Da/uy84lvgNfINCILAUxGIiAJi4fZVI9/w5fP7GnhaSoMo+T5RlIiIPCDKItSk6tuwpi9p0262YdQNDmcSEQVMJBuB0TBg1A3PP7cvgVfP1SFKTtITEVFwqAkVUkjypZfnS+A18g2EM2GuvaNAmZycxMTEBHbs2OF3KUS+EQQB4aEwGvmG5xtKC7bHqwBb1RbmD80juy0LJaZ4+amJekK5XEYymUSpVEIikfC7HCLPtfU2Zn83i/SWtKe3tjzv4dVzdUiKxLAjIgooWZWhRBXPT1DwNPBs20az0ORkFSKigNNSGpqlpqdbjXkaeHpJh2Va3EqMiCjg1KQK27I93WrM08Br5BsIRUI81ZyIKOBC4RAkRfJ0WNOzwLMtG80ihzOJiMixMKzpFc8CT6/osG0bWpJbiRERkRN4Zsv0bBG6d4FX1iEpEocziYgIAKDEFIiS6Fkvz9PAY++OiIgWCIIANaF6dh/Pk8AzWybazTbUBLcSIyKic7SUBqNuwDTMrn8uTwKvWWpCEAQocS42JyKic9SECkEQPOnleRJ4ellHKBriUUBERLSEKItQYgr0kt79z9XtT2DbzsJCDmcSEdHFqEkVelnv+mbSXQ88o2bAMi0GHhERXZSW0mDbNvRyd3t5XQ88vaxDlEWEIqFufyoiIupDsipD1uSu38fzJPDUuMqz74iIaFlaUoNe6eMentW20Kq1eLI5ERGtSIkpMFsmzFb3lid0NfAWxmN5/46IiFYSijq3vVq17p2e0PXAC4VDkEJSNz8NERH1OSkkQVZltKp9HHjs3RER0WqEoiEYte5tJN21wGvrbZiGyd1ViIhoVZSYAqNudG09XtcCbyGllSgDj4iILk2JKrBtu2vHBXUt8Fq1FmRVhihzOzEiAJicnMTExAR27NjhdylEPUkOyxBEoWsTVwTbtrvSd5w/NA9JkZDeku7G5Yn6VrlcRjKZRKlUQiKR8Lscop6SO5yDIAnIXJNx/dpd6X4tdEm5uwoREa1FNyeudCXw2o02bMvm/TsiIloTJabANEy09bbr1+5K4LVqLQiCwB4eERGtyUJHqRu9vK4EnlE3Fm8+EhERrZYoi5C17ixA71rghcLs3RER0dopUaUrMzVdDzzbttFutDmcSUREl6VbC9BdD7x2sw3bthl4RER0WeSwDMDJEze5H3gNp8CFgomIiNZC1vok8Iy6AUmRIErcYYWIiNZOlESIsuj60gT3A6/BBedERHRlZE3ujx4eZ2gSEdGV6PnAs0wLVttaHH8lIiK6HLIqw9RNV6/pauCZLac4SeUJ50REdPlkTYZlWjAN90LP3cDrpLGkMPCIiOjyLYwUutnLc72HJwgCpBADj4iILt/CSKGb9/FcDzz27oiI6EoJggBZlV1dmuBq4LX1Nu/fERGRKyRVYg+PiIgGn9tLExh4RETUkxaWJti2O5tIuxZ4tmU7a/BUrsEjIqIrJ2sybNteXPJ2pVwLvIUbi+zhERGRG9xemuBa4C0uOmfgEV3U5OQkJiYmsGPHDr9LIeoLYsiJKLcWnwu2S4OjtbkayifLGLtpDIIguHFJooFULpeRTCZRKpWQSCT8Loeop02/No34ujhio7ErvparPTxJkRh2RETkGlEWYZs9NmnF1DlDk4iI3CXKIqy25c61XLkKuCSBiIjc15OBZ7UtiDJPOSciIveIUg8Gnm3ZECTevyMiIvf0Zg/PtCCIDDwiInKPKIuwzB4LPNuyIUoc0iQiIvf0XA9vIX3ZwyMiIjeJsgjbsmFbV740wZXAWyiE9/CIiMhNC7niRi/P3cBjD4+IiFy0MPvfjft47gReZxU87+EREZGbFgOPPTwiIhpkPRd4i5NWeA+PiIhctDBy2DOBxx4eERF1iyAKgAv7R/MeHhER9TRBEODGSXau9fDYuyMioq7plR6eZVrs3RERUXe41J9iD4+IiHpezwxpwoZrCUxERPRuguBOwMireZBt26hUKsv+faVWQaPcgFbWXCmKaJDoug5d1xffX/hZmj0xi0a84VdZRH2jXCmjrbZhR1fu5cXj8RXDUbBX0U8sl8tIJpNrr5KILjAG4H8A8E0AMz7XQjRISqUSEonEsn+/qsC7VA+veraKD+z+AF77r69dXpWXsGPHDrz66qt9de1uXbdcLmPDhg04efLkiv+xl6vfvh79cO3ze3j5/f8P3rvvz3Dk6f8DIx+5/UpLvMAt778Fr/zyFdev26/X7tZ1y+UyNm7aiKkTU135Wey3r4cX177Uz+KlenirGtIUBGHF/1CxIUISpK78pwOAJPXftbtZMwAkEom+qrsf/w+7dW0xmwUARKNRJFPuj5xIktSV6/brtbt13YWJeslUsms/M/309fDi2lf6dXZl0oogCtj73+1141IXdd999/XdtbtZczf149ejX6/dLf369ejH771u6sevR69fe1VDmpdSn6+jeKKIdTevu+KCaGUL91MvNVZNvevMD/8R6+66E2d+8COs+6OP+l0OXSb+LPYfd5YldIZM3VgnQStTVRVf+tKXoKqq36XQZZKyWTzdaal/8Wex/7jSw2sUGii8U8D4TeNcgE50CbmZHLLjWcxPz2NobMjvcogCw517eJ1ZMQunJhDRCppNXNNpicg7HNIk8ph0+BDe7rRE5B1Xe3hu7GZNRETUDa4tSwA4pElERL2LQ5pERBQIHNLsA4VCAXv27EEymUQymcSePXtQLBZXfM7evXshCMKSt507d3pTMNGAePzxx7FlyxZomoabb74Zv/jFL5Z97IEDBy74mRMEAW+99ZaHFdNKVrW12KUIkhN4Vtty43J0nk9+8pM4deoU9u/fDwD4i7/4C+zZswc/+tGPVnze7bffju985zuL7yuK0tU6aXXM9/4zCADm3/vP/C6FVvC9730Pn//85/H444/j937v9/DNb34Td9xxBw4ePIiNGzcu+7xDhw4tWYg+PDzsRbm0Cq4EnhSSAACmYbpxOXqXN998E/v378fLL7+M97///QCAb33rW9i1axcOHTqEbdu2LftcVVUxNjbmValEA+XRRx/Fpz/9adx7770AgMceeww/+clP8I1vfAOPPPLIss8bGRlBKpXyqEpaC9cmrYiSCMtgD89tL730EpLJ5GLYAcDOnTuRTCbx4osvrvjcAwcOYGRkBFu3bsW+ffswOzvb7XJpFcS3j+DFTku9qdVq4de//jV279695OO7d+++5M/dTTfdhPHxcdx222342c9+1s0yaY3cmbQCQAyJ7OF1wczMDEZGRi74+MjICGZmlj9N7Y477sAzzzyDn/70p/ja176GV199FR/60IeWHFND/hDqdezqtNSb5ufnYZomRkdHl3x8dHR02Z+78fFxPPHEE3j22Wfx/e9/H9u2bcNtt92Gn//8516UTKvgypAm4Axrsoe3eg8//DC+/OUvr/iYhbOfLna+k23bK577dPfddy/++cYbb8T73vc+bNq0Cc899xz++I//+DKrJgqW83/GVvq527Zt25JbDLt27cLJkyfx1a9+FR/84Ae7WietjmuBJ4ZEmC328Fbr/vvvxz333LPiYzZv3ozf/OY3OHv27AV/Nzc3d8Grz5WMj49j06ZNOHKEw2hEl5LNZiFJ0gW9udnZ2TX93O3cuRNPP/202+XRZXKvh6dIMGqGW5cbeNlsFtlV7Ja/a9culEolvPLKK7jlllsAAL/85S9RKpVw6623rvrz5XI5nDx5EuPj45ddM1FQKIqCm2++GS+88ALuuuuuxY+/8MIL+PjHP77q67z22mv8meshrt3Dk0IS7+F1wfXXX4/bb78d+/btw8svv4yXX34Z+/btw0c/+tElwyfbt2/HD37wAwBAtVrFgw8+iJdeegnHjx/HgQMHcOeddyKbzS754SV/WBs24l93WupdDzzwAJ588kl8+9vfxptvvokvfOELmJqawmc+8xkAwEMPPYRPfepTi49/7LHH8MMf/hBHjhzBG2+8gYceegjPPvss7r//fr/+CXQeV4c0bcuGZVoQJddylAA888wz+NznPrc4Y+xjH/sY/vZv/3bJYw4dOoRSqQQAkCQJv/3tb/Hd734XxWIR4+Pj+IM/+AN873vfQzwe97x+WspOp/EMgP89nfa7FFrB3XffjVwuh6985SuYnp7GjTfeiOeffx6bNm0CAExPT2Nqamrx8a1WCw8++CBOnz6NcDiMG264Ac899xw+/OEP+/VPoPO4ch4eALSqLcwfmsfIDSOQNddylGjg5H93CP/ze7bjf/ntW8jcuPw6SiJyl6vLEgAuPidazuTkJCYmJvDnH/5DTAIQz5z2uySiQHH1Hh4ALk0gWsZ9992HgwcP4qknn/K7FKJAci3wFnZbYQ+PiIh6kauzS8QQtxcjIqLe5GrgcWkC0aXZsRh+0mmJyDvuBp4qod1su3lJooFjXX0Nbu+0ROQdVwMvFA6h3Wjz5HOilZgm4p2WiLzjauDJYRm2bcPU+YNMtBzpjd+h3GmJyDuu9/AAwGhwT00iIuot7s7SlEWIssj7eERE1HNc3/Ry4T4eERFRL3E98OSwzCFNIiLqOd3p4TXbsC3O1CS6GPP6CQx3WiLyTld6eAB4H49oOaEQ5jstEXnH/cDrHA3EYU2iixOPH8P/2WmJyDuuB54oiZAU7rhCtByhXMbHOi0ReacrR5NzpiYREfWargQeZ2oSEVGv6VoPz2yZsEweFURERL2haz08ABzWJLoIa3wcD3RaIvJOdwJPkyFKIvSK3o3LE/U1e3gEX++0ROSdrgSeIAhQYgpa1VY3Lk/U14RiEX/SaYnIO10JPABQ4k7g8Ww8oqXEqRP4L52WiLzTtcBT4ypsy4ZR42xNIgCYnJzExMQE9t671+9SiAKpa4Enh3kfj+jd7rvvPhw8eBBPPfmU36UQBVLXAo/38YiIqJd0LfAA3scjuigtjH/qtETkna4GHu/jEV3I3LoVN3daIvJOVwOP9/GIiKhXdDXwFu/jVXgfj2iB9NvfoNlpicg7XQ08oHMfr8b7eESLbBtqpyUi73Q98Hgfj4iIekHXA4/38YiIqBd0PfB4H4+IiHpB1wMPANSEila1xfPxiACY123FDZ2WiLzjSeBpKQ22bUMvcViTCOEwDnZaIvKOJ4EnKRKUqIJGoeHFpyPqaeLJk/hWpyUi73gSeACgpTXoJZ3DmhR4QiGPezstEXnHs8ALp8Mc1iQiIt94Fngc1iQiIj95FniAM3lFL+uwLe4wQURE3vI28NIabMtGs9T08tMS9RR7eASPdFoi8o6ngSerMkKREJoFBh4FlzU+jv+p0xKRdzwNPMCZvNIsNTmsScFVreD3Oy0RecfzwOOwJgWd9M47ONBpicg7ngeerMoIhTmsScEzOTmJiYkJ7L13r9+lEAWS54EHOL08DmtS0Nx33304ePAgnnryKb9LIQokXwIvnA7DtmzoZS5CJyIib/gSeLLmDGs28lyETgEUCuFUpyUi7/gSeAAQyUbQLDZhtbm3JgWLef0ENnRaIvKOb4EXzjhHo9Tn636VQEREAeJb4ImyCC2tMfAocKQ3D+JkpyUi7/gWeAAQHY6irbehVzh5hQLEMHBVpyUi7/gaeEpMgazJqM+xl0dERN3la+ABTi+Pk1eIiKjbfA88Tl4hIiIv+B54oiwinAmjNleDbXPnFRp85tVX4190WiLyju+BBwDRkSjMlolmkftrUgDE4vj/Oi0ReacnAi8UCUGJKZy8QoEgTk/jf+20ROSdngg8wOnl6RUdRoNTtWmwCXOzeKjTEpF3eibwtJQGKSShNlvzuxQiIhpAPRN4giAgOhJFI9/gEgUiInJdzwQe4GwoDRuozbGXR0RE7uqpwBNlEZHhCGpna7BM9vJoMNnpDJ7stETknZ4KPACIjcVg2zZqZ9nLo8FkbdiAfZ2WiLzTc4EnhSREshHUZtnLowHVaGCi0xKRd3ou8AD28miwSUcO441OS0Te6cnAk0ISosNR9vKIiMg1PRl4ABAdjbKXRwNlcnISExMT2HvvXr9LIQokwe7hHZvLp8qoz9cxcuMIRLlns5loTYr/90+R+pe3ofiT/xep3R/yuxyiwOjpFFm4l1c9W/W7FCL3CAL0TktE3unpwBNlEdGRzr087r5CA8J8z3uhdVoi8k5PBx4AxEZjAMBeHhERXZGeDzz28mjQSIcP49edloi80/OBB7CXRwOm2cB/02mJyDt9EXjv7uWZhul3OURE1If6IvAAp5cniAIqpyt+l0JERH2obwJPlEUk1idQz9XRqrb8LoeIiPpM3wQe4JyXp0QVlKZK6OH18kQrsjZuwic6LRF5p68CDwCSG5MwGgbqc3W/SyG6LHYqhX/otETknb4LvFAkhOhwFJUzFU5gob4kzM3iC52WiLzTd4EHAPF1cUAAJ7BQXxKnp/FopyUi7/Rl4HECCxERrVVfBh7ACSxERLQ2fRt4AJDYkOAEFiIiWpW+DjwlqiCSjXACC/UVO5HA/9Vpicg7fR14AJBYn+AEFuor1uYt+HinJSLv9H3gcQILXYlCoYA9e/YgmUwimUxiz549KBaLKz5n7969EARhydvOnTtX/0kNA9lOS0Te6fvAAziBhS7fJz/5Sbz++uvYv38/9u/fj9dffx179uy55PNuv/12TE9PL749//zzq/6c0psHMddpicg7st8FuCW5MYn5t+ZROVNxhjmJLuHNN9/E/v378fLLL+P9738/AOBb3/oWdu3ahUOHDmHbtm3LPldVVYyNjXlVKhG5YCB6eICzA0t8XRzVmSr0iu53OdQHXnrpJSSTycWwA4CdO3cimUzixRdfXPG5Bw4cwMjICLZu3Yp9+/Zhdnb5XVN0XUe5XF58q9c5q5jIDwMTeAAQHY1CjasoHi/CMnk6Oq1sZmYGIyMjF3x8ZGQEMzMzyz7vjjvuwDPPPIOf/vSn+NrXvoZXX30VH/rQh6DrF3+h9cgjjyzeI0wmk/joXXe69m8gotUbqMATBAGpzSnYpo3SVMnvcsgnDz/88AWTSs5/+9WvfgXA+Z45n23bF/34grvvvhsf+chHcOONN+LOO+/Ej3/8Yxw+fBjPPffcRR//0EMPoVQqLb794w9+5M4/lIjWZGDu4S2QFAnJjUkUjhWgJTWEM2G/SyKP3X///bjnnntWfMzmzZvxm9/8BmfPnr3g7+bm5jA6Orrqzzc+Po5NmzbhyJEjF/17VVWhquri+8aOW5AAcOyGG1f9OYjoyg1c4AFAOBNGs9REaaoEJaZAUiS/SyIPZbNZZLPZSz5u165dKJVKeOWVV3DLLbcAAH75y1+iVCrh1ltvXfXny+VyOHnyJMbHx1f3BElCpdMSkXcGakjz3ZIbkxAkAYVjBS5VoIu6/vrrcfvtt2Pfvn14+eWX8fLLL2Pfvn346Ec/umSG5vbt2/GDH/wAAFCtVvHggw/ipZdewvHjx3HgwAHceeedyGazuOuuu1b1ecV3jmJ/pyUi7wxs4ImSiNTmFFrVFmpna36XQz3qmWeewXve8x7s3r0bu3fvxnvf+1783d/93ZLHHDp0CKWSc09YkiT89re/xcc//nFs3boVf/qnf4qtW7fipZdeQjweX9XnFKpV/MtOS0TeGcghzQVqXEVsLIbKmQrUhIpQJOR3SdRjMpkMnn766RUf8+4RgnA4jJ/85CfdLouIumBge3gL4uvikMOyM7RpcWiTiCioBj7wBEFAeksaZstE+VTZ73KIiMgnAx94ACBrMhJXJVCbq6FZavpdDgWctW497uu0ROSdQAQeAESHo9CSGorHimg3236XQwFmZ7N4vNMSkXcCE3gAkNqSghgSkT+a59Zj5BuhUMB/32mJyDuBCjxREpG5NgPLsFB4h+vzyB/iySk83WmJyDuBCjwAkFUZ6avT0Ms6T0knIgqQwAUeAKgJFckNSVTPVlHP8agWIqIgGOiF5yuJjkRhNAyUTpQgazKUqOJ3SURE1EWB7OEtSG5MIhQNoXC0ALNl+l0OBYQdieClTktE3gl04AmCgPTVaUAA8kfz3ImFPGFdex1u7bRE5J1ABx4ASCEJmWsyaDfbKJ4o+l0OERF1SeADDwBCkRBSm1No5BuoznAHe+ou6Tf/FXanJSLvBHbSyvnC6TDa422UT5chh2VoSc3vkoiIyEXs4b1LfF0cWkpD4Z0CWrWW3+UQEZGLGHjnSW9JIxQOIf92HkbD8LscIiJyCQPvPIIoIHNdBpIiIX8kz42myTWTk5OYmJjA3nv3+l0KUSAJNjeUvCirbWH+0Dxsy0Z2WxaSIvldEg2I3PHTeP+Wq/DLY6cwtJlHBBF5hT28ZYiyiKGtQxAEAbnDOZgGF6aTSzQNRzstEXmHgbcCKSQhc10GtmUjfyQPq80jhejKiVMn8Hedloi8w8C7BFmVMbR1CKZhInckx3P06IoJxSL+daclIu8w8FZB1jqhp5vIv80tyIiI+hEDb5VC4RAy12Vg1A3uu0lE1IcYeGugRBVkrs2gVW2hcIwnphMR9RMG3hqpcdU5Mb2ko3i8yNCjNbNGR/FwpyUi7zDwLoOW1JDa4mw2XTzG0KO1sUfH8OVOS0TeYeBdpnA6jMw1GTSLTeTfznP2Jq2aUKlgd6clIu8w8K6AltKciSw1A7nDOa7To1URj72Dn3RaIvIOA+8KqXEVQ9uGYBkW5t+aR1vn3ptERL2IgeeCUDiEoW1DAIDcoRyMOk9ZICLqNQw8l8iqjOz2LMSQiNzhHPSK7ndJRET0Lgw8Fy1sOB2KhpA/kkej0PC7JOpFqoq3Oy0ReYeB5zJREpG5NgMt7ZycXpur+V0S9Rhz23Zc12mJyDuy3wUMIkEQkN6ShiiLKE2VYLUtxMfjfpdFRBRo7OF1UXJDEon1CVTOVFCaKnGBOgEApINvYLbTEpF32MPrsthYDKIsoniiCLNlIrUlBVHi64xAa7cxDKDY5hIWIi/xN68HItkIhq4bQqvawvyb8zAaXLZAROQ1Bp5H1ISK7PVZCJKA+bfmUc/V/S6JiChQGHgeklUZ2W1ZhNNhFI8XeV8vYCYnJzExMYG99+71uxSiQBJs/sb1RX2+jtJUCaFICOmr05AUye+SyCO5oyfw0Ws34x/fPo6hazb5XQ5RYLCH55NINoKhbUMwDRNzb85xZ5YgicbwcqclIu8w8HykRBUMXz+MUCSE3OEcqjNVv0siD4hnzuBrnZaIvMPA85koOzuzxMZiKJ8uI3+UZ+sNOmF+Dg90WiLyDgOvBwiCgMT6BDLXZNCqtJxjhppco0VE5CYGXg/RUpqzdEEQMPfmHJcuEBG5iIHXYxaOGVpYupA/mudJ6kRELmDg9SBBFJDanHKGOKstzL4xi2ax6XdZ5BI7M4TJTktE3uE6vB5nGiZKUyU0i01EhiJIbEhwL84+l5vJITuexfz0PIbGGHpEXuFvzh4nhSRkrskgtTmFZrGJuYNz0Mtcs9fX6nXc1GmJyDsMvD4RGYpgeGIYsiojdySH4okily/0KentI/inTktE3mHg9RFJkTC0dQjJjUk0C03MvTHHe3tERKvEwOtD0eEohiecHVryR/MovFPgTE4ioktg4PUpSZGQuTaD9JY09IqO2Tdm0cg3/C6LiKhn8cTzPhfOhKEmVJSmSigcK6A2V0PiqgSUqOJ3abQcUUS50xKRd7gsYYDoFR3lk2UYDQPhTBiJ9QkeO9SDuCyByB98iTlA1LhzqnpqUwqtirNgvXKmAtviaxoiIgbegBEEAZFsBMM3DCM6EkV1porZ381yX84eIh06hN91WiLyDgNvQImSiMT6BEZuHIESU1A8XuRBs71Cb+KGTktE3mHgDThJkZC+Oo3sducUhtzhHPJH82jrPH6IiIKFgRcQSlRBdnsW6S1pGHUDc2/MoXyqzN1aiCgwuCwhYMKZMLSUhurZKqozVdRzdURHooiORLkpNRENNAZeAAmigPh4HJFsBNVpJ/hqZ2uIDEcQHYlCCnEpQzdMTk5icnISYV3HegBPbdrsd0lEgcJ1eATTMFGbraE+V4dt2YhkI4iORiGrfD3UDVyHR+QP/kYjSCEJifUJxMZiqM/VUT1bRW2uhnAmjNhYDKFwyO8SB4owexb/Y6cFA4/IM7xpQ4tESURsLIbR94wiuTGJVrWFuYNzyL+dR6va8ru8gSHOzOCRTktE3mEPjy4giAKiw1FEshE08g1UZ6qYPzQPJaYgNhaDltT8LpGIaM0YeLQsQRAQGYogMhRBs9hEdaaK/Nt5hMIhREejCKfDEETB7zKJiFaFgUeroqU0aCkNekVHdaaK4vEiyifLCGfCiGQjCEV4n4+IehsDj9ZEjatQ4yraehv1+ToauQZqczWEwiFEshGEM2GIMm8Nr8ROJvFfAPy3yaTfpRAFCpcl0BWxbRt6WUd9vg695OzTqaU0RLIRKHEFgsAhz/NxWQKRP9jDoysiCAK0pAYtqcFqW6jn6qjP15E7koOkSM49wGyE5/K9W6uF9Z2WiLzDwCPXiLKI2GgMsdEYWrUW6vPOmr7KdAVqXEUkG4GW0gI/0UV6602cAlB8601g47jf5RAFBgOPukKJKlCiCpIbkmgUGqjP11E4VoAgClATqjMJJqnxfh8ReYaBR10liOeWNrSbbTSLTTSLTRSPFwEASkxxhkRTGmSN345E1D38DUOekTUZsbEYYmMxmIYJvaSjWWqiMl1B+XQZsipDS2lQkyqUGCe8EJG7GHjkCykkIZJ1JrTYlg29oqNZbDo7u5ytQpREqEln6FNNqDy6iIiuGJclUM8x6oYz9FlqwqgbEAQBSkxZfAtFQ30dgLkzc1i3fgRnTs9iaN2w3+UQBQZ7eNRzQpEQQpEQ4uviMFsmmqUm9JKO2lwNlemK85hwCKFoyAnBqNJf9/9EEa1OS0Te6aPfEhREkiIhOhxFdDgKAGg322jVWmhVWzBqBurzdQDOkgglqiwJwV5d/iAefRs/67Q8HojIOww86iuyJkPWZESGIgAAy7Rg1IzFEKydraFyptMLjIScEIyEFp/XC8sghFoN/wJAsVbzuxSiQGHgUV8TJRFqQoWaUBc/ZjQMJwSrLegVZyj03Y+XNRmSKjkhqMqL7/fzfUEiujQGHg2cUDi0uJk1ANiWjbbehqmbaDfbzpvehl7WYbWtxedJIWlpGGoypJAEURYhymLPDpES0eow8GjgCaKwGILns0wL7ebSMDTqBhr5Bmxr6QRmQRAWw0+URQjS0vdFWYQoiUsDUjj33IU/c2I0kT8YeBRoouRMdkH0wr8zDROWYcFqn/dmnvuzUTdgtS3Ypg3LtC68yEXUm1H8FYD/bcMGd/8xRLQirsMjcolt2+fCr23BtuxzvTkb+Ma3voEnnnwCpmniyNEjKJVKSCQS/hZNFCAMPCKPVY4dwxeuvhpff+cdxLds8bscosDgtDQijwknT+LJTktE3mHgERFRIDDwiIgoEBh4RB5Tr7sOeOoppyUiz3DSChERBQJ7eEREFAgMPCIiCgQGHhERBQIDj4iIAoGBR0REgcDAIyKiQGDgERFRIDDwiIgoEBh4RL2i1QJqNaDdXv4xjYbzGO4XQbRmDDwiN9g2sH8/YBjAK68AX/868PzzawumRx8FvvpV4MSJ5R/zwx8Cf/RHTjgS0Zow8IiuVL0O/OVfAskk8PDDwFtvAc0mcNddwL//92u71l/+JXDNNUClAhw4AOzbB/z935/7+3/1r4APftDN6okCQ/a7ACI//c3f/A2ee+45vP7661AUBcVi8ZLPsW0blmVBkiSnR/fnfw488ggwOQncdhtwxx3OA5tN4IkngM99bu2F1WqApgHV6spDnES0auzhUaC1Wi184hOfwGc/+9lVP+fpp592wg5whi4/8AHAsoAzZ4Dbbz/3wA99CJiZubzCxsaAnTsBma9JidzCnyYKtC9/+csAgKeeempVj7dtG7/73e+cdywL+OY3gX/6JyCRAP7zf176YMsCxsddrJaIrgQDj2gNjh07hlAo5Lzz+uuAaTr37i7m178G/vAP3/1k4Be/ACQJuOce4Kc/dXqFv//7wObN3S6dKPA4pEm0BjMzM5AXhhnzeWfCysVmYloW8A//APybf3PuY1u2OPfmjhxxQm96Gti1i2FH5BEGHg2chx9+GIIgrPj2q1/96rKv31pYEnDttcDcHHCxa/27fwd89rPApk1LP75nD/D97wM//rFzj2/r1suug4jWhkOaNHDuv/9+3HPPPSs+ZvNl9qrGxsZw6tSphYsAd94JfOpTTm/uhhuAUgl48kkn6P7kTy68QCwG/PN/Dpw9e24253JsGyiXnSUKRHTFGHg0cLLZLLLZbFeuvWXLFpx498Lwp54C/u2/BT7xCSfkbrnF6dmNjV38Aq++CnzsY04o7t27/Cc6eRL4D/8BuOoq4M03gb/5G+CLX3Tzn0IUOAw8CrSpqSnk83lMTU3BNE28/vrrAIBrr70WsVgMALB9+3Y88sgjuOuuuyAIAj7ykY+gWq06f5/JODukPPropT/Z668DguD0/L75TeDw4eWHNDdscEKOiFzDe3gUaH/913+Nm266CV/60pdQrVZx00034aabblpyj+/QoUMolUqL7//VX/0VpqenV/9JbNvp0X3968C2bc7ElRtvdBalr+U6RHRFBNvmLrREPeE//SdnW7I/+7Ple37f+Q7wwgvOUKqieFoeUb9j4BERUSBwSJOIiAKBgUfksXK5DEEQUC6X/S6FKFA4pEnkMdu2UalUEI/HIQiC3+UQBQYDj4iIAoFDmkREFAgMPCIiCgQGHhERBQIDj4iIAoGBR0REgcDAIyKiQGDgERFRIPz/IZy+pu+eAW4AAAAASUVORK5CYII=\n",
      "text/plain": [
Luke Naylor's avatar
Luke Naylor committed
       "Graphics object consisting of 4 graphics primitives"
Luke Naylor's avatar
Luke Naylor committed
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
Luke Naylor's avatar
Luke Naylor committed
    "plot_central_charge(-O, r\"$\\mathcal{O}_X[1]$\", argument=3*pi/2, radius = 0.9)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "8149044b",
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAb4AAAHWCAYAAAAIKuzeAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAA9hAAAPYQGoP6dpAABLNklEQVR4nO3deZScd33n+/dTT+17VVd3q7W1FkuyJC/ItrDlJCbBibEJhEs2GDJOSMAJxNwsXOae49kSyEy451wg5A4ihGQCQ3AyJDHkhhgcPBfMEmSwjeVNsnaptbV6rX2v57l/lLptyVpaUlU9tXxe5/TpRVVPfbvVXZ/67YZt2zYiIiIDwuV0ASIiIp2k4BMRkYGi4BMRkYGi4BMRkYGi4BMRkYGi4BMRkYGi4BMRkYGi4BMRkYGi4BNxgG3bZLNZtH+ESOcZ2rlFREQGiVp8IiIyUBR8Iq00Pd18Pz8PL78M9fqV3f/f/Tv4L//l4v9eLMJ/+k+Xvo2IXJLb6QJE+kI+D//1v8JP/zQ8+SQcPw5f/zokEvCd70A0urTr+P3wH/9j8+OZGfibv4GpKVi9Gt7zHggG4Y/+6JXbiMgVU4tPBt53vvMd3vrWt7J8+XIMw+Af//EfL3ufF1988ZVPcjl43/vgd34H/sf/gJ//efjMZ+CZZ+DECfjrv77yomwb/uRP4P3vbwbdt78NH/zglV9HRF5DwScDr1AocPPNN/OpT31qSbc/cuQI//RP//TKFz7wAfjt34Zdu2BkBDZvbn49lYKf/VnYu/fKi5qchI9+FPbsAcOAX/1V+Nznrvw6IvIa6uqUgXffffdx3333Lfn2n/nMZ165/dGjsG8f3Hkn7N8Pb3vbuTeORiEcvvKili1rdpVu2dL8PJttBqmIXDO1+ESu0K5du9iwYUPzk7/7O7jjjubHGzeCaZ574wMHYMeOVz6fm2u25ubnwbKaH09ONj9+NcOAN70JPJ5mt+enP90cQxSRa6bgE7lCk5OT+P3+5if79oHXe+Ebzs42W4T33vvK1yYm4Md+DHbvbs74/NCH4Ic/bIbbxXz84/DOd8K/+Tet+hZEBpq6OkWuQqPRaH6wbBn8679e+EZ/9Efw3/7bucH4utfBgw/C8883J8V88pOX7sL8+7+HDRuaXah79rzS9SkiV00tPpErtGzZMrLZbPOTX/iF5ozLv/3bV27QaMCf/mmzZfczP/PaC7z73fCJT8C2bZcOvW98A06ehDVr4Ec/gi99qZXfhsjAUotP5Art2LGDXbt2cd1118EttzRD7P3vb47DrV8Py5fDr/96s6V2IS+9BDfe2JwMs2rVhW+TTjcXs+fzzVYjwBvf2JbvR2TQKPhk4OXzeQ4ePLj4+ZEjR9i9ezfJZJLVq1fz0EMPcfLkSb7whS8A8L73vY8//uM/5v7772/e4fd/H37rt5oLzZctay5Cv5innoJkEj7ykebuK3fffeHbxePw3HMt+g5F5NXU1SkD7+mnn2bbtm1s27YNgA9+8INs27aN//yf/zMAp0+fZmJiYvH2a9eu5YEHHjj3IsFgs0vyYqG3sC7vy1+GrVshFoPvfx/+5/9sx7ckIpeg0xlEOqFeh0oFXC4IBKBWg2q1+W+h0Cu3+5Vfaa4JfPDBC1+nVGruCrNvX/O9iFwxBZ9IN8lkmpNjkskL/7tlNcf/PB6IRDpamki/UPCJiMhA0RifiIgMFAWfiIgMFAWfiIgMFAWfiIgMFAWfiIgMFAWfiIgMFAWfiIgMFAWfiIgMFAWfiANqc3Oc+Zu/oTYz43QpIgNHwSfigMqLLzL6K79CZc8ep0sRGTgKPhERGSgKPhERGSgKPhERGSg6gV2kDayGhW3Z2A27+d5qHoKycBhKKT7M7wH/qRaivH8Gl8sFBhguA9Nn4vF7MH0mbr8bl6nXpyKtpOATuYxGrYFVsy78vn5uwC0E3sV8/u8+z+f//vPUyjWOAj/39ZOEfOlLPr7LdGF6TVweF26vG5fXhdvvxhvy4ov68Ea8+GP+V95HvZhus6U/A5F+ovP4ZOBZdYt6uU69UqdertOoNKhX6lg1qxls5/2JuNwuTE8ziFxuFy7TheEyMEwDw2U0P3/1x67mxwCcfXfw/32Sj73jJ/gvT71AZHxtMyztZktx4fHrlbO1lJu11Mt1GtUGjWrza9V8lWq+Sq1Qw7Ksc2r0BDx4Q148YQ/BoSCBoQChkRChkRD+hL/ZwhQZUGrxyUCwbZtGpUGtVHsl3MrNoLMar4SG6W12L3qCnsVwO/+9YRjXXI9v7gx/CeSLc4SGb7ima9m2Ta1Qo5wpU8lWFt+quSqVbIX00TSTz04uhqPLdOFP+AkkAgRTQYKpIKGRENGVUUyvWorS/xR80nds26ZerlMr1s55W+iCdLmbXYVuvxt/3L84lub2uV9pmfUQwzDwhr14w15YceHbWJZFabZEcbpIYbrQ/HimyNSLU1RylcXrBIYChJeFiSyPEF0ZJbYqpjCUvqPgk55XrzS7/WrFGrVCjVrplZBz+9x4Qh78cT+eoAdPwIPLPXjdfC6Xi9BwiNBwiGGGz/m3RrVBfjJP9mSW3Kkc+VN5ZvbMYFlWMwyTAULLQkRXNIMwNh7ThBvpaQo+6Tn1Sr3ZjZdrduc1ag2AxS5Kf+JsyAU9eoJeAtNrElsdI7Y6tvg1q2GRO5Uje+KVMJzdN4tVt3C5XURXRkmsTZBYnyC2WkEovUXBJ13vYkHnDXkJJAN4I81uvl568rUCQZ4AbguHnS7lglymq9m6W3VuGGYmMswfnid9NM3E9yY48q0jmB6T6Moo8TVxBaH0BM3qlK5jWzaVXIVyujlZo1F9Jei8YW9PBt359n5zL1vu3kImkyEajTpdzlVZDMJDzSDMHM9g1S1Mj0lsPEbq+hTDW4fxR/1OlypyDrX4pCs0ag0qmcrizETbsnH7mpNPfFFfzwfda1gW3rPve5XLdDW7O9cmgLNBeCzD3KE55g7OceDRA+z/5/1ExiIMbRxieOsw0RW9GfLSX9TiE8fUSrVm2KXLVAtVALzh5iJsf9yP29+/r8uOfOYR1r7/F8l/+9uE77rL6XLaopqvMvXSFLP7Zpk7NIdVt/BFfc0Q3DJM8rpkf72YkZ7Rv88s0pXqlTql2RKluRL1Sh3DZeCP+YkPx/HH/AM547JfecNeVt6+kpW3r6RRbzC3f47pvdPM7Jvh1NOncPvcpK5PMfq6UYY2DLVkfaTIUij4pO2sukVprhl21UK1uYA67ie6Koov4uvJtXNyZUy3yfCWYYa3NJdSZI5nmHphiqkXp5h8bhJv2MvIDSOM3TKm7lBpOwWftIVt2ZTTZUpzJSrZ5gJpX8xHYl0Cf8yvsBtwCzNGN7x5A/OH5zn97GnOPHeGE0+eIJQKMXJTMwQDiYDTpUofUvBJS1VyFYozRcrpMrZl4w17ia6KEkgE1I0pF5RYlyCxLoHVsJjZO8Pk7kmOfecYR755hNjqGMtet4yxW8YwPdpBRlpDk1vkmll1i+JskeJMkXq5jtvvJpAMEEgGcPv02upC9v7Lc/zMva9jz/Q00VTK6XK6Tq1U48zzZzjz3BnSx9J4/B5Gbh5h5e0rCY9259pH6R0KPrlq1UKV4nSR0lwJAH/CTzAVxBfxOVxZ9+uHdXydUpgucGLXCSafm6RerhNfE2fl7SsZvmFYp0zIVVHwyRWxbZvSXHOz42qhiuk1CQ2HCKaC6sq8Age/+A2evf9N3Lt7N5Gbb3a6nJ7QqDWY3D3JyR+eJHcqhy/iY+zWMVbesVIvtuSKqB9KlsSqWxSmChSmC4vrsZLXJfFFfZqGfhXMfI5fAvKZjNOl9AzTY7Ji+wpWbF9BeiLNiV0nmPjeBBPfmWBo8xBr3rBGM0JlSRR8ckmNamMx8ACCQ82z2/p5cbl0v/jqOPHVcaqFKid/eJKTPzjJU59+iuR1Sda8YQ2JdQmnS5QupmcvuaB6uU7+TJ7SbAnDZRAeDRMaCak7U7qKN+Rl7U+tZfyucU49fYqJ703wo//+I6Iro6z5yTUMbx6+/EVk4Cj45By1Yo38ZJ7SfAnTYxJZEWmO32lrKeliLtPFyttXsuL1K5h6YYpj3z7G8198ntBoiPGfGGf05lFNhJFFmtwiQHNfxdzpHJVsBbfPTWg0RHAoqIXmbbL/H77D537pDfz7/fuJbNjgdDl9aWbfDMe+fYz0sTSBRIDVP76a5duX60WcKPgGXa1UI3cyRzlTxhPwEF4Wxp/wa8JKm2k5Q+ekj6Y5+sRRZg/MEkgGWHv3WsZeN+Z0WeIgBd+AalQb5E7lKM4WcfvcRJZHCCS1PVS77dy5k507d+LNFVl94hhfPHaM6OrVTpc1EDITGQ7+y0HSR9OEl4VZ/6b1pDZq84BBpOAbMFbdInc6R3G6iMvtIjwWJpgKqoXXYYNwLFG3mtk/w+FvHCZ3Okd8TZzr3nQdsdWxy99R+oYmtwwI27LJn8lTONNclhAeCxMeDWsMTwZOamOK1MYUk89Ncvh/HebpP3+a4euHWX/vekLDIafLkw5Q8PU527YpzhTJncphN2xCIyHCy8JaliADb9nNyxi5YYSTPzzJ0SeO8sP/54cse90y1t+7Hm/I63R50kYKvj5WLVTJTGSoFWsEh4JElkcwvdrhXmSBy3Sxascqlt+6nGPfO8bEdyeY3jvNup9ex4rbV2gIoE8p+PqQVbfInsxSnCniCXpIbUrhDesVbDexvF5eAsb9fqdLEcD0mqx74zqW37qcA48eYN9X93Hq6VNs+rlNGv/rQ5rc0kcWuzVP5gAWF5/rVWv30XKG7jZ3cI79X91PYabA2C1jXHfvder+7CNq8fWJ87s1oyujGscTuUrJ65K8/ndez8T3Jjj27WPM7J1h7d1rWXnHSr2Q7AN6ZuxxVsMifSzNzMszAKSuTxFfE1fodTnfwb1kANfzzztdilyEy3Sx5g1ruP13byd5XZL9/7yfpz71FJnjOlGj1+nZsYeVM2WmX5qmPF8mtjpG6vqUumN6hGHZRAEsy+lS5DL8MT83vPMGtv3GNqy6xY8++yMOPnYQq6H/u16lrs4eZDUsssezFGeL+GN+YuMxTI9ma4q0U3J9s/vz0DcOMfG9CWZenmHzL2wmtkqTX3qNWnw9ZrGVly4TXxMneV1SoSfSIS7TxYb7NnDrb94KNmr99Si1+HrEa1p5q2NakyfikNjqmFp/PUwtvh5wTitv/GwrT6HX0yqr13ILYG3c6HQpcpXU+utdWsfXxWzLJnM8Q3GmiC/qIz4eV+D1Ca3j6y9Ww+LQNw5x/F+PExwJcuM7byQ0on0/u5VafF2qXq4z8/IMpbkS8fE4QxuGFHp9xH3mFJ8CjOPHnS5FWuDVrb9GpcFTn36KU0+fcrosuQgFXxcqzhaZ3juNbdukrk8RTAWdLklazJ1J8yBgzM46XYq0UGx1jNf/768nuSHJ3q/s5aW/e4lGteF0WXIeTW7pIlbDIjORoTRXIpgKElsV07FBIj3G4/dw06/cxIknT3Dg6wfInsiy9Z1biS5Xl3a3UIuvS9SKNWZenqGcLpNYmyA+HlfoifSwlXes5Lb33YZhGDzz589wfJe6tbuFgq8LFKYLzLw8g+EyGN48TCAZcLokEWmByFiE7Q9uZ2TrCPv/eT/PP/w8tVLN6bIGnoLPQbZlM39knsxEhmAqSGpTCrdfvc+DoJ5I8gnAHh52uhRpM9NrsvWXt7L55zczd2COp3Y+RX4y73RZA03B55BGtcHMvrNdm+sSxFZrPG+Q1IeX8X8A9ooVTpciHbL81uVsf3A7hsvgmc8+w9RLU06XNLAUfA6o5qvMvDyDVbdIbUoRSKhrc9AYpQJ3AOT1yn+QhIZD3Pb+24iujPLi377IkW8dcbqkgaTg67DibJHZ/bOYPpPhzcN4gh6nS5IO2rlzJ1u2bOH/fO//xi7AdfCg0yVJh3kCHm5+982svGMlh//XYV742xdo1LTkoZO0c0uH2LZN7mSO/Jl8c6nC6pgOtBxgRz7zCGvf/4vkv/1twnfd5XQ54pBTT59i31f3ERoOcdP9N+GP+Z0uaSCoxdcBVsNi7uAchakCsdWx5lIFhZ7IwFt+23K2/cY2KrkKT+18ivkj806XNBAUfG1WL9eZ2TtDrVgjuSFJaFj794nIK+Ljcbb/9nZ8UR+7P7dbW511gIKvjaqFKjP7ZsCA1PUpfBGf0yVJl7BNk2kAt5avSPOU91t/61ZSm1Ps/cpeTXppMwVfm1SyFWb3z+L2u5vr83x6gpNXVNZvYgSwbrjB6VKkS5gekxv/zY2s2rGKw//rMPu+ug9NwWgPPRu3QWmuRPpoGl/UR2JdQuvzRGTJNr5lI96wl0OPH6JWqLHll7bgMtVGaSX9NFssfybP/JF5AskAifUKPbkw79EDHABce/c6XYp0oTU/uYbNb9/M9EvT7P78bp3w0GIKvhbKnsiSPZElvCxMfI1mbsrFuao1rgOoVJwuRbrU8tuWc8O7biAzkeGZzz5DNV91uqS+oeBrAdu2SR9Nkz+TJ7YqRnSFjh8RkWs3vHmYbb+xjXK6zNN//jSluZLTJfUFBd81si2b+UPzlOZKJNYmCI1ouYKItE58PM4tv3kLVt3i6c88rQ2uW0DBdw1sy2bu0ByVXIXkdUkdJyQibREeCXPb+27DHXDz7H9/ltzpnNMl9TQF31WyLZu5g3NU81WS1yXxRbVGT5auumI1bwKsdeucLkV6hD/m59YHbsUT9PDsf3+W7Mms0yX1LAXfVVgMvUKVoQ1DWpguV8wKhfkGQFTjwbJ03rCXWx64BW/Ey+7P7SZzIuN0ST1JwXeFzg89b9jrdEnSg9yz0/wBYExOOl2K9Bhv2Mst770Fb9TLc59/Ti2/q6DguwIKPWkV9+w0f4iCT66ON/RK+O3+3G6F3xVS8C3RwkQWhZ6IdANv0Mstv3EL3rCX3Z/frdmeV0DBtwS2bTN/eH5xIotCT0S6wWK3Z9DLs3/1LPkphd9SKPiWIH00TSVbIbk+qYksItJVvGEv2967DdNn8tznn6OcKTtdUtdT8F1G5niG0lyJ+Nq4lixIyzQiUb4I2PG406VIH/BFfGz79W1YDYvdf7WbWqnmdEldTcF3CfnJ/OKp6YGEFqdL69TGVnI/YK9Z43Qp0icCyQA3/9rNVHIVnvsfz9GoaWPri1HwXURxtkj2ZJbIWESnpkvLGdUK6wHK6paS1okuj3Ljr9xI7nSOF//2RSzLcrqkrqTgu4BypkzmWIZgKkhkecTpcqQP+Y4e4iDgevllp0uRPpNcn2TLL2xhdv8sL39Fv18XooNoz1MtVJk/PI8v5iO2OuZ0OdJndu7cyc6dO7luPss/OV2M9K3Rm0apFqrs/+f9eMNernvTdU6X1FUMW2fbL6qX68y8PIM74GZow5AOkZW2OfKZR1j7/l8k/+1vE77rLqfLkT516PFDHH3iKBvevIHVP7ba6XK6hlp8Z1l1i9kDs5hek+R1SYWeiPS89T+znmq+ysGvH8Qf9zOydcTpkrqCxvhoLlCfOzSHbdkkr0viMvVjEZH+sOltm0isT7D3kb3a3eUsPcMDmYkMtUKN5Pokptd0uhwZAOWNWzAA63Wvc7oU6XMul4sb3nkD3rCX57/4PNVi1emSHDfwwVeYLlCcKRJbHdNWZCLSlzwBDzfdfxO1Yo0XHn4BqzHYyxwGOvgquQrZ41lCIyGCqaDT5cgA8R4/wvcB14EDTpciAyI0HGLrO7aSOZbhwKOD/Xs3sMFXr9SZPzyPN+IlulKHgUpnuUoldgAUCk6XIgMktSnFup9Zx4kfnODED044XY5jBnJWp9WwmD80j8t0kViXwDA0g1NEBsOaN6yhcKbAgX8+QGgkRGJtwumSOm4gW3zpo2ka1YZmcIrIQLr+568nvDzMi3/7IqX5ktPldNzAPevnJ/OU02Xia+O4/QPZ4BWRAWe6TW76tzdhuAxe+NvBm+wyUMFXzVfJncoRGYvgj/mdLkcGWHXZcv4tYK3WbhriDF/Ex9Zf3kr+VJ6DXz/odDkdNTDBZ9Ut5o/M4wl5CI+FnS5HBpwVjfMwQDLpdCkywBLrEqz5yTUc33Wc6b3TTpfTMQMTfOljaWzLJrFWk1nEeWZ6jt8GjJkZp0uRAbf27rXE18TZ+8jegTm9fSCCrzBVaI7rrYlrZxbpCp6pSXYCxonBnVIu3cEwDG545w0YLoMX/+dgnOHX98FXK9bInsgSHg1rXE9E5AJ8ER+bf3EzmYkMhx8/7HQ5bdfXwWc1LOYPz+MJeois0IGyIiIXk9qYYvWPr2biuxPMHph1upy26uvgyxzLYNUt4mvjGtcTEbmM9fesJ7oyyp6/30M137+bWfdt8BVni5TmS8TGY7h9Wq8nl/bpT3+atWvX4vf7ufXWW/nud7970ds+8cQTGIbxmreXX355yY/XCAb5F4CIeiKke7hMF1vfuRXbstnzD3ucLqdt+jL4GrUG2eNZAskAgUTA6XKky33pS1/i937v9/gP/+E/8Oyzz/ITP/ET3HfffUxMTFzyfvv27eP06dOLbxs2bFjyY9ZWruFewFq//hqrF2mtQDzAhp/dwOyBWU49c8rpctqiL4MvcyyD4TKIrYo5XYr0gE984hO85z3v4b3vfS+bN2/mk5/8JKtWreLP/uzPLnm/kZERli1btvhmmlcwY7jRIHL2vUi3Gds2Rur6FAe/fpBytv+WOPRd8BVni5QzZWKrY7jcffftSYtVq1WeeeYZ7rnnnnO+fs899/D973//kvfdtm0bY2Nj3H333XzrW9+6osf1H9pHFnC98MKVlizSEZvfvhmAl7+89C78XtFXyfDqLk5/XEsX5PJmZmZoNBqMjo6e8/XR0VEmJycveJ+xsTE++9nP8sgjj/DlL3+ZTZs2cffdd/Od73znoo9TqVTIZrOLb8VysaXfh0irecPevu3y7KtZH+rilKt1/qxf27YvOhN406ZNbNq0afHzHTt2cPz4cT72sY9x1113XfA+H/3oR/nwhz+8+Pk24EfXXrZIW41tG2PqxSkOfv0gyQ1J/NH+aFD0TYtvsYtzXF2csnSpVArTNF/TupuamnpNK/BS7rjjDg5c4jT1hx56iEwms/j213/yhauuWaST+rHLsy8SYqGLMzgU1O4sckW8Xi+33norjz/++Dlff/zxx7nzzjuXfJ1nn32WsbGxi/67z+cjGo0uvgX9wauuWaSTvGEvG9+6sa+6PPuiq3OhizO6Kup0KdKDPvjBD3L//fdz2223sWPHDj772c8yMTHB+973PqDZWjt58iRf+EKzlfbJT36SNWvWsHXrVqrVKl/84hd55JFHeOSRR5b8mOV1GxgGDm3d2o5vSaSllt28jDPPn+HAowdIbUrhDXudLuma9HzwlTNlypkyiXUJnaYuV+Ud73gHs7OzfOQjH+H06dPccMMNfO1rX2N8fByA06dPn7Omr1qt8qEPfYiTJ08SCATYunUrjz76KG9+85uX/qBuDzMAHk9rvxmRNtn89s3s+sQuDnztAFt/ubdfsBm2bdtOF3G1bMtmes80ps9kaMOQ0+WILNnBL/4Le+6/l5969lkir3ud0+WILMmx7x3j4NcPcst7byGxNuF0OVetp5tI+TN5GtWGZnFKzzHzeX4OMLJZp0sRWbJVd64iNBpi31f39fTxRT0bfPVKnfzpPKHREG5/z/fYioh0PZfLxaa3bKJwpsCJXb17lmTPBl/2eBaXx0VkTJv8ioh0SmJdgtGbRjny/x3p2RMcejL4Fia0RFdGMVw6bkhEpJM2vLm5IfuBr1187Wo367ngsy2b7PEsvqhPJy9Iz6qlRvggYC9f7nQpIlfMF/Gx5qfWMPncJOmjaafLuWI9F3ya0CL9oJFM8SeAPTLidCkiV2XVnasIjfTmRJeeCr5GrUF+UhNapPe5chl+EWB+3ulSRK6Ky3Sx8a0byU/mOfVUb+3o0lPBlz+dx3AZhJeFnS5F5Jp4T5/k7wHXsWNOlyJy1ZLrkqSuT3H0W0dp1HrnbMmeCb56uU5xpkh4WVg7tIiIdIn1P7Oear7KxL9OXP7GXaJnEiR3KofL4yI0HHK6FBEROSu8LMzoTaMc/95xaqWa0+UsSU8EX61YozRfIrI8ouULIiJdZt0962jUGhx94qjTpSxJTwRf9kQWt99NIKnlC9IfLJ+veRBtQL/T0vsC8QDLb1vOiSdPUM6WnS7nsro++CrZCpVcheiK6EVPxBbpNdXx9dwKWK86yV2kl629ey0u08WRx484XcpldX3wZU9m8Ya9+OM6YFZEpFt5g15W3bmKyd2TFKYLTpdzSV0dfKX5ErVijegKHTAr/cV3YA9lwPXcc06XItIy43eN4w64OfSNQ06XckldHXy5Uzn8MX/Pn/Yrcj7DBh9A7x6HKfIaptdk/A3jTO+ZJj+Zd7qci+ra4Cuny9TLdS1Wl76yc+dOtmzZwu9++HedLkWkLVbcvgJvyNvVMzy7Nvjyk3m8Ya9ae9JXHnzwQfbs2cOf/sGfOl2KSFuYbpOVO1Yy/dI0pbmS0+VcUFcGXyVXoVqo6qw9EZEetOrOVbi8Lo5++6jTpVxQVwZffjKPJ+jBF/U5XYpIW1TG17EVLWeQ/uT2uVm+fTmTuye78rDargu+WrFGJVvR2J70NdvnZw9oAbv0rfEfHwfg2He7byP2rgu+/GQet8+tdXvS19xnTvIXgDHROxv7ilwJb9jL2C1jnHrqFLVyd+3h2VXBVy/XKc2XCC8La5cW6WvuTIb3AsbcnNOliLTN+BvGsWoWJ3adcLqUc3RV8OUn85gek8CQun9ERHpdIB5g+IZhTuw60VXn9XVN8DVqDUpzJUKjIbX2RET6xJqfXEO1UOX0j047Xcqirgm+4kwRDAimgk6XIiIiLRIeDZNYm+DkD086Xcqirgg+27YpThcJDgV1uroMhHpiiI8C9siI06WItN3KO1aSn8yTPpJ2uhSgS4KvnC7TqDUIDqu1J4OhPjzKvwfs5cudLkWk7Ya3DuOL+jj+g+NOlwJ0SfAVpgr4Ij48AY/TpYh0hKtY4A0AuZzTpYi0nWEYLL9tOTN7ZqjkK06X43zw1Uo1qvmqWnsyULwnjvEE4DrU3ce3iLTKittXAHTFWJ/jwVecLmJ6TC1YFxHpY76wj9TmFKefOY1lWY7W4mjw2ZZNaa5EYCigJQwiIn1uxe0rKKfLzLw842gdjgZfaa6E1bAIDYecLENERDoguS5JaDjEqR+ccrQOR4OvOFPEH/Njek0nyxDpONvt5gSARxO6ZLAs376cuUNzlOadO6vPseCrl+tUC1VtTyYDqbJuI6sAa+tWp0sR6aixW8YwTIPJ3ZOO1eBY8JXmSrhMF/6YJrWIiAwKT8BD8rokZ54741gNjgafP+7HcGlSiwwe3+H9HAdcL73kdCkiHbfs5mUUpgtkT2UdeXxHgq9aqFKv1NXNKQPLqNdZCVDrrnPKRDohtSWF2+9m8kfOdHc6EnyluRKmx8Qb9jrx8CIi4iDTbTK8eZipF6ewbbvjj9/x4LPts2v3klq7J4Nn586dbNmyhd/98O86XYqIo5ZtW0YlV2HuUOcPY+548FWyFay6RSCpbk4ZPA8++CB79uzhT//gT50uRcRRiXUJfFEfk892vruz48FXmivhCXjwBLV+SQZXdeU4PwlY69c7XYqIIwzDYPTGUWb2znT8dPaOBp9t2ZTTZbX2ZOBZwRDfBohEnC5FxDHLti2jXqkzs7ezW5h1NPjKmTK2ZeNPaO2eDDb39Bn+GDBOObt1k4iTImMRQqkQZ17s7Jq+jgZfJVPBE/Dg9rk7+bAiXcc9P8tDgDE15XQpIo4aun6I+YPzWI3OndjQseCzbZtypqzjh0REZNHwlmHqlTrzh+c79pgdC75aoYZVt/DFfJ16SBER6XKx1TG8IW9Hx/k6FnzldLm5aD2kResiItJkGAbJ65LM7p/t2GN2LvgyZbX2RM6qx2L8JWAnk06XIuK44S3DlOZL5CfzHXm8jgRfvVynXq5rfE/krProCh4A7NWrnS5FxHFDG4dwuV1MvdSZyV4dCb5ypozhMvBF1OITATAqZbYAlJw7jFOkW5hek8TaRMe6OzsTfOkyvohPRxCJnOU7dpiXANe+fU6XItIVUtenyJ3MUc1X2/5YbQ8+q2FRzVfVzSkiIheV2pLCtm2m90y3/bHaHnzVXDO9fVF1c4qIyIX5o35Cw6GOrOdre/BVchXcPjem12z3Q4mISA+LjcdIH0u3/XE60uLTgbMi57INqADoTEqRRcn1SSrZCqW59k76amvwWXWLWqmGN6LgE3m1yoYt+AHr5pudLkWkayTWJwCYPdDe2Z1tDb6F2TlaxiAiIpfjDXkJjYRIH0m39XHaGnwa3xO5MO+xQzyDljOInC82HiN9NN3Wx2hvi0/jeyIX5KpUuAW0gF3kPMl1SSq5CsXZYtseo23Bp/E9ERG5UonrmuN8cwfn2vYYbQs+je+JiMiV8ga9hEbbO87X1uAzvabG90ReZefOnWzZsoXf/fDvOl2KSNeKrYqROZ5p2/XbF3yFqs7eEznPgw8+yJ49e/i/P/MP/BJgjY87XZJI14mujFJOl6kW2rNvZ9uCr1as4Ql62nV5kZ5mRWL8A0Ai4XQpIl0ntioGQPZ4ti3Xb0vw1ct1bMtW8IlchDk3w+8DxlRnzh8T6SWh0RCm1yR7soeCr1asASj4RC7CMzPFJwDj1CmnSxHpOoZhEBoJkTuVa8v12xZ8bp8bl7sjx/2JiEifiSyPkJ/Mt+XabUmmaqGq1p6IiFy16Ir2TXBpzxhfqa7gExGRqxZdGQXaM8Gl5cFXL9exGpaCT+QSGuEw/wTY0ajTpYh0pdBoCNPTngkuLQ8+TWwRubza8tW8DbDXrXO6FJGuZBgGodH2THBpffCVapheUxNbRC6lXiMFUKs5XYlI1wovC1OYKrT8um3p6nT73a2+rEhf8R8+wDTgeuklp0sR6VrBVJBKuoLVsFp63ZYHX6PSwO1T8ImIyLUJpoJYlkVprrXHd7U0+GzbVotPRERaIjwaBmh5d2dLg8+qWdi2reATEZFr5k/4cbldFGdaeyhtS4OvXq4DYPp0FJGIiFwbwzAIxAPdH3yGYegMPpHLKK/fRBSwbrzR6VJEulpgKEBptovH+OqVOqbPxDCMVl5WpP+YJrmz70Xk4gLJQHdPbtHEFpGl8Zw4ymOA69Ahp0sR6WrB4SCVXIVGtdGya7Y0+BqVhoJPZAnMYpE3AeTac+yKSL8IDYcAKEy3bmZny4LPtm3qlbrW8ImISMsEh4MALR3na1nwWfXmynqXR1uViYhIa3jDXgzDoJwtt+yarQu+WjP4TI8G60VEpDUMw8AT9lDNte5cvpb1SzZqzYFHtfhELm7nzp3s3LmTULbI64H/a+VKp0sS6Xq+kK+lwdfyFp9OZRC5uAcffJA9e/bwhS98nU8DdirldEkiXc8T9lDNd2HwNWoNTI/W8IkshSub5lcA5uacLkWk6/kiPir5Ssuu19IWn7o5RZbGO3mKLwKuiQmnSxHpet6Il1qudWdXtrzFJyIi0kq+qI9qsdqyc/nU4hMRka7mi/gAqORa093ZuhZfVS0+ERFpPV/0bPBluiz4rLqlGZ0iS2QFAuwCCIWcLkWk63nDXoCWzexsSVLZto1t2ximZnSKLEV11VruBKwNG5wuRaTruQPNJeet2qi6NcHXsJsXM9XiExGR1lo447VR6abgs5rBZ7jU4hNZCv/+PdiAa/dup0sR6Xou04XL7aJeqbfmeq24yMIUU3V1iohIO5heUy0+EREZHKbX1BifiIgMDtNjdldXp1p8IiLSTqavy1p8GuMTuTKVNeu5DrCuv97pUkR6QleO8RmGoZMZRJbI9vo4BOD3O12KSE/ovjE+y1Y3p8gV8Jw+wV8DxtGjTpci0hNMj0mj3kXBhw0o90SWzMxl+beAkU47XYpITzBcBrTmcIbW7dUpIiLSC1q2V6fG90REpF0Ml4FltabJ517KjWzbJpfLXfTfc7kcpXyJQDbQkqJE+k2lUqFSeeVIlZn0DGuBM4fPMJfQKewilzM3M0elXCGbzV7ydpFI5LINMcO2bftyD5jNZonFYldWpYhc1DLgt4A/ByYdrkWkn2QyGaLR6CVvs6Tgu2yLbzLHXW+6i2efe/bKq1yC7du389RTT+naNF+ErFq1iuPHj1/2P/dq9NrPo93XbtV1z2/xnT59mte//vU8/+zzrFi+4pqvf74dP76DXd/b1fLrtvPavVhzNpdl7XVrOXLwCNFI6/8ee/Fn0raf9YksP/0LP82PnvvRJW+3lBbfkro6DcO45JOsUTAwXWZbnogBTFPXPl80Gm3LtXv159Gua7ezZoBEKkFyJNny65pusy3Xbee1e7Fmt7/5FJocTrbn968Hfybtuq6r5GpZzrRkcothGPz6L/96Ky51QQ8++KCu3SG9+vNo17V78f8Q9LPu5LXbqRd/Ju2s+Tfe9Rstuc6SujovJz+ZJ38mz7Kbl7WiJrmEhfHWpfRjS/c6ceLEYpf1ypUrnS5HrpL+HjsnfSxNvVQndX3qmq/VmnV8Bs1F7NJ2Pp+PP/iDP8Dn8zldilyDhf8//T/2Nv09dlALN0ppSYuvMF0gM5Fh+a3LW1GTSN9TS0HkyswdmsO2bIY2DF3ztVrS4ls4h2/heCIREZFWsi27ZWe+tmZyy9kNqhV8IiLSDnajdYchtCb4zp7Dt3Aun4iISCvZlt2yM1/V4hMRka5nNazuavEtjvE1FHztMD8/z/33308sFiMWi3H//feTvsxxNu9+97sXDwdeeLvjjjs6U7BIH/n0pz/N2rVr8fv93HrrrXz3u9+96G2feOKJ1/zdGYbByy+/3MGK+1Mrx/iWtHPL5SyksLo62+Nd73oXJ06c4LHHHgPgN3/zN7n//vv56le/esn73XvvvXzuc59b/Nzr9ba1TpF+86UvfYnf+73f49Of/jQ/9mM/xp//+Z9z3333sWfPHlavXn3R++3bt++c2brDw8OdKLevtXKMrzXBZ6qrs1327t3LY489xpNPPsntt98OwF/8xV+wY8cO9u3bx6ZNmy56X5/Px7Jl2lRA5Gp94hOf4D3veQ/vfe97AfjkJz/Jv/zLv/Bnf/ZnfPSjH73o/UZGRojH4x2qsv/Ztt08/q4rx/jU1dlyu3btIhaLLYYewB133EEsFuP73//+Je/7xBNPMDIywsaNG3nggQeYmppqd7kifaNarfLMM89wzz33nPP1e+6557J/e9u2bWNsbIy7776bb33rW+0scyAsZEtXjfEZhtE8JFBdnS03OTnJyMjIa74+MjLC5OTFD7S57777ePjhh/nmN7/Jxz/+cZ566ine+MY3nnNCgIhc3MzMDI1Gg9HR0XO+Pjo6etG/vbGxMT772c/yyCOP8OUvf5lNmzZx9913853vfKcTJfethWzpqjE+ANNjYtUVfEv1h3/4h3z4wx++5G0WjsO50BEblzv1/h3veMfixzfccAO33XYb4+PjPProo/z8z//8VVYtMnjO/zu71N/epk2bzhl+2LFjB8ePH+djH/sYd911V1vr7GcL2eLydFnwuTwurJqCb6k+8IEP8M53vvOSt1mzZg3PP/88Z86cec2/TU9Pv+aV6KWMjY0xPj7OgQMHrrhWkUGUSqUwTfM1rbupqakr+tu74447+OIXv9jq8gbKQraYHrMl12tpi69Ra7Tqcn0vlUqRSl1+l/EdO3aQyWT44Q9/yOtf/3oAfvCDH5DJZLjzzjuX/Hizs7McP36csbGxq65ZZJB4vV5uvfVWHn/8cd7+9rcvfv3xxx/nbW9725Kv8+yzz+rv7ho1ag0Mw8Dl7qIty0AtvnbZvHkz9957Lw888ABPPvkkTz75JA888ABvectbzulSuf766/nKV74CQD6f50Mf+hC7du3i6NGjPPHEE7z1rW8llUqd8wcsIpf2wQ9+kL/8y7/kr/7qr9i7dy+///u/z8TEBO973/sAeOihh/jVX/3Vxdt/8pOf5B//8R85cOAAL730Eg899BCPPPIIH/jAB5z6FvqCVbNa1s0JavH1hIcffpjf+Z3fWZxd9nM/93N86lOfOuc2+/btI5PJAM1Tw1944QW+8IUvkE6nGRsb46d+6qf40pe+RCQS6Xj9Ir3qHe94B7Ozs3zkIx/h9OnT3HDDDXzta19jfHwcgNOnTzMxMbF4+2q1yoc+9CFOnjxJIBBg69atPProo7z5zW926lvoC41ao2XdnNCiY4kAirNF0kfTjN0ydslJFyKiY4lErsTcwTkwILk+2ZLrtaztuJDG6u4UEZFWalRb2+Jr6RgfoO5OkUvYuXMnW7ZsYfv27U6XItIzGrVGS8f4WtbVadUtJp+bJLk+iT/ub8UlRfqWujpFlsa2bU7/6DTx8TjBVLAl12xdi8/twjAMGlW1+EREpDUWhs9a2eJr3ZUA02dSr9RbeUkRERlgC5ni9rVsEUJrg8/td1MvK/hERKQ16uU6hmFg+rpwcgs0E7lRUVeniIi0RqPSwPSaLV0m1/oWX6VOi+bLiIjIgKuX67j9revmhDYEH6BWn4iItETXB99CH6zG+URE5FrZtt1cvN7C8T1odfB5TFymS8EnIiLXrFFpYNt2d7f4QEsaRESkNdqxlAHaEHxa0iAiIq1QL9cxXAamt4u7OkHBJyIirdGOiS3QhuDzBDxYdUtbl4mIyDWpFWt4Ap6WX7f1wRdqFlkr1lp9aRERGRC2bVMv1RczpZVaP7nFY2J6TAWfiIhctXqpuRmKJ9gDwQfgCXqoFqrtuLSIiAyAaqGKYRi90dUJzeBTi09ERK5WrVjD7XdjuFq3R+eCtgWfJriIiMjVqhVrbenmhDYGH2iCi4iIXLnFiS29FHym18Tldin4RETkii1ObGnDjE5oU/ABeENeBZ+IiFyxhexox8QWaGPwaWaniIhcjWqhiifgacvEFmhniy/ixapb1Epq9Yks2LlzJ1u2bGH79u1OlyLStaq5Kt6It23XN+w2HZduWzaTuyeJrowSGgm14yFEelY2myUWi5HJZIhGo06XI9I1GtUGZ144Q3J9En/c35bHaFuLz3AZeEIeKrlKux5CRET6zEJmeMPta/G1LfgAfBEf1bzG+UREZGmq+eb4nsvdvnhqa/BpnE9ERK5Eu8f3oN3BF/JiGAbVnFp9IiJyaY1qg3qlji/ia+vjtDX4NM4nIiJL1YnxPWhz8IHG+UREZGmqufaP70EHgk/jfCIishTVfPvH96ATwRfyYrgMKhl1d4qIyIXVy/Xm+F60veN70IHgM1wGvqiPcqbc7ocSEZEeVc6Um3nR5okt0IHgA/DH/FTzVay61YmHExGRHlNOl/FFfG3bn/PVOhJ8vlgzwdXqExGR81l1i2q+2rYtys7XkeAzPSbekJdyWsEnIiLnWmgULTSS2q0jwQfgj/upZCvYVlv2xBYRkR5VyVTwhryYHrMjj9ex4PPFfNiWrTV9IiKyyLZtKtlKx1p70MHg8wQ8uH1udXeKiMiiaq6K1bA6Nr4HHQw+aLb6NMFFREQWlDNlTK+JJ+Dp2GN2NPj8cT+NaoNaUbu4iIhIcxmDP9a51h50OPi8YS8ut4vSXKmTDysiIl2omq/SqDYIJAMdfdyOBp9hGASSAUpzJWxbsztFRAZZcbaI6TXbfhrD+ToafACBZIBGraEz+kREBpht25Tnyx1v7YEDwecNeXH73OruFBEZYJVMBathERwKdvyxOx58AIGhAKX5khazy8DZuXMnW7ZsYfv27U6XIuKo0lwJT9CD2+/u+GMbtgODbfVKnakXp0isSxBIdL6ZK+K0bDZLLBYjk8kQjUadLkeko6yGxZnnzhBZESE8Gu744zvS4nP73HhDXkqz6u4UERk05fkytm07Mr4HDgUfNCe5VLIVHVUkIjJgSnMlfBFfx/bmPJ+jwQdQmlerT0RkUDSqDSq5CoEh54a5HAs+l9uFL+qjOFN0qgQREemw4mwRw2V0dG/O8zkWfADB4SC1Yo1qQWv6RET6nW3bFGeKBJIBXKZz8eNo8PmiPkyvqVafiMgAqGQrNKoNQsMhR+twNPgMwyCYClKaK2E1NMlFRKSfFaeLeIIePMHOncRwIY4GH0AwFQQbLW0QEeljjWqDcqbseGsPuiD4TI+JP+6nMF1wuhQREWmTwnQBl+lybO3eqzkefAChkRD1cp1KtuJ0KSIi0mK21ZzUEkwFMVyG0+V0R/B5w148AY9afSIifag0X8KqWwSHO78h9YV0RfBBs9VXTpdpVBtOlyIiIi1UmCrgj/lx+zq/IfWFdE3wLazrKEyp1Sci0i+q+Sq1Yq1rWnvQRcFnuAxCIyEK0wXt3yki0ifyk3k8AQ/+mHM7tZyva4IPmt2dgMb6RET6QK1Yo5wpE17W+aOHLqWrgs/ldhFMBSmcKeiQWhGRHpefzOP2ufEnuqe1B10WfADh0TC2ZavVJyLSw+qVOqX5EqHREIbh/BKGV+u64DO9JoFkoNnq6/zh8CIi0gL5yTymxyQ41D2TWhZ0XfABhJeFadQa2sZMRKQHLTx/h0ZCXbFg/XxdGXxuvxt/3E9+Mq9Wn4hIjymcKWC4jK5awvBqXRl8AJGxCPVKnXK67HQpIi2zc+dOtmzZwvbt250uRaQtrLpFYbpAaCTk6Jl7l2LYXdykmj0wi1WzGN4y7HQpIi2VzWaJxWJkMhmi0ajT5Yi0TO5UjvyZPKM3juJyd2fwdWdVZ0WWR6iVapTmNNYnItLtGrUG+TP5ZmuvS0MPujz4vCEv/rif7MmsxvpERLpc/nQewzC6bsH6+bo6+ACiK6I0qg2K00WnSxERkYuoV+oUZ4qEl4W7dmxvQXdXR3OGZzAVJHc6h9XQHp4iIt0odyqHy+1a3Hqym3V98EFzhqfdsHVyg4hIF6oVm3MxIssjXblu73w9EXym1yQ0EiI/mdfJDSIiXSZ7Movb7yYwFHC6lCXpieCD5m4uhmGQO51zuhQRETmrkqtQyVaIroh23Z6cF9MzwedyuwgvC1OcLlKv1J0uR0REgNzJ3OIM/F7RM8EHLK4NyZ1Uq09ExGmluRLVQpXIiojTpVyRngo+w2UQXRmlNF+ikq04XY6IyMCyGhbZE1kCiQC+iM/pcq5ITwUfQCAZwBv2kjme0aJ2ERGH5E/nsRoW0ZW9t+VezwUfQGx1jEalQeGMljeIiHRarVSjMFUgMhbB9JpOl3PFejL4PAEPweHmovZGteF0OSIiAyV7PIvpMwmNdv9i9QvpyeCD5gbWLtNF9kTW6VJERAZGaa5EJVchtirWM8sXztezwecyXZroIiLSQedMaIn21oSWV+vZ4ANNdBER6aRentDyaj0dfKCJLiIindDrE1pereeDzxPwEBoJkTudo17Wji5yZebn57n//vuJxWLEYjHuv/9+0un0Je/z7ne/G8Mwznm74447OlOwiANs2yZ9NN3TE1pereeDD5oTXUyPSfpoWl2eckXe9a53sXv3bh577DEee+wxdu/ezf3333/Z+917772cPn168e1rX/taB6oVcUbhTIFasUZ8TbxnJ7S8mtvpAlrBcBnE18SZ2TdDYapAeLS7T/+V7rB3714ee+wxnnzySW6//XYA/uIv/oIdO3awb98+Nm3adNH7+nw+li1b1qlSRRxTK9XIncoRXhbGG/I6XU5L9EWLD8Ab9hIeDZM7qS5PWZpdu3YRi8UWQw/gjjvuIBaL8f3vf/+S933iiScYGRlh48aNPPDAA0xNTV3y9pVKhWw2e86bSLd7dRdnZKy39uO8lL4JPjjb5elVl6cszeTkJCMjI6/5+sjICJOTkxe933333cfDDz/MN7/5TT7+8Y/z1FNP8cY3vpFK5eLLaj760Y8ujiPGYjFWrVrVku9BpJ3O6eLsgQNml6qvgm+hy7NaqOq09gH2h3/4h6+ZfHL+29NPPw1wwfEK27YvOY7xjne8g5/92Z/lhhtu4K1vfStf//rX2b9/P48++uhF7/PQQw+RyWQW344fP37t36hIG/VjF+eCvhjje7VXd3n6Y37c/r77FuUyPvCBD/DOd77zkrdZs2YNzz//PGfOnHnNv01PTzM6OrrkxxsbG2N8fJwDBw5c9DY+nw+fr3cX/Mpg6dcuzgV9mQqR5RHKmTLpo2mGNg31xSwkWbpUKkUqlbrs7Xbs2EEmk+GHP/whr3/96wH4wQ9+QCaT4c4771zy483OznL8+HHGxsauumaRbpKfzFMr1khdn+qrLs4FfdXVucBwGcTHm12e+dN5p8uRLrV582buvfdeHnjgAZ588kmefPJJHnjgAd7ylrecM6Pz+uuv5ytf+QoA+XyeD33oQ+zatYujR4/yxBNP8Na3vpVUKsXb3/52p74VkZZZeN7sxy7OBX0ZfNDs8owsj5A7naOS016ecmEPP/wwN954I/fccw/33HMPN910E3/91399zm327dtHJpMBwDRNXnjhBd72trexceNGfu3Xfo2NGzeya9cuIpH+6xKSwWI1LNJH0niCHiLL+/f32bD7ePqjbdvMHZijXq4zvGUYl7tvc156TDabJRaLkclkiEZ7e99D6R/zh+epZCukNqdw+/pyJAzo4xYfNGfsxdfGsW2b+SPzTpcjItK1CtMFSvMlYuOxvg496PPgAzA9Jom1CSrZCvlJjfeJiJyvVqqRPZ4lNBwikAg4XU7b9X3wAfiiPsLLwuRO5ajmq06XIyLSNWzLZv7wPO6Au+ePG1qqgQg+aC5x8IQ8zB+Zx6pbTpcjItIVMhMZGtUGibWJvly6cCEDE3yGYZBYm8Bu2KSPpZ0uR0TEccXZIsXZIvHx+EBt9jEwwQdgek3ia+KU02WN94nIQKsVa2QmMgRTQQLJ/h/Xe7WBCj4Af9xPZCxC9mSWcqbsdDkiIh3XqDWYOzSH2+8mtirmdDkdN3DBB83xPn/cT/pIWkcYichAse3mZBZsSF6XHJhxvVcbyOADiK+JY3pN5g7OabKLiAyMzLEMtUKNxPoEpsd0uhxHDGzwuUwXifUJrIbF/JF5nd8nIn2vMFWgOFskNh7r2304l2Jggw/A7XOTWJegmquSPaETsUWkf1WyFTLHM4RHwwSHgk6X46iBDj4AX8RHdFW0+Upopuh0OSIiLVev1Jk/PI8v6iOyon83n16qwVm4cQmh4RD1Up3MRAa33403PLhdACLSX6yGxdzBOVweF4l1CZ1Pilp8i6KronjDXuYOzmmmp4j0BduymT80j1WzSK5P4jL1lA8KvkWGYTRnOXlNZg/M0qg2nC5JROSqLZxKUy1USV6XHKidWS5HwfcqLtNFckMSgNkDs1rmICI9K3s8SzldJrE2oeGb8yj4zmN6TIY2DGHVLeYOzWFbWuYgrbNz5062bNnC9u3bnS5F+ljudI7CdIH4eBx/3O90OV2nr09gvxbVQpXZ/bP4Ij4S6zUgLK2lE9ilXQrTBTITGaIrooSXhZ0upyupxXcR3pCX5Ppkc+3LRMbpckRELqucLpOZyBAaCSn0LkHBdwm+qI/4mjjFmSLZk1rgLiLdq5KrMH94nkAiMJAbT18JTfO5jEAyQKPWIHsii8vtIjyqV1Ei0l1qxRrzh+bxhr3E18adLqfrKfiWIDwaxm7YZE9kMQyD0EjI6ZJERIBm6M3un8Xtd2s+whIp+JYosjyCbdtkjmfAaO72IiLipFqpxuyBZuglN2iB+lIp+K5AdEUUbBYnuyj8RMQptVKzpWd6TYXeFVLwXaHoymiz5afwExGHvDr0hjYMKfSukILvKizMmMpMZMBGY34i0jG1YrN7czH03Aq9K6Xgu0qxVTEMwyBzPINt25rtKSJt9+qJLOrevHoKvmsQXRkFA7InstiWTWRM51yJSHtU81XmDs4p9FpAwXeNoiuiGC6D3KkcVt3SwlERablypsz84Xm8IS+J9QmF3jVS8LVAZCyCy+0iM5HBqlnE18a1lkZEWqI4UyQzkcEX85FYm8Bw6bnlWullQ4uEhkMk1iUop8vMHZjDauhIIxG5NvnJPOljaYKpYPP0dIVeSyj4WiiQCJDckFwcgNZ5fiJytTLHM2RPZoksjxBbHVMvUgsp+FrMF/ExtGkIq2Yx8/IM9Urd6ZJEpIcsnJxemCoQWx3TpLk2UPC1gSfgYWjTEACz+2apFWsOVyQivcBqWMwdnKM8Xya5PqkNMtpEwdcmbp+b1PUpXB4Xs/tnKWfKTpckIl2sUW00XygXaiQ3JHVyehsp+NrI5XaR2pTCG/Eyd3CO/GTe6ZJEpAtV81Wm905jWzZDm4bwRXxOl9TXtJyhzQyXQXJ9ktypHNmTWWqlGvHxuGZniQjwynIFb9hLYl1CW5B1gIKvQyLLI7gDbtJH08yUZ0iuT2J6TafLEhGH2HbzjM/CVIHQcIjoqqhmbnaIXlp0UCARILUphVVvzvisFqpOlyQiDrDqFnMH5ihOF4mtjmm5Qocp+DrME/SQuj6F6TOZ3TdLcbbodEnSQTt37mTLli1s377d6VLEIbVSjZmXZ6iVmpNYNHOz8wzbtm2nixhEC2f6FWeKhEZCRFeqm2OQZLNZYrEYmUyGaDTqdDnSIeV0mfTRNKbXJLE+gdun0SYn6KfuEMMwiI/H8QQ8ZE9kqRVqxNfG9Ycg0odsyyZ7sjmeF0gEiI3HtNG0g/Qs67DQSAhPyEP6SJqZvTPExmMEEgGnyxKRFqlX6swfnqdeqhNbHVPXZhfQS44u4A15SW1O4Yv6mD88T2Yig22pB1qk15XmSkzvmcZu2KSuTyn0uoRafF3CZbpIrEvgnfaSPZ6lmq+SWJfA7dd/kUivsS2bzPHmGH4gGSC2Wl2b3UTPql0mNBzCG/Yyf3ie6b3TxMfjBJLq+hTpFbVSjfnD8zSqDeJr4gSHgk6XJOfRS5Au5Al4GN48TCARYP7IPOmjaZ3vJ9IDCtMFZl6ewTAMUtenFHpdSi2+LmW4DOJr4ngjza7PSq5CfDyOL6o9/ES6TaPaIH00TSVXIZgKElsV07aEXUzB1+WCQ0F8ER/pY2lmD8wSTAWJroxqvECkSxSmC2RPZHG5XQxtGNKL0x6g4OsBptdkaMNQczPb4xkqWbX+RJzWqDZIH0tTyVb0grTHKPh6SDAVxBdV60/EacWZItkTWQzTUCuvByn4esxC62+he0WtP5HOUSuvPyj4elRoOIQ/5id9tNn6CyQDRFdGMT066kik1WzbpjBVIHcqp7G8PqDg62Gm12Ro4xDF2Wa3y9SLU0SWRwiNhLThtUiLVHIVMhMZ6uU6oZEQkeURtfJ6nIKvDwSHgvjjfnInc2RPZCnNloiuiuKL6BWpyNVqVBvNv6f5Et6wl+Etw3gCHqfLkhZQ8PUJl+kitjpGMBUkM5Fhdr+6P0Wuhm3bFM4UyJ3O4TJd2n2lD+k8vj610P1pW7a6P7uQzuPrTpVshcxxdWv2O7X4+tRi9+epZvdncaZIdEUUf9zvdGkiXaderpM9maWcLqtbcwAo+PqYy3QRWxUjOBQkezLL3KE5vCEvkRURjf+J0BzHy53KUZwtNk9FX5vQpvADQME3ADxBD0MbhqjkKmRPZJndP4sv6iO6IoonqFe1MnisukV+Mk9hqoBhGs0XiMNBDQcMCI3xDaDSfIncyRz1Sp1AMkBkeQS3T6+BOkljfM6wLZv8mTyFMwVs2yY8GiY0GtI43oDRs90ACiQC+ON+ijNF8qfzTL80TXA4SHhZWDNApS/Ztr34+27VLYLDQSJjEVxuBd4gUvANKMMwCA2HCA4FKUwVyE/mKc4UCQ6dDUCvArAddu7cyc6dO2k0Gk6XMhBs62zgncnTqDbUwyGAujrlLKthUZgqUJgqYDdsAskAodGQZra1ibo620u/z3IpetkjQHMGaGQsQng0vPgKuThbxB/3E14WxhvyOl2iyGU1ag0KZwoUpgsA6sGQC1LwyTkMl0FoJERwOEhprkR+Ms/MyzP4Ij7Cy8LamFe6Ur1Sp3CmQHGmuPg7HBoJacxaLkjBJxdkGAbBoSCBZIByukx+Ms/sgVk8QQ+hkRCBRADDpanf4qxKrkJxukhpvoTL7SKyPEJwOKhZmnJJCj65JMMwCCQCBBIBKtkKhakC6aNpsieyBFNBQsMhdSNJR9mWTWmuRGGqQK1Uw+13N/epHQrqxZgsiYJPlswX9eGL+prdSlMFitNF8pN5/DE/weHm6fBaACztUi/XKUwXKM2WsBoW/rif6Mqout/limlWp1w127IpzhYpzhSpFWuYXpNgKkgwFdTYymVoVufS2LZNeb5MYbpANV/F5Xapp0GumYJPWqJaqC6OtWA3W4eBZHOhvLqfXkvBd2nVfJXibJHyfBmrYeGL+AgONzdeV6+CXCt1dUpLeENevCEv0VVRSnMlSnMl5o/MY7gM/HE/waEg3ohXT1pyUbVSbfF3p1FtNHsQhoMEh4K4/XqqktbRb5O0lMt0ERoOERoOUa/UF5/IZudmMT0m/kQzBLU5tkBz3d3C70itWMNluggkAwSSAbxhrR2V9lDwSdu4fW4iYxEiYxFqxVdezRemCrj9bgKJAL6YT4vjB0yj2qCcKVOeL1PJVTCMZq9AZCyCL6YJUtJ+Cj7pCE/QgyfoIbIiQjVXXQzA3OkcpsfEF/Phj/vxRXwaE+xDtWKNcrpMOVOmVqxhGAbeiJf4eBx/wq91d9JRCj7pKMMwFpdF2LZNrXD2CTFdXtx1wxf14Y/58cV8mh3ao2zLppKrUE6XqWQqNGoNXKYLX+yVHYAUduIUBZ84xjAMvGEv3rCX6Moo9XK92QWWLpM+lgbOTpqJePFFfHjDXrUGu5Rt29RLdSq5CtVclUqugm3ZuH1uAsmzXdphTW6S7qDgk67h9rsJ+8OER8NYdYtyptlaKM40F8obhoEn5MEbVhA67Zygy1ep5qpYDQvDZeANeYmMRfDH/ZqNKV1Jv5XSlVxuF8Gh5lR2aO7asdCauFgQeoIeHSzaJrZlUyvVqBVqi/8Prw660Gio+X8Q8qhVJ11PwSc9we134/a7CQ2HgAsHIYDpNfEEPXhD3sUJNQrDK7MYcsVm0NWKNerlOrZtL3ZPLwZd0KNWt/QcBZ/0pNcEYaW++CRdK9bIT+axGhbwShguvLl9bkyfqZYJzXV09XKdeqm++LN7dci5A248IQ/B4eDiz08/N+l1Cj7pC26fe3EixYJ65ZUn81qhRuFMYTEMDcPA9JqLAWr6zn7sc/fdHpBWw6JertOonA25cp16pfnetpo7Fp4fct6QF3fArZCTvqTgk761GIaJV8JwsYXzqiAoZ8o0Kg0Wtq01XM1QND0mLo8L02Niel/5eOG90118tm1j1SwatcZF3zeqDay6tXgf09MMeE/QQyAZeCXs1QKWAaLgk4FiepqB5ouce5SNbds0qq8KxerZ8Kg2qBVqNGqNxdbRApfpwuV2YbgMDNPAcBm4zIt8vhCSZ9+VciUAinNF3LWzf4Z2sw7bsrEbzfdWwzr347P/ZjWscwINmq02l9u1GMzekBdX3HVOuGntnIiCTwQ429V3toVI7MK3sRrWa1tW9XODyWpY2NXXBtVCa/Lzf/d5Pv/3n6dhNQDIHMtghc8NMJfpumCQml7z3M8957ZCXW6XWm0iS6BjiUQ6aOHPLZvJEk/ESafTi8cSKbREOkMtPpEOWgi3ha5PwzAUeCIdpg5/EREZKAo+EREZKAo+EREZKAo+EREZKAo+EREZKAo+EREZKAo+EREZKFrALuIA27bJ5XJEIhGt4xPpMAWfiIgMFHV1iojIQFHwiYjIQFHwiYjIQFHwiYjIQFHwiYjIQFHwiYjIQFHwiYjIQPn/AXezbS6l0N4FAAAAAElFTkSuQmCC\n",
Luke Naylor's avatar
Luke Naylor committed
      "text/plain": [
       "Graphics object consisting of 5 graphics primitives"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "plot_central_charge(O, r\"$\\mathcal{O}_X[2]$\", argument=2*pi, radius = 0.9) + \\\n",
    "disk(\n",
    "    (0,0),\n",
    "    0.9,\n",
    "    (0, float(pi/2)),\n",
    "    alpha=.4,\n",
    "    fill=False,\n",
    "    thickness=1,\n",
    "    rgbcolor=\"purple\"\n",
    ")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
Luke Naylor's avatar
Luke Naylor committed
   "id": "4af73988",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnQAAAGFCAYAAAB0T27uAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAA9hAAAPYQGoP6dpAABExElEQVR4nO3dd3hUZdrH8V8IJEASmnSkKU2KSlcQFUQQ1vKqoAgWsCAqsCogrrIWBGWXBUWaCGui7tqpSrAsSABpRmkhEAFBCGDokIQUkpn3jyc9M8kkmeTkJN/PdeUa5pwzc+6ZDHnuuZ9yfJxOp1MAAACwrQpWBwAAAICiIaEDAACwORI6AAAAmyOhAwAAsDkSOgAAAJuraHUAAACUlCNHjui9995TcnKyzp07p759++r++++3OiygyEjoAADlwvz58/XOO+/oq6++UocOHZSQkKBbbrlFO3fu1NSpUyVJmzdvVqVKldS5c2eLowUKhi5XAECZN2fOHI0dO1affPKJOnToIEmqUqWKXnjhBU2bNk1RUVGSpODgYDVu3NjKUIFCIaEDAJRpe/fu1YQJE/TYY4/lqrx169ZNDodDS5YsUWpqqk6fPq26detaFClQeCR0AIAybcaMGUpMTNSoUaNy7atTp44k6dChQ/rggw/00EMPlXR4gFf4cOkvAEBZVrNmTVWuXFnHjx93ud/Hx0f33nuv4uPjFRoaKh8fnxKOECg6JkUAAMqs06dP69y5cxo4cGCex4WGhuqnn34imYNt0eUKACizKlWqJEl5jovz8fHRHXfcoY4dO5ZUWIDXkdABAMqsatWq6dprr9WpU6dc7n///ffl5+cnh8MhSbp06VJJhgd4DQkdAKBMe/vttxUWFqbt27dnbAsPD9djjz2moKAgPfHEE/rtt9/kdDo1dOhQJSQkWBcsUEhMigAAlHnh4eGaN2+efH195evrq2uvvVbDhg1TUFCQYmNjNW7cOCUkJOjJJ5/UDTfcYHW4QIGR0AEAANgcXa4AAAA2R0IHAABgcyR0AAAANkdCBwAAYHMkdAAAADZHQgcAAGBzJHQAAAA2R0IHAABgcyR0AAAANkdCBwAAYHMkdAAAADZHQgcAAGBzJHQAAAA2R0IHACh3nE6nLly4IKfTaXUogFeQ0AEAyp3Y2FhVr15dsbGxVocCeAUJHQAAgM2R0AEAANgcCR0AAIDNkdABAADYHAkdAACAzVW0OgAAAIpbUlKSkpKSMu5fuHDBwmgA7yvRhO7YsZI8G4DyIjZWOv77MXXuLAUFWh0NSqM5M2ZoxsyZ2bY1qCH5JB6XKsZZExTspWpDqyPIk4+zBFdVJKEDUFwa1uAPDNxLSk5ScnJyxv3Y2Fh16dJVUVF7FRQYZGFksI1SntDR5QoAKPP8/fzl7+dvdRhAsWFSBADb27VLatjI3AJAeURCBwAAYHMkdAAAADZHQgcAAGBzJHQAAAA2xyxXALbXsqX0009SgwZWRwIA1iChA2B7lStLDetbHQUAWIcuVwC2d/iwNHqMdPiI1ZEAgDVI6ADY3vnz0pIl0vlzVkcCANYgoQMAALA5EjoAAACbY1IEAKDcCA4JVkhIiBypDqtDAbzKx+l0OkvqZMeOldSZAJQnMTHSmm+O6cGHpHp1rY4GdhAbF6vWrdsoKmqvggKDrA4HdlC1odUR5IkKHQDbq1dPGjfO6igAwDqMoQNge7Gx0towKTbO6kgAwBokdABs79AhaehQ6dBBqyMBAGuQ0AEAANgcCR0AAIDNkdABAADYHAkdANvz85OaNZP8/K2OBACswbIlAGyvdWtp409WRwEA1qFCBwAAYHMkdABsLzJSat9BitxjdSQAYA0SOgC2l5oqnTkjpaZYHQkAWIOEDgAAwOaYFAEAKDeCQ4IVEhIiR6rD6lAAr6JCBwAoN0YMH6GwtWEKXRVqdSiAV5HQAbC9K66QVqyQrrjS6kgAwBp0uQKwvYAAqWVnq6MAAOtQoQNge8eOSa+9Lh07bnUkAGANEjoAtnf6tPT++9LpU1ZHAgDWIKEDAACwORI6AAAAmyOhAwAAsDkSOgC2V6uWNHy4uQWA8ohlSwDYXqNG0ptTrY4CAKxDhQ6A7SUkSLt2mVsAKI9I6ADY3v79Uv/bzC0AlEckdAAAADbHGDoAQLkRHBKskJAQOVIdVocCeBUVOgBAuTFi+AiFrQ1T6KpQq0MBvIqEDoDt+fhIgYGSD3/RAJRTdLkCsL327aXfoqyOAgCsw/dZAAAAmyOhA2B7v/0m3dxb+m2f1ZEAgDVI6ADYXlKSSeqSEq2OBACsQUIHAABgcyR0AAAANkdCBwAAYHMkdABsr0kTKThYatLU6kgAwBqsQwfA9qpXl/r3szoKALAOFToAtnfihDR7jnTipNWRAIA1SOgA2F5MjPTWW1LMn1ZHAgDWoMsVAFBuBIcEKyQkRI5Uh9WhAF5FhQ4AUG6MGD5CYWvDFLoq1OpQAK8ioQMAALA5EjoAtletmnT77VK16lZHAgDWYAwdANtr2lR6f4HVUQCAdajQAbC95GTp2HEp+ZLVkQCANUjoANheVJTUpYsUtdfqSADAGiR0AAAANkdCBwAAYHMkdAAAADZHQgcAAGBzLFsCwPbatZMOHpQqVbI6EgCwBgkdANurUEHy97M6CgCwDl2uAGzvwAHp3kHSgd+tjgSlVVJykmLjYjN/YmOtDgnwKip0AGzv4kVp0ybpYrzVkaC0mjN7jmbMnGl1GECxIaEDAJR5o8eM1sgnR2bcj42NVZcuXS2MCPAuEjoAQJnn7+cvfz9/q8MAig1j6AAAAGyOhA6A7TVqJP3rX+YWAMojulwB2F6tWtLQB6yOAgCsQ4UOgO2dOSN98qm5BYDyiIQOgO0dPSqNH29uAaA8IqEDAACwORI6AAAAmyOhAwAAsDkSOgC2V7WqdP31UtUAqyMBAGuwbAkA27vySmnxV1ZHAQDWoUIHwPYcDikpWXI4rY4EAKxBQgfA9nbvlpo3l3ZHWB0JSrvgkGDddPNNGjhgoNWhAF5FQgcAKDdGDB+hsLVhCl0VanUogFeR0AEAANgcCR0AAIDNkdABAADYHMuWALC91q2l8HCpdm2rIwEAa5DQAbA9Pz+pYQ2rowAA69DlCsD2/vhDGvmk9MdhqyMBAGuQ0AGwvQsXpG++kS6ctzoSALAGCR0AAIDNkdABAADYHAkdAACAzZHQAbC9evWkv/1Nqlff6kgAwBosWwLA9urWlcaMtjoKALAOCR0A2zt/Xtq1VbruOql6NaujQWkWHBKskJAQOVIdVocCeJWP0+l0ltTJjh0rqTMBKE927ZJGDDmm776VOnSwOhrYQWxcrFq3bqOoqL0KCgyyOhzYQdWGVkeQJ8bQAQAA2BwJHQAAgM2R0AEAANgcCR0A2/P3l1q1kvwrWx0JAFiDWa4AbK9VK2ntj1ZHAQDWoUIHAABgcyR0AGwvIkJq1VqK2G11JABgDRI6ALbndEpxcZKTtWIBlFMkdAAAADZHQgcAAGBzJHQAAAA2R0IHwPZatJC++9bcAkB5xDp0AGyvShXpyg5WRwE7CA4JVkhIiBypzKBB2eLjdDqdJXWyY8dK6kwAypOjR6WlnxzTM09LjRpZHQ3sIDYuVq1bt1FU1F4FBQZZHQ7soGpDqyPIE12uAGzvzBkpJMTcAkB5REIHAABgcyR0AAAANkdCBwAAYHPMcgVge5ddJo0cKV1W2+pIIElKTJDWrZMcDqlVK+mKK62OCCjzqNAB5VhCgvTdd9KqVdKBA1ZHU3gNG0qvvSo1bGB1JF7Uv780bKj0W1T+xy54zxw7dUrxx5Wfi/HSkiXSrf2k2wZIr79udURAuUBCB5RT8Wntbr9+0gCbt7vx8VL4L1L8RasjycfHH0mD7jVrq/S9RRoz2vw887Q0eJDUrKn017Hm2LZXSf/9RGrV2tx3OqUN66V778n9vE+OMseeOuX9mE+dlF75uzR8uLm/das08gnp2mukNq2ld97OfvzKldIDQyUfH+nIYcmHZgYoCXS5AsXoo4+k5culzZulNm2kq64y2x0O6cQJ6eefpbvukt591zvnO3nSPNeRI2YZj61bpYULzW1iojRqlPTcc+bYlSuloWnt7uHDUgUbt7u//y6NGGKuFtGhNC8w/NDD5nIWgwZLM2ZI11ybfX/oSumXX3I/7usV0po1JnONji6RUDN8+620Nkzq2lV6c6rkX1l6Z5apxD32mDT9X9Kjj0rVqkvJSVJAgLRzhxQZKYWG2vubAmAjJHRAMXo4rf0enNZ+X3tt9v0r3bTfhfXtt1JYWts7dapUubI0a5bJAx57TPpXWttbubJpd3fQ7pa8LVukoEDXmed117tegf2OO83PF5+bX1phJSVKX3whbdps7qevK39LH5NkuvLgQ9I//ymt/p90++3SSy+b7VWrSnfeaVZ1Dqpmtm3bZj58R49K8pH+/DPzHACKFQkdUMy2bJEC3bTf17tpvyVTUfviC1PdkzLbxT59TILoykNpbe//0trel120vdWqmZjS210f2t2StWWL1L27VMHX3L940XSl9usvpaZILVsVz3mjo6UXJphvGXNmZ54/P7EXpNNnpBo1MpM5STp9Wvr+e+m998yHSJJiYqTu10l16krXdpT275d+/VW6kovsAsWNhK4ErFhhxjf7+1sdCayQ3n77Zmm/1683n4mUFDMJMKfoaGlCWts7e3bmY/Nz4YK5WkKNGpnJnJS77Y2Jka67TqpbV+qYpd3l4vbFzJFqSrLPP5+5LThY6tTR/LtOXemmusVz3uefMyXaJk0L9ti9aZMyHnrQfDOIiJC++lJKSpLeeUdqkOVySGfPZn9sRIQ05P4ihQ7AMyR0xeyrr6TYWJK58irVTfvdMa39rlvX/OR8zHNpbW/TAra9UWlt74NZ2t4vs7S9DdPaXlft7v02bnd9faVatSTf0v4XLSLCzNzYvFnas0fat8/c7t1bvOf99lvzDaKgyZxkYvaRmRQxYby0erUUc0Jq1FDq0cN0BUuZS5UMHSpV8pPWrJa6dyu+iiOAbEr7nz/bSkkxY6a2bZM++8zqaGCViAhTkStI+53e9hY0mUs/n2Ta3vFpbe+JEyaR69HDdLsmZGl3/fzMMd26ua4U2kXbtlLELquj8MDWrVLVKtIHH5guz0vJ0hMjzUSD4hQeLj31VOEeGxEhXXONVL+BmQCRvu2vY6Wnn5aaXyG1b2/O8corZkBmxYpmkObz47z3GrwkOCRYISEhcqQ6rA4F8CoSumKQlCQ98IDpapPMCgXe9p//SL17e/954V1bt0pV0tpvX18pOdksgFs5j/bbG21vgwamwpe+bWxa23vFFaY6l7PdHVf62t2yacsWU55NH79WyU+6rnvm/k8/MUt+eFtKivTKq+735zUpYusW6e67s29r395Mp+7XX/ppg7l/8qR0Qy+paTOvhV0cRgwfoRHDRyg2LlatW7exOhzAa0joioG/v/Txx9Jrr0l//CHNmeP9c9So4f3nhPelt9/pY+D8/Mx4unSffGIqZVmlpEiv5tH25jUpYksebW///tKGDaaLt1cvqVmzAr+cUisqShr0vOnObl2aK41bt+b+hY9Ky97PnpHOny+e87ZrZwZNFjTZio6Wfj8o3Xpr7n1Xpl39oUHaas5l6QMF2BAJXTEJCJCmT5cWLZLWrpXuu8/qiGAFV+13evXtjJv2O73tLWj7GB0tHfSg7W3cuGDPawfJydKhQ2YZtFLrwH7p1GmpSxfX++fMlR57tHjOfedd0pQ3pClTCrbQb9ha86HpcHXufZs2mex5wABzv2Mnr4QKoHBsvJSoPTz+uKm4JJXmhgbFYv9+M7vUXfs9d64Z05bTXXdJ779vFh8uiLVrTdt7tZu2t1Va29uJdtca6WMwcv4CHKnShyHmEl8N8xmf4XBIjkKsL1OlijRokDRxYsEWJl6xQpLTrHGzZLEZgHn+nLR0iRlH8N//mm5jAJajQlcCclZoUD64a79TU02XfFSU6/GVWdvev/5Vuvxyz863YoVZS65rV+nFF6WBA80XiTVrzCW+/vtf0+WLErZ9mzR/vpkZ41tBmjTJXJbD6TQzXiMizGKE7y9w/xw/rjH981u2mDXh7rrLrDEzY4bncXTsJNWrZ9auOXTIDKBM16+fNOSB7Mfv329m8Sx431wh4vnnpb//3QzSHDTIXAaFy3oBpQYJHeBl27K03xVytN8Xs7TfC/Jovzvl0/Y+4Kbtff99c4UIV22vnS/tZWvXdjRJUVH07mN+iqphI2nyG54d6+8n/fCDdFltaemyop8bQLEioQO8rGNHk1gVVaNG0hsetr1+aW1v7drSsmVFP7fdNGtmCljNmlsdSRnSuInVEQAoABI6oAxoUs7b3qAgqfVNVkcBANahEwaA7cXEmOFkMSesjgQArEFCB8D2TpyQZsyUTsRYHYkXRe6Rhg01s1/zs+A9c2zt2sUfF4BSiS5XACiNvvvO82OfHGV+AJRbVOgAAABsjgodAKDMS0pOUnJycsb92NhYC6MBvI+Ezgt27TLri2W9RmdZlpAgrVtnFq1v1SrzslKAVapXl+65R6pew+pIUFrNmT1HM2bOtDoMoNj4OJ3OQlxHpnCOHSupMxXemTPSkSMmabnuuvyPP3JE+vxzafz4op13yhTpuefMNWBLs/h4s87Z0KGSj4/08MNm0VrAag1r2OAPDCzjqkLXpUtXRUXtVVBgkIWRwTaqNrQ6gjyVqTF0p09LkydLgwdLo0eb1fKDg80K/U8+6dlz7NxpVtgfO9az4994Qxo50vW+xYulO+4wC8Q2aZL9OSdONNvbtJGGDzc/r77q2Tm96eRJ83qHDzf3t26VnnjCXGGgdWvp7bezH79yZWYyd/hw+bz6wK5dmZf1Ko1Ke3zFITFROnhISuSayXDD389fQYFBmT9BJHEoW8pMc/zrr1KfPubi5J9/Ls2ZI82cKV11lTRsmPTNN549z803SxMmeHbszz9Ll10mVavmev+990pffy21bSu1bCm9+27mvscek2691TS8ISHmep1XXCGtWuXZub3l22+lsDDzOqZOldavl2bNkv73P5PQ/etf0vnz5tikJFNB3LHDrMr/0ksmGSxPjhwxkw9La/e6N+KbMsVUYj3dXhrs2yf17Cnt+83qSADAGmViDN3Jk6bCNGKEqS5ldd115jqXYWGeP5+vr2fHffhhZmUrL0OHmut57thhKl+nTkmLFkkLF0qVKmUeN2yYeb4BAzyPtageekj65z9NAnf77dLLL5vtVatKd94pHT2ambBu22Yu/H70qKnQ/fmnqX5a7eJFKTLS3F5xhecXsy+MN94wSa47KSnSZ59J//2vtHev5O9vqrMtWpjLczVuLI0bZ018H30kLV9urjHbpo35siOZsZAnTpgvKHfdZarHr76a+3lGjHC9PacLF8yVG3x8CvcanE5p9mwz/CEw0FSCp0xx/8XJFnbtkhIuSt1K6TeB0h4fgHyViQrdu++ayoG7btX77pNq1vT+ebdula6+Ov/j7rlHqlzZNPIXL0pvvWUqW1mTOckM7K5WzSQnnhg3zjSeRXHhgmk4U1IykznJdF9//725OHx6wxwTI9Wta65V+sAD0k03mcqolZYvl7p1M4nIAw+YytSAAWac38mT0sGD0rRp3jlXfhXZixelIUNMMrNggfTVVyYBuvlmKSrKJCaeduUXR3wPP5yZTM6YYarYc+ZI8+aZWOfNM+vSNmrkulrsbntOmzaZLydLl0qpqQV/HSEh5jlee82MTW3ZUvrrXwv+PKVGdFrZtLQmS96Ib+oU6aKL8q277QC8zvYVOqfTNOqtW5uqkit+flLnzpn3HQ7TeJ08KdWoYSpOo0ebC3xntXWraVji4kyV7913TWMtSb//LtWqZZ47P9Wrm+rXsmXm3BMmZD5PTl26SD/+aLpp8xIfbxKxolYtotIWoX/wQfP+RURIX35pulffeUdqmGUM6Nmz2R8bESHdf3/Rzl8Ue/aY39stt0i9epnfyZIlZhzkM89kHvfBB945X34V2alTzazfu+4y95s0MZ+7bduk1atNVTZnEu9NnlSMt2wxVa8OHXLvu/76zIlL7qrFnlSR+/eXeveWvvjCVIAHDjRfqjz5vyKZ/5t/+1vm/UGDTFJ+7Fj2z2ORXLwo7YmU4ktBWTfVRVm3cVpZ1z+trPu8RWXdj9PKups2S1d5UNad7qasm3O7tzmd0py0sm5AoHQkrawbZOeyLlAwtk/ozp411aSuXfM+bs6czH+/9JJpbF9/3dw/etRU0UJDTYVDMgnTnj2ZlYFBg0yXVXqi8OefBbvKzrBhpgrSpo1Uv77749q2NZMp8vPWW9Ljj3t+fnciIszt8OGmGrJ6tfk73bCh1KOH6XaVMpcqGTrUNMyrV5vKWKtWRY+hsH78UerUyVR00r30kmkbv/1WunTJJBL9+3vnfFu3mnGZ7jz3XO7PRHS09MIL0r//bb4AFKf84pNMQte9e+awgosXzbjJ/v1NlTb995m1Wpz1y4W77Tn5+ZkvCQ88IK1YYdr1G280CZ67L16SdOCASdyyfq4aNDBfgDZsML/PIlux3JSjz2T5hnLN1abE37On+Qb3+efSxBeLfq7wtLKpu8Ti4kXpkYfNG7NggfmW+dtvprS8erV5w6dPL3ochY3voYdNYjlosCnrXnNt9v2hK6VffpEaNpKaXyF9u0q6LUu27257YRzYL507J3XuknvfhyHm2/cnn5r7c2ZLY/9qZsUB5USZ6HKV8h+vk14R27tX+s9/TGKSrlEj0zjNnZu57dw5k8Sla9jQVIDSnT5dsOrYzp0mkVu2LO/jatSQDh1yvz8x0XzxbNjQtD1FFRFhxvU1aGC+pG/bZnpfgoKkp5/OTPjCw6VXXjFJ7zffmNmtxTkWzBPNm5v3I2u3XkCAmXAyfrxJnpOSMl9DSoqp6D34oHnMihWmO/T++00b/ttv0lNPmSrpsmUmiU3nSUU2ZzJ35oxJuqdPz1399TZP4ktNNW3v9ddnbgsOzvy/Ubeu6UZPl14tzsnddld8faW77zb/55o3N+/vzJnm/5crf/xhbnNWsAMDs///y6lDB+nYUdeVx2z2ppV1O3eR3pgsjRguVa8m7dgpPf2MSVh63mDGFXjDhx+a2VHuvJlW1r3zLlOV69RZGvJAZgl00SKpYjGXdfOKTzLfAoLclHWvu9788ZCkB4dJ77+f+xh32z115LBZsuBfM8yH1JV588zyBukGDTJjRo6zlA3KD9tX6GrWNF8w81vj7o8/pKZNTTWiYkWTOGVVr57000/ZnzfrmnAVKmRPHBwOz2Ncvtw09k89ZXof9uzJ7LnIqXp1yd0C5k6nSUQjIkw341NPeR7D9debMVQ5bdliGtys2rc33cv9+5uqSPv2pnDQq1fxJyYFMWCAqRQNHWralPQxdJGRppIzYYKpDCWlLWUxc6a0e7dJ3MaONVXZnj1NhW/YMNMt3q2bWZrl+eelTz/NnC3qaUU2Ls4sY3P0qEn6JTNruLgLBZ7EFxFhCkKbN5vP4L595nbvXtfHu6sWe1pFzsrHR+rXz/xs2GCW8/nb30x3bFbpM6pzVvGqVs3cVyRrfpQ6elDW7VeCZd3L8ijr1ixFZd0KWcq6G9ab9yg1RWqZVk6tlla+3RMpXZWlfOtue36OHTV/iM6ckZ59VmrbzvVxvx+QjuYo69ZvIFVLK+sO9kZZFyj9bJ/Q+fhI//d/mRMO3HXnfPGFaeBTU02lJjU1+2zWxESz3VO1annWwGzcaLown3jCVCXefNNUK6ZOdX28w+F+nJWPj+kJmjHD/H18+mnP43UlOtr07Nx6a+596Vd/SP/yXZoSuaxcVQl37DC/y927zX1/f3M7YYL5nc+caRK4Hj3M0KmQEFM9Sp8hPX++uc269IenFdnAQDOR5IknzO+6Sxf3CVNhuLtKhyfxbd0qValixhT6+krJySb5rFzZ9fHuqsX5VZHdSa+KfvGF+XLRu3fuY9LXNcw50zwlJe//n/v3S09ONusmtsjryiVXpJV1HamZCUrVAOnRx8w3l6VLzWJ2ERHmfmqKSSbOnTPVrJUrzQeobj3pg3+bb5Iz3zZjySZNkm7rL1WuYp73YFrZtFIeZdOcydzZM9ITaWXdps3yeCFe4El8jrSy7vPPZ24LDpY6pVUw69SVbspSNevSxSTNORM3d9tdORFjpjkfPWre+6uvyft4d2XdgHzKukAZUya6XMeONQ3pggWu9585k1lR69bNVLpiYrIfc/hw/uPwsqpXL/ckgZyiokzjm54o1KhhKhJLlmTvzsvq/PnMcXyuVKokvfiiWfpk40bP43Vl7VqTsLmaqbtpk0kY0nt+OnUq2rlKUnpv2YoV0po1mdtddcvnNZ4rK08rsufOmd/3+PGm+ufvb7q0vSE+3nx2+vUzv5f0MaCexrdli3lv0pMlP7/sSesnn2Q/3l21OK8qsiuXLplq54MPmgrmhx+a96hKldzHpn/2c76eixfzTlgTEkzekXAxn2BuGyD17yc9MFR67llpa9oKzHvSyrqXX27Kuq1amu3pZd01P5o/NAEBUo+epvIzbJhZD6l7N+nkCWnc82ZsRTpPy7rxcdKwodLNN0m9bpSOREvvzMr/cUXlaVk3Pq2s+9ex0sAB0r+mu/+D0LZt5jcpT7Zndfq06QYfP950QX8QnH8yJ7kv6wZ4q6wL2EOZSOhq1zYTFj76yCR1WRuDEyfMDLn0CQSdOpnunqxdRgcPmsHYY8aY+05n7vXVcm5r0cJ0Q7qrGmzbZs77wgvZtw8bZpYK+fpr1487cyazKpaXCRPMOnZFsWKFeU1du5r3IyHBJCRLlpgqzn//6/nMRCuEh5uB9u3bZ7/8WJs2mTNNJ00qWPe4O55UZM+ckUaNMuPtc3452LvXJDJZ/e9/ZpjP8uXmfn4V4ryu0uFJfFu3Zp/tLWV22585k/vx7qrFeVWRs0pMNMWchx82wxw+/tiM/c/rM9Wkibk9eTL7+S5cyNxXZM+PM6Xut9/JXKpjxw7pUooUkZZ0+KWVdcdPkP5yu/n3sGFSn1vMJAbJlHUff0IaPsIMkK1dO/vSH6dPezbLMiCtrFuvvvT+AunXX0yy6S2JCdL335mJCb8fKFh8W7dKVdPKurPeNR/Wm26W/PMo6/5xyPPt6b5eId03WOrbV/ro44J9u07/j1ChgGVdoIwpEwmdZKogq1ebhuD++81s1IkTTSPy2mvZq16zZ5vG5u9/N2Pa5s41SUyjRuZ20iST2AwdahK9Z581w2u++y6z56FCBZMc5lwzbtUq0wV8xx3mS23WBY3DwzPXeps0SXr0UdOQZrVjh2fXkK1Sxbymwq5Dt3+/GUP13nvmdT//vGnsn3rKJHkffeTFJSKKwblzZoLGyy+b92DSJHOVhHTpV/v44w/p+PHM7RfTKjjpFdKc9yXzu82ZBOZXkT1+3EySfPXV3BXPpCQzXjtrj9DLL5tKWb9+pmfp22/NZBN3yWd+V+nIL779+0373cXFBEHJ/B9In9Gczl21OL8qcny8GaP+2GNSnTrmi8HgwSapy0+TJqZ7/0CWvOPAAfP6vTEJyK2sZd0fS7ise/6cSQwnjDfVPz9/zypTnriYVta9tZ+pTha2rJueLFXyk67LkrR+6qKse8FNWdfV9nQ9epoZOYsWSdu35R9XVukfRqeLsi7LlqAcsf0Yuqxq1TINe34qVXJ/ea977jE/Wb3zjvnJ6f77TdUkawM+YID7Nbq6dMneBejKL7+YxWk9kd+K/Xnx85N++MEUFfKbeVsaBQSYrmc/P9NT9umnZmxW+pi65s3NMX5+JqmQTMITHm7+/d13ZmLJDz+Y+1FRmeuc1atnKqwvvmg+J5ddlr0imzMx2bfPJP8dOpgkc/NmU2WtV8+0Y//7n5lQkv65Cg83Fbv0cZStW5svAjfc4P7auPldpSOv+KTMa7vm7ClLTTVfeqKiTGKflbtqcX5V5HXrTHd9Ycd4Dh5slvi5Pa0w9vnnJvFt0aJwz5fNL+Hm28uZM6Z8/lBata11Wll3+XLzR2TDBsmniN93a9WSLuRTNj17RnrqaWnSy1KHHN8EovaaX9zDj2RuW/0/k4glJJhuydQUyTePP+MrV5ruZR8fM1vUJ0dZN7/4tm7NviSAJI16KjN2V2VdPzdlXVfb0112mfTKq1LMn+Yb9/z5pou7Xfu845Oyl3XTxyQ608q6TXOXdYNDghUSEiJHqhdK90ApUqYSupLWp4/piUhIcD0eqKCOHTNf/r3WtZSHkjhHccra5TdhgllG44MPzJf8tm1NkuJwmEQ8vYuvZk0zmfGf/zSVNMmsYBEfb8Z0pa8TN3asSeYefDDzy3/WimzWBH7/fjP54dZbzbiyrVtNte7ECXO/Th1TiZ09O/MxW7dmr8JWqmQSttdec/96Y2LMY9Kv1LF/v7lKR3qS4y6+bdtM27h5szlm0iRz63SaAkZEhPncuRp/6q5anF8VuaiXrnvmGfOevvyyGTcXE2MmO+Tl8svNe3x54zwOOn/OZNwvvSyNfMK8Gb17Zz7ohQkmoTv0h3TseGaG62lZ91KO7r38yqZ/HjcftsmTpTY5pr0np5V1s1bUJr1suiRv6SvNn2fKuomJ0l13uk4+k9PKujt3mA9GaGj258svvgP7pVN5lHXnzJUeezT7tnNuyrfutudUr740Zar55vLuLCk2zozda93G/WMaN5GaNzOl3PT38cABM7nFRVl3xPARGjF8hGLjYtU6r+cFbIaErohefNEsWuyu4lcQH3yQOY4PnqtXz7RVH31kug4rVDDj6DZuzL1sVZUqmclcuokTs9+/5RYzaTEnVxXZFi0KviRJu3YmwZJMRW3fPlOlS0zMvHbp5Zdnr9Z5cpUOV/F17Fj4JcDcVYsLUkUujEqVcv+O8lOzpnTvPfkclF7WrZRW1v0kraybfhWGZs2lwACzv25aWffcWVPVk0xZt0+Osu7xY1KDLGXdv71oxt1lLeu6qqLtz1LWXRlqrsTQoIFUP0dZ9+60F/VLWll3So6ybs8b3FcSs5Z1lUdZ112Vz11Z15FW1v0tyiwcnFVhy7o5NWok/eOfZtzdrFlmZs2zz0pXuinTDh4sfflV5njHz9LKuu6OB8ogEroiat/ezDjdsMF0mRVWRIRZRsJOs0lLk3r1vJNU58VbFdmbbjLJW/r4yoceMt2kkZFmksejj5rxcemL/Hp6lQ5vVozdVYs9qSIvWeJ+0k9WgwfnXoeusE6flr5bJt15Rx6FoIouyrr/TivrXtVW+k9aWXfWO5lLedTIUtZ9JS3LfCZLWbdmlrLuxBelYVnKuj5ZyqZZu1MP7Jem5ijr/llMZd3u15mlRa7NUtZNT3Lcxbc9S1nXN0dZNz5LWfd9F2XdnW7Kt+6256dpM7MszP79pkx7yy2ZSW5W6WXdSWll3RMx0jv5lHWBMsbH6cw5n7P45Lf4b3k2a5YZc1Sc1/pE0UVEmMJIcSaPly6ZKlh6+7d+vbmc5/btJvELCHC9hps345syxSRbOb9guNtutV27pBFDjum7bz24WkS6EzGmrLs70iQubdpIjzxiEiBv+XGNqXS9+Lf8j81LWJi0ZbP0wkRTUfvwQyn8F7PMx2W1zfi4yy/PXq37MER6ZHjm/fvvl6a8kbkQsDfjy3qOf0033aCebLdQepdrVNReBQW6ubg2kFXVUjxTUFToSo30a8aidPNWRTYvP/5oetvSFeQqHd6Iz121uMxVkevWM92jxal3CZR1e6WVdSdlKesmZinrVvKT1qw26+W1zFHW9VZ8kul+rlo1d9LmbjsAr6JCB5Qy8fHZLzv3668lm0S5qxaX5ipyoSp0JSUiwqwBV5zJY0paWbd7Wll3g4uy7s15lHW9Ed/UtPJtx06ebbcYFToUGBU6AAWRNZmTSr4i5q5aTBW5kNq3l06fkn7aYCYxFAdXZd0benl2+TBvxLc7rXybM2lztx2A15HQAbC9gADT2xgQaHUkbtx0c/E+f8+emdeQlQp+8eWixrd6dfaVrvPbDsDr6HIFUCY0rMEfGHiOLlcUWCnvci0zl/4CUH6lppo1aFn8H0B5RUIHwPYiI81au5G7rY4EAKxBQgcAAGBzJHQAAAA2xyxXAEC5ERwSrJCQEDkYcIkyhgodAKDcGDF8hMLWhil0VajVoQBeRYUOgO21aSPt2ilVq251JABgDRI6ALZXqZJ0WQ2rowAA69DlCsD2Dh2Shg+XDv1hdSQAYA0SOgC2Fxsrff+DFHvB6kgAwBokdAAAADZHQgcAAGBzJHQAAAA2R0IHwPbq15defVWq38DqSADAGixbAsD26tSRnhxpdRQAYB0qdABs79w56etvpHPnrY4EAKxBQgfA9o4ckZ58Ujpy2OpIAMAaJHQAAAA2xxg6AEC5ERwSrJCQEDlSHVaHAngVFToAQLkxYvgIha0NU+iqUKtDAbyKhA6A7VWuLLVvL1WuYnUkAGANulwB2F7LltL331kdBQBYhwodAACAzZHQAbC9iAipWXNpV4TVkQCANUjoANie0yklJ0tyWh0JAFiDhA4AAMDmSOgAAABsjoQOAADA5li2BIDttWgh/fij1KSJ1ZEAgDVI6ADYXpUqUsMGVkcBANahyxWA7UVHS+PGm1sAKI9I6ADY3tmz0qefmlsAKI/ocgUAlBvBIcEKCQmRI9VhdSiAV1GhAwCUGyOGj1DY2jCFrgq1OhTAq0joAAAAbI6EDoDt1a4tjR4t1a5jdSQAYA3G0AGwvQYNpJf+ZnUUAGAdKnQAbC8uTtq4SYqLtzoSALAGCR0A2zt4UBo0SDr4u9WRAIA1SOgAAABsjoQOAADA5kjoAAAAbI6EDoDtVaxoZrpWrGR1JABgDZYtAWB7V10l/RJudRQAYB0qdAAAADZHQgfA9vbskTp3kfbstToSlFZJyUmKjYvN/ImNtTokwKvocgVgeykp0vHjUsolqyNBaTVn9hzNmDnT6jCAYkNCBwAo80aPGa2RT47MuB8bG6suXbpaGBHgXSR0AIAyz9/PX/5+/laHARQbxtABAADYHAkdANtr3lz66iup+RVWRwIA1qDLFYDtBQZKra63OgoAsA4VOgC2d/y49OZb0vE/rY4EAKxBQgfA9k6dkubMkU6dtDoSALAGCR0AAIDNkdABAADYHAkdAACAzZHQAbC9mjWlBx4wtwBQHrFsCQDbu/xyaca/rI4CAKxDhQ6A7SUkSFG/SQmJVkcCANYgoQNge/v3S717S/v3WR0JSrvgkGDddPNNGjhgoNWhAF5FQgcAKDdGDB+hsLVhCl0VanUogFeR0AEAANgcCR0AAIDNkdABsD0fH8nPT5KP1ZEAgDVYtgSA7bVvLx06aHUUAGAdKnQAAAA2R0IHwPb27ZP69Zf27bc6EgCwBgkdANtLTJQiIqTEBKsjAQBrkNABAADYHAkdAACAzZHQAQAA2BwJHQDba9xYWrBAatzE6kgAwBqsQwfA9mrUkO643eooAMA6VOgA2N7Jk9KC96WTp6yOBKVdcEiwbrr5Jg0cMNDqUACv8nE6nc6SOtmxYyV1JgDlya5d0oghx/Tdt1KHDlZHAzuIjYtV69ZtFBW1V0GBQVaHAzuo2tDqCPJEhQ4AAMDmSOgAAABsjoQOAADA5kjoANheUJDU71YpqJrVkQCANVi2BIDtNWsmhYRYHQUAWIcKHQDbu3RJOn1aupRidSQAYA0SOgC2t3ev1OFqae8eqyMBAGuQ0AEAANgcCR0AAIDNkdABAADYHAkdAACAzbFsCQDba9tWioqSqla1OhIAsAYJHQDb8/WVggKtjgJ2EBwSrJCQEDlSHVaHAngVXa4AbO/336UHhkq/H7Q6EpR2I4aPUNjaMIWuCrU6FMCrSOgA2F58vBQWJsXHWR0JAFiDhA4AAMDmSOgAAABsjoQOAADA5kjoANhew4bS1KnmFgDKI5YtAWB7l10mjRhudRQAYB0qdABs7+xZafES6ew5qyMBAGuQ0AGwvehoacwYKfqI1ZEAgDVI6AAAAGyOhA4AAMDmSOgAAABsjlmuAGyvShWpc2epSlWrIwEA74qNjdXChQsVGRmp5ORktW3bVuPGjVOlSpWyHefjdDqdJRXUsWMldSYA5U3DGvyBQf6CQ4IVEhIiR6pDB37/XVFRexUUGGR1WLCDqiW/0GV4eLhmz56t8ePHq0OHDnI6nRo8eLACAgL04YcfZjuWhA5AmUBCh4KIjYtV69ZtSOjguRJO6Hbt2qW5c+dqzpw5qlgxs0N1/vz5Gj16tM6dO6egoMzPLmPoANjerl1Sw0bmFgDsLikpSW+88YbeeeedbMmcJEVHR8vhcCg5OTnbdhI6AACAUmTevHkaPXq0KleunGvf119/rY4dO+qyyy7Ltp2EDgAAoBT5+eefdeONN+baHhoaqsjISM2fPz/XPma5AgAAlBI7duxQp06dJEmHDx/W1KlT5XA4lJSUpB9//FHLly9X9+7dcz2OCh0AAEApsXbtWvXp00eSFB8fr/S5qxUrVtTJkydzjZ1LxyxXALaXmCj5JB5TgwZSZX+ro4EdMMsVBVZCs1xHjBihRYsWydfXN9e+IUOGaPPmzTp06FCufVToANhe5cpS82YkcwDsz+FwuEzmJOmKK67QH3/8od9//z3XPhI6ALZ3+LA0eox0+IjVkQBA4UVGRqpNmzZu958+fVqSFBcXl2sfCR0A2zt/XlqyRDp/zupIAKDw1q1bpxtuuMHt/u3bt6tixYpq1qxZrn0kdAAAAKXAunXrdOWVV7rcd+rUKf3yyy/q06ePqlWrlms/CR0AAEApsH79ep09e9blvkWLFik1NVUvv/yyy/1FWofO6XQqNjbW4+NjY32KcjoAcCk+XpJiFX9Ris09tARQclKykpKTMu7Hmw9NgdowlG/OlAsFOj4oKEg+Pp7nPb///rvatGmjJUuWqF27dtn2RUdHa9q0aRo3bpzLBYelIi5bcuHCBVWvXr2wDwcAACiTzp8/77Jr1J2QkBA1atRIoaGhuvrqqzV8+HD5+PgoMjJS9913nwYOHKh//OMfbpPEIiV0BanQXbhwQY0bN9aRI0cK9AIlqWvXrvr5558LHJ8dHsf74hrvi2u8L66Vh/elsI8tD++NJ49LSkpSUlJmhe748ePq1q2bIiMj1ahRo1ITp9WP4/PiWmHel4JW6J588knNnDlTVatW1SeffKKVK1fK19dXgYGBevzxx9W5c+c8H1+kLlcfH58C/8KrVatW4Mf4+voW+DF2epzE++IO74trvC+uleX3paiPLcvvTVHel6CgIN4XF/i8uFaY98VTFy9eVEBAgCRp2LBhGjZsWIEeb4tJEc8880yZflxh2eX18b5493GFZZfXx/tSPI8tyfPZ5XGFZZfXx/vi3ccVp2PHjhW4UpxTiV36K328XUH7lMs63hfXeF9c431xjffFPd4b16KjozO60C6//HKrwyk1+Ly4Vtzvy6effqoaNWpowIABhX6OEqvQ+fv769VXX5W/P9fmyYr3xTXeF9d4X1zjfXGP98a19PeD9yU7Pi+uFff7cvLkSfXs2bNIz1FiFToAAEoLKlEoa2wxhg4AAADukdABAADYHAkdAACAzZHQAQAA2BwJHQAAgM2R0AEAANgcCR2AcufXX3/V+vXrrQ4DALyGhA5AsTp58qRmzJihqVOnaubMmVaHo0OHDmn58uXq1atXoZ9j4sSJiouL82JUKClz585V27Zt1bVrV6tDAbyKhYUBFJvdu3dr8eLFGj9+vKpWrapJkyapbt26Gjt2rGUxDR48WIsWLVL16tWzbX/vvff02WefKSwsTB06dNDVV18tSXI4HDp+/Lh++uknDRkyRB999JGOHDmiyZMna+HChVa8BHgBCwujrKlodQAAyqakpCR9+eWXeu211zK21a9fX6GhoZYldBs3blSdOnVyJXOSNGrUKLVp00a9e/fWv//971wVnMWLF2vTpk2SpMaNG6tly5ZaunSp7r777hKJHQDyQpcrgGIxf/58PfPMM9m2RUZGqmnTphZFJM2bN08PPfSQ2/3r169XtWrV1KlTp1z7brrppmwXcR85cqTefvvtYokTAAqKhA5AsUhISFCdOnUy7iclJWnZsmV6/PHHLYtpw4YN6ty5s9v969evV69eveTr6ytJio+P14oVKyRJKSkpatu2bcaxNWrUUPXq1bVz587iDRoAPEBCB8Drtm3bpnbt2kmS3nrrLb300ku64447NH78+DwTqsKKjY1V06ZNderUqYxtwcHBGjVqVMb9ffv2qXbt2vLz83P5HKmpqdq0aZNuvvnmjG1z5szJGF9Vv3599evXL9tjevTooVWrVnnxlQBA4TCGDoDXpU8gcDqdiouLk4+Pj5KTk3X48OFiOd+GDRvkcDhUu3btjG1fffWVrr/++oz7R48eVd26dd0+x7Zt2xQXF6ewsDDt3LlTe/bs0c6dO3X+/Hm3j7nmmmv08ccfe+dFAEARkNAB8LqYmJiM5Grq1KmSzNpvPXr00LRp01S5cmVJpip27733KjExMd/nrF69uj7//HOX+9auXasbb7wx477D4dDGjRs1YcKEjG0nT550ORki3YYNGxQQEKBly5bJ19dXycnJGjRoUEasrtSqVUsHDhzIN3YAKG4kdAC8zuFw5Nrm4+OjpKQkHT9+XM2bN5ck+fr6atmyZUU+X1hYmIYPH55xPyIiQvHx8erevXueMWW1fv16devWLWP8nJ+fX7YkcdGiRbnG/9WsWTPPCh4AlBTG0AHwqmPHjrns2jxw4IAqVqyYZ7dnYcTHx+uXX35Rz549M7atX79enTt3VpUqVTK21a5dW+fOnXP7PBs2bMjWRStJ48ePlySdPn3a5WNTU1PdjskDgJJEhQ6AV23cuFH169fPtX3ZsmW69dZbFRAQkLHNG12uP/30k1JSUnTllVdmbFuzZo169uyppKQkzZ8/X88++6waNGig06dPu3zuqKgonThxQj169HC5f9q0aS7Xzjt79my2mbwAYBUSOgBe9csvv+RaeT8mJkbffPONvv/++2zbvdHlGhYWJslMemjZsqVWrFihnTt3qk+fPlq5cmXGzNQ2bdooJiZGKSkpqlgx+5++devWSZKuu+66bNtTU1O1YMEC7d69W40bN8517lOnTmVbmw4ArEJCB8CrUlJSlJCQoN27d6tdu3aKi4vTI488on/+85/q1q2b18+XPn7ugw8+UJ06ddStWzeFhoZq4cKFCgwMzFg7rkKFCurevbt27NiRsXTK1q1bNX36dK1bt06+vr4aM2aMKlSokDE7d9u2bTpy5Ii+/PJLl+cODw/PNs4OAKzCtVwBeE1iYqL+8Y9/6OWXX9Z7772nmJgYJSUlafDgwcVyMfSEhATVqFFDGzdu9Gh9u1WrVmn9+vV68803vXL+vn37auHChRmTPGAfXMsVZQ0VOgBeEx4eru7du6tixYoaPXp0sZ9v48aNCggIUMeOHT06fsCAAZo9e7YuXryoqlWrFunc0dHRCggIIJmzmblz52ru3LlKTU21OhTAq5jlCsBrtm7d6nZiQXHYt2+f7rnnHlWo4PmfsqlTp2ratGlFPvfs2bP10ksvFfl5ULKeeeYZRUZG6ueff7Y6FMCrqNAB8JoGDRqUaPdV1kt7eapjx446ceKE1qxZoz59+hTqvNu3b1dSUlK2de4AwEqMoQOAApo6dapeeOEFVapUyepQUEiMoUNZQ0IHACh3SOhQ1jCGDgAAwOZI6AAAAGyOhA4AAMDmSOgAAABsjoQOAADA5kjoAAAAbI6EDoDX/Prrr1q/fr3VYbhV2uMDgMIioQPgFYcOHdLy5cvVq1cvq0NxyRvxTZw4UXFxcR5vB4CSwsLCQBkWHx+vnTt3Ki4uTq1atVLTpk2L7VyDBw/WokWLVL16dZf7U1JS9MEHH2jhwoXatWuXKleurObNm6tNmzby9/dX8+bN9eqrr1oS33vvvafPPvtMYWFh6tChg66++mpJksPh0PHjx/XTTz9pyJAhmjp1qiZPnqyFCxdme/yRI0dcbvc2p9Opt956S6dOnVJQUJAOHjyo2bNnu33P4R4LC6Os4VquQBn1+eefa/To0Tp16lTGti5dumjcuHHq3bu3YmNjFRwcrKlTpxb5XBs3blSdOnXcJhbx8fG6/fbbNWrUKH3xxReKiYlRZGSk9u3bp5UrV6pmzZrFmgzlF9+oUaPUpk0b9e7dW//+97/VtWvXbPsXL16sTZs2qXHjxmrZsqWWLl2qu+++O2O/u+3eNm/ePIWFhem7776TJL311lt6+OGHtXz58mI7JwCbcAIoc3bu3On09fV13nHHHc5Zs2Y5R48e7axRo4ZTUraf5cuXe+V8w4YNc27cuNHt/tGjRzuffvrpXNu3bNni7NChg/PUqVNeicOd/OJzOp3OyZMnO6tVq+ZMSUnJte/kyZPOt99+2+l0Op1nz5519urVK9cx7rZ7U5MmTZz/+c9/Mu5HR0c7JTmPHDlSrOctixwOh/P8+fNOh8NhdSiAVzCGDiiDVq1ape7du2vFihUaO3asZs+erSNHjmjWrFnq3bu3brjhBi1cuFB33nmnV863YcMGde7c2e3+v//975o7d262bX/88YdGjhyppUuX6rLLLvNKHIWNT5LWr1+vXr16ydfXV5KpKq5YsUKS6S5u27atJKlGjRqqXr26du7cme3x7rZ7y2+//abDhw+rXbt2GdsaNWqk6tWra/Xq1cVyzrLMx8dH1apVk4+Pj9WhAF5BQgeUQS1btlRCQoJSU1MztgUGBmrs2LF6/fXX1b59eyUmJmrbtm2STMLy4IMPauDAgUpNTdUXX3yhdu3aqW/fvoqNjVVkZKSGDBmixo0b69NPP1VCQkLG8+7bt0+1a9eWn5+f23jq1q2b7f7p06d17733auHChbryyiu9/Oqz8yS+1NRUbdq0STfffHPGtjlz5mSMrapfv7769euXsa9Hjx5atWpVrudxt90bDhw4IEm5xnsFBQXpt99+K5ZzArAPEjqgDLr77rt11113qX///ho+fHjGUh07d+7ULbfcoqZNm2r06NEZVafXX39d27dv16pVq/Twww8rMDBQvXv31urVqzVgwAB9//336tWrl/788089+uijCg8PzzjX0aNHcyVsrsTGxuq2225T27Zt1bp1ax06dEhTpkwpnjcgC0/i27Ztm+Li4hQWFqaHH35YXbt21SuvvKLrrrvO5fHXXHONtm/f7vF2bzh79qwkKSAgINv2wMDAjH0Ayi8mRQBllKsZo+Hh4bp06VJG0uHv7y9Jmjx5snx9fbV792498cQTuvnmm9W0aVPNnTtXLVq00LPPPitJmj59uiRlW/rj5MmTHs2yDAoK0ueff657771Xc+fOVY8ePbRr164ivspMCQkJ+uGHH5Samqp27dqpVatWHse3YcMGBQQEaNmyZfL19VVycrIGDRqkypUruzy+Vq1aGRUzT7ZLpgp47733KjExMd/XUr16dX3++efZtqV3Baffprt06ZJSUlLyfU4AZRsJHVCOdOvWTZKZAfvII49owIABkuRyHFHOSpA7DofDo+POnj2rwYMHa/LkyerRo4ckM+vWG+Li4vTpp5/q8ccfl4+Pj26//XZ98803Hse3fv16devWLSNZ8vPz04033pixf9GiRXr88ccz7tesWVPnz5/P9TzutksmEVu2bFlBXlY2derUkZT79cTHx7NsCQC6XIGyaNOmTbrqqqtUp04dvffeexnb27dvryFDhkiSxowZ43EylpfatWvr3LlzeR5z+vRp3XffffrHP/6Rkcyli4iI0Pz587NtW7lypb7++uuMKlV+FajFixdnJHMHDx5UhQqZf9o8iW/Dhg26/vrrs20bP358Ruw5H5+amupyTJ677d7QvHlzSVJMTEzGNofDoXPnzumKK64olnMCsA8qdEAZc/bsWS1evFjTpk3ToEGDNGbMGN12221q1qyZJOmNN97QZ599pgMHDig6OlpNmjSRZCo9knTx4kWX9yWpQoUKunTpUrbzNWjQQKdPn3Ybz9GjR/XQQw9p1qxZ6tChQ7Z9SUlJ6tOnj95+++2MbWPGjNFf/vIX/eUvf9H06dO1bNkyJSQk6P7778+WqGV9jsDAQIWHh2vHjh1asmSJ3nnnHY/ji4qK0okTJ3IlmummTZumsWPHZtt29uzZjIqZJ9ulone5Nm/eXC1atFBUVFTG+xgVFaXExET16dMn3+cEULaR0AFlTGBgoN588035+flp+PDhWrRokT788MOMMXUtWrRQUFCQ/Pz8VL9+fUnSmTNntGnTJknS8uXLNXDgQH399deSpN27dys6OlqXX365GjZsqC1btuipp57S5MmTVadOHbVp00YxMTFKSUlRxYrZ/6Ts2bNHt912mzp16qTFixcrLCws43lq1qypb775Rn379tWwYcMkmcrirl27NHv2bElSu3bttHTpUvXp08dlMidJW7ZsUc+ePXX48GH5+Pjo6NGjcma5AE5e8UnSunXrJCnXBIjU1FQtWLBAu3fvVuPGjbPtO3XqlC6//PJcz+Vuu1T0LldJeuSRR/TRRx9p0KBBkqTg4GDdeeedat26dZGeF4D9kdABZUylSpUy/j158mStWrVK7777rvr166err75aCxYskMPh0IcffpjRPVirVi398MMPmjRpkmbMmCFJevHFFxUXF6fnnnsuY524l156SaNGjdLIkSMzKlEVKlRQ9+7dtWPHjmxrvUVFRenFF1/U7bffrgsXLmjDhg06evSojh8/rgsXLqhevXoaM2aMPv7444zHbNiwIdvYNT8/P23ZskUzZ850+3qPHz+uG2+8UfXr11e3bt20d+9ebd68OSPJcRff1q1bNX36dK1bt06+vr4aM2aMKlSoIKfTqbi4OG3btk1HjhzRl19+meuc4eHh2eLMb7u3TJw4US+++KLGjBmj6tWr6/jx4woJCSm28wGwD67lCpRxx48f1/z587Vjxw75+vqqffv2evrppzOqc96watUqrV+/Xm+++WaRnuf777/XunXrNGXKFKWkpGjevHnatGmTZs2apbp16+rgwYNq2rRptmrdvHnz9PTTT2fc79u3r959992MJVm8GV/WcyxcuDBjXFt+2wGguFGhA8q4Bg0aaPLkycV6jgEDBmj27Nm6ePGiqlatWujn6devnxISEvT9999LMtdYrVSpknbu3Km+ffvq//7v/zR9+vSMRX7Tlyp5/PHH5efnp9DQUPXq1StbMufN+CQpOjpaAQEBuZI2d9sBoCRQoQPgFdu2bdPSpUuLNXm8dOmSNm3alNGtuXr1ajVr1kw///yzKlasqMDAQN12223FGt/EiRN1zz33qHv37h5tB4CSQIUOgFd07NhRJ06c0Jo1a4pt1uW3336rvn37Ztz/888/dcstt3h0+TBvxLd9+3YlJSXlStrcbQeAkkKFDoBtxMXFKTAwMOP+li1bSjSJmjp1ql544YVsE0/y2g4AJYWEDgAAwOa4UgQAAIDNkdABAADYHAkdAACAzZHQAQAA2BwJHQAAgM2R0AEAANgcCR0AAIDNkdABAADYHAkdAACAzZHQAQAA2BwJHQAAgM39P1HK2q4grQybAAAAAElFTkSuQmCC\n",
Luke Naylor's avatar
Luke Naylor committed
      "text/plain": [
       "Graphics object consisting of 7 graphics primitives"
Luke Naylor's avatar
Luke Naylor committed
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "ymin = 0\n",
Luke Naylor's avatar
Luke Naylor committed
    "ymax = 2\n",
    "xmin = -3\n",
    "xmax = 1\n",
    "\n",
    "mu = -1\n",
Luke Naylor's avatar
Luke Naylor committed
    "\n",
    "p = line(\n",
    "    [(mu,0), (mu,ymax)],\n",
Luke Naylor's avatar
Luke Naylor committed
    "    linestyle = \"dashed\",\n",
    ") + text(\n",
    "r\"\"\"$\\beta = \\mu(E)$\n",
    "$\\mathfrak{Im}(\\mathcal{Z}_{\\alpha,\\beta}(E)) = 0$\"\"\",\n",
    "    (mu, -0.25),\n",
    "    fontsize = \"large\",\n",
    "    rgbcolor = \"black\"\n",
    ") + text(\n",
    "r\"\"\"$\\longleftarrow$\n",
    "$\\mathrm{Coh}(X) \\leftarrow \\mathcal{B}^\\beta$\"\"\",\n",
    "    (-2.5,1),\n",
    "    fontsize = \"large\"\n",
    ") + text(\n",
    "r\"\"\"$E \\in \\mathcal{B}^\\beta$\n",
    "$\\mathfrak{Im}(\\mathcal{Z}_{\\alpha,\\beta}(E)) \\geq 0$\"\"\",\n",
    "    (-1.5,1),\n",
    "    fontsize = \"large\"\n",
    ") + text(\n",
    "r\"\"\"$E[1] \\in \\mathcal{B}^\\beta$\n",
    "$\\mathfrak{Im}(\\mathcal{Z}_{\\alpha,\\beta}(E)) < 0$\"\"\",\n",
    "    (-0.5,1),\n",
    "    fontsize = \"large\",\n",
    "    rgbcolor = \"red\"\n",
    ") + polygon(\n",
    "    [\n",
    "        (mu,0),\n",
    "        (xmin,0),\n",
    "        (xmin,ymax),\n",
    "        (mu,ymax)\n",
    "    ],\n",
    "    rgbcolor = \"blue\",\n",
    "    alpha = 0.1,\n",
    "    zorder = 102\n",
    ") + polygon(\n",
    "    [\n",
    "        (xmax,0),\n",
    "        (mu,0),\n",
    "        (mu,ymax),\n",
    "        (xmax,ymax)\n",
    "    ],\n",
    "    rgbcolor = \"orange\",\n",
    "    alpha = 0.1,\n",
    "    zorder = 102\n",
Luke Naylor's avatar
Luke Naylor committed
    ")\n",
    "\n",
    "p.axes_labels([r\"$\\beta$\",r\"$\\alpha$\"])\n",
    "p.tick_label_color(\"white\")\n",
Luke Naylor's avatar
Luke Naylor committed
    "p.show()"
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "d52b9473",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "SageMath 9.7",
   "language": "sage",
   "name": "sagemath"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.10.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}