From 768724ecc72ab521b47a0f4779bce23d175fb63f Mon Sep 17 00:00:00 2001 From: Logan Ward Date: Wed, 4 Oct 2023 13:15:39 -0400 Subject: [PATCH] Add example of using this new data --- .../1_fit-forcefield-using-mbtr.ipynb | 64 +++++++++++-------- 1 file changed, 38 insertions(+), 26 deletions(-) diff --git a/notebooks/2_testing-fitting-strategies/1_fit-forcefield-using-mbtr.ipynb b/notebooks/2_testing-fitting-strategies/1_fit-forcefield-using-mbtr.ipynb index 30793bc..ffa3e2b 100644 --- a/notebooks/2_testing-fitting-strategies/1_fit-forcefield-using-mbtr.ipynb +++ b/notebooks/2_testing-fitting-strategies/1_fit-forcefield-using-mbtr.ipynb @@ -51,7 +51,7 @@ }, "outputs": [], "source": [ - "db_path = '../1_explore-sampling-methods/data/random-dir-variable-dist/water_b3lyp_def2-svpd_d=5.00e-03.db'" + "db_path = '../1_explore-sampling-methods/data/along-axes/caffeine_pm7_None_d=5.00e-03-N=2.db'" ] }, { @@ -98,7 +98,7 @@ "name": "stdout", "output_type": "stream", "text": [ - "Loaded 55 structures\n" + "Loaded 5185 structures\n" ] } ], @@ -171,8 +171,8 @@ "name": "stdout", "output_type": "stream", "text": [ - "CPU times: user 18.3 s, sys: 444 ms, total: 18.8 s\n", - "Wall time: 1.59 s\n" + "CPU times: user 3min 38s, sys: 9.7 s, total: 3min 48s\n", + "Wall time: 20.1 s\n" ] } ], @@ -180,8 +180,7 @@ "%%time\n", "mbtr = MBTRCalculator(\n", " model=GridSearchCV(\n", - " KernelRidge(kernel='rbf', alpha=1e-6), {\n", - " 'alpha': np.logspace(-10, -7, 8),\n", + " KernelRidge(kernel='rbf', alpha=1e-10), {\n", " 'gamma': np.logspace(-5, 5, 32)\n", " }),\n", " descriptor=MBTR(\n", @@ -238,7 +237,7 @@ }, { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAUwAAAFQCAYAAADUeke/AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAAzWklEQVR4nO3dd1hT9+I/8HcIUxmCoFXBTWkRF6ItapWhoiJ2OX+IOG+tQh3X24qzapVWW1u1ah2gthW1irsqooKjbsZ1tU6siKh1EJas5Pz+8JKvCGoICSfj/XqePLc5Cee881x4+zkj5yMRBEEAERG9lonYAYiI9AULk4hIRSxMIiIVsTCJiFTEwiQiUhELk4hIRSxMIiIVsTCJiFRkKnYATVIoFLh79y5sbGwgkUjEjkNEekIQBOTk5KB+/fowMXn5ONKgCvPu3btwcXEROwYR6an09HQ4Ozu/9HWDKkwbGxsAzz60ra2tyGmISJf9+uuvGDduHAAgJCQEv/zyi7JDXkZiSN8lz87Ohp2dHWQyGQuTiF4qOjoao0aNgiAIGDduHObNm4datWq9tjt40oeIjEpUVBRGjhwJQRAQFhaGpUuXqnzOg4VJREbF3NwcEokEn332GZYsWVKpE8QGdQyTiOh1QkJC4Obmhvbt21f6ahqOMInI4G3YsAGZmZnK5x06dFDr0kMWJhEZtOXLl2PIkCHw8/ODTCar0rpYmERksJYtW6a8dKhPnz5VvnqGhUlEBmnp0qUICwsDAPznP//BggULqvwNQJ70ISK9UiJX4PbjfDR0qAFTqUmFr+3cEIVJEycAAL744gtERkZq5OvSLEwi0hslcgU+Wn4C5zNkaNXADtvGdlSWZulrJ/ZtxaN9iwEAERERmDdvnsbuLcFdciLSG7cf5+N8xrMTN+czZLj9OL/ca5aN28K01hv4dMJkjZYlwMIkIj3S0KEGWjWwAwC0crZDQ4ca5V4ztXVEt4hoLF74tcbvWsZdciLSG6ZSE2wb27HcMczvvvsOLi4u2Da230uPb2pk+xpfIxGRFplKTdDUyVr5fMGCBfjiiy8glUrx3/+2QIsWLbS2be6SE5He+vrrr/HFF18AAGbMmKHVsgRYmESkp+bPn4+IiAgAwJw5czBr1iytb5OFSUR656uvvsK0adMAAHPnzsWMGTOqZbs8hklEemX//v3Kgpw3bx6mTp1abdtmYRKRXgkICMDYsWPh4uKCKVOmVOu2WZhEpBfkcjmkUikkEgl+/PFHUWaG5TFMItJpgiBg1qxZ6NevH4qKigBAtGm0dbYwS78sP2HCBLGjEJFISstyzpw52LFjB/bt2ydqHp3cJT979ixWrVqFVq1aiR2FiEQiCAJmzJiBefPmAQAWLVqE999/X9RMOjfCzM3NRXBwMFavXg17e3ux4xCRCARBwLRp05Rl+f3332PixIkip9LBwhw3bhwCAwPRrVu31763sLAQ2dnZZR5EpN8EQcDUqVMRGRkJAFi8eLHOHJrTqV3yTZs2ITk5GWfPnlXp/ZGRkZg9e7aWUxFRdUpLS8OSJUsAAEuWLEF4eLjIif6PzhRmeno6xo8fjwMHDsDS0lKln4mIiMCkSZOUz7Ozs+Hi4qKtiERUDZo2bYq9e/fi8uXL+PTTT8WOU4ZEEARB7BAAsGPHDnz44YeQSqXKZXK5HBKJBCYmJigsLCzzWkWys7NhZ2cHmUxW5cmOiKj6CIKAzMxM1K9fX5Ttq9odOnMM09/fHxcuXEBqaqry4eXlheDgYKSmpr62LIlIPwmCgH//+99o3bo1Ll68KHacV9KZXXIbGxt4eHiUWVazZk3Url273HIi0j8VTV4mCAImTpyIxYufzcFz9uxZnf5715nCJCLDVdHkZVKTZ19MKT3Bs3LlSgwfPlzkpK+m04WZmJgodgQi0oAXJy/7+1Eefpg7FT/++CMAYPXq1Rg1apSYEVWiM8cwichwPT95WcsGtvhu9hTlDTSioqL0oiwBHR9hEpFheH7yMicrCXr1jFCWpa7vhj+PhUlE1eL5ycv27duHxMREBAUFiZyqcrhLTkRap1AoytxpyMbGRu/KEmBhEpGWKRQKjBkzBr1798bXX38tdpwq4S45EWmNQqHAv/71L0RFRcHExATOzs5iR6oSFiYRaYVCocDo0aMRHR0NExMT/PzzzwgODhY7VpWwMIlI4+RyOUaNGoV169bBxMQEv/76KwYPHix2rCpjYRKRRgmCUKYsN2zYgEGDBokdSyN40oeINEoikaBNmzaQSqWIiYkxmLIEdOj2bprA27sR6Y7r16+jefPmYsdQid7d3o2I9FdJSQlmz56NrKws5TJ9KcvKYGESUZWUlJRg6NCh+PLLLxEYGAiFQiF2JK3hSR8iUltJSQlCQkKwadMmmJqa4j//+Q9MTAx3HMbCJCK1lJSUYMiQIdi8eTPMzMywZcsW0ecN1zYWJhFVWklJCYKDg/Hbb7/BzMwMW7duRd++fcWOpXUsTCKqtAkTJuC3336Dubk5YmNj0adPH7EjVQvDPdhARFoTHh6Oxo0bY9u2bUZTlgBHmESkBjc3N/z111+wsLAQO0q14giTiF6rqKgIgwcPRlxcnHKZsZUlwMIkotcoKipC//79sWnTJgwaNAgymUzsSKJhYRLRSxUWFqJfv37YtWsXLC0tsXnzZtjZ2YkdSzQ8hklEFSotyz179sDS0hK7du1C9+7dxY4lKhYmEZVTUFCAjz/+GHv37oWlpSV2796Nbt26iR1LdNwlJ6Jyli9fjr1798LKygp79uxhWf4PR5hEVM5nn32GS5cuITg4GH5+fmLH0RksTCIC8Gw33MzMDFKpFKampoiKihI7ks7hLjkR4enTp+jbty9GjRpl0LdnqyoWJpGRy8/PR9++fREfH48tW7bgypUrYkfSWdwlJzJipWV56NAh1KxZE/v27cPbb78tdiydxcIkMlL5+fkICgrC4cOHYW1tjX379qFz585ix9JpLEwiI5SXl4egoCAkJCTA2toa+/fvR6dOncSOpfN4DJPICCUlJeH48eOwsbFBXFwcy1JFHGESGaEuXbogNjYWjo6O8Pb2FjuO3mBhEhmJ3NxcPHr0CI0aNQIABAUFiZxI/3CXnMgI5OTkoFevXnjvvfdw8+ZNsePoLRYmkYErLcvjx48jOzsbjx8/FjuS3mJhEhmw7Oxs9OzZE3/88Qdq1aqFgwcPwsvLS+xYekvtY5jFxcW4d+8e8vPz4eTkBAcHB03mIqIqkslk6NmzJ06dOqUsy3bt2okdS69VaoSZm5uLlStXwsfHB3Z2dmjcuDHc3d3h5OSERo0aYfTo0Th79qy2shKRimQyGQICAnDq1CnY29vj0KFDLEsNULkwv//+ezRu3BirV6+Gn58ftm3bhtTUVFy5cgUnT57ErFmzUFJSgu7du6Nnz564du2aNnMT0SuUlJSgoKAADg4OOHToEDw9PcWOZBAkgiAIqryxf//+mDlzJlq2bPnK9xUWFiIqKgrm5uYYNWqURkKqKjs7G3Z2dpDJZLC1ta3WbRPpmocPH+Lu3bto1aqV2FF0nqrdoXJhAkBqairatGmjiXxawcIkY5aVlYX4+Hj0799f7Ch6R9XuqNQxTE9PT7Rr1w4rVqww6qk2iXTNkydP0K1bNwwYMABr164VO47BqlRh/vHHH/D09MSUKVNQr149DBkyBAkJCdrKRkQqePz4Mbp164akpCQ4OjrysiEtqlRhent7Y/Xq1bh37x5WrFiBO3fuoFu3bmjWrBnmzZuHO3fuaCsnEVWgtCyTk5Ph5OSEhISE155nIPWpdeG6lZUVQkNDkZiYiKtXr2Lw4MFYuXIlmjRpgt69e6sVZMWKFWjVqhVsbW1ha2sLb29v7Nu3T611ERmDR48ewd/fHykpKahTpw4SEhLg4eEhdiyDVqmTPi+Tm5uLDRs2YOrUqcjKyoJcLq/0Onbv3g2pVIrmzZsDANavX4+FCxciJSUFLVq0UGkdPOlDxiI/Px+dOnVCamoq6tati8OHD8Pd3V3sWHpL1e6o0t2Kjhw5gujoaMTGxkIqlWLAgAEYOXKkWut68c4p8+bNw4oVK3Dq1CmVC5PIWFhZWSEoKAiZmZlISEjgtBLVpNKFmZ6ejnXr1mHdunVIS0tDx44dsXTpUgwYMAA1a9bUSCi5XI4tW7YgLy+P9+ojqoBEIsHs2bMRHh4OJycnseMYjUoVZvfu3ZGQkAAnJycMHToUI0aMgJubm8bCXLhwAd7e3igoKIC1tTW2b9/+yt2MwsJCFBYWKp9nZ2drLAuRrnnw4AFmzpyJRYsWoUaNGpBIJCzLalapwrSyskJsbCz69OkDqVSq8TBubm5ITU1FVlYWYmNjERoaiiNHjry0NCMjIzF79myN5yDSNffv34efnx8uX76MnJwcbNiwQexIRkkjJ320pfSSpZUrV1b4ekUjTBcXF570IYPyfFk2aNAACQkJcHV1FTuWQdHKN32ed+zYMQwZMgTe3t7IyMgAAPzyyy84fvy4uqssRxCEMoX4IgsLC+VlSKUPIkNy7949+Pr6KssyMTGRZSkitQozNjYWAQEBsLKyQkpKirLUcnJyMH/+fLWCTJ06FceOHcOtW7dw4cIFTJs2DYmJiQgODlZrfUT6LjMzE76+vvjzzz/h7OyMxMRE5WV3JA61CvOrr77CTz/9hNWrV8PMzEy5vGPHjkhOTlYryP379xESEgI3Nzf4+/vj9OnT2L9/P7p3767W+oj0mSAI6NevH/766y+4uLiwLHWEWtdhXrlyBV26dCm33NbWFllZWWoFiYqKUuvniAyRRCLB0qVLMXLkSMTGxqJp06ZiRyKoOcKsV68erl+/Xm758ePH+X8sURU8fw7W09MTSUlJ/JvSIWoV5ieffILx48fj9OnTkEgkuHv3LjZs2IDJkydj7Nixms5IZBTu3LkDLy8vnDx5UrnMxITzFOoStXbJP//8c8hkMvj6+qKgoABdunSBhYUFJk+ejLCwME1nJDJ46enp8PX1xY0bNzBmzBikpKSwLHVQla7DzM/Px+XLl6FQKODu7g5ra2tNZqs03nyD9NHt27fh6+uLmzdvonHjxkhMTESjRo3EjmVUtH7zjYKCAly8eBEPHjyAQqHAvXv3lK/17dtX3dUSGZXbt2/Dx8cHaWlpaNKkCRITE9GwYUOxY9FLqFWY+/fvR0hICB49elTuNYlEotbt3YiMzd9//w1fX1+kpaWhadOmSExMhIuLi9ix6BXUOkgSFhaGAQMGIDMzEwqFosyDZUmkmvnz5yMtLQ3NmjVjWeoJtUaYDx48wKRJk1C3bl1N5yEyGosXL4aJiQmmTZsGZ2dnseOQCtQaYfbr1w+JiYkajkJk+B4+fKi81tLS0hIrVqxgWeoRtc6S5+fno3///nByckLLli3LfD0SAD777DONBawMniUnXXbz5k34+vqiX79++PbbbyGRSMSORP+janeoVZhr1qzBmDFjYGVlhdq1a5f5P14ikeDmzZvqpa4iFibpqhs3bsDX1xfp6elwc3PDqVOnUKtWLbFj0f9o9bKi6dOnY86cOZgyZQovriV6jRs3bsDHxwd37tzBW2+9hcOHD7Ms9ZRahVlUVISBAweyLIle4/r16/Dx8UFGRgbeeustJCQk4I033hA7FqlJrcYLDQ3F5s2bNZ2FyKBcu3YNXbt2RUZGBtzd3ZGYmMiy1HNqjTDlcjkWLFiAuLg4tGrVqtxJn0WLFmkkHJE+S05ORmZmJtzd3XH48GFehmcA1CrMCxcuoG3btgCAixcvlnmNZ/6Inhk4cCDMzc3RqVMn1KlTR+w4pAE6PQlaZfEsOYntypUrsLW1Rb169cSOQpWg9UnQiKisP//8Ez4+PvD19UVmZqbYcUgLVC7M27dvV2rFpTNJEhmDy5cvw9fXF/fu3YOlpWW54/pkGFQuzPbt22P06NE4c+bMS98jk8mwevVqeHh4YNu2bRoJSKTrLl26BF9fX9y/fx9t2rTBoUOH4OjoKHYs0gKVT/r8+eefmD9/Pnr27AkzMzN4eXmhfv36sLS0xJMnT3D58mVcunQJXl5eWLhwIXr16qXN3EQ64eLFi/Dz88M///yDtm3bIj4+HrVr1xY7FmlJpU/6FBQUYO/evco5xJ8+fQpHR0e0bdsWAQEB8PDw0FbW1+JJH6pOpSPLf/75B56enoiPj4eDg4PYsUgNWv0uua5iYVJ1unPnDnx8fFCrVi0cOHCAZanHtD5FBZGxc3Z2RmJiImrWrAl7e3ux41A14GVFRJWQmpqKrVu3Kp87OzuzLI0IC5NIRSkpKfDz88OgQYMQFxcndhwSAQuTSAXJycnw9/fHkydP4OXlhXfffVfsSCQCtQvz2LFjGDJkCLy9vZUXqf/yyy84fvy4xsIR6YKkpCR069YNT548wbvvvosDBw7Azs5O7FgkArUKMzY2FgEBAbCyskJKSgoKCwsBADk5OZg/f75GAxKJ6dy5c8qy9Pb2RlxcHK/AMGJqFeZXX32Fn376CatXry7zFbCOHTsiOTlZY+GIxHTr1i10794dWVlZ6NSpE8uS1Lus6MqVK+jSpUu55ba2tsjKyqpqJiKd0KhRIwQHB+O///0v9u7dCxsbG7EjkcjUKsx69erh+vXraNy4cZnlx48fR9OmTTWRi0h0EokES5cuRUFBAaysrMSOQzpArV3yTz75BOPHj8fp06chkUhw9+5dbNiwAZMnT8bYsWM1nZGo2pw6dQpDhw5FUVERgGelybKkUmqNMD///HPIZDL4+vqioKAAXbp0gYWFBSZPnoywsDBNZySqFidPnkRAQABycnLQrFkzzJo1S+xIpGOq9F3y/Px8XL58GQqFAu7u7rC2ttZktkrjd8lJXX/88Qd69uyJ3Nxc+Pj4YM+ePahZs6bYsaiaVMt3yWvUqAEvL6+qrIJIdMePH0evXr2Qm5sLPz8/7N69GzVq1BA7FukgtQpzzpw5r3x95syZaoUhqm7Hjh1Dr169kJeXB39/f+zatYtlSS+l1i556YyRpYqLi5GWlgZTU1M0a9ZMtGsxuUtOlZGfn4+mTZvi/v376NatG3bu3MmyNFJa3SVPSUmpcIPDhg3Dhx9+qM4qiapdjRo1sHnzZvzwww+IiYnh2XB6LY3eQPjixYvo06cPbt26palVVgpHmKSKoqIimJubix2DdIgo0+xmZWVBJpNpcpVEGpWQkIA333wTFy5cEDsK6SG1dsmXLFlS5rkgCMjMzMQvv/yCnj17aiQYkaYdPnwYffr0wdOnT7Fw4UL8/PPPYkciPaNWYX7//fdlnpuYmMDJyQmhoaGIiIjQSDAiTTp48CCCgoJQUFCA3r17Y/Xq1WJHIj2kVmGmpaVpOgeR1sTHx6Nv374oKChAYGAgYmNjYWFhIXYs0kOVPoZZXFwMX19fXL16VRt5iDTqwIEDyrIMCgpiWVKVVLowzczMcPHiRUgkEm3kIdIYQRCwcOFCFBQUoG/fvti6dSvLkqpErbPkQ4cORVRUlKazEGmURCLBtm3bMH36dGzZsoWXElGVqXUMs6ioCGvWrEF8fDy8vLzK3aRg0aJFlV5nZGQktm3bhr/++gtWVlbo2LEjvvnmG7i5uakTkYzYjRs30KxZMwCAjY0N5s6dK3IiMhRqjTAvXrwIT09P2Nra4urVq0hJSVE+UlNT1Qpy5MgRjBs3DqdOnUJ8fDxKSkrQo0cP5OXlqbU+Mk6///473N3dObcUaUWlRpgjRozA4sWLkZCQoPEg+/fvL/N87dq1qFOnDpKSkiqcDoPoRXv27MFHH32E4uJipKSkQKFQwMSEM0mT5lTqt2n9+vV4+vSptrKUUfqNIQcHh5e+p7CwENnZ2WUeZJx2796tLMv+/fsjJiaGZUkaV6nfKA1+7fy125k0aRI6d+4MDw+Pl74vMjISdnZ2yoeLi0u15CPdsnPnTnz88ccoLi7GgAEDEBMTU2Y2UyJNqfQ/wdVxOVFYWBjOnz+PjRs3vvJ9ERERkMlkykd6errWs5Fu2bFjB/r164fi4mIMGjQIGzZsgKlple6LTfRSlf7NevPNN19bmo8fP1Y7UHh4OHbt2oWjR4/C2dn5le+1sLDgdXVGLj09HSUlJRg8eDB+/vlnliVpVaV/u2bPng07OzuNBxEEAeHh4di+fTsSExPRpEkTjW+DDE94eDiaN2+O7t27syxJ6yp1P0wTExPcu3cPderU0XiQsWPHIiYmBjt37ixz7aWdnZ3KN3bl/TCNQ1xcHDp06AB7e3uxo5CB0Mr9MLV5/HLFihWQyWTw8fFBvXr1lI/NmzdrbZukf7Zs2YLAwEB0796dV0VQtavUPow2z5JX1xl40l+bN29GcHAw5HI5WrRowWlwqdpVqjAVCoW2chC90qZNmxAcHAyFQoHQ0FBERUVBKpWKHYuMDK/sJZ0XExOjLMvhw4ezLEk0LEzSaVu2bEFISAgUCgVGjBiBNWvWsCxJNLwOg3Ra69atUbduXQQGBmLlypX8uiOJioVJOu3NN9/E2bNnUa9ePZYliY6/gaRzfv75Z8TFxSmfN2jQgGVJOkHjv4UmJibw8/NDUlKSpldNRmDt2rUYNmwY3n//fVy6dEnsOERlaLwwo6Oj0bVrV3z22WeaXjUZuOjoaIwcORKCIGDUqFFwd3cXOxJRGZX6aqSu41cj9deaNWswevRoAM/uVrVkyRJOtEfVRitfjSTShtWrVyvL8rPPPmNZks6q8lnyrKwsxMXFISMjAxKJBPXq1UNAQABvjEAqOXToEP71r38BAMaPH4/vv/+eZUk6q0ojzKioKHTo0AGnTp2CQqGAXC7HqVOn8O6773IaXlJJ165dMXjwYEycOJFlSTqvSscw3dzckJSUBGtr6zLLc3Jy0K5dO1y9erXKASuDxzD1hyAIynKUy+UwMTFhWZJoquUYpkQiQW5ubrnlubm5/OWnl1q2bBlGjhwJuVwOAJBKpfx9Ib1QpWOY3377Lbp27QoPDw80aNAAAHDnzh1cunQJ3333nUYCkmFZunSp8pKzXr16oX///iInIlJdlQqzT58+6NWrF86cOYO7d+9CEAQ0aNAAHTp04A0SqJzFixdjwoQJAIAvvvgC/fr1EzcQUSVVqjBDQkKwcuVK1KhRQ7lMKpXC29tb48HIsPzwww+YOHEigGezfc6bN4+74aR3KnUMMyYmpswxy08++QRPnjwp857i4mLNJCOD8f333yvLctq0aSxL0luVKswXT6hv3LixTGHev38fNjY2mklGBiEtLQ1TpkwBAEyfPh1z585lWZLeqtIxzIquSCoqKqrKKsnANGnSBFu3bkVKSgpmzJjBsiS9pvH7YfIPggBAJpMp568PCgpCUFCQyImIqq7S12HGxMQgOTlZeaySBUkv+vrrr9GyZUvcvHlT7ChEGlWpEWbnzp0xa9Ys5OTkwMzMDCUlJZg6dSo6d+4MT09PODk5aSsn6Yn58+dj2rRpAIDdu3dj/PjxIici0hy1vhp57do1JCUlITk5GUlJSUhJSUFWVlaZr7qJgV+NFNe8efMwffp0AMDcuXOV/02k61TtDrWOYbq6usLV1RWDBg1SLktLS8O5c+eQkpKizipJz82dOxczZ84E8Kw4p06dKnIiIs3jDYSpymbPno0vv/wSABAZGam8jIhIX2h1hElU6unTp9i+fTuAZyd7vvjiC5ETEWkPC5OqxMrKCgcPHsTu3bsxfPhwseMQaVWVp6hITk7mxepGRhAEnDlzRvnc0dGRZUlGocqF2b59e9y6dUsDUUgfCIKAGTNm4J133sFPP/0kdhyialXlXXIDOmdEryEIAqZNm4bIyEgAz45fEhkTHsMklQiCgIiICHzzzTcAnt2ujRelk7FhYdJrCYKAKVOmYMGCBQCe3Qi49K7pRMaEhUmvJAgCPv/8c3z77bcAnk0xERYWJnIqInGwMOmVJBKJ8g77P/74I8aNGydyIiLxsDDptb788ksEBgaiQ4cOYkchElWVLysiw1NcIsfcBT9Alp0D4Nkok2VJpIHCnDVrFhwdHTWRhXRAcYkcbt0GY+YXE9G8vQ+KikvEjkSkMzRSmA4ODprIQiITBAGjPg1D2pEtAABF0064k1Ugcioi3cFjmATgWVmGh4fj5zU/ARIJavcMR8fe/dHQocbrf5jISLAwCYIgICwsDMuXL4dEIsGq1avh13cgGjrUgKmUh7mJSlXqryEkJAT5+fnaykIimTJlirIso6OjMWrkSDR1smZZEr2gUn8RMTExyM3NVT7/5JNPysxLDkA5ORrpj0GDBsHR0RFr167FsGHDxI5DpLMqVZgv3mhj48aNZQrz/v37sLGx0UwyqjZt27bF9evXERoaKnYUIp1WpX2uiu5UxHtj6j6FQoHx48fjxIkTymWlc4gT0ctp/CAV5ynXbQqFAqNGjcKSJUvQp08fZGVliR2JSG9UujBjYmKQnJysPFapyYI8evQogoKCUL9+fUgkEuzYsUNj66Zn0x+PHDkSa9euhYmJCZYvX45atWqJHYtIb1TqsqLOnTtj1qxZyMnJgZmZGUpKSjB16lR07twZnp6ecHJyqlKYvLw8tG7dGsOHD8fHH39cpXVRWaVluX79ekilUmzYsAEDBw4UOxaRXlFrmt1r164hKSkJycnJSEpKQkpKCrKyspSjTblcXvVgEgm2b9+ODz74QOWf4TS7FZPL5Rg+fDh++eUXSKVSxMTEYMCAAWLHItIZWp1m19XVFa6urhg0aJByWVpaGs6dO4eUlBR1VklatHTpUmVZbtq0Cf369RM7EpFe0tg3fZo0aYImTZqgf//+mlrlaxUWFqKwsFD5PDs7u9q2rU/GjBmDhIQEDB06lIc6iKpAr78aGRkZidmzZ4sdQyfJ5XKYmJhAIpHA0tISO3bs4BUMRFWk1999i4iIgEwmUz7S09PFjqQTSkpKMGTIEPz73/9WXivLsiSqOr0eYVpYWMDCwkLsGDqlpKQEwcHB+O2332BmZobhw4ejZcuWYsciMgg6VZi5ubm4fv268nlaWhpSU1Ph4OCAhg0biphMPxQXFyM4OBhbtmyBmZkZtm7dyrIk0iCdKsxz587B19dX+XzSpEkAgNDQUKxbt06kVPqhuLgYgwcPRmxsLMzMzBAbG4ugoCCxYxEZFJ0qTB8fnwq/n06vVlxcjEGDBmHbtm0wNzdHbGws+vTpI3YsIoOjU4VJ6jl+/Dh27NgBc3NzbN++Hb179xY7EpFBYmEaAF9fX6xfvx61a9dGr169xI5DZLBYmHqqqKgIWVlZqFOnDgBgyJAhIiciMnx6fR2msSosLES/fv3QpUsXZGZmih2HyGiwMPVMaVnu3r0bf//9N65cuSJ2JCKjwV1yPVJQUICPP/4Ye/fuhaWlJXbv3g0fHx+xYxEZDRamnigoKMBHH32Effv2wcrKCrt374a/v7/YsYiMCgtTDxQUFODDDz/E/v37YWVlhT179sDPz0/sWERGh4Wp40rkCvz3egauX7+OGjVq4Pfff+duOJFIWJg6rESuwEfLT+B8hgyuId9gdVdH+HTtInYsIqPFwtRRT58+xda9h3A+49lt2a7lW6Chu6fIqYiMGy8r0kH5+fno27cvhg14H473kwAArZzt0NChhsjJiIwbR5g6Jj8/H0FBQTh8+DCsra3x7TAfNHT3REOHGjCV8t83IjGxMHVIXl4egoKCkJCQAGtra+zfvx+dOnUSOxYR/Q8LU0fk5eWhT58+SExMhI2NDfbv34+OHTuKHYuInsPC1AEFBQUIDAzEkSNHYGNjg7i4OHh7e4sdi4hewINiOsDCwgIeHh6wtbXFgQMHWJZEOkoiGNAtzrOzs2FnZweZTAZbW1ux41SKIAhIS0tD06ZNxY5CZHRU7Q6OMEWSnZ2NGTNmoKioCMCzaXBZlkS6jccwRZCdnY2ePXvi5MmTSE9P5wRvRHqCI8xqJpPJEBAQgJMnT8Le3h7h4eFiRyIiFXGEWY1Ky/L06dOwt7fHwYMH4enJrzsS6QuOMKtJVlYWevTogdOnT8PBwQGHDh1iWRLpGY4wq4EgCPj4449x5swZZVm2adNG7FhEVEkcYVYDiUSCWbNmoVGjRjh8+DDLkkhPcYRZTbp06YKrV6/C3Nxc7ChEpCaOMLXk8ePHCAgIwPnz55XLWJZE+o0jTC24/+Af+Pl3w+WL53Hnzh1cuHABJib8t4lI37EwNez+g3/wpmcnZGdcg7mNPWI2bmJZEhkI/iVr0MOHD+Hr54/sjGswqVkLtQfMg029JmLHIiIN4QhTQ/755x/4+/vjz0sXYGHjgNoD5qFdGw9OK0FkQFiYGjJt2jRcuHABb7zxBuIPHkKNOg05rQSRgWFhVlGJXIHbj/PxzYKFePz4Mb766iu89dZbYsciIi1gYVZBdk4uhqxLxfkMGVo1sMO2zb9xRElkwPjXrab79+/Dq317HNu6GgBwPkOG24/zRU5FRNrEwlTDvXv34Ovri2tX/kLBf/dBXpDLecOJjAB3ySsp/U4GfP38cOPaVTg7OyP+4CGYO9TnCR4iI8DCrIS0v9Ph0b4T8v9Jh6V9HRw8dBhub7qKHYuIqgmHRCq6nX5HWZZSWyc4DJgPM/t6YsciomrEwlRBiVyBFRu2K8uy7uBIeLg15zFLIiPDXfLXKJEr8NHyEzif1RRvBE6AibMHPNxcsXNcRx6zJDIyLMxXyMjIwN2cEpzPkAEALDy6Yf3w9ujU3JFlSWSE+Ff/Eunp6ejatSu6d+8BeUEuAKBlAzuWJZER419+BW6m3ULn97rixo0byMuRQSh6CgBYNKA1y5LIiHGX/AVXrt1A23ffw9PHmahRuz68wxfj+lMrtHK2QxPHmmLHIyIRsTCfc+XqdbR59z0UPLkH01r1YD9gPlaMCYDURMIL04mIhVnq0l9X0c67Cwqz7sPUvh7qDopEC9fGaOJYk0VJRAB08Bjm8uXL0aRJE1haWqJdu3Y4duyY1reZlVeA3j8koqSkGKb29VF3cCRcmzTEznGdWJZEpKRTbbB582ZMmDAB06ZNQ0pKCt577z306tULt2/f1sr2SuQKJN96hHZfHYKkVgPUHRyJuoMj0djFGfvGvwdLcw7Aiej/SARBEMQOUeqdd96Bp6cnVqxYoVz29ttv44MPPkBkZORrfz47Oxt2dnaQyWSwtbV95XtL5Ar4TI/BtZu3YNWkrXK5mYkEKTO6wdqKU+ISGQtVu0NnRphFRUVISkpCjx49yizv0aMHTpw4UeHPFBYWIjs7u8xDVbuOnMOpZRPxIHYOCv5+Nne41AQ4O82PZUlEFdKZwnz48CHkcjnq1q1bZnndunVx7969Cn8mMjISdnZ2yoeLi4tK27p+/TpC+wdBnvMQpnZ1YVb72c+djvBFrZqWVfsgRGSwdKYwS0kkkjLPBUEot6xUREQEZDKZ8pGenv7a9V+7dg1du3ZF7uMHMKvtgjf+XySk1vY4+O/OcLThzTSI6OV05qyGo6MjpFJpudHkgwcPyo06S1lYWMDCwkLlbVy9ehU+Pj7IzMyEu7s7Nu3YjXXJ2fhPz+Z4w866SvmJyPDpzAjT3Nwc7dq1Q3x8fJnl8fHx6NixY5XXn56erizLFi1aICEhAS1dm+K7gW1YlkSkEp0ZYQLApEmTEBISAi8vL3h7e2PVqlW4ffs2xowZU+V1169fH76+vrhw4QIOHToEJycnDSQmImOiU4U5cOBAPHr0CHPmzEFmZiY8PDywd+9eNGrUqMrrlkqlWL9+PXJycmBvb6+BtERkbHTqOsyqevFaqsuXL2PNmjVYuHAhpFKp2PGISEepeh2mTo0wNenSpUvw9fXFP//8A3t7e8yYMUPsSESk53TmpI8mXb58WVmWbdu2xdixY8WOREQGwCBHmIGBgXj8+DE8PT0RHx8PBwcHsSMRkQEwyBHm48ePlZcosSyJSFMMaoRZev6qVatWiI2NhampaaW+X05Exqm0J153DtygzpLfuXNH5e+TExG9KD09Hc7Ozi993aAKU6FQ4O7du7CxsXnp98+fl52dDRcXF6Snp7/2dnD6jJ/TcBjDZwSq/3MKgoCcnBzUr18fJiYvP1JpULvkJiYmr/zX4WVsbW0N+pevFD+n4TCGzwhU7+e0s7N77XsM8qQPEZE2sDCJiFRk1IVpYWGBWbNmVeoWcfqIn9NwGMNnBHT3cxrUSR8iIm0y6hEmEVFlsDCJiFTEwiQiUhELk4hIRUZbmMuXL0eTJk1gaWmJdu3a4dixY2JH0rijR48iKCgI9evXh0QiwY4dO8SOpHGRkZFo3749bGxsUKdOHXzwwQe4cuWK2LE0bsWKFWjVqpXyQm5vb2/s27dP7FhaFRkZCYlEggkTJogdRckoC3Pz5s2YMGECpk2bhpSUFLz33nvo1asXbt++LXY0jcrLy0Pr1q3x448/ih1Fa44cOYJx48bh1KlTiI+PR0lJCXr06IG8vDyxo2mUs7Mzvv76a5w7dw7nzp2Dn58f3n//fVy6dEnsaFpx9uxZrFq1Cq1atRI7SlmCEerQoYMwZsyYMsveeustYcqUKSIl0j4Awvbt28WOoXUPHjwQAAhHjhwRO4rW2dvbC2vWrBE7hsbl5OQIrq6uQnx8vNC1a1dh/PjxYkdSMroRZlFREZKSktCjR48yy3v06IETJ06IlIo0RSaTAYBB3wdVLpdj06ZNyMvLg7e3t9hxNG7cuHEIDAxEt27dxI5SjkHdfEMVDx8+hFwuR926dcssr1u3Lu7duydSKtIEQRAwadIkdO7cGR4eHmLH0bgLFy7A29sbBQUFsLa2xvbt2+Hu7i52LI3atGkTkpOTcfbsWbGjVMjoCrPUi7d/EwRBpVvCke4KCwvD+fPncfz4cbGjaIWbmxtSU1ORlZWF2NhYhIaG4siRIwZTmunp6Rg/fjwOHDgAS0tLseNUyOgK09HREVKptNxo8sGDB+VGnaQ/wsPDsWvXLhw9elStW/zpA3NzczRv3hwA4OXlhbNnz2Lx4sVYuXKlyMk0IykpCQ8ePEC7du2Uy+RyOY4ePYoff/wRhYWFok+XbXTHMM3NzZXz/TwvPj4eHTt2FCkVqUsQBISFhWHbtm04fPgwmjRpInakaiMIAgoLC8WOoTH+/v64cOECUlNTlQ8vLy8EBwcjNTVV9LIEjHCECQCTJk1CSEgIvLy84O3tjVWrVuH27dsYM2aM2NE0Kjc3F9evX1c+T0tLQ2pqKhwcHNCwYUMRk2nOuHHjEBMTg507d8LGxka552BnZwcrKyuR02nO1KlT0atXL7i4uCAnJwebNm1CYmIi9u/fL3Y0jbGxsSl37LlmzZqoXbu27hyTFvckvXiWLVsmNGrUSDA3Nxc8PT0N8jKUhIQEAUC5R2hoqNjRNKaizwdAWLt2rdjRNGrEiBHK31cnJyfB399fOHDggNixtE7XLivi7d2IiFRkdMcwiYjUxcIkIlIRC5OISEUsTCIiFbEwiYhUxMIkIlIRC5OISEUsTCIiFbEwiYhUxMIkUsOjR49Qp04d3Lp1S+woAIB+/fph0aJFYscweCxMKqNLly6QSCTlHsHBwXqzjWHDhkEikVR4M5WxY8dCIpFg2LBhyvd+8MEHld5GZGQkgoKC0Lhx46qFfYWgoKCX3nX85MmTkEgkSE5OBgDMnDkT8+bNQ3Z2ttbyEAuTniMIAlJTU/Htt98iMzOzzENT91ysjm0AgIuLCzZt2oSnT58qlxUUFGDjxo1VvlPT06dPERUVhVGjRlU15iuNHDkShw8fxt9//13utejoaLRp0waenp4AgFatWqFx48bYsGGDVjMZOxYmKV27dg05OTno0qUL3njjjTIPa2trvdkGAHh6eqJhw4bYtm2bctm2bdvg4uKCtm3bVmnd+/btg6mpaZn5dHx8fBAeHo4JEybA3t4edevWxapVq5CXl4fhw4fDxsYGzZo1Kzc1riAIWLBgAZo2bQorKyu0bt0aW7duBQD06dMHderUwbp168r8TH5+PjZv3oyRI0eWWd63b19s3LixSp+NXo2FSUpJSUkwNTXV6tSm1bGNUsOHD8fatWuVz6OjozFixIgqr/fo0aPw8vIqt3z9+vVwdHTEmTNnEB4ejk8//RT9+/dHx44dkZycjICAAISEhCA/P1/5M9OnT8fatWuxYsUKXLp0CRMnTsSQIUNw5MgRmJqaYujQoVi3bh2ev6nYli1bUFRUVO4QRocOHXDmzBmDuqmwrmFhklJycjLkcjlq164Na2tr5WP06NHVuo09e/bAzc0Nrq6uWLNmjdrbCgkJwfHjx3Hr1i38/fff+OOPPzBkyJAqf4Zbt26hfv365Za3bt0a06dPh6urKyIiImBlZQVHR0eMHj0arq6umDlzJh49eoTz588DeDZv/KJFixAdHY2AgAA0bdoUw4YNw5AhQ5SHJ0aMGIFbt24hMTFRuZ3o6Gh89NFHsLe3L7P9Bg0aoLCwkJP5aZFR3nGdKpaUlIT+/ftj3rx5ZZa/+IcJAF9++SVmz579yvWdPXu23EjsddsoKSnBpEmTkJCQAFtbW3h6euKjjz5Sa9pcR0dHBAYGYv369RAEAYGBgXB0dKz0el709OnTCifpen7ULJVKUbt2bbRs2VK5rHTOqAcPHgAALl++jIKCAnTv3r3MeoqKipSHDd566y107NgR0dHR8PX1xY0bN3Ds2DEcOHCg3PZL7zD//AiWNIuFSUopKSmYM2eOcqKtVwkLC8OgQYNe+Z6KziC/bhtnzpxBixYt0KBBAwBA7969ERcXh8GDB7/+A1RgxIgRCAsLAwAsW7ZMrXW8yNHREU+ePCm33MzMrMxziURSZlnprKQKhaLM//7+++/Kz1vKwsJC+d8jR45EWFgYli1bhrVr16JRo0bw9/cvt/3Hjx8DAJycnNT5WKQCFiYBAG7evImsrCyVT4g4OjpWerSmyjbu3r1bpjycnZ2RkZFRqe08r2fPnigqKgIABAQEqL2e57Vt2xa//vprldfj7u4OCwsL3L59G127dn3p+wYMGIDx48cjJiYG69evx+jRoyucEvrixYtwdnbWyCiaKsbCJADPdpWBZ7uNLx4Dq1OnDkxMqn64W5VtVDRjSlXmi5dKpfjzzz+V/10RmUyG1NTUMsteNVFcQEAAIiIi8OTJkwoPV6jKxsYGkydPxsSJE6FQKNC5c2dkZ2fjxIkTsLa2RmhoKADA2toaAwcOxNSpUyGTyZTXkL7o2LFj6NGjh9p56PVYmAQAygug33zzzTLLzczMkJOTU2YXUZvbaNCgQZkR5Z07d/DOO+8on69btw7Dhw+vsFhfxtbW9pWvJyYmlhv1hoaGlrucp1TLli3h5eWF3377DZ988onKOSoyd+5c1KlTB5GRkbh58yZq1aoFT09PTJ06tcz7Ro4ciaioKPTo0aPCIi8oKMD27dsRFxdXpTz0apwEjXRKSUkJ3n77bSQmJipP+pw6dQq1a9cG8OxkU2JiYpmzxmLYu3cvJk+ejIsXL2pk9F1Vy5Ytw86dOys8GUSawxEm6RRTU1N899138PX1hUKhwOeff64sSwCIi4vD4sWLRUz4TO/evXHt2jVkZGTAxcVF7DgwMzPD0qVLxY5h8DjCJCJSkfj7EkREeoKFSUSkIhYmEZGKWJhERCpiYRIRqYiFSUSkIhYmEZGKWJhERCpiYRIRqYiFSUSkIhYmEZGKWJhERCr6/wkuKz21tE6NAAAAAElFTkSuQmCC", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAUwAAAFQCAYAAADUeke/AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAA6W0lEQVR4nO3dd3RU1drH8e+ZSe89EBKk1xAgBDQQKdKkSREQBAQpooYmokCkXECI16t4pYhIVynqpfmqgCgECIqUEGlKEZCSRkvPZJjy/hEZibSUSWaSPJ+1snTOHOY8k0x+2Wfvc/ZWjEajESGEEI+ksnQBQghRVkhgCiFEAUlgCiFEAUlgCiFEAUlgCiFEAUlgCiFEAUlgCiFEAUlgCiFEAdlYugBzMhgMJCQk4OrqiqIoli5HCFFGGI1GMjIyCAgIQKV6cDuyXAVmQkICQUFBli5DCFFGXb58mcDAwAc+X64C09XVFch7025ubhauRghhzT7//HMiIyMBGDJkCJ999pkpQx5EKU/3kqenp+Pu7k5aWpoEphDigVauXMnIkSMxGo1ERkYyd+5cPDw8HpkdMugjhKhQVqxYwYgRIzAajYwZM4aFCxcWeMxDAlMIUaHY2dmhKArjxo1jwYIFhRogLld9mEII8ShDhgyhbt26NG/evNBX00gLUwhR7q1du5bExETT4xYtWhTp0kMJTCFEufbRRx8xePBgnnrqKdLS0or1WhKYQohya/HixaZLh7p3717sq2ckMIUQ5dLChQsZM2YMAG+88Qbvvvtuse8AlMAUQpQrOr2BGfP+w7hx4wCYPHky//73v81yu7SMkgshyg2d3kDYC2/x67p3AJg8eQrR0fPMNreEBKYQoty4dDOb6+71sPGohFP91oyaGGXWiXgkMIUQ5UZVLydC69ckfuh/aVwzgMe8nc36+hKYQogy7/333ycoKIj+/fuz6dWWXLqZTVUvJ2zU5h2mkcAUQpRp7777LpMnT0atVtOwYUMaNmxIDV+XEjmWjJILIcqsd955h8mTJwMwffp0GjZsWKLHk8AUQpRJ8+bNY+rUqQDMnj2bmTNnlvgxJTCFEGXO22+/zVtvvQXAnDlzmD59eqkcV/owhRBlyvbt200BOXfuXKKiokrt2BKYQogypXPnzrz66qsEBQUxZcqUUj22BKYQokzQ6/Wo1WoURWHRokUWWRlW+jCFEFbNaDQyc+ZM+vbti1arBbDYMtrSwhRCWK07YTlnzhwAtm3bRs+ePS1WjwSmEMIqGY1Gpk+fzty5cwGYP3++RcMSJDCFEFbIaDTy1ltvER0dDcAHH3zAhAkTLFsUEphCCCtjNBqJiorinXfypmj78MMPTXNbWpoEphDCqly4cIEFCxYAsGDBAsaOHWvhiv4mgSmEsCo1atTgu+++49SpU7zyyiuWLicfxWg0Gi1dhLmkp6fj7u5OWlpasRc7EkKUHqPRSGJiIgEBARY5fkGzQ67DFEJYlNFo5PXXX6dx48acOHHC0uU8lASmEMJijEYjr732Gh988AHXr1/n0KFDli7poaQPUwhhEUajkQkTJpgGeJYuXcqLL75o4aoeTgJTCFHqjEYj48aNY9GiRQAsW7aMkSNHWriqR5PAFEKUKqPRyNixY1m8eDGKorB8+XKGDx9u6bIKRAJTCFGqNBoN8fHxKIrCihUrrP40/G4SmEKIUuXo6Mi2bduIiYmhR48eli6nUGSUXAhR4gwGA9u2bTM9dnV1LXNhCVYcmNHR0SiKYhU33Ashis5gMPDyyy/TtWtX0/3hZZVVnpIfOnSITz75hJCQEEuXIoQoBoPBwEsvvcSKFStQqVQEBgZauqRisboWZmZmJoMGDWLZsmV4enpauhwhRBEZDAZGjRplCstPP/2UwYMHW7qsYrG6wIyMjKRbt2506NDhkfvm5uaSnp6e70sIYXl6vZ4RI0awcuVKVCoVn3/+OYMGDbJ0WcVmVafkGzZsIC4ursC3R0VHRzNr1qwSrkoIURhGo5GRI0eyevVqVCoVa9euZcCAAZYuyyyspoV5+fJlxo8fz+eff46Dg0OB/s3UqVNJS0szfV2+fLmEqxRCPIqiKDRp0gS1Ws26devKTViCFU3vtmXLFnr37o1arTZt0+v1KIqCSqUiNzc333P3I9O7CWE9zp07R61atSxdRoGUuend2rdvz/Hjx4mPjzd9hYWFMWjQIOLj4x8ZlkIIy9HpdMyaNYvU1FTTtrISloVhNX2Yrq6uBAcH59vm7OyMt7f3PduFENZDp9PxwgsvsH79er7//nv27duHSmU1bTGzsprAFEKUPZpcLf0GDuKbzf/DxsaGN954o9yGJVh5YMbExFi6BCHEA2hytdRs2YWEuF0oahu++PJLevXqZemySlT5/VMghCgxOp2OZ58bSELcLlDZ4NNzKk0iOlq6rBJn1S1MIYT10ekNDB8dyXdbN6FS2+LdaypPtO1IVS8nS5dW4iQwhRAFptMb6PPRTxxRheLoVZm1Kz6mcaunqOrlhI26/J+wSmAKIQrs0s1sjl1Nw9Y7EN8XP6Jxq6eo4eti6bJKTfn/kyCEKDatVsvAgQP5/dA+Qqq4A9C4mk+FOA2/m7QwhRAPpdVq6devH19//TXbt2/n3B/nSdPbVpjT8LtJYAohHig3N5e+ffvyzTff4ODgwBdffIG3lyfeli7MQiQwhRD3lZuby7PPPsu3336LYmNHk+HzeKr9o6ddLM8qVntaCFEgGo2GPn36mMLS99kZJLrW4dLNbEuXZlHSwhRC3OOjjz7iu+++w9HRkcbD55HoUpuQQPcKN8jzTxKYQoh7jBs3jpMnTzJo0CBat2nLpZvZFXKQ55+sZj5Mc5D5MIUousysbBIzblPd17XCBWOZmw9TCGE5t9IyqBbahqYd+9J7USw6vcHSJVklCUwhKrj0jExqtWjHjTOHyT4dS9zxUxV+cOdBJDCFqMCys7Pp0q0HN88cQbF1wK/fvwhuWL/CD+48iASmEBVUekYmHZ7uyk/7YlDbO+LXfxaNw8LZGtmqwvVhFpSMkgtRAaWlZ1AzrC03zsahtnfk++3bqNawmYyEP4IEphAV0LaY/dz44xiKnSM+fWdTrWGzCjXrUFFJYApRAfXt1on5I97mco6a5o8/IX2WBSSBKUQFkZmZyY0bN3jsscewUav4ackbckF6Icl3SYgKICMjg6ef7kJ4qwjOnD0HgI1aRQ1fFwnLQpDvlBDl3K3UNNp26MT+/bEkXb/FyE9i5ML0IpLAFKIcu3krlZqhEcQdPIDK3hn/597mkrqyXJheREXuw7x9+zZJSUlkZ2fj6+uLl5eXOesSQhRTWloaHTp24taFE6jsnfEbMBf7SrVk1qFiKFRgZmZmsnbtWtavX8/BgwfJzc01PRcYGEinTp146aWXaN68udkLFUIUXPK1G7Tr2Inffo3D1skV735zCGsWyvz+Taju4yz9lkVU4MD84IMPmDt3LtWqVeOZZ55hypQpVKlSBUdHR27evMmJEyfYt28fHTt25IknnmDhwoXUrl27JGsXQtxHZo6WVtE/cCnxFipHV37Y+T2BtYNlNNwMCjy9W79+/ZgxYwaNGjV66H65ubmsWLECOzs7Ro4caZYiC0qmdxMVlU5v4ML1LHJv6xj9WRxX0zTos9PQZ96kbr2G7Hy9jYTlQxQ0Owo1H2Z8fDxNmjQxR30lQgJTVEQ6vYHei/fz6/kEci7G41wv4p59dr3eRu7keYgSmQ8zNDSUZs2asWTJEtLS0opdpBCi+C5czyL+fALJX0zj+tZ3yDy2EwBbtQIggzxmVKjA3L9/P6GhoUyZMoXKlSszePBgdu/eXVK1CSEeQac38OqKPaRseAtt0jlUjm7YVa7NY54OHJ3WgV2vt2HTKy3ldNxMCvVdDA8PZ9myZSQlJbFkyRKuXLlChw4dqFmzJnPnzuXKlSslVacQ4j4On77E3gUT0Cb/gcrJHf+B87DzrcbHQ8JwcbSTO3nMrEjfSUdHR4YOHUpMTAxnzpxh4MCBLF26lOrVq9O1a1dz1yiEuI8/ryTSoUPHv8LSA/8BeWEZHOBGbX9XS5dXLpllEbQ712dGRUWRmpqKXq83R22FJoM+oiLQ6Q0cu5hMy1YR5CafR+WcF5a169Rl6ZAwavlJq7KwCpodxZqtaM+ePaxcuZKNGzeiVqvp378/I0aMKM5LCiEe4s6I+LGraTjUbMHtzJtUGhBN9Vp12DGhNQ52MgFZSSr0d/fy5cusXr2a1atXc+HCBVq2bMnChQvp378/zs7OJVGjEOIvvyemczwhHUVRcI8YhGuzHji5ebB9fISEZSko1He4Y8eO7N69G19fX1544QWGDx9O3bp1S6o2IcRfdHoD+4+f45kXx+P+1EhUtg4oioLayZ2tkRG4ONpZusQKoVCB6ejoyMaNG+nevTtqtbqkahJC3EWnN9D1nW+I+WAst29cQqfNwbfHGwA0quJGLT+5IL20FCowv/7665KqQwhxHzq9ga0/nzKFpdrFG49Wz2Orgs2vtqJeZTcZ4ClFRf5O79u3j8GDBxMeHs7Vq1cB+Oyzz4iNjTVbcUJUZBqtjvZvb2Fgr66msPR/PpratWpxdHpHggM9JCxLWZG+2xs3bqRz5844Ojpy9OhR0zRvGRkZzJs3z6wFClER6fQGOs/byk8LxnP7xmXUrj74Px/NyrHd+eH1ttJnaSFFCsy3336bjz/+mGXLlmFra2va3rJlS+Li4sxWnBAVUWaOlsW7zvLL8unobl5B7eqL/8BobD0DsLNRSavSgop0HcLp06dp3br1Pdvd3NxITU0tbk1CVFjXM7IJm5s3P4Nnh9Hc2LYA395R2HpUwtFWTXgNbwtXWLEVKTArV67MuXPnqFatWr7tsbGx1KhRwxx1CVHhZOZoeXzuLiBvliH7SrWoPOy/BFfx4PWOdWhZ01uutbSwIrXtR48ezfjx4/nll19QFIWEhATWrl3LpEmTePXVV81doxDlnk5voMOcLVxZPYHcq7+Zts9+Jpitka14qr6/hKUVKNJP4M033yQtLY127dqh0Who3bo19vb2TJo0iTFjxpi7RiHKNZ3ewOZ9v3JkyQR0qYnc2LGYyi8uwEZR8VzzIOmztCLFmnwjOzubU6dOYTAYaNCgAS4ulr2AVibfEGWNRquj09ub+HnheHSpSajd/ak0MBobdz++jmxFSJCHpUusEEpkxvW7aTQaTpw4QVJSEklJSezatYuvv/66yBe3L1myhJCQENzc3HBzcyM8PJxt27YVtTwhrF5mjpYnZ3zJzwvGoUtNwsbdn2avzMfG3Y/gAFcaBMgffWtTpFPy7du3M2TIEG7cuHHPc4qiFGl6t8DAQN555x1q1aoFwJo1a+jZsydHjx6lYcOGRSlTCKul0epo9MY6rnw+BV1aMjYelag0MJod0/pwM1snKzxaqSL9RMaMGUP//v1JTEzEYDDk+yrqXJg9evSga9eu1KlThzp16jB37lxcXFw4cOBAkV5PCGv2428pXI/94q+wrIz/wGjUbr78eiVdZkm3YkVqYaakpDBx4kT8/f3NXQ8Aer2er776iqysLMLDwx+4X25urukuI8jrhxDCmun0BuL+vE7k+qN4dXgJFAX38OewcfPB0VbN49W9LF2ieIgiBWbfvn2JiYmhZs2aZi3m+PHjhIeHo9FocHFxYfPmzTRo0OCB+0dHRzNr1iyz1iBESdHpDXT597ecSVNQFAXFxg7vzpEARPcOpnfTKnLpkJUr0ih5dnY2/fr1w9fXl0aNGuW7PRJg3LhxRSpGq9Vy6dIlUlNT2bhxI8uXL2fPnj0PDM37tTCDgoJklFxYpR8PHqNzxw441W2FZ7sRKEreBeoNKjnz9djWchpuQQUdJS9SYC5fvpyXX34ZR0dHvL29TT94yBv0OX/+fNGq/oc7K1IuXbq0QPvLZUXCGmm0Otb/cJBXnu9Fbto1bLwCqTzkPVQOLvy3f2O6Nw6QsLSwEl3TZ9q0acyePZspU6agUpXcD9poNOZrQQpR1mTmaGk06XMufzYZfcZ1bLwC8R84D5WDC7V8nCQsy5giBaZWq+W5554za1hGRUXRpUsXgoKCyMjIYMOGDcTExLB9+3azHUOI0qTTG3hq5pdc/vRN9Jk3sPEKpNLAaNQuntTxc+HrMa0kLMuYIgXm0KFD+eKLL4iKijJbIcnJyQwZMoTExETc3d0JCQlh+/btdOzY0WzHEKK0ZOZomffFbuKWjEefeRNb76r4D5yLvYsnWyJlpvSyqkiBqdfreffdd9mxYwchISH3DPrMnz+/0K+5YsWKopQihNXJzNESPGsnWb/tRZ95yxSWgZX82TmxjUz+W4YVKTCPHz9O06ZNAThx4kS+5+4eABKiotHpDSzedQ4A5/qtUdQ22FdpQPUq/nw/sa1cNlTGFemnt3v3bnPXIUSZl5mjpe30dSRp1Ni45F2A7lSnJQDTewRLWJYD8hMUwgyuZ2TT5LVPSdoQhWLvjP/AeabQdLBReLK2j4UrFOZQ4F7nS5cuFeqF76wkKUR5l5mjJeS1NSRumIo+6xaKjS2KSg3AwgGNiZ/RSVqX5USBA7N58+aMGjWKgwcPPnCftLQ0li1bRnBwMJs2bTJLgUJYM53ewL8+3UHy+igMWanY+tXAf8BcOjatyYmZHenRJFDCshwp8E/yt99+Y968eTz99NPY2toSFhZGQEAADg4O3Lp1i1OnTnHy5EnCwsL4z3/+Q5cuXUqybiEsSqc38HtiGs9Hr+fk6igM2WnY+dfE77k5qB3dmNChroyGl0OFvjVSo9Hw3XffsW/fPi5evEhOTg4+Pj40bdqUzp07ExwcXFK1PpLcGilKg05voOfi/Rz99TjJG+4Oy7dRO7pSz8+Jb8a3kessy5ASuzXSwcGBPn360KdPn2IVKERZdTY5g5MJ6ajsnVHZOWHj5otf/zmoHV35eFAoHRr4S1iWU9K5IkQB6fQGLlzP4uVP8/rxbdx88B8YjWLngNrBhbp+ThKW5ZwEphAFoNMb6LkolqPxv3L7VgLO9SKAvNAE+M+zDekdWlXCspyTwBTiEXR6A9uOJ3I0Pp7kDW9hyM1GZeeIY41mADSo7CphWUFIYArxEHdalnFHj5LyxTQMmkzsKtfFvko9AD4b0YzwGn4SlhVEkX/K+/btY/DgwYSHh5suUv/ss8+IjY01W3FCWNqF61n5wzKgLv7PzUFl78y6EY/zZO1KEpYVSJF+0hs3bqRz5844Ojpy9OhR0yS/GRkZzJs3z6wFCmEpGq2OlZt2krLhLQyaTOwD6uHffw4qeydqeNnTooYsWFbRFCkw3377bT7++GOWLVuWb2q3li1bEhcXZ7bihLCU1CwN9SZ8yruvvYAhNwv7Kg3w6z8blb0TAIsHN5eWZQVUpD7M06dP07p163u2u7m5kZqaWtyahLCo1CwNTeb8iNHVD+eGbdGmXMCv779MYdmgkgu1/V0tXKWwhCIFZuXKlTl37hzVqlXLtz02NpYaNWqYoy4hLOJOWELe3K6eHV7GqNOisrUH4H+jw2lS1UNalxVUkX7qo0ePZvz48fzyyy8oikJCQgJr165l0qRJvPrqq+auUYhSkZqlof7Li7n+zfsY9beBvNBU2drjYa8QP709YdW9JCwrsCK1MN98803S0tJo164dGo2G1q1bY29vz6RJkxgzZoy5axSixCWlZdJ0zMckfzkDozYHG4/KeEQ8D4Cfk8LeKTJFmyjiuuR3ZGdnc+rUKQwGAw0aNMDFxcWctRWaTL4hCkunN/DLH9d5dvYaUr6aiVGbg33VRvg9OxOVnQMuNnB4RmcJy3KuRNclv8PJyYmwsLDivIQQFqPTG+j24R5+jTtEylf/wqjNweGxEHyfnYHK1gGAmMntJCyFSZE+CbNnz37o8zNmzChSMUKUFp3ewNdHr/DrkYN5YXlbg8NjjfF9djoqWwda1vDko0GheDg7WLpUYUWKFJibN2/O9/j27dtcuHABGxsbatasKYEprJpOb6D7h3s4dfUm17a+81dYNsH32WmmluUnQ8JkAmBxjyIF5tGjR+/Zlp6ezrBhw+jdu3exixKiJP3yx3V+T8lGZeuA7zOTST+8FZ8eb5guHTr8VjsJS3FfxRr0+acTJ07QvXt3Ll68aK6XLBQZ9BEPo9Hq2HjkKlGbjqKobe953s/Fll2vt5WwrIBKZdDnn1JTU0lLSzPnSwphFhqtjsazdpB2/hjXv/svfn1nYOdbzfS8l6OKvW8+JQM84qGK9OlYsGBBvsdGo5HExEQ+++wznn76abMUJoQ5ffNrIqnnf+Xa/2Zj1OWS/stGfLq/bnp+1yQZDRePVqRPyAcffJDvsUqlwtfXl6FDhzJ16lSzFCaEuVzPyCby/U+5tnE2Rp0WxxpheD89DoBREdUY3762nIaLAilSYF64cMHcdQhRIq7cyqDFyx9ybdOcvLCs2RzfXlEoNrZU87Rjcpf6cqujKLBCf1Ju375Nu3btOHPmTEnUI4TZnE5OpdnL/yXpTljWamEKy0rOKra/1k7CUhRKoVuYtra2nDhxAkVRSqIeIczi3LU0Os2PJf2XTX+F5eP49pqCorbFx1FNzOQO0mcpCq1If15feOEFVqxYYe5ahDCLizfS6fB+LIqi4Ns7Cvfw50xhCRAjo+GiiIr0qdFqtSxfvpydO3cSFhaGs7Nzvufnz59vluKEKKwrtzJoFfUltp6VAVDZO+HReggAlV3U7Hz9KRngEUVWpMA8ceIEoaGhAPf0ZcqpurCUpLRMQl96n2ub5+LR6nncw/ubnls8oDGdGwVIn6UolkIF5vDhw/nwww/ZvXt3SdUjRKFptDp2/ZbMsDmfcG3zPDDo0Cb/gdFoQFFUzO8XTLcmgZYuU5QDhfpzu2bNGnJyckqqFiEKTaPV0WjmDobNWWoKS6e6Efj0eANFyft4NwqU1R2FeRSqhWnG286FKLbULA3jNhwh7cwBrm15Jy8s6z2JT49JKCo1AMEBblT3cX7EKwlRMIXuw5Q+SmEN7ixWln3mZ65tfQcMepzqt8an++umsPxmTCvqVXaTfkthNoUOzDp16jwyNG/evFnkgoR4FI1Wx4QNeVMM6jKu/xWWbfDpPhFFpWZkeBCTujSQS4eE2RX6EzVr1izc3d1LohYhHkmj1VFvxg7TY7dmPbD1qIxD9aYoKjXVPe2Z0j1YWpWiRBQ6MAcMGICfn19J1CLEI63ee4Gc80ewC6iL2iFv0T3HmnnrSvk527LttbYSlqLEFOqTJf2XwlIyc7R8+N1vTF+0kpT/zSLli2kYcrPz7bNq+ONyGi5KlIySC6uXmaMleNZOsn7by/X/ew+MBmx9qqL8taQEQF1fR+pWcrVglaIiKFRgGgyGkqpDiHvo9AYuXM/ine2/knVqD9e/eR+MBpyD2+PdZZxpNHzdiMdpUcNLTsVFiZPzF2GVNFodzyyK5UxKFlmnYrj+zfy8sGzUAe+nx5rCMnZyawI9pWUpSocEprA6Or3h77D8PfausOyId5expjt4JCxFabOac5jo6GiaN2+Oq6srfn5+9OrVi9OnT1u6LGEB51IyOZOSBYCdX3XUzh64hHQyhaW7nUL89PYSlqLUWU1g7tmzh8jISA4cOMDOnTvR6XR06tSJrKwsS5cmSlFmjpbnFu8zPbb1qkKlF+bj9fQYFEVFgIuaX6Z1wsPZwYJViorKrOuSm9O1a9fw8/Njz549tG7dukD/RtYlL9s0Wh3BM3eQevxH1E4eONZolu95bycb9k9pL5cOCbMraHaYvYWpUql46qmnOHLkSLFe58765l5eD55pJjc3l/T09HxfomzSaHVEf3uS1F93cuPb/5Ky6W201/40Pa8C9rwhS+EKyzJ7YK5cuZI2bdowbty4Ir+G0Whk4sSJREREEBwc/MD9oqOjcXd3N30FBQUV+ZjCcq5nZFNvxg4WL1vJjW0LACOujTth61MVgKHNq3BqdmeZKV1YnFWekkdGRvLtt98SGxtLYOCDJ37Nzc0lNzfX9Dg9PZ2goCA5JS9D7lyUnvHrDm5uXwiAa2h3PDuMRlEU+jTyZf6gFhauUpR3BT0lt7rzm7Fjx/L111+zd+/eh4YlgL29Pfb29g/dR1gvnd7Af3b8Tkb8dm7uWASAa7MeeLZ/yXQb7vA2dSxZohD5FDswU1NT2bFjB1evXkVRFCpXrkznzp3x9PQs1OsYjUbGjh3L5s2biYmJoXr16sUtTVgxnd5Ap/kxnDzy811h+Qye7UeZwrJhZVfqVZYzBWE9itWHuWLFClq0aMGBAwcwGAzo9XoOHDjAE088UehleCMjI/n8889Zt24drq6uJCUlkZSUJEtilFPbTyZw/kYODlUb4VS/Da5hPfOF5bZxEWwdEyG3OwqrUqw+zLp163LkyBFcXFzybc/IyKBZs2b3rCj50EIeMBPSqlWrGDZsWIFeQy4rsn6ZOVo+++lP3vn+tOlnbjToQVGhKAqRraoytnN9GQ0XpapU+jAVRSEzM/OewMzMzCz0VHBWOPYkzMw0wBP3DblJ50z3hN+5L3xA88q80aORhasU4sGKFZjvvfcebdq0ITg4mCpVqgBw5coVTp48yfvvv2+WAkX5oNMbmL7lGOlH/o9bPywFwLFGGM71Ikz7DHmipqXKE6JAihWY3bt3p0uXLhw8eJCEhASMRiNVqlShRYsWqNVqc9UoyrjrGdlEzN1NyuGt3PpxGQBuj/fFqW4r0z71/Z1kPkth9QoVmEOGDGHp0qU4OTmZtqnVasLDw81emCgfUrM0hM3dTfqhrdza9VdYPtEPj9YvoCgKLz1ZjT6hQdTyc5EBHmH1CvUJXbduHZmZmabHo0eP5tatW/n2uX37tnkqE2VeZo6WyLVHSD+05e+wDH/OFJY1fBx48+n6shSuKDMK9Sn958DM+vXr8wVmcnIyrq5yWiXyTsODZ+0kJu53bu1ZDYB7+HN4PDkYRVFoW8eL719rJ0EpypRi9WHeb2Rbq9UW5yVFGafTGzidlEG3hbEA2HpUwrfXVLTJ53FvOcB09cR/n2sqYSnKHLNf7CYrS1ZcOr2Bngv3cTIpE0NuFip7ZwCcaj2OU63HAXAE9r3VTuazFGVSof/Er1u3jri4OFNfpQSkuONUQjonkzJJO/AVCSvGcDs1Kd/z3UJ8OT63Cz6uTg94BSGsW6FamBEREcycOZOMjAxsbW3R6XRERUURERFBaGgovr6+JVWnsHKZOVqeWbyftJ+/JHXvpwDknPsF27Cepn3+3buJnIaLMq1It0aePXuWI0eOEBcXx5EjRzh69Cipqamm1qZerzd7oQUht0aWvjt9lt0XxpL60xek7vsMAPcnB+PRcgAAzmrYH9VeTsOF1SrRWyNr165N7dq1GTBggGnbhQsXOHz4MEePHi3KS4oySKc30HNRLCcTM0jdv5602LUAeLR+Affw/gC0r+vB0hfCpWUpygWzDfpUr16d6tWr069fP3O9pLByZ5Mz8sIydh1p+9cB4NFmKO5P/P0ZWDzocQlLUW7IlDCiSDJztPRcEIvhdi7ZZ38GwKPNMNyf6AvAyJZBTHq6gcw6JMoV+TSLQtNodQTP2gmAytYe/+feJufcQVxCOgIQ4AxTugVLy1KUO8X+RMfFxcnF6hWITm9g1jfHyE04bdqmdnI3hWUlV3t2Te4sYSnKpWJ/qps3b87FixfNUIqwdqlZGprN2sZH771L0mevk3H0u3zPB7jZEvNGWzkNF+VWsT/ZMvFvxZCapaHx7B9I3fsp6Qe+AsCo+/vMYmR4EFO6y2m4KN+kKSAeSac3ED7nB1L3rCH9l/8B4Nl+FG53XZQ+qUsDCUtR7klgiofSaHVM2xxPQsxq0n/ZCIBn+5dwC3vGtM+BqW3kNFxUCPIpF/el0xs4lZBOr0Wx3IhZRfrBTQB4dhiNW7Mepv0Ov9VO7g0XFYYEpriHTm/gmYX7OJWUCYqCYmMPgFfHl3EN7W7aL3663O4oKhYJTHGPk1fT8sLyL+4Rz+NYMwz7gLoAOChwYJqEpah4pJde5JOUlskzi/eTcfQ7DFoNkDeF352wtFHB4RkdJSxFhVTswJw5cyY+Pj7mqEVYWGqWhsfnxXDrx2Xc/P4jrm2ajdFoMD3vqob46R1xcbSzYJVCWE6xT8lnzpxpjjqEhen0Bl759Gdu/fgJGUf+DwCn+m1QlL//pv48TcJSVGzShynQaHU0mb6dhB8+JiPuW0DBu8tYXEI6AeBuC/ujJCyFkMCs4O7cwXNz58dkHr0TluNM94Z7Oaj4KaqjXGcpBIXswxwyZAjZ2dklVYsoZalZGprM+ZHUPav/Dsuu401hqQB7J7eXsBTiL4UKzHXr1pGZ+fflJqNHj863LjlgWhxNWLc7YQngXL81Kkc3vLtOwKVRB9M+28Y/KafhQtylUIH5z4k21q9fny8wk5OTcXV1NU9lokRotDq2xF0xhSWAnX9Nqoxehkuj9qZtjaq4UcvPxRIlCmG1inWudb+ZimRuTOul0eoInrGD20YDt35chlO91jgE1gcwrSEOsHRwKO3r+8tkGkL8g9l/I2Sdcuuk0xtY+sNZbhsN3Ni2gIwj/8e1jbMwaDLz7dewsouEpRAPUOgW5rp162jdujWNGjUCJCDLAo1WR5cPYjh/I5sb2xaSdeIHUFR4dXwFlUPeabebDXwZ+SS1/FwkLIV4gEIFZkREBDNnziQjIwNbW1t0Oh1RUVFEREQQGhqKr69vSdUpikinN9Bj0f6/wnIBWSd+BEWFT49JONdvbdrvp7fkOkshHqVQgbl3714Azp49y5EjR4iLi+PIkSNMnz6d1NRUaW1aoWOXb3EmKY0b3/2XrJO7/wrLN3Cu/6Rpn/jp7SUshSiAIg361K5dm9q1azNgwADTtgsXLnD48GGOHj1qtuJE8VzPyKbPxwfIOPLN32H5zJs414sw7SNTtAlRcIqxHC3Kk56ejru7O2lpabi5uVm6HIvR6Q0cunCNgcsPA3lr71zb+m+cg5/CuW4r034HprahkrtcOiREQbNDbuEoZzRaHV0X7OOPlAxQVCiKgmJjh2+faaYuk5BKTnw6KlxalkIUkgyHliM6vYGuH8TwR0oG1795n1u7lpuulb27f3nd6FYSlkIUgbQwy5EjF2/xx41srv/fe2T/vg9UNriEdMTOt5ppnwNT28gAjxBFJIFZTqRmaej/cWxeWJ6OBZUNvr2mmsKyawMf3u3XVMJSiGKQwCwHdHoDreft4PrX75J95qe8sOw9Fadaj5v2mT+gmcw6JEQxSR9mGafR6piz9Sh//O+vsFTb4Ns7yhSWvRv78PvszhKWQpiB/BaVYdczsgmbuxvNn8fIPnsA1Db49X4Lx5rNAXAC/tO/udzqKISZSGCWUUlpmTwRvQcAh8dC8O72GmoHVxxrhpn2+WBQUwlLIcxIArOM0ekNxF9K5dmP9mLQZKF29gDApWG7fPs52qpoU9fPAhUKUX5ZVfNj79699OjRg4CAABRFYcuWLZYuyaro9Aa6fbCbZxfv5dqWaJLWTUGXefOe/ZYOasrR6bIOjxDmZlWBmZWVRePGjVm0aJGlS7FKx6+k8XtSBte2zCPn3EH06Snobl7Jt8/ht9rRuVGAhKUQJcCqfqu6dOlCly5dLF2GVcrM0dJrYQzXNs8j5/zhvNsdn52BQ9UQAMa1qcmr7WtJUApRgsr0b1dubi65ubmmx+np6RaspuRk5mhpOP1bUjbPRXP+CIqNPb7PTsexWhMAHvO0Z1ynOjLAI0QJK9O/YdHR0bi7u5u+goKCLF2S2ej0Bs5fy0Sj1dFm7rekbLorLPvOMIXlG0/V4MdJT0lYClEKynQLc+rUqUycONH0OD09vVyEpk5voNfi/ZxISMfP2ZaU1Cx0qQkotvb49Z1pOg33dVQY3b6uhKUQpaRMB6a9vT329vaWLsPszqVkciIhr3shJes2ahdP/AdEo0tPxiGwIQAKsPvNDhKWQpQi+W2zQpduZGO4nUvOxXjTNhs3H1NY+jgqHJ8pa/AIUdqsqoWZmZnJuXPnTI8vXLhAfHw8Xl5eVK1a1YKVlZ7ULA1jV+3n2sa30Vw6ds/6O9W8nNg+4UkZDRfCAqzqt+7w4cO0a/f3HSt3+ieHDh3K6tWrLVRV6UlKy6TF7B1c2zgbzZ/HUOwcUbt65dtn2dAwCUshLMSqfvPatm1LOVpiqMA0Wh07TyUR+ekvpGycTe6lvLD06zcLh8AGpv1CAt2p7uNswUqFqNisKjArIo1WR5M5O8nOyiZl4yxyLx3/Kyxn4xBY37TfqmFhPFnbVwZ5hLAgCUwL+/n8DXJyNKT871/kXj6BYueIf//Z2Ff5OyyDA9wkLIWwAhKYFpSapWH9L+dBbYud72Nok8//FZb1AJjepT6t6/lS3cdZwlIIKyCBWcp0egPnUjI5lZDKxK+OA3krOnp2eBnX5r2x9agEQE1PW4ZGVJOgFMKKSGCWorvv4DHkZpP+y0bcWw1AUduiKIopLDeMfJyw6l4SlkJYGQnMUnTpZrYpLFO+nEFuwu/oMq7j0+010z7x09vLmuFCWClpwpSiql5OeKs1JH85ndyE31E5uODarIfp+ff7hkhYCmHFpIVZCnR6A78nprPlp9OcWPUW2sTTqBxc8Hvubewr1QLATg3dQipbuFIhxMNIYJYwnd5AjwV7OflnMslfzkCbeAaVgyv+A97Gzr8mAO/2acQzTWSWdCGsnfyGliCd3sD2E0mcSsrk2pZ5d4XlXOz8a5j2C6vuJWEpRBkgv6UlRKc30HXBXs4kZ6EoCu6tnuf2rST8np2Gnd/fYVnX35WqXk4WrFQIUVASmCXk5/MpnEnOMj12CAqmyktLUdS2pm11/V3YGtlSLh8SooyQwDSzvAGeNJ5fFMP1r9/Fs91w7PyqA5jCMnZya7Q6hapeThKWQpQhEphmpNMb6Ll4P8f+uELyhmncTjnP9a//Q+URi1CUvGCs7W1PJTe51VGIskh+a83obHJGvrBUOXng03OyKSwBzt7I5dLNbAtWKYQoKglMM9DpDZy4cosZ62NJ3vBWXlg6e+A/cB52vo8BUM0zb+2hkEB3GeQRooySU/Ji0ukNdF+wl5MXEvLC8tpFVM4eVBoQja1P3gqWtbzs+GZCWxLSNNJvKUQZJoFZDDq9gW9+vcLvyVmk7v2U29cuonb2xH/gPGy988Jy1QvNeLKuHzZqFTV8XSxcsRCiOCQwi0ij1dFj0X7OpmQC4NluBIacDDxaDzGF5drhzWlVx8+SZQohzEgCswh0egM9FsVy+uoNVLZ5k2Wo7J3w7R1l2qeahx2P1/SxVIlCiBIgnWmFoNMbOH8tk5NXU/n9whWS1kwk7ecv79mvqoc92ye2k75KIcoZaWEWkE5voM9HP3Hsahr6zFskb4ji9o3LZMR9i0vTrqgd8vonFz3fmKcbBkhYClEOSWAW0IXrWRy7moYu8ybJ66PQ3byC2tUH/4HzTGEZ88aTVPN2s3ClQoiSIoH5CDq9gQvXs4j8/NA/wtI3bzTcM28Oy69eelzCUohyTgLzIXR6A70X7+d4Qjq6jBskb4hCd/Mqajdf/AdG/71gmZcdTR/zsnC1QoiSJoH5EBeuZ3E8IR0AzcX4+4bl3ddZCiHKNwnMh9DpDab/d2nUHqNBj8NjIaaw3PFaK+r6e1ioOiFEaZPAfACNVseXPxxFn5OB2tEVANfGnfLtY6uSb58QFYn8xt9HZo6WBq+vI3F9FCp7J/wGzDWNhN8hk2gIUfFIYP7lzmh47m0dgz/4PxLXT0WXmoSNuz9GbQ78FZhbIp/AzcFeJtEQogKSwOQfo+FpKSSvn4ouLRkbj0r4D4zGxs0XgB9ej6CWr7uFqxVCWIoEJnDpZvZfYZlM0voo9GnJ2HhU/issfVABcdPb4+HsYOlShRAWVOEDU6c3kHtbjzEtmaR1U9Gnp2DjWRn/AXlhWdndgZ0TnsTF0c7SpQohLKxCB+bd94fr9Dow6LDxDMB/4DxsXPNmGprTs6GEpRACqOCBeef+cABbryr4D4xGsXXAxtU7b5sCEbVkijYhRJ4KG5g6vYERC74m58IFHKs3BfJCE8DdwYY5PYPp1NAfB7sK+y0SQvxDhU2DH3/5lZ8WTkCfnYZ/v1k4PBZiem7dqCdoWEVGw4UQ+VXICwnPnTvH4D7d0Gdcx8bd37SkBEDDyq7UreRqweqEENaqwrUwz549S9u2bbmenIidTxB+A+Zh7+LJulEt8HR2oLqPs1yQLoS4rwoVmGfOnKFt27YkJibSoEEDvtv+PRezbXm8upf0VQohHqnCpMTly5dNYdmwYUN27dqFn58fj1m6MCFEmVFhzj0DAgJo164djRo1Yvfu3fj5yfK3QojCqTAtTLVazZo1a8jIyMDT09PS5QghyqBy3cI8deoUEydORK/XA2BjYyNhKYQosnLbwjx58iTt2rXj2rVreHp6Mn36dEuXJIQo48plC/PUqVOmsGzatCmvvvqqpUsSQpQD5bKF2a1bN27evEloaCg7d+7Ey0tWdBRCFF+5bGHevHmTZs2aSVgKIcyqXLUwjUYjACEhIWzcuBEbGxvS09MtXJUQwtrdyYk7GfIgivFRe5QhV65cISgo6NE7CiHEfVy+fJnAwMAHPl+uAtNgMJCQkICrqyuKoli6HLNLT08nKCiIy5cv4+bmZulySkRFeI9QMd5nWXqPRqORjIwMAgICUKke3FNZrk7JVSrVQ/86lBdubm5W/wEsrorwHqFivM+y8h7d3R89pWO5HPQRQoiSIIEphBAFJIFZhtjb2zNz5kzs7e0tXUqJqQjvESrG+yyP77FcDfoIIURJkhamEEIUkASmEEIUkASmEEIUkASmEEIUkARmGbB371569OhBQEAAiqKwZcsWS5dkdtHR0TRv3hxXV1f8/Pzo1asXp0+ftnRZZrVkyRJCQkJMF3KHh4ezbds2S5dVoqKjo1EUhQkTJli6FLOQwCwDsrKyaNy4MYsWLbJ0KSVmz549REZGcuDAAXbu3IlOp6NTp05kZWVZujSzCQwM5J133uHw4cMcPnyYp556ip49e3Ly5ElLl1YiDh06xCeffEJISIilSzEfoyhTAOPmzZstXUaJS0lJMQLGPXv2WLqUEuXp6Wlcvny5pcswu4yMDGPt2rWNO3fuNLZp08Y4fvx4S5dkFtLCFFYpLS0NoNzOZ6rX69mwYQNZWVmEh4dbuhyzi4yMpFu3bnTo0MHSpZhVuZp8Q5QPRqORiRMnEhERQXBwsKXLMavjx48THh6ORqPBxcWFzZs306BBA0uXZVYbNmwgLi6OQ4cOWboUs5PAFFZnzJgxHDt2jNjYWEuXYnZ169YlPj6e1NRUNm7cyNChQ9mzZ0+5Cc3Lly8zfvx4vv/+exwcHCxdjtnJrZFljKIobN68mV69elm6lBIxduxYtmzZwt69e6levbqlyylxHTp0oGbNmixdutTSpZjFli1b6N27N2q12rRNr9ejKAoqlYrc3Nx8z5U10sIUVsFoNDJ27Fg2b95MTExMhQhLyHvfubm5li7DbNq3b8/x48fzbXvxxRepV68ekydPLtNhCRKYZUJmZibnzp0zPb5w4QLx8fF4eXlRtWpVC1ZmPpGRkaxbt46tW7fi6upKUlISkDepq6Ojo4WrM4+oqCi6dOlCUFAQGRkZbNiwgZiYGLZv327p0szG1dX1nn5nZ2dnvL29y0V/tARmGXD48GHatWtnejxx4kQAhg4dyurVqy1UlXktWbIEgLZt2+bbvmrVKoYNG1b6BZWA5ORkhgwZQmJiIu7u7oSEhLB9+3Y6duxo6dJEAUkfphBCFJBchymEEAUkgSmEEAUkgSmEEAUkgSmEEAUkgSmEEAUkgSmEEAUkgSmEEAUkgSmEEAUkgSmEEAUkgSlEEdy4cQM/Pz8uXrxo6VIA6Nu3L/Pnz7d0GeWeBKbIp3Xr1iiKcs/XoEGDyswxhg0bhqIovPzyy/c89+qrr6Ioiun+9GHDhhVpqrzo6Gh69OhBtWrVilfsQ/To0eOBM5b//PPPKIpCXFwcADNmzGDu3Lmkp6eXWD1CAlPcxWg0Eh8fz3vvvUdiYmK+L3PN11gaxwAICgpiw4YN5OTkmLZpNBrWr19f7BmecnJyWLFiBSNHjixumQ81YsQIdu3axZ9//nnPcytXrqRJkyaEhoYCEBISQrVq1Vi7dm2J1lTRSWAKk7Nnz5KRkUHr1q2pVKlSvi8XF5cycwyA0NBQqlatyqZNm0zbNm3aRFBQEE2bNi3Wa2/btg0bG5t8a/G0bduWsWPHMmHCBDw9PfH39+eTTz4hKyuLF198EVdXV2rWrHnPsrpGo5F3332XGjVq4OjoSOPGjfnf//4HQPfu3fHz87tnRqrs7Gy++OILRowYkW/7M888w/r164v13sTDSWAKkyNHjmBjY1Oiy6KWxjHuePHFF1m1apXp8cqVKxk+fHixX3fv3r2EhYXds33NmjX4+Phw8OBBxo4dyyuvvEK/fv1o2bIlcXFxdO7cmSFDhpCdnW36N9OmTWPVqlUsWbKEkydP8tprrzF48GD27NmDjY0NL7zwAqtXr+buScW++uortFrtPV0YLVq04ODBg+VqQmJrI4EpTOLi4tDr9Xh7e+Pi4mL6GjVqVKke45tvvqFu3brUrl2b5cuXF/lYQ4YMITY2losXL/Lnn3+yf/9+Bg8eXOz3cPHiRQICAu7Z3rhxY6ZNm0bt2rWZOnUqjo6O+Pj4MGrUKGrXrs2MGTO4ceMGx44dA/LWm58/fz4rV66kc+fO1KhRg2HDhjF48GBT98Tw4cO5ePEiMTExpuOsXLmSPn364Onpme/4VapUITc31zT5sjA/mUBYmBw5coR+/foxd+7cfNv/+YsJ8K9//YtZs2Y99PUOHTp0T0vsUcfQ6XRMnDiR3bt34+bmRmhoKH369CnScrs+Pj5069aNNWvWYDQa6datGz4+PoV+nX/Kycm57wJfd7ea1Wo13t7eNGrUyLTN398fgJSUFABOnTqFRqO5ZwJhrVZr6jaoV68eLVu2ZOXKlbRr144//viDffv28f33399z/Dsz09/dghXmJYEpTI4ePcrs2bOpVavWI/cdM2YMAwYMeOg+9xtBftQxDh48SMOGDalSpQoAXbt2ZceOHQwcOPDRb+A+hg8fzpgxYwBYvHhxkV7jn3x8fLh169Y9221tbfM9VhQl3zZFUQAwGAz5/vvtt9+a3u8d9vb2pv8fMWIEY8aMYfHixaxatYrHHnuM9u3b33P8mzdvAuDr61uUtyUKQAJTAHD+/HlSU1MLPCDi4+NT6NZaQY6RkJCQLzwCAwO5evVqoY5zt6effhqtVgtA586di/w6d2vatCmff/55sV+nQYMG2Nvbc+nSJdq0afPA/fr378/48eNZt24da9asYdSoUabwvduJEycIDAw0Syta3J8EpgDyTpUh77Txn31gfn5+qFTF7+4uyDHut2LK/cKhoNRqNb/99pvp/+8nLS2N+Pj4fNsetsBc586dmTp1Krdu3bpvd0VBubq6MmnSJF577TUMBgMRERGkp6fz008/4eLiwtChQwFwcXHhueeeIyoqirS0tAeucbRv3z46depU5HrEo0lgCgDTBdB16tTJt93W1paMjIx8p4gleYwqVarka1FeuXKFxx9/3PR49erVvPjii/cN1gdxc3N76PMxMTH3tHoftsBco0aNCAsL48svv2T06NEFruN+5syZg5+fH9HR0Zw/fx4PDw9CQ0OJiorKt9+IESNYsWIFnTp1um+QazQaNm/ezI4dO4pVj3g4WQRNWBWdTkf9+vWJiYkxDfocOHAAb29vIG+wKSYmJt+osSV89913TJo0iRMnTpil9V1cixcvZuvWrfcdDBLmIy1MYVVsbGx4//33adeuHQaDgTfffNMUlgA7duzgww8/tGCFebp27crZs2e5evUqQUFBli4HW1tbFi5caOkyyj1pYQohRAFZ/lxCCCHKCAlMIYQoIAlMIYQoIAlMIYQoIAlMIYQoIAlMIYQoIAlMIYQoIAlMIYQoIAlMIYQoIAlMIYQoIAlMIYQoIAlMIYQooP8HDqiRiwMUbpEAAAAASUVORK5CYII=", "text/plain": [ "
" ] @@ -276,8 +275,18 @@ "metadata": { "tags": [] }, - "outputs": [], + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "CPU times: user 14min 25s, sys: 43.9 s, total: 15min 9s\n", + "Wall time: 1min 23s\n" + ] + } + ], "source": [ + "%%time\n", "model = MBTREnergyModel(reference=data[0], calc=mbtr)\n", "with warnings.catch_warnings():\n", " warnings.simplefilter(\"ignore\")\n", @@ -293,7 +302,7 @@ { "data": { "text/plain": [ - "{'alpha': 1e-10, 'gamma': 0.03534981105030106}" + "{'gamma': 9.284145445194744e-05}" ] }, "execution_count": 11, @@ -349,7 +358,7 @@ "name": "stdout", "output_type": "stream", "text": [ - "Maximum force: 1.46e-03 eV/Angstrom\n" + "Maximum force: 1.51e-02 eV/Angstrom\n" ] } ], @@ -367,7 +376,7 @@ "outputs": [ { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAC+CAYAAADa6ROSAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAAOTUlEQVR4nO3cfWxTdfvH8U+3seNW1ipMnbBujTzMeIcnswQhMUBQWcaDghgwGAYSE8CIMVESE8MIiovo2KKJaCIDIUZFBaIYiDwIIooPIEGjQfQPHgQEDFkHxo5u398f/rjuuzJj2Vq2te9X0pBzenrO9S3X2afnfLf6nHNOAABIyursAgAAXQehAAAwhAIAwBAKAABDKAAADKEAADCEAgDAEAoAAEMoAAAMoZDGfD6fNm7c2NllIIOtXr1a1157bWeXcUW6Y83JRCgkyeeff67s7GxVVFRc0evC4bDq6+tTUxRwhWbNmiWfz3fZI5G+bquXp02bpp9++ilF1f5Xpv8gT6aczi4gXTQ0NOjRRx/V66+/rqNHj6qkpKSzSwLapaKiQqtWrYpb53leu/aVl5envLy8ZJSFq4QrhSS4cOGC1q1bp3nz5mnChAlavXp13PMffPCBysvLdc0116iwsFBTpkyRJI0ePVpHjhzR448/bp/IJGnx4sUaOnRo3D7q6+sVDodt+euvv9Zdd92lwsJCBYNBjRo1Svv370/lMJEhPM9TUVFR3OO6666T9FdvlpSUyPM89enTRwsWLJD0z73890/wl3q7oaFBJSUl6tmzp+bNm6eWlhYtW7ZMRUVFuuGGG7R06dK4mpYvX65BgwbJ7/crFApp/vz5On/+vCRp586dmj17thobG+3YixcvliQ1Nzdr4cKF6tu3r/x+v4YPH66dO3fG7Xv16tUqKSlRfn6+Jk+erN9//z0F72r3QSgkwTvvvKOysjKVlZXpwQcf1KpVq3Tpy2c/+ugjTZkyRePHj9e3336r7du3q7y8XJK0fv16FRcXa8mSJTp58qROnjyZ8DGbmppUVVWl3bt3a+/evRowYIAqKyvV1NSUkjEC7733nurq6vTaa6/p8OHD2rhxowYNGiTpynr5l19+0ebNm7Vlyxa99dZbamho0Pjx43X8+HHt2rVLzz//vJ5++mnt3bvXXpOVlaWXXnpJ33//vd544w3t2LFDCxculCSNHDlS9fX1CgQCduwnnnhCkjR79mzt2bNHb7/9tg4ePKj7779fFRUVOnz4sCTpyy+/1EMPPaT58+frwIEDGjNmjJ599tlUvYXdg0OHjRw50tXX1zvnnLt48aIrLCx0W7dudc45N2LECDdjxox/fG1paamrq6uLW1ddXe2GDBkSt66urs6Vlpb+435isZgrKChwH374oa2T5DZs2HBFY0Fmq6qqctnZ2c7v98c9lixZ4mpra93AgQNdc3Nzm69tq5dXrVrlgsGgLVdXV7v8/HwXiURs3bhx41w4HHYtLS22rqyszNXU1PxjnevWrXO9e/f+x+M459zPP//sfD6f+/XXX+PWjx071j311FPOOeceeOABV1FREff8tGnTLttXJmFOoYMOHTqkr776SuvXr5ck5eTkaNq0aWpoaNCdd96pAwcO6OGHH076cU+fPq1FixZpx44d+u2339TS0qI//vhDR48eTfqxkFnGjBmjFStWxK3r1auXLly4oPr6et18882qqKhQZWWlJk6cqJycK/sxEg6HVVBQYMs33nijsrOzlZWVFbfu9OnTtvzJJ5/oueee0w8//KBIJKJYLKY///xTFy5ckN/vb/M4+/fvl3NOAwcOjFsfjUbVu3dvSdKPP/6oyZMnxz0/YsQIbdmy5YrGlE4IhQ5auXKlYrGY+vbta+ucc+rRo4fOnTvXrkm2rKwsu/10ycWLF+OWZ82apTNnzqi+vl6lpaXyPE8jRoxQc3Nz+wYC/D+/36/+/ftftr5Xr146dOiQtm7dqm3btmn+/Pl64YUXtGvXLvXo0SPh/f99W5/P1+a61tZWSdKRI0dUWVmpuXPn6plnnlGvXr302Wefac6cOZedF/+rtbVV2dnZ2rdvn7Kzs+Oe69mzpyRddp6BUOiQWCymNWvWqLa2VnfffXfcc/fdd5/efPNNDR48WNu3b9fs2bPb3Edubq5aWlri1l1//fU6deqUnHM2YXfgwIG4bXbv3q1XXnlFlZWVkqRjx47p7NmzSRoZ0La8vDxNmjRJkyZN0iOPPKJbbrlF3333nW677bY2ezkZvvnmG8ViMdXW1trVxLp16+K2aevYw4YNU0tLi06fPq077rijzX3feuutcXMXki5bzjSEQgds2rRJ586d05w5cxQMBuOemzp1qlauXKm6ujqNHTtW/fr10/Tp0xWLxbR582abJAuHw/r00081ffp0eZ6nwsJCjR49WmfOnNGyZcs0depUbdmyRZs3b1YgELD99+/fX2vXrlV5ebkikYiefPJJfvUPSRGNRnXq1Km4dTk5Odq0aZNaWlo0fPhw5efna+3atcrLy1Npaamktns5Gfr166dYLKaXX35ZEydO1J49e/Tqq6/GbRMOh3X+/Hlt375dQ4YMUX5+vgYOHKgZM2Zo5syZqq2t1bBhw3T27Fnt2LFDgwYNUmVlpRYsWKCRI0dq2bJluvfee/Xxxx9n9K0jSUw0d8SECRNcZWVlm8/t27fPSXL79u1z77//vhs6dKjLzc11hYWFbsqUKbbdF1984QYPHuw8z3P/+9+xYsUKFwqFnN/vdzNnznRLly6Nm2jev3+/Ky8vd57nuQEDBrh33333sok+MdGMK1RVVeUkXfYoKytzGzZscMOHD3eBQMD5/X53++23u23bttlr2+rltiaa//5LFFVVVe6ee+6JWzdq1Cj32GOP2fLy5cvdTTfd5PLy8ty4cePcmjVrnCR37tw522bu3Lmud+/eTpKrrq52zjnX3NzsFi1a5MLhsOvRo4crKipykydPdgcPHrTXrVy50hUXF7u8vDw3ceJE9+KLL2b0RLPPOW6qAQD+wt8pAAAMoQAAMN0qFKLRqBYvXqxoNNrZpaRMJoxRypxxJlOmvGeZMM6uPMZuNacQiUQUDAbV2NgY95s46SQTxihlzjiTKVPes0wYZ1ceY7e6UgAApBahAAAw7f7jtdbWVp04cUIFBQX2V7epFolE4v5NR5kwRunqj9M5p6amJvXp0yfuO3bai/5PnUwYZ2eMMdFzoN1zCsePH1coFGp3gUBnOHbsmIqLizu8H/of3dW/nQPtvlK49C2HPxw6HPeNh+ko19fa2SWknPOl953EpqYmDRgwIGm9mlH9r1hnl5ByLiv9v/En0XOg3e/EpUvmgoKCLjd7nmyEQvpI1q2ejOp/QiGt/Ns5kBk/CQAACSEUAACGUAAAGEIBAGAIBQCAIRQAAIZQAAAYQgEAYAgFAIAhFAAAhlAAABhCAQBgCAUAgCEUAACGUAAAGEIBAGAIBQCAIRQAAIZQAAAYQgEAYAgFAIDJSXTDaDSqaDRqy5FIJCUFAV0R/Y9MkfCVQk1NjYLBoD1CoVAq6wK6FPofmcLnnHOJbNjWJ6VQKKRjJ04pEAikrMCuINfX2tklpJzzpfedxEgkoqKiIjU2NrarXzO6/xXr7BJSzmUlfNOk20r0HEj4nfA8T57nJaU4oLuh/5Ep0vvjIQDgihAKAABDKAAADKEAADCEAgDAEAoAAEMoAAAMoQAAMIQCAMAQCgAAQygAAAyhAAAwhAIAwBAKAABDKAAADKEAADCEAgDAEAoAAEMoAAAMoQAAMIQCAMDkdHQHPbL+eqS3tB8g2ikz+r/DPybQjaR9OwMAEkcoAAAMoQAAMIQCAMAQCgAAQygAAAyhAAAwhAIAwBAKAABDKAAADKEAADCEAgDAEAoAAEMoAAAMoQAAMIQCAMAQCgAAQygAAAyhAAAwhAIAwBAKAABDKAAATE6iG0ajUUWjUVuORCIpKQjoiuh/ZIqErxRqamoUDAbtEQqFUlkX0KXQ/8gUPuecS2TDtj4phUIhnTp1SoFAIGUFdgW+zi7gKkioCbqxSCSioqIiNTY2tqtf6f/0lu79LyV+DiR8+8jzPHmel5TigO6G/kemYKIZAGAIBQCAIRQAAIZQAAAYQgEAYAgFAIAhFAAAhlAAABhCAQBgCAUAgCEUAACGUAAAGEIBAGAIBQCAIRQAAIZQAAAYQgEAYAgFAIAhFAAAhlAAABhCAQBgcjq6A19Ls3wtzcmopcuKZeV2dgkpl+3r7Aq6J/o/PdD//8WVAgDAEAoAAEMoAAAMoQAAMIQCAMAQCgAAQygAAAyhAAAwhAIAwBAKAABDKAAADKEAADCEAgDAEAoAAEMoAAAMoQAAMIQCAMAQCgAAQygAAAyhAAAwhAIAwBAKAABDKAAATE6iG0ajUUWjUVuORCIpKQjoiuh/ZIqErxRqamoUDAbtEQqFUlkX0KXQ/8gUPuecS2TDtj4phUIh/fbrUQUCgZQV2BXEsnI7u4SUy/Z1dgWpFYlEVFRUpMbGxnb1K/2f3tK9/6XEz4GEbx95nifP85JSHNDd0P/IFEw0AwAMoQAAMIQCAMAQCgAAQygAAAyhAAAwhAIAwBAKAABDKAAADKEAADCEAgDAEAoAAEMoAAAMoQAAMIQCAMAQCgAAQygAAAyhAAAwhAIAwBAKAACT094XOuckSU1NTUkrpquKZeV2dgkpl+3r7ApS61KfXurbjqL/00u697+U+DnQ7lC4dID+t/ynvbsArrqmpiYFg8Gk7Eei/9H9/Ns54HPt/OjU2tqqEydOqKCgQD7f1YnZSCSiUCikY8eOKRAIXJVjXm2ZMEbp6o/TOaempib16dNHWVkdv2tK/6dOJoyzM8aY6DnQ7iuFrKwsFRcXt/flHRIIBNK2WS7JhDFKV3ecybhCuIT+T71MGOfVHmMi5wATzQAAQygAAEy3CgXP81RdXS3P8zq7lJTJhDFKmTPOZMqU9ywTxtmVx9juiWYAQPrpVlcKAIDUIhQAAIZQAAAYQgEAYAgFAIAhFAAAhlAAABhCAQBg/g+IUzrL4FJNvgAAAABJRU5ErkJggg==", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAC+CAYAAADa6ROSAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAAS3ElEQVR4nO3df4xcVd3H8c/M7O7t7uzMlHaR7bKz3dAfa0zKD58mtU0MEFQ2C0VbMMVgKIWQQI0YEyUxMbRBcWO1dKOJaGJ/CDFqVUDFtBFaQeSHSEsfIPCUyvOkLf1BW+yzd/trdmfmPH/02aNDd9O7M3M7O+e+X8lkM3dm7nzP7LnzveeeM+fEjDFGAABIitc6AADA5EFSAABYJAUAgEVSAABYJAUAgEVSAABYJAUAgEVSAABYJAUAgEVScFgsFtOTTz5Z6zAQYZs2bdLUqVNrHcaE1GPM1URSqJIXX3xRiURCvb29E3pdd3e3BgYGwgkKmKA77rhDsVjsnFuQej1WXV62bJneeeedkKL9l6h/kVdTQ60DcMWGDRv05S9/WT/96U+1b98+dXV11TokoCy9vb3auHFjyTbP88raV3Nzs5qbm6sRFi4QWgpVcPLkSW3evFn33nuvbrzxRm3atKnk8d///veaP3++pkyZora2Ni1dulSSdM0112jv3r366le/as/IJGn16tW68sorS/YxMDCg7u5ue//vf/+7Pv3pT6utrU2ZTEZXX321du7cGWYxERGe56m9vb3kdtFFF0k6Wze7urrkeZ46Ojp03333SRq/Ln/4DH60bm/YsEFdXV1qbW3Vvffeq0KhoDVr1qi9vV0f+chH9NBDD5XE9PDDD2vevHlKJpPKZrNauXKlTpw4IUl69tlntWLFCg0ODtr3Xr16tSRpeHhY999/vy699FIlk0ktWLBAzz77bMm+N23apK6uLrW0tGjJkiX64IMPQvhU6wdJoQp+9atfqaenRz09PfriF7+ojRs3anTy2T/+8Y9aunSpbrjhBr322mvatm2b5s+fL0l6/PHH1dnZqQcffFCHDh3SoUOHAr/n0NCQli9frueff14vv/yy5syZo76+Pg0NDYVSRuA3v/mN1q1bp5/85Cfas2ePnnzySc2bN0/SxOryu+++qy1btmjr1q36xS9+oQ0bNuiGG27Qe++9p+eee07f/e539c1vflMvv/yyfU08HtcPfvADvfnmm/rZz36m7du36/7775ckLVq0SAMDA0qn0/a9v/a1r0mSVqxYoRdeeEG//OUv9frrr+vzn/+8ent7tWfPHknS3/72N915551auXKldu3apWuvvVbf/va3w/oI64NBxRYtWmQGBgaMMcaMjIyYtrY28/TTTxtjjFm4cKG57bbbxn3tzJkzzbp160q2rVq1ylxxxRUl29atW2dmzpw57n7y+bxJpVLmD3/4g90myTzxxBMTKguibfny5SaRSJhkMllye/DBB83atWvN3LlzzfDw8JivHasub9y40WQyGXt/1apVpqWlxfi+b7ddf/31pru72xQKBbutp6fH9Pf3jxvn5s2bzfTp08d9H2OM+cc//mFisZg5cOBAyfbrrrvOfOMb3zDGGPOFL3zB9Pb2ljy+bNmyc/YVJfQpVGj37t165ZVX9Pjjj0uSGhoatGzZMm3YsEGf+tSntGvXLt19991Vf98jR47ogQce0Pbt2/X++++rUCjo1KlT2rdvX9XfC9Fy7bXX6pFHHinZNm3aNJ08eVIDAwO67LLL1Nvbq76+Pi1evFgNDRP7Gunu7lYqlbL3L7nkEiUSCcXj8ZJtR44csff//Oc/6zvf+Y7eeust+b6vfD6vM2fO6OTJk0omk2O+z86dO2WM0dy5c0u253I5TZ8+XZL09ttva8mSJSWPL1y4UFu3bp1QmVxCUqjQ+vXrlc/ndemll9ptxhg1Njbq+PHjZXWyxeNxe/lp1MjISMn9O+64Q0ePHtXAwIBmzpwpz/O0cOFCDQ8Pl1cQ4P8lk0nNnj37nO3Tpk3T7t279fTTT+uZZ57RypUr9b3vfU/PPfecGhsbA+//w8+NxWJjbisWi5KkvXv3qq+vT/fcc4++9a1vadq0afrrX/+qu+6665zj4t8Vi0UlEgnt2LFDiUSi5LHW1lZJOuc4A0mhIvl8Xo8++qjWrl2rz3zmMyWP3Xzzzfr5z3+uyy+/XNu2bdOKFSvG3EdTU5MKhULJtosvvliHDx+WMcZ22O3atavkOc8//7x+9KMfqa+vT5K0f/9+HTt2rEolA8bW3Nysm266STfddJO+9KUv6aMf/ajeeOMNffzjHx+zLlfDq6++qnw+r7Vr19rWxObNm0ueM9Z7X3XVVSoUCjpy5Ig++clPjrnvj33sYyV9F5LOuR81JIUKPPXUUzp+/LjuuusuZTKZksduueUWrV+/XuvWrdN1112nWbNm6dZbb1U+n9eWLVtsJ1l3d7f+8pe/6NZbb5XneWpra9M111yjo0ePas2aNbrlllu0detWbdmyRel02u5/9uzZeuyxxzR//nz5vq+vf/3rDP1DVeRyOR0+fLhkW0NDg5566ikVCgUtWLBALS0teuyxx9Tc3KyZM2dKGrsuV8OsWbOUz+f1wx/+UIsXL9YLL7ygH//4xyXP6e7u1okTJ7Rt2zZdccUVamlp0dy5c3Xbbbfp9ttv19q1a3XVVVfp2LFj2r59u+bNm6e+vj7dd999WrRokdasWaPPfe5z+tOf/hTpS0eS6GiuxI033mj6+vrGfGzHjh1GktmxY4f57W9/a6688krT1NRk2trazNKlS+3zXnrpJXP55Zcbz/PMv/87HnnkEZPNZk0ymTS33367eeihh0o6mnfu3Gnmz59vPM8zc+bMMb/+9a/P6egTHc2YoOXLlxtJ59x6enrME088YRYsWGDS6bRJJpPmE5/4hHnmmWfsa8eqy2N1NH94EMXy5cvNZz/72ZJtV199tfnKV75i7z/88MNmxowZprm52Vx//fXm0UcfNZLM8ePH7XPuueceM336dCPJrFq1yhhjzPDwsHnggQdMd3e3aWxsNO3t7WbJkiXm9ddft69bv3696ezsNM3NzWbx4sXm+9//fqQ7mmPGcFENAHAWv1MAAFgkBQCAVVdJIZfLafXq1crlcrUOJTRRKKMUnXJWU1Q+syiUczKXsa76FHzfVyaT0eDgYMlIHJdEoYxSdMpZTVH5zKJQzslcxrpqKQAAwkVSAABYZf94rVgs6uDBg0qlUvZXt2Hzfb/kr4uiUEbpwpfTGKOhoSF1dHSUzLFTLup/eKJQzlqUMegxUHafwnvvvadsNlt2gEAt7N+/X52dnRXvh/qPenW+Y6DslsLoLIf/+dbukhkPXZQv1jqC8GU8t68kDg0NafacOVWrq6P72b7zLbW2ul7/62YsStkum9pU6xBCNzQ0pDkBjoHASSGXy5UMnxpdzCWVSik1yXrPqy0KSSHteFIYVe6lnvHqf2trSq0p1+u/+0khnXY/KYw63zEQ+Jugv79fmUzG3mg6I0qo/4iKwH0KHz5T8n1f2WxW/73/IC0FB0x1vKXg+74uaW8ve1z4ePX/lXf201JwwOyL3G8p+L6v9gDHQODLR57nyfO8qgQH1BvqP6Ki4vUUkvG8WuP5asQyaeXiwVeVqleunwuGVb6udKPSabfrx/un3D6+JSlm3L8cELSMbl8zAABMCEkBAGCRFAAAFkkBAGCRFAAAVsWjj2Km6HzPfUPiwkx4hvpzaqSoxIjb9X/GG7+rdQihKy64udYhhM7EgrUBaCkAACySAgDAIikAACySAgDAIikAAKyKRx9FQaI4XOsQQmcS7s8SGYYT+aJijo8+mtY5u9YhhC53/qdEBi0FAIBFUgAAWCQFAIBFUgAAWCQFAIBV8egjE2+Qibs9iCkfc7t8kpSodQB1akbhn0oXRmodRqjyF3XVOoTQuT5/m8TKawCAMpAUAAAWSQEAYJEUAAAWSQEAYFU8rCZnEsoZt8eueMrXOoTQmQiMsApDsXmqii3pWocRqhNF9+tGa60DmERoKQAALJICAMAK3C7M5XLK5f41wazv+6EEBExG1H9EReCWQn9/vzKZjL1ls9kw4wImFeo/oiJmjDFBnjjWmVI2m9W+A4eUTrvd0ebFCrUOIXSuT1Xi+77a29s1ODhYVn0dr/6/f2Cf8/U/Eh3N7hdRvu/rkhkd5z0GAn8UnufJ87xzt2tEntye+4WuF4xX/+NvbFM82VKDiC6c5H8srnUIoSsEOjWub4WA32N82wEALJICAMAiKQAALJICAMAiKQAALJICAMCqeHSuX2iQKbg9yLe10f3c6X4Jw5HomKVEyu3p1P6Zc3+pyozn/hEQjwV8XrhhAADqCUkBAGCRFAAAFkkBAGCRFAAAVsXDhqYWfKXzjs8mNez+cpyF1otrHUJdKqTaVUinah1GqKaNfFDrEEJX8Kj/o2gpAAAskgIAwCIpAAAskgIAwCIpAACsikcf/W8irUKD22vUTkkEnDSkjjXVOoA6tSc3Ra1nmmsdRqjaW90eXSVJbs9eNTG0FAAAFkkBAGCRFAAAFkkBAGCRFAAAVsWjj9KJvNIJx+cGSrg/Nsfx2atCM1dHldbpWocRqmLM/XmBjHH/GI+ZYCvo0VIAAFiBWwq5XE65XM7e930/lICAyYj6j6gI3FLo7+9XJpOxt2w2G2ZcwKRC/UdUxIwxgS4nj3WmlM1m9f6BfUqn3f5FM30K9c/3fbW3t2twcLCs+jpe/T/2XzuUTrn9e9hiKgJ9ChE4xn3f1yUzOs57DAS+fOR5njzPq0pwQL2h/iMqKh59ZBJNzmdZ92c+cr+MYZVvMDlDxVa3W8pe3P3xKI2uHwCSTCzY/9H9/zYAIDCSAgDAIikAACySAgDAIikAAKyKRx/F5P7IleFgU4bUtSZOD8qSbIyrtdHtD+9MwfVfsUhNzn+LBf+edrs2AwAmhKQAALBICgAAi6QAALBICgAAq/K5j+T+DJuODy5BBWLFvGJFt1ceNCZR6xDCF3BVsrrGymsAgIkiKQAALJICAMAiKQAALJICAMCqePRR0Zy9uazx9PFahxC6YstFtQ6hLp0qJtRQdHt0TktspNYhhM7E3F49UmLlNQBAGUgKAACLpAAAsEgKAACLpAAAsJj7KIDhZvdH5sQc/yeGtXhYInb25rJTprHWIYTOc7z+S8FHidJSAABYgVsKuVxOuVzO3vd9P5SAgMmI+o+oCNxS6O/vVyaTsbdsNhtmXMCkQv1HVMSMMYGuNI11ppTNZnXg0GGl0+nQAsSF4fhlcfm+r44Z7RocHCyrvo5X//cdOOR8/Q+rP2Yy8VzvGFLwYyDw5SPP8+R5XlWCA+oN9R9RUfHooyhoMG6vrCVJJu52VYiHdCLoxY28uNun0oUIjEcJq35MJkHL6P5/GwAQGEkBAGCRFAAAFkkBAGCRFAAAFiuvBZCPuT0yR5Iaim6PsIqFVL6C4s6PzjmdL9Y6hNC1Nrr9P5SC/xbJ/U8CABAYSQEAYJEUAAAWSQEAYJEUAABWxcNq4jH35w2JwASKMo6PsAprbqejp/M60+D2yK1Uk/vnjo4PoJQUvIzu/7cBAIGRFAAAFkkBAGCRFAAAFkkBAGCRFAAAVsXj9BpHTqtxxO3hjMWmllqHELp4/kytQwhVWOWbET+tdMLt+j/SMP4i764Ia8LEySRoGWkpAAAskgIAwCIpAAAskgIAwCIpAACsyodNmOLZm8NcX25UkkbiU2odQqhy8eFwdhxvOHtzWBQmhCxU4atwsisEnPSSlgIAwAqcHnO5nHK5nL3v+34oAQGTEfUfURG4pdDf369MJmNv2Ww2zLiASYX6j6iIGWMCXTEf60wpm83qyN53lU6nQgtwMsg3tdY6hNAVHO838X1f2Y52DQ4OKp2e+C90x63/+/7H+fofhV/0R6Hf0Pd9dcw4/zEQ+PKR53nyPK8qwQH1hvqPqKi4yz3flHT+THo4AqcRUxwfYtIY0pCKk/EpSsSbw9n5JDEy7PboQklKR2DJ0aDLJrv/SQAAAiMpAAAskgIAwCIpAAAskgIAwKp49FHC5JUwbq9a1OJ4+SQpb9ye+yisAWReIibP8ZFbqVOHax1C6E43zKh1CKEbCTiIjJYCAMAiKQAALJICAMAiKQAALJICAMCqfLmhWPzsDXXN9RFkYZUvEYvGymSuazIhrcw3iQQtI9/mAACLpAAAsEgKAACLpAAAsMruaB5dxXNoaKhqwUxWsYL7nVAmXvmYg8lstJ4GXH32vKJU/xMRKGOh0FjrEEIX9Bgo+5tg9A1mz+0pdxfABTc0NKRMJlOV/UjS7DlzKt4XcCGd7xiImTJPnYrFog4ePKhUKqVY7MKMyRtdLH3//v1lLb5eD6JQRunCl9MYo6GhIXV0dCger/yqKfU/PFEoZy3KGPQYKLulEI/H1dnZWe7LK5JOp52tLKOiUEbpwpazGi2EUdT/8EWhnBe6jEGOATqaAQAWSQEAYNVVUvA8T6tWrZLnebUOJTRRKKMUnXJWU1Q+syiUczKXseyOZgCAe+qqpQAACBdJAQBgkRQAABZJAQBgkRQAABZJAQBgkRQAABZJAQBg/R9mSKmFESumkAAAAABJRU5ErkJggg==", "text/plain": [ "
" ] @@ -408,8 +417,8 @@ "name": "stdout", "output_type": "stream", "text": [ - "CPU times: user 204 ms, sys: 3 µs, total: 204 ms\n", - "Wall time: 203 ms\n" + "CPU times: user 4min 8s, sys: 3.23 s, total: 4min 12s\n", + "Wall time: 1min 24s\n" ] } ], @@ -437,9 +446,9 @@ { "data": { "text/plain": [ - "array([[ 6.69634496e+01, -6.52252610e-03, 1.55513525e-08],\n", - " [-6.52252610e-03, 4.33040064e+01, 1.54795588e-06],\n", - " [ 1.55513525e-08, 1.54795588e-06, 5.49449873e-03]])" + "array([[ 1.42193910e+01, 2.19101005e+01, -2.39586674e-04],\n", + " [ 2.19101005e+01, 7.86328494e+01, -3.46912831e-05],\n", + " [-2.39586674e-04, -3.46912831e-05, 2.52655314e+00]])" ] }, "execution_count": 17, @@ -462,9 +471,9 @@ { "data": { "text/plain": [ - "array([[ 6.69440158e+01, -6.20564720e-03, 0.00000000e+00],\n", - " [-6.20564720e-03, 4.27768994e+01, 0.00000000e+00],\n", - " [ 0.00000000e+00, 0.00000000e+00, 4.91597802e-03]])" + "array([[ 1.50773132e+01, 2.24398727e+01, 4.01418338e-03],\n", + " [ 2.24398727e+01, 8.10175928e+01, -5.86781535e-03],\n", + " [ 4.01418338e-03, -5.86781535e-03, 2.09022244e+00]])" ] }, "execution_count": 18, @@ -486,7 +495,7 @@ "outputs": [ { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXYAAADJCAYAAAA3tRlxAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAAS5klEQVR4nO3dbWwU5frH8Wtb2rVdti2lKgK1ldMCBSwghQDCISgqCogJCkcJUHxIxEjwMZEogvLCF1rFF0ZRsWrkrxgBE1CDRkHRUrBAQbRSnilCwykP2wKyben1f0E6x4UCQ2fazt58P8kGmJ29957pNb8O83T7VFUFAGCMmLbuAADAXQQ7ABiGYAcAwxDsAGAYgh0ADEOwA4BhCHYAMAzBDgCGIdgBwDAEOwBPyM/Pl3vuuaetu2EEgv0y5Ofni8/nO+81evToVvn+efPmSb9+/VrluxC9ioqKJDY2ttXq0i1vvvmmfPjhhy3+PVfCL5B2bd2BaDN69GgpLCyMmOb3+9uoN8D5PvjgA5k5c6a8//77sn//frn++utb9Pvq6uokLi7OcTvJycku9AYiIqKwbdq0aTp+/Pgm31u9erXGxcXpTz/9ZE177bXXtGPHjnrw4EFVVf3mm2/05ptv1uTkZE1NTdUxY8bozp07I9qpqKjQSZMmaYcOHTQxMVEHDBigxcXFWlhYqCIS8SosLGypRUWUOnHihAaDQf3zzz910qRJ+tJLL1nvrV69WkVEV65cqbm5uer3+3XQoEG6detWa57CwkJNTk7W5cuXa3Z2tvr9fh01apTu37/fmmfu3Lnat29fXbRokd5www3q8/m0oaFB9+3bp3fffbcGAgENBoN63333aWVlpaqqlpWVaUJCgi5evNhqZ+nSper3+63vP3f7GjFihD7++OM6a9YsTUlJ0WuuuUYXLlyoJ06c0Pz8fG3fvr1269ZNv/76a+sz9fX1+uCDD2pmZqZeddVV2r17d12wYEFE38/djlavXq2qqgcOHNCJEydqSkqKpqam6t1336179uxx5efS2gj2y3CxYFdVffbZZzUjI0OPHz+upaWl6vf7ddmyZdb7X3zxhS5dulTLy8t18+bNOm7cOL3xxhv1zJkzqqpaU1Oj3bp10+HDh+vatWt1x44dumTJEi0qKtJTp07p008/rb1799ZDhw7poUOH9NSpUy29yIgyixYt0ry8PFVVXbFihWZmZmpDQ4Oq/i/Yc3Jy9Ntvv9WtW7fq2LFjNTMzU2tra1X1bLDHxcVpXl6eFhUVaUlJiQ4aNEiHDh1qfcfcuXM1EAjoHXfcoZs2bdItW7ZoQ0OD9u/fX4cNG6YlJSVaXFysN910k44YMcL63FtvvaXJycm6d+9e/euvvzQ1NVXfeOMN6/2mgj0YDOr8+fO1vLxc58+frzExMXrnnXfqu+++q+Xl5Tpjxgzt2LGjnjx5UlVVa2tr9cUXX9QNGzbo7t279ZNPPtHExERdsmSJqp7dxiZOnKijR4+2tqNwOKwnT57U7OxsffDBB3Xr1q36xx9/6AMPPKA9evTQcDjcEj+qFkWwX4Zp06ZpbGysBgKBiNfLL7+sqqrhcFj79++vEydO1N69e+vDDz980fYOHz6sIqK//fabqqouXLhQg8GgHjlypMn5G/eUgAsZOnSotYdaV1enaWlp+t1336nq/4L9s88+s+Y/cuSIJiQkWMHX+D/D4uJia56ysjIVEV2/fr2qnq3DuLg4PXz4sDXPt99+q7GxsRF79r///ruKiG7YsMGaNmbMGB0+fLjeeuutetttt1m/dFSbDvZhw4ZZ/66vr9dAIKBTpkyxph06dEhFRNetW3fBdfLYY4/phAkTLvg9qmd/Ifbo0SOiP+FwWBMSEnTVqlUXbNurOMZ+mUaOHClvv/12xLTU1FQREYmPj5dPPvlEcnNzJSMjQxYsWBAx365du2TOnDlSXFwsVVVV0tDQICIi+/fvlz59+khpaan079/fag+4HNu3b5cNGzbIsmXLRESkXbt2MmnSJPnggw9k1KhR1nxDhgyx/p6amio9evSQsrIya1q7du0kLy/P+nfPnj0lJSVFysrKZNCgQSIikpGRIVdffbU1T1lZmaSnp0t6ero1rVevXtbnBg4cKCJnj/93795dYmJiZNu2beLz+S66TLm5udbfY2NjpWPHjnLjjTda06699loRETl8+LA17Z133pH3339f9u3bJ3///bfU1tZe8qKDjRs3ys6dOyUYDEZMP336tOzateuin/Uigv0yBQIBycrKuuD7RUVFIiJy9OhROXr0qAQCAeu9cePGSXp6urz33nvSuXNnaWhokD59+khtba2IiCQkJLRs52G0RYsWSX19vXTp0sWapqoSFxcnx44du+hnzw3YpgL3n9P+WdeN39PUZ86dvmXLFjl58qTExMRIZWWldO7c+aL9OvekrM/ni5jW2HbjTtLnn38uTz75pBQUFMiQIUMkGAzKq6++KuvXr7/o9zQ0NMiAAQNk8eLF5733z19g0YLLHV20a9cuefLJJ+W9996TwYMHy9SpU62CO3LkiJSVlckLL7wgt956q+Tk5Jy3seXm5kppaakcPXq0yfbj4+PlzJkzLb4ciD719fXy8ccfS0FBgZSWllqvLVu2SEZGRkRgFRcXW38/duyYlJeXS8+ePSPaKikpsf69fft2OX78eMQ85+rVq5fs379fKioqrGl//PGHhEIhycnJEZGzOzv5+fny/PPPy/Tp02Xy5Mny999/u7L8jdauXStDhw6Vxx57TPr37y9ZWVnn7XE3tR3ddNNNsmPHDrnmmmskKysr4hWNV+sQ7JcpHA5LZWVlxKuqqkrOnDkjU6ZMkdtvv12mT58uhYWFsm3bNikoKBARkQ4dOkjHjh3l3XfflZ07d8oPP/wgTz31VETb999/v3Tq1Enuuece+eWXX2T37t2ydOlSWbdunYiIZGZmyp49e6S0tFSqqqokHA63+vLDm1auXCnHjh2Thx56SPr06RPxuvfee2XRokXWvC+//LJ8//33sm3bNsnPz5e0tLSI67rj4uJk5syZsn79etm0aZNMnz5dBg8ebB2GacqoUaMkNzdXJk+eLJs2bZINGzbI1KlTZcSIEdZhnUcffVTS09PlhRdekNdff11UVZ555hlX10NWVpaUlJTIqlWrpLy8XObMmSO//vprxDyZmZmydetW2b59u1RVVUldXZ1MnjxZ0tLSZPz48bJ27VrZs2eP/PjjjzJr1iw5cOCAq31sFW17iD+6TJs27bxLpUREe/TooS+99JJed911WlVVZc3/5Zdfanx8vG7evFlVVb/77jvNyclRv9+vubm5umbNGhURXb58ufWZvXv36oQJEzQpKUkTExM1Ly/POml1+vRpnTBhgqakpHC5IyKMHTtW77rrribf27hxo4qIFhQUqIjoihUrtHfv3hofH68DBw7U0tJSa97Gyx2XLl2q3bp10/j4eL3lllt079691jwXOol/scsdP/roIw0EAlpeXm7NX1JSovHx8frVV1+patMnT2fNmhXxHRkZGRFX0qhqxDZ0+vRpzc/P1+TkZE1JSdEZM2boc889F9Hfw4cP62233abt27ePuNzx0KFDOnXqVE1LS1O/36/dunXTRx55REOhUJPr1ct8qgxmDVwJ1qxZIyNHjpRjx45JSkpKk/N8+OGH8sQTT8jx48dbtW9wF4diAMAwBDsAGKZFgz0cDsu8efMcneRzow0v9cW05fFaX7zEa+t3zZo1cvr06QsehhE5+4CsCx2G8dry0JeLaMkD+KFQSEXE0ckHN9rwUl9MWx6v9cVLTFu/pi2PiX1pxKEYADAMwQ4Ahmn2IwUaGhrk4MGDEgwGL/i8h+rq6og/m8ONNrzUF9OWpzX7oqpSU1MjnTt3lpiYtt0naa36d6sdr7RBX5y1Y3cbaPZ17AcOHIh44A/QWioqKqRr165t2gfqH23pUttAs/fYG5+C1u+ZxRLrT2xuMyIi8vXMIZee6RLcusvK9/OnjtuITety6ZlsqO/5b8dtuLE8IiK//mu84zbyrnNWJzU1NZKdnX3eE/jaAvV/YdR/05zWv4j9baDZwd74389Yf6K0uypwibkvLikpydHnRVws7ETnT1iMbe/8BygiUu/CenFjeUREAi6EaVKSszppdKlHvbYG6v/CqP+muVX/IpfeBjh5CgCGIdgBwDAEOwAYhmAHAMMQ7ABgGNtXxYTD4YiH0zi9GB+IJtQ/oontPfZXXnlFkpOTrRc3Z+BKQv0jmtgO9tmzZ0soFLJe/xy0FjAd9Y9oYvtQjN/vF7/f35J9ATyL+kc04eQpABiGYAcAwxDsAGAYgh0ADEOwA4Bhmv3Y3kZfzxzi+LGjBUXOLx17aqg71xXHZ/d13IbW1brQE3eU9/2PK+0MD21x3EatDHShJ97iRv0vKD7guB+zBrsz8Aj137Roq3/22AHAMAQ7ABiGYAcAwxDsAGAYgh0ADEOwA4BhCHYAMAwDbQA2UP+IJgy0AdhA/SOaMNAGYAP1j2jCQBuADdQ/ogknTwHAMAQ7ABiGYAcAwxDsAGAYgh0ADONTVW3OB6urqyU5OVkqKysdDzTgJb627oDLmvXDbYIX1kt1dbVc26mThEKhNq85U+vfS9yoObfq3w1uLI/dbYA9dgAwDMEOAIYh2AHAMAQ7ABiGYAcAwxDsAGAYgh0ADEOwA4BhGEEJsIH6RzRhBCXABuof0cT2IwWa2mNJT0837pZqL9w67yYeKeCOK6X+vYRHCpzP7jbACEqADdQ/ogknTwHAMAQ7ABiGYAcAwxDsAGAYgh0ADGP7qpgL8f38qfgSExy1EZ/d12k3pLaL8zZERP5v238dt9Hr6vYu9ESk37XO1quIiO9MrQs9cUlsfFv3wHVeqf+wS/X/qYfqv69h9a8u1L/dyzfZYwcAwxDsAGAYgh0ADEOwA4BhCHYAMAzBDgCGIdgBwDAMtAHYQP0jmjDQBmAD9Y9oYjvYZ8+eLaFQyHpVVFS0ZL8AT6H+EU0YaAOwgfpHNOHkKQAYhmAHAMMQ7ABgGIIdAAxDsAOAYRwPtBGb1kVi2yc6akPrvPMwfDcGCfC3887vy1qfO4NbxPpcaMN5E55D/Z+P+m9aa64V7/wEAACuINgBwDAEOwAYhmAHAMMQ7ABgGIIdAAxDsAOAYQh2ADAMIygBNlD/iCaMoATYQP0jmjCCEmAD9Y9owghKgA3UP6IJJ08BwDAEOwAYhmAHAMMQ7ABgGIIdAAzjeASl+p7/lvqkJDf64ogLA5yIiEi/axNcask5daGNOJd+dbuxfp0ujxvrw22m1X9fD9W/G9yq/2hzhS42AJiLYAcAwxDsAGAYgh0ADEOwA4BhCHYAMAzBDgCGYaANwAbqH9GEgTYAG6h/RBOfqtq6oa+pPZb09HSprKyUJIPuvPMSL91p6YU7T6urq6VTp04SCoVaveaulPr3Us3hfHa3AQbaAGyg/hFNOHkKAIYh2AHAMAQ7ABiGYAcAwxDsAGAYxwNt+H7+VHyJzh7OX973P067Idkd4h23ISLiO1PruI1anzt9cWOQgE+3/dd5IyIy9aodjtuoyxrqQk+8hfo/H/XftNasf/bYAcAwBDsAGIZgBwDDEOwAYBiCHQAMQ7ADgGEIdgAwDANtADZQ/4gmDLQB2ED9I5rYDvbZs2dLKBSyXhUVFS3ZL8BTqH9EEwbaAGyg/hFNOHkKAIYh2AHAMAQ7ABiGYAcAwxDsAGAYgh0ADON4BKVf/zVeAsGgozaGh7Y47YbUdRjouA23xPrcaceNZtwY+UVE5LcOAxy3kePw8y6tVld5pf5rDax/N3ip/nu60A+72GMHAMMQ7ABgGIIdAAxDsAOAYQh2ADAMwQ4AhiHYAcAwjKAE2ED9I5owghJgA/WPaMIISoAN1D+iCSMoATZQ/4gmnDwFAMMQ7ABgGIIdAAxDsAOAYQh2ADCM44E28q5LlKSkgKM2asX5IAGuPds/Nt55Ey50Q0REXWijLmuoC604HyRDxPnyuLE+3GZa/asL9e+lvUW36r81B8lwg5d+BgAAFxDsAGAYgh0ADEOwA4BhCHYAMAzBDgCGIdgBwDAMtAHYQP0jmjDQBmAD9Y9o4lNVWzf0NbXHkp6eLpWVlZKUlNRiHbTLtTtPPcRLd1q6sX6dLk91dbV06tRJQqFQq9fclVL/Xqo5nM/uNsBAG4AN1D+iCSdPAcAwBDsAGIZgBwDDEOwAYBiCHQAM0+yBNhqvkqypqXGtM05wuWPL8sLljo21ZvMK3RZlav23/ZrFxdjdBpod7I1fkJ2d3dwmgGapqamR5OTkNu+DCPWPtnGpbcD2DUrnamhokIMHD0owGBSfr+n9hcabOCoqKpp9E4cbbXipL6YtT2v2RVWlpqZGOnfuLDExbXsUsbXq3612vNIGfXHWjt1toNl77DExMdK1a1db8yYlJTm+O8+NNrzUF9OWp7X60tZ76o1au/7dascrbdCX5rdjZxvg5CkAGIZgBwDDtGiw+/1+mTt3rqNnbLjRhpf6YtryeK0vXmLa+jVteUzsS6NmnzwFAHgTh2IAwDAEOwAYhmAHAMMQ7ABgGIIdAAxDsAOAYQh2ADAMwQ4Ahvl/Il/FRVH3GVEAAAAASUVORK5CYII=", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXYAAADJCAYAAAA3tRlxAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAABAj0lEQVR4nO2de5xdRZXvf+fVJ92dTifpBEJCTMwNBBIIBgICwiAjOvgAncsIowwQdOZz1ZGLjnpH7sjwcD7DZ+4I4p3xqiiP8YJXnUG9g44X0AFhjAEDJjEQCa+QEBJCXt2dpNPd55y6f+xa+9SuU1W7ap/TSffJ+n4+/enuvavWrr332nuvWrVqVU4IIcAwDMO0DfnD3QCGYRimtfCLnWEYps3gFzvDMEybwS92hmGYNoNf7AzDMG0Gv9gZhmHaDH6xMwzDtBn8YmcYhmkz+MXOMAzTZvCLnWGYccGKFSvwgQ984HA3oy3gF3sAK1asQC6Xa/i58MILD8nxb7zxRrzlLW85JMdiJi4rV65EoVA4ZHrZKr7yla/gnnvuGfPjHAkfkOLhbsBE48ILL8Tdd9+d2FYulw9TaximkbvuugvXXHMNvvWtb2Hz5s1405veNKbHGx0dRalUalpOb29vC1rDAAAE481VV10l3v/+9xv3PfLII6JUKonHHnss3valL31J9PX1iddee00IIcRPf/pT8ba3vU309vaK6dOni/e+973ihRdeSMjZsmWLuOyyy8S0adNEV1eXOO2008SqVavE3XffLQAkfu6+++6xOlVmgrJv3z7R09Mjfve734nLLrtM3HTTTfG+Rx55RAAQP/7xj8XSpUtFuVwWZ5xxhli3bl1c5u677xa9vb3ihz/8oTjuuONEuVwWF1xwgdi8eXNc5oYbbhCnnHKKuPPOO8Wb3/xmkcvlRK1WE6+88oq4+OKLRXd3t+jp6REf/OAHxfbt24UQQmzYsEF0dnaK++67L5Zz//33i3K5HB9ff77OO+888clPflJce+21YurUqeKoo44S3/jGN8S+ffvEihUrxOTJk8WCBQvEv/3bv8V1KpWK+MhHPiLmz58vJk2aJI4//nhx++23J9quP0ePPPKIEEKIV199VVx66aVi6tSpYvr06eLiiy8WL7/8ckvuy6GGX+wBuF7sQgjxuc99TsybN0/s3btXrFmzRpTLZfGDH/wg3v8v//Iv4v777xcbN24Uv/nNb8RFF10kTj75ZFGtVoUQQgwODooFCxaIc889Vzz++OPi+eefF9/73vfEypUrxYEDB8RnPvMZsWTJErFt2zaxbds2ceDAgbE+ZWaCceedd4rly5cLIYR44IEHxPz580WtVhNC1F/sJ554onjooYfEunXrxPve9z4xf/58MTIyIoSIXuylUkksX75crFy5UqxevVqcccYZ4uyzz46PccMNN4ju7m7xB3/wB+Lpp58Wa9euFbVaTSxbtkycc845YvXq1WLVqlXi1FNPFeedd15c76tf/aro7e0VmzZtElu3bhXTp08XX/7yl+P9phd7T0+P+OIXvyg2btwovvjFL4p8Pi/e/e53izvuuENs3LhRfPzjHxd9fX1i//79QgghRkZGxF//9V+LJ598Urz00kvi3nvvFV1dXeJ73/ueECJ6xi699FJx4YUXxs/R8PCw2L9/vzjuuOPERz7yEbFu3Trx7LPPig9/+MNi0aJFYnh4eCxu1ZjCL/YArrrqKlEoFER3d3fi5+abbxZCCDE8PCyWLVsmLr30UrFkyRLxp3/6p055O3bsEADEb3/7WyGEEN/4xjdET0+P2LVrl7E8WUoMY+Pss8+OLdTR0VExY8YM8fDDDwsh6i/27373u3H5Xbt2ic7OzvjFRz3DVatWxWU2bNggAIgnnnhCCBHpYalUEjt27IjLPPTQQ6JQKCQs+2eeeUYAEE8++WS87b3vfa8499xzxTve8Q7xzne+M/7oCGF+sZ9zzjnx/5VKRXR3d4srrrgi3rZt2zYBQPzqV7+yXpNPfOIT4pJLLrEeR4jog7ho0aJEe4aHh0VnZ6d48MEHrbLHK+xjD+T888/H1772tcS26dOnAwA6Ojpw7733YunSpZg3bx5uv/32RLkXX3wR119/PVatWoWdO3eiVqsBADZv3oyTTjoJa9aswbJly2J5DBPCc889hyeffBI/+MEPAADFYhGXXXYZ7rrrLlxwwQVxubPOOiv+e/r06Vi0aBE2bNgQbysWi1i+fHn8/wknnICpU6diw4YNOOOMMwAA8+bNw8yZM+MyGzZswNy5czF37tx42+LFi+N6p59+OoDI/3/88ccjn89j/fr1yOVyznNaunRp/HehUEBfXx9OPvnkeNvRRx8NANixY0e87etf/zq+9a1v4ZVXXsHQ0BBGRkZSgw6eeuopvPDCC+jp6UlsP3jwIF588UVn3fEIv9gD6e7uxsKFC637V65cCQDYvXs3du/eje7u7njfRRddhLlz5+Kb3/wmZs+ejVqthpNOOgkjIyMAgM7OzrFtPNPW3HnnnahUKpgzZ068TQiBUqmEPXv2OOvqL1jTC1fdpuo1HcdUR9++du1a7N+/H/l8Htu3b8fs2bOd7dIHZXO5XGIbySYj6fvf/z4+/elP49Zbb8VZZ52Fnp4e/P3f/z2eeOIJ53FqtRpOO+003HfffQ371A/YRIHDHVvIiy++iE9/+tP45je/iTPPPBNXXnllrHC7du3Chg0b8IUvfAHveMc7cOKJJzY8bEuXLsWaNWuwe/duo/yOjg5Uq9UxPw9m4lGpVPDtb38bt956K9asWRP/rF27FvPmzUu8sFatWhX/vWfPHmzcuBEnnHBCQtbq1avj/5977jns3bs3UUZn8eLF2Lx5M7Zs2RJve/bZZ9Hf348TTzwRQGTsrFixAn/1V3+Fq6++GpdffjmGhoZacv7E448/jrPPPhuf+MQnsGzZMixcuLDB4jY9R6eeeiqef/55HHXUUVi4cGHiZyJG6/CLPZDh4WFs37498bNz505Uq1VcccUVeNe73oWrr74ad999N9avX49bb70VADBt2jT09fXhjjvuwAsvvIB///d/x1/8xV8kZH/oQx/CrFmz8IEPfAC//OUv8dJLL+H+++/Hr371KwDA/Pnz8fLLL2PNmjXYuXMnhoeHD/n5M+OTH//4x9izZw8++tGP4qSTTkr8/NEf/RHuvPPOuOzNN9+Mn//851i/fj1WrFiBGTNmJOK6S6USrrnmGjzxxBN4+umncfXVV+PMM8+M3TAmLrjgAixduhSXX345nn76aTz55JO48sorcd5558VunY997GOYO3cuvvCFL+C2226DEAKf/exnW3odFi5ciNWrV+PBBx/Exo0bcf311+PXv/51osz8+fOxbt06PPfcc9i5cydGR0dx+eWXY8aMGXj/+9+Pxx9/HC+//DJ+8Ytf4Nprr8Wrr77a0jYeEg6vi39icdVVVzWESgEQixYtEjfddJM45phjxM6dO+PyP/rRj0RHR4f4zW9+I4QQ4uGHHxYnnniiKJfLYunSpeLRRx8VAMQPf/jDuM6mTZvEJZdcIqZMmSK6urrE8uXL40GrgwcPiksuuURMnTqVwx2ZBO973/vEe97zHuO+p556SgAQt956qwAgHnjgAbFkyRLR0dEhTj/9dLFmzZq4LIU73n///WLBggWio6ND/P7v/77YtGlTXMY2iO8Kd/ynf/on0d3dLTZu3BiXX716tejo6BA/+clPhBDmwdNrr702cYx58+YlImmEEIln6ODBg2LFihWit7dXTJ06VXz84x8Xn//85xPt3bFjh3jnO98pJk+enAh33LZtm7jyyivFjBkzRLlcFgsWLBB/9md/Jvr7+43XdTyTE4IXs2aYI4FHH30U559/Pvbs2YOpU6cay9xzzz341Kc+hb179x7StjGthV0xDMMwbQa/2BmGYdqMpl7sw8PDuPHGG1s+iDeR5E6ktrLc1jLRzvmss87CDTfc4AyrXbFiRbAbZqJdhyPhPdCUj31gYAC9vb3o7+/HlClTsoqZ0HInUltZbmuZaOfMcsdO7nhrK7tiGIZh2gx+sTMMw7QZmVMK1Go1bN26FUDUXWglJG8iyJ1IbZ3ocoUQGBwcxOzZs5HPH16bhPWf5Y61TJNc32fA28c+PDyccOBv3boVixcvbqbNDJOJLVu24Nhjjz2kx2T9Z8YTac+At8V+yy234KabbmrY/q+r1uHaO9YCAFb97TszNDGavgkAOe1/ojBU/wqKQpQAqNaRHNnX0w9VFCEdQ1HulVrXdKN8tW7xjRcAAGuuvQ4AcNJ3/jlRp6ZULuT82m9qHxEycm2Sr8v1kWer49P+Zs7VdUwfBgYHcdxxxzVk4DsUuPT/z/8xmrK++kvvtdYPOV+9bF7V/3wh+l3uTpSFVqeq7ChJ/a9K/ScdzhsaU9L0f4nU/5D7nVUn9Pqt0DnX8xDy7NjaYqo/VmUBYNDzGfB+sV933XWJ3CYDAwOYO3cuuif3IN/RBQCZR4NTX+yl+pb6i70rUcb5Yi+ORnW6pxjlJ17sBycDACYXo0tD58Qv9vCyLrK82OM6Kalex4Jm9b+pF7uq//GLfXKiLLQ6iRe71P+q1H/ni92i//xiN7fFVH8sX+xxe1KegczhjhSG8/r27fHNH6nV95c094/rBGwn5WoYKWchQGMqWh2fF2XuYGQt1SZNadhva5/PjUyTYWufTb4Nn5e0a3vaA5P1Pqadg6nO4MAAjp41a1yEPpL+b9++Hb2yLaox4dIxGyH3lV7c+svZdQ9tdVz188P7AAA1+SFx3W+XnqbpcBYDx9WGkGfQVs517Fbpf+izPeD5DHBUDMMwTJvBL3aGYZg2o6UrKKnul0UfiwZd/uO2/wwAmNkV+QYL+94AAOzvnBGX7dT8KXqXhtwhKrlJyW6IT1dOd9u4ulMx8jhrL3o3AOC0++6Kd1WnHGM8jqstzfhafc6x4do52qXLNcnP4qdrxnXkqjMe05Cq7VT1i/T/sdv+EAAwqyt61PL7dwEA9k+qL39I+m+7bqr+58hz2tmbqOO6xrFbxVLIpSPkyyf9P/U798T7aj1HG+W49KkVuuB0nXrIS3OZ+MgIcVu65IW6VX2fAbbYGYZh2owxW/OULPWC9ukYkpb6JMOop89AiChNMtYhfCw9lwVgG5QlS6VWSl+X1PU1DxlwtZF1cDNtgCmk3c1GKIT0QtLkjhfU9pGlnkfSGj/Y2QcAmKQ8FyE9HD3MV5cR0lMzytf0n+pST7XWUV/rNIsepdV1lXURMohqK5PF2s86kGuTq7dRl+P7DLDFzjAM02a01GJXvybkUydGZShk2eBPTA3dmtQY1lO1WNbNfCVNbaAyuj8RqPs+qX0hVpNPW5oJd/SxqNMsLlcZr/EJh9xmyow3TOdNPnXaR/rakW+sY4srj8so+h/rY0roYog+qZA8vQyNJ6nbKRSS/PCuY9v2ZdEjEyFWcivGjZolpOee5Zlgi51hGKbNaNpiF5a/KfqFfOpkqS/9bw8CAH59y7visrQPtQoAoJqTEQRy84gyjY5mXJHlQ1bzcEdk1egTowAgr00yEj+LVmzPn/PBxHYAyI8cAABUSl2JNhR3bQIA7J3yprjsZK0nkRNRt0Tk8olrAACiWI5+63Vof2Ozg6wEqwVQHWloX63oHqcIke9TtpVRAeMNAfM9pOgX8qmTvi677iEAwKq/qaff0PW/IvWfNg9XG8+8Q+7MW/Q/Mbai6T889H9U6n/Bof/dhklLKgn9p7Epy4xZEy3RDXlNASAv/07Tfx+/ueu4Pr52275WRRKxxc4wDNNmtDSlgMqQtDL06Beyvs//u8fibSv/+3kA/L5ehaF+AEBVxvGS754OE0/llv4/ABAyl0eDJS0t2KriPyfjqKj5GncfrAIAjtm1Pi47POeURNv0dueVNhS2RvUqC85M1GnWJxjXk+eSGz0Y/UvRQ0LJ8yDPn65DTqtT0aw0VxtCLIws52iymsZjSgFV/9VzOJii/xd9/cl420N//tZEfVdUiN771I150tuctLyBuv7HMhRLGgBqk2fGf1c0/SdI/4/e/Uz9XGYvTbRXv4eq/he3bQAAjM4/3Vg265hYXJ90WfZQRaFDVmrU//i3VqdSiJ4ZV8oFvb1Z9T/tXWcryykFGIZhjlD4xc4wDNNmtMQV02PoEujdEVt4ItDotmlm8kqmUCllgEXkiwl5OolJF1QvX2zYZ5WhdxVN7bEc02sQJaVNxiZZut+hx05z1/i4YlyTW3y7oYeCNFeMfh6ubKSk/2mpNUzyyG1AZXwyOOrkFbeNng5bJzEAWIlcePpgpPNYljqm5mYaYPXQf/3aNaP/WQf0QyZhqbo0MDCAWeyKYRiGOfJoSbijjyUWWxDyi0qWMVC31E+VoWAPXv8OAMBRcpKTIdqrbvko1nZ0IGlxP/rteNPwuVcAqIeV0cBSbnQokt/dV6+eEj6pWvd0rPzg6wCA/V1HAWi0vADlyy8t9dKOjVHz5SCwOoBr/XrLwZ7hWl1+w4SvfPot1a1JvbnqwBsKUp5sd3weWmhnoj4NWNEgrXY8U1uyTMMeD5isdH074Nb/Ti0U+Gc3XgDArf+xpU7hrKSLdM0V/R+S+t+p638lWuqv2jWtLlcOeA6XorDEhvBhRf/J6i5I/R+S+l82dUu0Oi79Jxp0QurViKg3qiG8U9N/U0v0Xo6ePiEfoP/w0H/CFBoeov85y3YTbLEzDMO0GU1b7Dn4hfMQNPmogEbIUn/bp6KUp8/f8cdRWcdnKn9gT3QcmtxBU5/ffmVcplMLkYy/juV6QiPiYClaS7Bs+eSplhbJeb0YTcKaufI7UZlzL0/sTxxT/q5JK6nY/1rUNoPFotetyu9wGfVJRwIdxrIhEx0a/H1KojWhWd10nel3tVAvWxDJfTSOYLTqtXa6fI7j2Xq36T8s22jykenBI0v9vM//FADw7P98HwC3v5z0nyA9UvW/i1JfFJI+2ZrUf1X8gWK0rdOQ+gCoP79A/RneXoj0/yip/5D677qHpP+FgW2Jdrsg/e8Qdv338VWnPQet1v9YvvIcjPVkO+8Xu75K+8BAY450hmlXWP+ZiYR3VMyNN95oXKV9u2WCUsiUWb0OsXc48lVNNZjPWQL8yXdGI/96JAHQ6H+zHdd1bNvIvwuqo9bLcg11XNEGIdE3ejQG/Z+IJCDfIv12RCakHXu8RcXY9N82QU8ny7T1foP++45BmK5fTtN/wrRAuw+2duvPmQ9jpf8uQiJTbJE06vWK/zSMo5iOZzq2S/+BMZigdN1116G/vz/+2bJli29VhpnwsP4zEwlvV0y5XEa5XG7Ynh/ej8JQ9H0hHzaAOPFQ/AWSU6BpSnWHwTTQv4JkqZz9t7+Iyzx23XlRw7WRbF2a6nvcX54KAJhEljpFEkgfmOrDLP7u0agtx50dyZdliuQLVCJoaF+D/1yzOFztpO1Cse5LO18CAOyaMh8AMEVmkMoSD+6zT7dUXBa2q0cT+yNzyfY6kfdiNBddyzg1RK3uRx3JGaaHH2Js+p8b3o+81P+aQf9jpP4PG/TfFoNO+v/2W/8jLvvQp88BYE52p8qi1BsAsF9GeHnp/8bHAQCV//TWRBnyhdeUCBpInW2wNh2Wuk131d5tSSYc29UTJRzT9d8kz6ZrrkglfbtLlj7Oob9/EnI8ItNi5L2o5pPvElX/K3JfETXvZ4CjYhiGYdoMfrEzDMO0Gc2HO9aqEIWSdb++RinlU1e7iqoLxwS5XwDggtuibukjnznHWWe0s95lpNbFvSktu5s6kaBywtuj9sl82uR6ETLzYfHFJ+pljz/X2QbXQKvLnUIumL7+yCVTmbmwQZ5+jGYm+oQMItmO75JLGQNLo41T12uym9mhNYK6puq+IUPY5OEmV6tC5BuDd3MyJoHWKG0YYFdcNTXDCmFqHXK/AMDF34iyQv6/j5+RKKPf99FJ9WeqQf8dg3qjUqcL0pVZ65Ihq5T5c9PT9bILz26or7bBdAwfvdw5OXLBzND0H4ZJcWlulGYHWrPI1dtk0n9yV8VuFu2CVBT9r7t98sbQYRPj70lhGIZhmqLpJGBquKMrbEofGM0SAqdu2ycTsfeU0gfqmhlgoTAsYQhd1JMg+UwOyhKqmHO0gXCF04XiMzkoZFKaSS4RJ2CC1ntSpq7TpJj9g+M7CVjVQ/9Nk41s65e6Qt9I/7tL5rQNpntgewZ99N8Uumt7NkIG8L2se22w1yRPr2vaH6LLNvk6Ic+ZKqshbNhD/ws5zsfOMAxzxNK0j13FtfJOvLKR9C2OdNS/NnHolvxKUagirexi+mqSpf7Bb68BAHz3irckjjOiRAUNVaJ/emXYVMfWtQCAbX0nAQD6JtV9pPEpWNYHzSnhRnESpKceiDYsilZHqsp2m77mIZMu9FBImvhhmvJsq2vCZ1KXLifEV2+Ta7IM6zqjJUzKN05db9XElFYioEzMySW3q9A+8q2Pluv6Hz8b0jKN9d8wzZ7kkqV+qdT/71/5lkQ5k/5T2GDptXUAgK1TlwCoJxsD7Pqvb1f3xfp/QuRzV0OCbbj0v2EfhWWS/hcVy10bL7CNOaUd0/S/Xj+tbNo+VVbju1LzPBj0PwS22BmGYdqMlljsIb4pSodrnIpL/9OUXA/5ZKlf+X8iK+S+D0frMKoTOCaNRhZQtSOyJF6YshgAsOC1XwMARrV1SAEkEvaobTBZyIVjj4v2SUurqqwhqbfXB+skDookUYTprcniW/fx/7ssIV95Ln+nq/cQ0ms41NjaZLsGpP95Q9mQlMsEWep/8p1I/++V+t/h0P+Xe04EAMzfHum/vg5vdKCk717fru7Lz10U/UE9DYokc7TbRcP7QP6O1y4OkOvTc3WV9el92uRm6RG45Kq9wzTYYmcYhmkzWrI0nml0Ns2f5TPV11TWJp/oH0n609Wy+hJYsW/UI6bYFaFgi1FWZaX56nwSMflEMaR98U1lQiIVXNja0mxcM5UZHOdL4/lgOl9bxEyWKI39MlpmcqlR//Xl+WL9NyyN14ooLlMkjY/cLPMqsvjCbTKy9ixttOq5BXhpPIZhmCOWliyNV9EsYRMNviS5BBcAiHK0DBct50WLBMQ+PCWhF80oJeuGRv/JQCFL/cybH4nr/PqaaOZabtpcAEBhT5SZjxKTqTNfi/1bAQDbSrMAKBEDMhpg1LAsF1noA4XoPHqkXOcXWkZA0FJ5IelSE3Is201WQpo/2FUuzdp3yTftt/USfC2X8UJVuGPUbeenLj+YJytZ0/84ckiZpU0zSvXoL+qgkqV+/m2/jOs89tH5AIBK75yortRx0RHpq0n/Xy1GETnHdMtXhNT/imILUht0/Z9csifBI/TYdFdUDP02XecQS923bEgEmEmOtVeiRBTpifJ8es0hsMXOMAzTZvCLnWEYps1o2hVTFUDHgShhVs0wMcE2sEjuF5VhWk1dW6OU8qkDjQmNaPJFHNIl20DuFwBYdNNvAQDP3h65YuIJDht/Ff0+5cK47BvlyAVz7K71UZu6TpEHlOuNVusrvYh81OU8INdJ7R2Jzq1a6k20UUXfZlqJxWeCQ1pZ3/1p8m3yQiaAZO3ahrTrcCEAlAz6r3exCZqgJwyD9ENS/7u0MvvVyXx6HYv+k/sFABbdsBYA8Mztc2QjpD0n9T+n6P/uzmjN4HlvRM/MSHek/xTmW1JWOqKJc7H+D+8FANRK0xrOzYZrYFHf55oASeh6GRLCO+ZhtR7Jy1zXIxfQNrbYGYZh2oyWhjtWFElZBgP1r5dp7UTb19UW0qiy40AVADCzKzlJ15RCOB64lftodZxmQ6GylMlCq62OkF6ErW5IKKbpOh/ONU91VP3vMSTBS0vo5TMBy7RGqWsgUW2DaSD3DYv+q+G+ek9C72H4hPuGEDJY3kwIpkmuT1hlmjUfEooZIt/2nvF9BrxdMbxKO3Mkw/rPTCS8X+y33HKLcZV2QAnLcoR7xWV/dmdU9syL421x0ixpoVCdmr5GI1D3U8np15TQi9IEvKkn8kJSSCNQ96nPlEmVjn7PzQCAl/7vFwAAk5Vwr4PSFOqUJ0OWOlnu/cWeuCyFlhV3vxKdx/R50fHiE7O322UJZLG2s/gHW+ETzxpWFhKmOR7CHF3677KSG66T1P+8Q/9juSb9z1MK12RCL0oTcCzpvwxbjAREZWdOkf7zD30NAPDst/8LAKBbsdKHZLe7S3Z5aZJdwaD/lIivKJ+1igwnji+DknqWEB76D22bq5fbTKhiK/TdR0+z9ARscnyfBW8fO6/SzhzJsP4zEwlvi922SrvNCtO/LLHf/JwPRn+MDjXIysltotydlKcm2NdW6abUu3FCr54ooVHCV6hFv5Cl/odymbGHP1lPglSSZldei8wZkRNDysqJUfsO9EZLeY1oi3+oKX5FSnrRkIk+LjkhNDO5w3VcH398WhnT/sNpudv0HzBb6jZi/a8MN+zLyW01qf8xjgUmKPUuJfSK9b9DiTqL9T+y2MlS/8DXomUeH/6vZ8VFOwpu/e80tOHAlGOjMrr+KxY7pd5tGEdoOLNGWh0x1creraktWcYAfHoCHBXDMAxzBNOSpfF6pzQ3Yh6n6dWsWh/fJRH3CAwJuIiCFuFCbNtftyxoCrWX786wGLapbkKOtiBvq6JisozWt5o0qz6rlUNl+j0TIB0KTEtDqoRETOnRL67IlzSLzRVJk7M8G9sd+p/l3oVESplI0wmfiKJmLGqXfJ92u9qQVjdtjImTgDEMwxyh8IudYRimzWg6pUBh/04UD24HAFSOOj7eTt0H3Z1Ck45GS/WuYkG6YMiNclBOUaaBnOLvHo3LVk54e0K+3v2hbmZRCfeiNAHTpAuGQhppoDTOYAfgyu9GU6m/ekk0KEUDQeR2ye97o34sfU1KWUZffUltr+62adiP+urvlYJc61SepGug1WcyR5auqA2XDJ/BqVSXQqC8w0Vh3w6UmtV/6TbJy4ynB4rR4CmF3BY3Ph6XHT3+3MTxG/RfylL1n9IE9Mpng0Ia6flS9f/q70epNP7hD6Pw4claIIBJ/+M2ONyMQW4Jqf9VD/3XCdmXJeQypA2tdA9RWR48ZRiGOUJpevB099pH8fznookbS3/0k3i/7atlSnpF6Lml65XqEzTyBwejTTLZUU4bwKTj0vRpADh2T2SFjMw5JSFWD+kCgH0yZGvOeddE57nyq5F4aUUIZVUY/Rz1cwuxrFUqmhyfgZUsVnaWEMaQwTSfQcAQueMxpYCq/ycr+k/o52VKeUH7Rn30f3h/tKlrWiTfMoCvBgTM2/sMAD/9pxWYZjeh/7qFbcKla7r+6/jop749bV+afNvxTPXGMjTS9xlgi51hGKbNaGkSsJwjmVBIeJPPlzXvWFdRl0H19XBHn4k4mwdHAQBze/SEqXZryddfZqPhOmi+yxCyWDdZJxKF+PvT5Jna4BvqdSgwhTvmDauCtWKylnGfpv9e/mJLuKPrOE79txwzRI9c2Nrl89yOdRuyyHS1hcMdGYZhGCdNR8UA9S+KaqWvvejdAIBTv3MPgPoIenHXJgDA691z47J9k5JpRBsaObCtfiwZTaBPMmr4sqqpB6Slqyf0omnSiRS/FNki/ZpzeyKLaNrvfQoAsOXnt8VF1ZXgfQnxtdd9lskogxCLxeVjTJPhU8dXThr6Ofmc43ggJ2pxu2rK4jGk/6fddxcAoCoTcJH+71D0f1qK/hcM+q+nl25ATWdBPT0todeItn4qnQ/QqP+9Z/85AGDrL/4hLttTMvdUs4zDmOTYUhP79HZN/6f1ALL4vVWase599Z+jYhiGYY5Qmvaxv7at7mNUv/xkZdQ6ophcsuYp6mT6jt/GZYflaL0ttYCatrT4YpS4qELxvBb/86hisJdrSX8ktaFciOqoUQj5wdejsrKHQReH6lx139q47A9WLEuU0dMFqAsvFKvmMQGXlUAxz0KZHq6TJTbXNzrGRBb5tuO45JgsrfGYUmDrtu2Y1ivHmJT9pEdC6j9Z8xR1Mu2N9XHZ0dlLoz+k/ldzlJo32pxTlqMrbHo6qrPw7GifJSqsqlzAUtWs/53FfOI4AFBI0f+P/XO93fd+OGo36bkexaK2oVCLnmGhJTRz6YS+0E5Ij1WX5Sqjlw2Jismq/yFjVyrsY2cYhjlC4Rc7wzBMmzHm2R3j/11rJ1Le5rx5LNckVw/30vf7TEzQwwmjjckurSsb46syFGyODAXzCuvzyO7YgKW77aqfdQDIFw53bFzzVEU/HwqFVAdYCT101+dapNUxtUUnZEKavh0AXtsXPbdzJjefEdJEiHslS8hlFr30baNKlrZwuCPDMAyToOlwxxz8vr76hKUElARMDty8XpwBAJjRWYiPoUOWSuGpB6Lfxx4HABg5+oRIljJZ6oBMKhYnVZJrlNLKRzSIqhIf0zH5iCz1ZX/5EABg3Y1vjdrmmADVMMnIMskp2Ri5vuW2aOBqz4zF8S495NI1MJploFWv29ATM5TV5bqskBB54xGBlIlEVE5a6k5dlvq/vRDp/1FdhfgYDfI0/c/PXQQAGJWJyEz6P4n0X65RSisfdSijnjntt64zakAAWerLr/8ZAOCpvzojKuN61jVMk/xsemPS/+6Su+ebpUfpo28+vZxmeiW2Mr7PgveLnVdpZ45kWP+ZiYS3j/3GG280rtKuphTw8ms7vtCUTrdz5Xei7edebm+4/E1pRPMH9gCoWyz7lHjH3opMHKZN6hjU1miMDhptqyIZCtbQfuUcaMLTh34Y9QTuk2FgXv5Ow/iC7YbQOU2p1qeuk3XkM/EpzW+axV+YNfTMVtbVlsOZBKxZ/dcxXYsRqf9lTf9dsvL7d0W/pf5XZi4EUNdtAOgdjfSTEocRpE9qry+Xpv+Gc6DewYp/3QwAuOvSkxraaT0HQ4izzYKmc5pSO1BvgzZmkUU/08J/XWS951kmQ+UwBknAeJV25kiG9Z+ZSLQkKiZOAqbs9xlVb2iMViYt0ZeJvBZ94zpmSASNC70OpUyd1Z1uhbh6Aj7X0HaNWu2XznIdQnyO7RYVQzTjb80ZJqil6bIpWizLWEpa21x1KGX2zC53qgSb/JB3R96QTjitTkj0VjP3L6Qn4DvhiaNiGIZhjlBakgTMhf4lKkifuCh11stofjIqq8foqvL0soTuczaVoYUL4phyw9fe9oVWowIoOZFeliz1r6x6Nd52zVuPNdYxRd/Yjm20ZmXbqacy3BGdf1V2xHq3r4vLDs0+JdEGvU65kp7CwMf6ayaG3qcXNZ6wtVOP7CBM+q/Htsd1Pe5DQx3DM9OA9GvHKTwM1r0OyTfpvw5Z6v/4ZH15vk+eMSe9XZZjutpWC9B/Sl2ijw3o+g+PHpKrTVni5dNkhMIWO8MwTJvBL3aGYZg2oyUTlHxCi+pdxTIAoLC1niWusuDMRBnXIERaV9ErzIkGJw0pDNLkUJZGIH31JnK/AMCXfhmFgv3lOW9KyDWdc5ZByBHZnZw03J/crqxzWdQG1oa1Ov3FaCKL6hiwXW+Xy8Rn0CtkoCkkHPNQ46v/BLlgits2xNtG5p9urBMSStrwnFnKJdqiZVpU69naQFkaTfX1OuR+AYD/8R9m/dfr6n+bypqIdfngHikkesZV/dcDDfQ6/aUoHLox6UO2NrnKWl1phv1Z3DJssTMMw7QZTVvsoYNHFIZIVjqA+mCmtAB8rDfTPhsN1oee790RaqjLSISRWRJ6xYMzSuPIUvnO+mjw7EMnzTS3LRCqT4NFIhdtKezbGW1XV5qigSZpuZRrUfsHpKU+rX8TAGB0xoK4Spq17HNvsoSMTcTBUxWrhScHSkcVKz3L+r02fELqqKcacl3j+6JY6Wn11f02/dfbljhWyv8qsf7no4Fbk/7HgQYyjLQstw9IS92k/yG6GxJWHSIrS3222BmGYdqMlvjYQ7abICugtGMjAKAmpz5XJ89sKOvrS1RXXYqRx2mQ4UjA5fpa2uqZ6tDfZKnc8tgrAIDPnjMPQHIVp4bj2HfF+yjBWa0Q9YhciZj08M5u+bsiLZWObfXxD0q4RNPOadJMRa69WRCVuGy+MpLYR22LVwLyPKeG9maoc6hwhtXC0GM1lZH3I0T/0+T76L+zTRnK+PQWSP//5tFNAIC//L35AADX8sE+bSH9F1L/Kwb9j3vd2opM5FMfDdD/UanjRYP+0z5CX1nKB9sYk68ottgZhmHajJb42EMsKZcFTOlui/2vAahbLC751ggSxRdIay8W9DIe7U07XiSwlvht8mHqxyRL/csrk9ECalmfNrhG09Nk2XyY/TOXxGV6ZAK1WklaQtIa6TgQJZ/aU5oWl+3W9tF4BVmgLr95lrGTiUCIf5T0n9YLNlnsROr9VvQ/xFcNbZ+rbLxP03/bgjkqZKn/w6oo585nzp7bkrbodVzybHVbrf+UfC1kPLBZ2GJnGIZpM8Y8pYDNJ26yOqtyZXT67UrwY/O7mlZM131cLgvAZhU4ozXI107RMYbUxHp98qnr0QIA8OGTzJZayMi8LaZcxXaOXcX6nloxslToXuSkNUaWSt/AprjsgekymoDSGctUsnq6WFd7s5YZD4T0LE3Y9N+VBM92TFfvqOXXkSLJHMs96vpPPnWy1L/7TF3//3iJOWLGx1J34RvBZNR/6VvPFaOekEn/D/Zp+i/Tedcc+q+3Ie1cfc+XLXaGYZg2g1/sDMMwbcaYhTtmkdHgVqEwsJ0vxdt2TZkPAJiixQeSDFNoUcjEGer+jhYmWeXZiLtRhklONjFURnW/3C6zQn7qzGMNNZKkyQ2p4yqruwN65dicOpmjRGVl15N+Fw4kp3kD9YkjlO2w1t0HANgro/Smj+6Ky9K+/IHdAS0/dPhc6yyDw7H+79oUb9s5OXLdTS375+wPmWSWs+h/lkF40zFsbVDdL/9rdRQ88Ynlsx1HNctvpqyr3ZRpk8pMkfpfUfSfgjPI9VjNoP8Do9H+aSMG/R/qR2Fo0NL6JGyxMwzDtBktXczaZ3p5ltBIstIBoK8/st5pbcdWT1qpaJYKtYFCJoujynqLlnzZzbaJLPVbVyZDwbJcw2ZJs8Zc91z/Pz/4evw3WSw1GdJHlmJPR3T9a+W6BRevCtQ51b/hLcal/yH3xRk2aKmzq6ceCjtD6v+o1H8fKzwkJLaaQf9t8n0mNZkgS/3vLInDTJPCfAdGTXiFdGaQq2PUf2mN04SyySWp/6W++rHkvlpnL2qjfkf2tthvueUW9Pb2xj9z59rjThmm3WD9ZyYS3muemiyWuXPnJtY8TQjW/m/VJKbYV6eFgvmEy4W0LS5jSfSl7ktbo9R2DFebVP73b3cAAP7k5KOs9cZqgk+IxZ4mQ4XGMvbWImclTdku74osUjVEjCbuHNi0Hn1Lzjwsa57a9N+25mnI9crybFDiPGFJkxEyiSfLGIHtGK7jmeQ5w4glFApsCwM2HbPV7wOdkPBcY09Dhk/uRfT+8tH/4q5NGBjc5/UMeLtiyuUyyuVyekGGaUNY/5mJREsmKJm+upmsEGkBV6WHyDQiH/+tpeCsBazT6dOWuKxj0gXtK8mkQbtlwqCekn9UDOGybshSX/v6EADg1J76Yh9p05Vb1YYQy9PHKqNJZ73aPoqyMfmIK33zDUc9/GQZ+zBeaw/9j3uJ0lIn/dfHe6pKJVMSNhtZfNa6/scJswxlbeMwrrJkqa9x6L+NkEl9hGvcKESuC3pfkc1Nx6QoG9NzOzpjAUY7Gsc2TXBUDMMwTJvREos9y6iy6Ys3XItKlWH2HxqPTaPIhlQCWbBZmabzoG2U2nNKdV/UFpkwyOfLH+ILJEtld6E33jY1QH4zfv4sPQKXjzHEz3k4ooF8Uad5h1xrEyMisrM6RFL/E8fTF3UpmhenyRtumE+ER0ivi4j1vyZ7z8bF5dxyXWUJ0v838nXfch/MhPj3XTTzOnEdZ6x1mi12hmGYNoNf7AzDMG1GS1IKZOnCm7aXaRUUNHZBiYauIWVSs+xXt9lkmdrlU4fK0mARuWBcWSltLh5TGR0aKJqqbKvQxBHTQJtFflqbfMqaSKsfMpDVrFvjcBAyUK1CZSnjoa7/rslG0HL/+1w/H/eKD1S/Ow4WiFwwPlkpdXyeA9J/1f1C+m9zwYboUZZQ0RA5PueY5r7x1Su22BmGYdqMMc/HrhN/kShsC2iY4OOiFYMOIYO9IdYNlSFLJX+wHpo00hFZ87SaerxGY4b2qnXIUt83Gl3PyfloDcacXH9RDYOjMDqqb8tTH9IWEyFWiF7GRasmXbWaTJNVDPpvu24hA8lZrlFI0INroDXebtD/4QD9DxnkJB3e76H/NZgnErrCQbMM8vt4KVw99izbddhiZxiGaTPGLNwxNbRutD7JgEK2bP5CH99lFp9mFny+6gRZKQAwabg/KpuLJNQKyZDIrD0QqkeWyhsjUfLQQq4TADBNSRVKf6X5/kLCEn33hTJRwh1t10q33hquo0wJACBeI5f033WeWfzjIRZwmn9Y3Z5WJqH/B/dEZfKRfoqCf0oIn3Nspf5nJSSMMgvsY2cYhjmCadpiT/P1WS2AkjJinkt+X3wsyJDoFZMcW1mfMmmQjKohv1ph304AQG1S0mLJaiXHZaRPkSyVvcNVAMC0SYV6Ga1Ow70hK9IwMSYEW/td92Q8W+Uu0nTHqv/qNc5lt6+a0WlT2Wb2NUSxqPovz5H0vzKpcRKfDa/oEov+95br+m9bNCSWX4usfuo5mY7Z7JiYTpbxKB/YYmcYhmkzWro0no+PMfYtKlaKsFgsrqm4WduVhv7F9IkPtrWpd/u6+O+ROacAqCfY95WRVja2JOToP/kUyVJfvW1/XPbUWd0AlCgAPSWxwVJvVZSFjbTrPd7JIT32OAFZhQrGlNAWGbrl6HO9sliBNvmusQ+9DT767yIkDlzXf7LU17xeXxjk1FldSbm6/hveQ2nX2fWecdUd6wgvttgZhmHaDH6xMwzDtBljNnhqneBA3R/DAEVcxuMYaYNGWbs6NtdLyOQaKjs0+5R4W1FOs0bANOu0tqnQ5CM9pIvcLwAwVJGTOOQUcL1OiLspS9c/q3tsPCMsf9vIS1eMKXMjkebiMJWx7ffBJzjBZ2BYb8vwnLr+62k2mgk5NLW3pukyuRvJ/QIAB2WS+klyp17H5RLzcU252pkm12d7yD1li51hGKbN8LbYbau0pwXN6186V2KgLKE/Pl/FEGwDID7ox1ZzYsdpBhwJwmxyfLBdKzVtAFnqw9JyKeeikDCRM08MM8nX25g1VE6X46ozHqx7l/4TPpavaY1eWkuA9CWrta0ex0SWyU0++PTqapYVz3wGnH2eSZ+eBlnqMvsAOiB7T7n012AzPRifcN9WD7R6W+y8SjtzJMP6z0wkckIYZtEYcK3STqtlu74yDRafkgQpJNxRx2cyTzNf0rTjucrmDEmQypXkGpXNhKIZj6mVUa9zvJamiCyVfdXof5qOfUCUAABdxcbWZJl0pNf1sWhNMmjf4MAAjp41K3WF9rHApv/bt29Hr6MtVh+4Qf9b0VMLudYu69DnuUprkykJmK7/eh3TsXzaQmu86pOQmtX/tGNntazTxkZs+wcGBjDL4xnwdsXwKu3MkQzrPzORaGkSMNe2BmtBSQJWKUVfb1v6TJ+voo8PLE2GD67JInoZUxKwgWIPAKBbK9uK9qhtia+HIQkS+RQny7U1KXHS0VWZ7mDyTKu8kLb40GwU0+Emrd3W+6MkAasUIv+zLY1yFl+465kJifDSt4f0ukxJwAZK0QSlyVqdLFE3Kg2WOtXJoP9VRf9DaIUut2rcg6NiGIZh2oyWWOw+X1sdstKB+mIRNhk+x3YdL0vMbIiP0baP/IkA0C8t9Wn9mwAAlRkLEscJGU/wiTOnhF5qmgC9LPkUyVJ5oRpZU29WxIT4XLNYGyFW6XiIjjERdF8kZKUDjcsahvQ+9V5js/cl7X6H9M4T+i8tddL/Uan/Jpp9JgDEqRtyjtQlPvpPZInYyzJO0aqIJbbYGYZh2oyWzDwNiZ0lTIvP5mSMK6X0zRmWzEv7koX42LPM7gr5kqoj/+RTJEulY9t6AED/zCUAzJEoWYilOBJ60TnQMcmnTpbK09vrltbyaXKmcHlyQoYehaBC+/LavqxzAiaKHz5kLEi/NkA9xjtOae1YMjItakvdn2YFunq5rh5Bavx6i/Q/pDceY5jZrpelY1Y1/X/itXrivLP6oqPXSP8pysYUyaclFWvYbdgW2hPwfQbYYmcYhmkz+MXOMAzTZrRk8JQICYVSifcV5JqPHhM2WjXIkEar3AB6e/fMWAwA6KkMAgBqxSnWsvp2Fz4TfWx1qCy5XwBgIB91p3vk/7aJIEB9arzunglx1ZlCSMezCyZU/5yDj5r+m0ibZJRlkD9raKRv6LFpXxb9zzJAnOX9QO4XwKD/NMnJIL+m7XO1xXbfTOeubvM9H7bYGYZh2owxC3ckgiwyOeCXxQoJ6SHY6qoEDZIGHJvKUEKuWimyVCg5GGBOFGU6TihpA2Lx9vLkuA5ZKnrisKpUHdU6j8/bNcAU2FZTe8cTOYT1Sl0Wma7/pkHoZq6Bj5WfJdSQ8LHUaV+D/o/UB+yFliDMq7dv2R7iRaDtNYP+j1DisFz0R1XIlL/G7oO2MlOGdptQdS0NttgZhmHajJamFGjWFxrLkZMLanLqr+mrmOYL9PEtu/Dxc/taAEb50kKhiVo5x8roLkJ6N76oMshqJEt9SEQexE5pElSURtat92SYXogll9ae8Ybq9/SxzLyseqn/ecdiNDoh19Gr9zDGNOh/sXEiXZA87X+XZyA1pNCQOIws9aFaVKszH/1fEXXbOJbn0VPNEs7IPnaGYZgjmKYt9jQfY4g1S1/KXLx8WNFaNsg61uqE+Natky8M+0IiO8hS6TiwCwCwpzQt3terzS1qJuIhRI7u2wXqVjj51MlS33EgsuBnyenYAFDrOdp4nCx6ETIp53Bi0389ysFLlqb/FdljVefu+EaKZPGRm7a59L+ZY4069H9KR2uPrddPjQ5SLG6KcCGfOlnqO4ai0rOqr8dl0/Tf5UXwOdcQ2GJnGIZpM8bcx27bZ/R9yS9lVSZIMs0y9vUlu6ynkCiWVvnY9bKU7J8slb6BTXFZmnZts/qa/Zqn1TdFutA28qmTpf5ibXpc9s26IEMiMl1uiFU2HqNibL2ykPOrC0jqv2mugOlYPtuztimkZxxStujQ/4ojQZjvsU1tSBvvcL2/8rH+R/eILHUf/Tel9wghi/6zxc4wDNNm8IudYRimzfB2xdhWaVe7B1kmEpj2kZsiDpszhA+ldStDJl00O+HHZ7vencpXom5atxxEOjC93v0saWWbGYRx4eOSojQBukuGBorU7ueQHHXtjJfCSt43n7C6EHfeocSm/0Dr3Uex/svB0ywyQiZLZXHfhJbRy+r6f7Cvrv/6lHwfubbr7Rogtt4jJdyxptm+VMek/wekn7KrICVqE5Ra5VLzwdti51XamSMZ1n9mIpETQngZBLZV2l/fvt24WnYmi5rCveSgg6AJGo6JGj5WrW9ddZtOlq+ua3IEWcIlGe6lWre1rnroVxpZrrPeJn1/TSkYD6RqaQJc14FCIQvylKZPsqVFCsN3hfaxIFT/bTSr/63Qd52s+h8SwKBDIbUdBv2vpuh/s4PUWQIvvPRfJEMhdf0PCdu2vTt8nwFvVwyv0s4cybD+MxOJlqyg1Czxl4y+hhQeZPg6plnHWXy0Lvk+5xcSPtlQRp5j/sCeuAxZ7M1YZ17HtmBMbKTdi7gIhXQpZQpSwMBwZMG4LPbxHMroi4/uwVGmvjOp/6axJV/9D7meh1L/GyD9H+qPN+kWu0+odINYR9lMYZ85s6/dpf/7ZOawPqn/PmMxaWOEvs8+R8UwDMO0GS1JKUBk/ZI27HMk0aEJMjR5qRUTXExtda3pmSbXp9cQR5dI60S1UgqK9a7uc51HluiYEGtMlxufj2HyxfRJ9DuyVPYcJJ97vVW9HZoF5Ep1Sr5nj/YeDlzWYRYL2IWu/zpBvUUHrSjj8oXHqSqkbtc89N8kN61trgi1EKz3UdF/KkP6T5b6boP+T0nRf9O+kDTYbLEzDMO0GWPuY7d9OX38kSZKo1G6z1pKMn7Xl9rHJ0gfV5fcLL5FrzrkdxyMpi1XLT532zabXJ97YG1TE/vIUnl1oO6P7J0hzRqbNaLEEjezYMdYE2qBu/Q+vpeOBRrS9L8higNwRnLY2qK3N0s8e+aepab/riixLONnae8ktU5IL992/qT/L++pR1W95ehOWWhs9H/8PjEMwzBMJvjFzjAM02a0ZPB0rAcqVKgLSmRxRfgM4Lrk6PiUsck1Dtx29iZ+03qoQq6FmnUATq/XzECeTzgZQQOlsfsF9TVUi/l8opJxsFrtlo4z0vSfcLkpfELeYjlS/2uW66W7c/RjmeSGBBxknfgT8h5o0P8Rs/vJp10+98RFiPvWBg2Uxu4X1NdQLeSS+p+2hqrvcdliZxiGaTNaMnga8lX0+aq7JsHU8h2JMlaLxXEc23FNcny2tyKkTa1T2PcGAKA2eSYAYG8tOufeAHlZeiUhE2x8rnPDPsXyJkt911AUCnZUVxQaVhjYBgAY6TkmLmu0YsYRLitcJ1PvVtH/St6c29tm9Zvk+rQhRIdDdMKHWP+7+wAAexH19GgCvU8vMWRwU8cWwp1G6nVQ9J8sdZv+V6YcA52cSaYFttgZhmHajJZOUAohxAc+mqtbKR3xTvr6mae6+3zVfdqil/E53yw9ALUOWSo56VufXJpkLWuTb5sg4zp2q3sjsXxTSJfcSZZKnPJXpkMd71Y6IRAWAprlWlcVK71o0X+X/JBeVzPttcmwbbPJivVf9lRI/1ulElmscN/tpuP46D/53Duk/jcLW+wMwzBtRkvWPPUhxLeubzdGSqQE7Tc7QakVk42yRgntlS7Vno7IUinveglAfS1UUxvS2uZTxyVHJ8ivbLhXdE/Jp0iW+ksD0SITb+41+JLVhEvjBF+/p4/vOqjXZUvKZvg/ZBKPb5tc9bJEraj7B0aj32Spk/6b1kIN6YXb2hfS289yr036n9f0nyz1lwcin/v8XsP7rTri/Qywxc4wDNNmtDSlgI/1Ztpv+2LGkS+1+leK/I20fBgtRhDk8zK0V29XWptMx3T5o0Os+umj0eIDtXIUFVNzJAFL84FniTv2sfYy+SkN8egU/UIWDFnqb8jFOoC6H3LnSGsW7DiUtOIaq/pPUTFFWj7PsghNVt2ztSkk2sanDa5nZ9qI1P9S5GvX01ib6oVExYT0crOMy/lExRB69AtZ6ib9f2OkgEHPZ4AtdoZhmDaDX+wMwzBthrcrxrZKeyvCHdNcGiOGcMeqbDp1THxcDj7dtpAB0bQudIhc9X893LEmp1aHhKe58O2SqvJCBukaBZrCvaJttrBG6n4C9VDIaS1aOzULNv1Xwx19JsWEDDrT9ooh3FHXf/14rja4ytrquvb5yPWdHAQ0hjuS/rcq3NGnDXrZkP0Nz4xD/22o+k+hkNMnFVBstSuGV2lnjmRY/5mJRE4I4TUOZlulfbuySntWazYV5eumr7ATYjU0U8bnK9/MRB91f36/HDySg0bFXZsA1MO9QgaRXITcG5+JK1kG6RqQVtobimVCKzE9s3Un3nr8m1JXaB8LWqX/JlLLOvRfl2EiS3ifT9m0Oi5cdWgFJbLUSf9N4b5ZwhFNx7S1LeT5d8mxHbuhjEP/Nw+MYt/gAM44fm7qM+DtitFXaafvweDgYLyNX+zZj514sR+IrmmtEt3Q4uA+AEClo979Dzm2jfH8YldH/6n7uW9fdF08bZGW0ir9N8Evdu3FPiT1fzQqRfo/KvXfJKfdXuxG/R8c9X4GMoc7kkIfd9xxWUUwTCYGBwfR29t72NsAsP4zh4e0Z8DbFaNTq9Xw3HPPYfHixdiyZUtLu8bUzZ0IcidSWye63J6eHgwODmL27NnI5w9vQBfrP8sda5kmuUIIr2cgs8Wez+cxZ84cAMCUKVPGxOc5keROpLZOZLmH21InWP9Z7qGSqcv1eQY4jp1hGKbN4Bc7wzBMm9HUi71cLuOGG25IRAu0gokkdyK1leW2lol2zix37OSOt7ZmHjxlGIZhxifsimEYhmkz+MXOMAzTZvCLnWEYps3gFzvDMEybwS92hmGYNoNf7AzDMG0Gv9gZhmHaDH6xMwzDtBn/Hxvoe3BLdPC3AAAAAElFTkSuQmCC", "text/plain": [ "
" ] @@ -538,7 +547,7 @@ { "data": { "text/plain": [ - "0.5785552396972465" + "4.684709110728037" ] }, "execution_count": 21, @@ -561,7 +570,7 @@ { "data": { "text/plain": [ - "0.5812977845446157" + "4.746888516975277" ] }, "execution_count": 22, @@ -625,7 +634,7 @@ "name": "stdout", "output_type": "stream", "text": [ - "Plotting at 16 steps: 5, 8, 11, 15, 18, ...\n" + "Plotting at 16 steps: 5, 350, 695, 1041, 1386, ...\n" ] } ], @@ -646,7 +655,7 @@ "name": "stderr", "output_type": "stream", "text": [ - "100%|█████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 16/16 [00:34<00:00, 2.16s/it]\n" + "100%|█████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 16/16 [18:07<00:00, 67.94s/it]\n" ] } ], @@ -680,7 +689,7 @@ "outputs": [ { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAVwAAAC+CAYAAAB57bcRAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAArlUlEQVR4nO3deVhTV/oH8G9IyEJYXNgSBcQFlIAbqODOoKhUW/WxtXa0dlrbOtVWH8aZ1rEt6thqbetYx2rVOrjNVDuD+rNqW9GCu6W1MLIo4gpiAAEhbEkgOb8/IrdGCCQQSKDv53nymNx7cu8bNC/H95x7Lo8xxkAIIaTNOdg6AEII+a2ghEsIIe2EEi4hhLQTSriEENJOKOESQkg7oYRLCCHthBIuIYS0E0q4hBDSTgS2DqCj0uv1uH//PlxcXMDj8WwdDiHEhhhjqKiogFwuh4OD6X4sJdwWun//Pnx8fGwdBiHEjuTl5aFnz54m91PCbSEXFxcAhh+wq6urjaMhhNiSSqWCj48PlxdMoYTbQvVlBFdXV0q4hBAAaLa8SINmhBDSTqiH20pVVVXg8/kNtvP5fIjFYqN2pjg4OEAikbSobXV1NUwt+Mbj8eDk5NSitjU1NdDr9SbjkEqlLWqrVquh0+ms0tbJyYnrUWg0GtTV1VmlrUQi4QY+tFotamtrrdJWLBZz/1YsaVtbWwutVmuyrUgkgkAgsLhtXV0dNBqNybZCoRCOjo4Wt9XpdFCr1SbbOjo6QigUWtxWr9ejpqbGKm0FAgFEIhEAw4BXdXV1q9o29Z01wkiLlJeXMwAmHzExMUbtnZycTLYdN26cUVt3d3eTbcPCwoza+vn5mWwbFBRk1DYoKMhkWz8/P6O2YWFhJtu6u7sbtR03bpzJtk5OTkZtY2Jimvy5PW7WrFlNtq2srOTazp8/v8m2RUVFXNs33nijyba3b9/m2i5btqzJthkZGVzbuLi4JtumpKRwbdevX99k26SkJK7t5s2bm2x79OhRrm18fHyTbb/++muu7ddff91k2/j4eK7t0aNHm2y7efNmrm1SUlKTbdevX8+1TUlJabJtXFwc1zYjI6PJtsuWLePa3r59u8m2b7zxBte2qKioybbz58/n2lZWVjbZtry8nDWFSgqEENJOeIzRAuQtoVKp4Obmhvv37zc6aEYlhcbbUkmhYUmhSlMHJyHfaMCFSgoN29pzSUGlUkEul6O8vLzJQXRKuC1Un3Cb+wET8iS9nuFKfjm+zyzAicwC3HxQBYXcFfPC/fD0YDmchJ1/aIUxhppaHSrUdVDV1EKlroVKXQcwgMcDHHi8Rw9DZ6CxP+vbtOa6I1PnevzYDg6/bjPVvrKiAnLP7pRw2wolXGKJWp0eP94qxYmsApzILESBqvFenYtYgGdDfTA33Be9PZzbOUrLlVRqUFSheZQ061ChruWeq2pqDQlV/Sih1jza/2hfnb7zpB69php5G59rNh90/l+lhNhItbYOZ64X40RmAU5dK0J5za/lBqmQj/GBnohWeGGob1d8m6HEvku5yC2txj/P38Y/z9/G6L7umBvuhwkDPCHg28dwC2MMN4oqcSKrECcyC/C/e+WtOp4DD3CVOMJV7AhnkQB8Bx70jEHPDOdiDI9eP/4cYGDQ6w1tdK3oMzLGjdZCX398vfG59E/E0ZrfEzbv4W7ZsgUff/wxlEolFAoFNm7ciDFjxjTaNjk5GZGRkQ22X716Ff379wcAjB8/HqdPn27QJiYmBseOHQMArFy5EqtWrTLa7+XlhYKCArPjph6uben1DEUVGuSWViO3tBqlVRpIhAI4i/iQCgWQigwPZxGfey4VGr7QbamsWouTV4vwfWYBzuY8gLr219p2d6kQEwZ4YVKwF0b2cYfY0Xg6oV7PcCbnAfZduotT14pQ/830dhXjhRG+eH64DzxdxGhvej1Dal4ZTmQVIDGzELeKfx1j4PGA7lIRXCUCuIgd4SoWPEqgAriKHbnnLmJHuEoM2x5//mTtuqMwStCMobxcBS/3bvbdwz1w4ACWLl2KLVu2YNSoUdi2bRumTJmCrKws+Pr6mnxfdna20Yfy8PDgnh88eNBo0KCkpASDBg3Cs88+a3QMhUKBkydPcq8bm0tLbKtCXYvc0mrkldYg71FizXto+PPewxpo60wP1JkiceQ3SMTOjydnoQAiRwej3o6eATo9475kOmZ4rnu0z9CGoUClxk93HkL3WBeoZ1cJJim8MUnhjVC/rk0mfAcHHsYHemJ8oCfySqvx75RcHPgpDwUqNTYkXsemUzmYHOyNeeF+GO7frU0TlaZOh4s3S3AiqxCJWYV4UPHrgJmQ74BRfbsjWuGNqAGeNvklYGs8Hg98HsCH4e/gyV+eJt9nyx7uiBEjMHToUGzdupXbNmDAAEyfPh1r165t0L6+h/vw4UN06dLFrHNs3LgR77//PpRKJTcCvnLlShw+fBhpaWktjp16uNalqdPh4C/5OHejGHml1cgrrcbDatMj/gDAd+ChRxcJfLpJ4OEsQk2tDlUaHaq0dajS1KFKo0OlxvC8PeuF/b1dMEnhjWiFF4Jkrq1KjJo6Hb5NL8Cei3fwS24Ztz3Ayxnzwv0wY2hPOIus02+qUNciOfsBvs8sQHL2A1Rqfp3J4SISILK/oQQyLsADLmJHq5yzszA3H9ish6vVanH58mW88847Rtujo6Nx4cKFJt87ZMgQqNVqBAUF4d133220zFBv586deP75542mGwFATk4O5HI5RCIRRowYgQ8//BC9e/du+QciLaKu1eHAT3n44vRNKMsbDiR1kwrh080Jvt2c4NNVAt/6592cIHMTm1XbZIxBU6c3TsLaOlRq6lCt0aFKU8cl5kptHTS1ejjweOA7PBqRdjAeua7fx3vsuQOPB4mQj9F93eHXXdpsTOYSCfiYPqQHpg/pgYz8cvzrx7s4nHof1wsr8d7/ZWLdt9cwc2hPzA33Q6B30wun6PQMJVUaFKk0KKpQo0ilQeGj53dLqvHj7RLU6n79xeTpIsLEIC9EK7wR0bs7hAL7qCN3ZDZLuMXFxdDpdPDy8jLa3lQtVSaTYfv27QgNDYVGo8HevXsRFRWF5ORkjB07tkH7lJQUZGRkYOfOnUbbR4wYgT179iAgIACFhYVYs2YNRo4ciczMTHTv3r3Rc2s0GqN5iCqVytKPTB5To9XhXz/exbYzt7j/rnq7ijE33Bf9vFy4pGqN3huPx4PYkQ+xIx/d7X/g36TgHm5YO3Mg3pkyAAd/uYe9l+7i1oMq7L10F3sv3cVw/26YNbQndIyhUKVGUYUGRY/+LFSpUVypNSp3NKa3uxTRCm9MUnhhUM8ucGjjmvdvjc1KCvfv30ePHj1w4cIFREREcNs/+OAD7N27F9euXTPrONOmTQOPx8ORI0ca7Hv99ddx4cIFpKenN3mMqqoq9OnTB3/5y18QGxvbaJvGBtoAUEnBQlWaOuy9dBdfnr2F4kpDrb1HFwn+OL4Png3rCZGAaunmYozhws0S7L14F4lXC5tNpsCvg1xeriJ4uYrh6SKCp6sYXq4ijPDvjr6eHfg3kg3ZfUnB3d0dfD6/QW+2qKioQa+3KeHh4di3b1+D7dXV1di/fz9Wr17d7DGkUilCQkKQk5Njss3y5cuNknH9+pfEPCp1LfZcuIOd525ztVnfbk5YFNkHM4b0pP+utgCPx8Oovu4Y1dcdyvIafJWSh/M3iuEmcYSXqwgeLoZE6vnYn+7OQruZYvZbZLOEKxQKERoaisTERMyYMYPbnpiYiGeeecbs46SmpkImkzXY/vXXX0Oj0WDu3LnNHkOj0eDq1asmp6MBhksi6y/vI+Yrr67FP8/fRvz524YriQD4u0uxOLIvnhkspy+/lcjcJIidGIDYiQG2DoU0wabTwmJjYzFv3jyEhYUhIiIC27dvR25uLhYuXAjA0KvMz8/Hnj17ABhmHPTq1QsKhQJarRb79u1DQkICEhISGhx7586dmD59eqM12WXLlmHatGnw9fVFUVER1qxZA5VKhfnz57ftB/4NKa3SYue5W9h94S432t3X0xlv/q4vpg6Ut/l8WELskU0T7uzZs1FSUoLVq1dDqVQiODgYx48fh5+fHwBAqVQiNzeXa6/VarFs2TLk5+dDIpFAoVDg2LFjiImJMTru9evXce7cOZw4caLR8967dw9z5sxBcXExPDw8EB4ejkuXLnHnJa3z5dlb2JB4HdVaw8Iz/b1d8Obv+mFKsDcNwpDfNJtfadZR0Tzcxt18UImoTw1X+gX3cMWbv+uHiQO8KNGSTs3uB81I53TsihIAMKafO/a8PLxDXrZJSFuhEQtiVfUJ9+lBckq2hDyBEi6xmpzCCmQXVsCRz0N0kLetwyHE7lDCJVZzLN3Qux3bzwNuTnStPSFPooRLrKa+nPDUwIbzogkhlHCJlVwvrEBOUSWEfAdMCDL/SkFCfkso4RKrOPqodzs2wAOutHQfIY2ihEtajTGGY1fuAwCmUjmBEJMo4ZJWyy6swM0HVRAKHBA1wNPW4RBityjhklarHywbT3cCIKRJlHBJqxjKCTQ7gRBzUMIlrXJVWYFbxVUQCRwQNYBmJxDSFEq4pFWOpRsGyyIDPa12M0NCOiuLviHl5eU4dOgQzp49izt37qC6uhoeHh4YMmQIJk2ahJEjR7ZVnMQOUTmBEMuY1cNVKpV49dVXIZPJsHr1alRVVWHw4MGIiopCz549kZSUhIkTJyIoKAgHDhxo65iJnci8r8KdkmqIHR3wu/40O4GQ5pjVwx00aBBefPFFpKSkIDg4uNE2NTU1OHz4MDZs2IC8vDwsW7bMqoES+1O/dsLv+ntCSuUEQppl1rckMzMTHh4eTbaRSCSYM2cO5syZgwcPHlglOGK/GGM4+uhih6dC5DaOhpCOwaySgoeHBw4fPgydTmfWQZtLzqTjS88vR15pDSSOfET2p79vQsxh9iyFWbNmoUePHnj77bdx7dq1toyJdAD1g2W/G+AJJyGVEwgxh9kJNzc3F2+++SYOHToEhUKB0aNHIz4+HlVVVW0ZH7FDhnKCIeFODaHZCYSYy+yEK5fLsWLFCly/fh0//PAD+vTpg7feegsymQwLFizAxYsX2zJOYkf+d68c+WU1cBLyMT6QZicQYq4WXfgwbtw47N69G0qlEhs2bMDVq1cxevRoKBQKa8dH7FD9ymBRA7wgEfJtHA0hHUerim/Ozs6IjIzEnTt3cO3aNVy/ft1acRE7ZXSxA5UTCLFIi3q41dXV2L17N8aNG4eAgAAcOHAAsbGxuHPnjpXDI/YmNa8M98vVkAr5GB9IsxMIsYRFPdzz58/jn//8J/7zn/+grq4OM2fOxMmTJxEZGdlW8RE7U9+7nRDkBbEjlRMIsYTZCTcgIAA3b97EkCFD8NFHH+GFF16Am5tbW8ZGWkBda5gr3RbJUK9nOJ5O5QRCWsrshDt58mS88sorGDRoUFvGQ1qhSlOHyE+S0dVJiCNvjoJIYN2km5r3EMpyNZxFAowNoHICIZYyu4a7adMmLtnW1dXh5MmT2LZtGyoqKgAA9+/fR2VlZdtESczyv3tlKKrQILuwAgmX861+/Pq5txOpnEBIi1g8aHb37l2EhITgmWeewaJFi7h1E9avX9+iBWu2bNkCf39/iMVihIaG4uzZsybbJicng8fjNXg8fuXbrl27Gm2jVqtbfN6OIjNfxT3/POkGanV6qx2bygmEtJ7FCXfJkiUICwvDw4cPIZFIuO0zZszAqVOnLDrWgQMHsHTpUqxYsQKpqakYM2YMpkyZgtzc3Cbfl52dDaVSyT369etntN/V1dVov1KphFgsbvV57V16fjn3PL+sBod+sV4v93LuQxSqNHARCTAmwN1qxyXkt8TihHvu3Dm8++67EAqFRtv9/PyQn2/ZF3zDhg145ZVXsGDBAgwYMAAbN26Ej48Ptm7d2uT7PD094e3tzT34fOP/3vJ4PKP93t7eVjmvvcu4b0i4Y/oZEuJmK/Zy62cnTFR4Wb02TMhvhcUJV6/XN7pq2L179+Di4mL2cbRaLS5fvozo6Gij7dHR0bhw4UKT7x0yZAhkMhmioqKQlJTUYH9lZSX8/PzQs2dPTJ06Fampqa0+r0ajgUqlMnrYk0pNHW4XG9a1+HBGCLpLhcgtrcb/pd1v9bF1j5UTptKdHQhpMYsT7sSJE7Fx40buNY/HQ2VlJeLi4hATE2P2cYqLi6HT6eDlZXzjQS8vLxQUFDT6HplMhu3btyMhIQEHDx5EYGAgoqKicObMGa5N//79sWvXLhw5cgRfffUVxGIxRo0ahZycnBafFwDWrl0LNzc37uHj42P2Z20PWfdVYAyQuYnh080Jr43tDcBQy61rZS/35zulKKrQwEUswOi+NDuBkJay+NLev//974iMjERQUBDUajVeeOEF5OTkwN3dHV999ZXFAfB4PKPXjLEG2+oFBgYiMDCQex0REYG8vDx88sknGDt2LAAgPDwc4eHhXJtRo0Zh6NCh+Mc//oFNmza16LwAsHz5csTGxnKvVSqVXSXdjEf12+AehrnRc8P98MXpm7hdXIVvrtzHjCE9W3zs+js7TFJ4Qyig+44S0lIWf3vkcjnS0tLw5z//Ga+//jqGDBmCdevWITU1FZ6e5q8c5e7uDj6f36BXWVRU1KD32ZTw8HCu99oYBwcHDBs2jGvT0vOKRCK4uroaPexJff02WG5IuFKRAAvGGHq5//jhBnR61qLjGsoJhp8VlRMIaZ0WdVckEgn+8Ic/YPPmzdiyZQsWLFhgNGPBHEKhEKGhoUhMTDTanpiYaNHdf1NTUyGTmU4EjDGkpaVxbax1Xnvzaw/3118E80f2QhcnR9x6UMX1Ui2VcrsUxZUauEkcMaovzU4gpDXMKilcvHgRERERZh2wqqoKd+7cMWupxtjYWMybNw9hYWGIiIjA9u3bkZubi4ULFwIw/Dc+Pz8fe/bsAQBs3LgRvXr1gkKhgFarxb59+5CQkICEhATumKtWrUJ4eDj69esHlUqFTZs2IS0tDZ9//rnZ5+1oarQ63CgyXHRSX1IAAGeRAK+M8senidfxj1M5mBoig4OD6bJJY46lGwbdJiu84cincgIhrWFWwn3xxRfRq1cvvPrqq4iJiYGzs3ODNllZWdi3bx/i4+Oxfv16sxLu7NmzUVJSgtWrV0OpVCI4OBjHjx+Hn58fAMPt2R+fG6vVarFs2TLk5+dDIpFAoVDg2LFjRoN1ZWVleO2111BQUAA3NzcMGTIEZ86cwfDhw80+b0dztUAFPQM8XETwchUb7Zs/qhd2nL2FnKJKfJtRgKcsKAvU6fT4LsNQTrDkfYSQxvEYY80W92pra7Ft2zZs3rwZN2/eREBAAORyOcRiMR4+fIhr166hqqoKM2fOxPLly03eSr0zUalUcHNzQ3l5uc3ruXsu3sH7/5eJyEAPxP9heIP9G09ex8aTOejv7YLjb40xu5d74UYxXvjyR3R1ckTKignUwyXEBHPzgVnfIEdHRyxevBjXrl3Djz/+iNdeew3BwcHo0aMHxo8fj23btiE/Px//+te/fhPJ1t48OUPhSX8Y6Q8XkQDXCipwIqvQ7OMefVT3nRxM5QRCrMHiaWFDhw7F0KFD2yIW0kIZj9ZQMJVw3Zwc8dKoXvjHDzew6VQOJim8mpwCBzxRTgiRWzdgQn6jqNvSwalrdbheaFixzVTCBYBXRvtDKuQjS6nCyatFzR730q1SlFZp0U0qRHjvblaLl5DfMkq4Hdz1wgrU6Rm6OjlC7iY22a6LkxDzR/YCAGw6lYPmSvdHH90ocnKwNwRUTiDEKuib1ME9Xk5orkywYExvOAn5SM8vR1K26V5urU6P7zIfXexASzESYjWUcDu49GYGzB7XTSrEvHDD1LfPTt0w2cu9cLMEZdW1cHcWYrg/lRMIsRZKuB1c5hOX9Dbn1bG9IXZ0wP/yynAmp7jRNseonEBImzD72xQTE4Py8l8XuP7ggw9QVlbGvS4pKUFQUJBVgyNNq9XpcU1pGDALMaOHCwDuziLMHfGol3vyeoNerrZOj+8zDVPHaHYCIdZldsL9/vvvodFouNcfffQRSktLudd1dXXIzs62bnSkSdcLK6DV6eEiFsCnm/lrWbw2rjdEAgf8kluG8zdKjPadv1mM8ppauDuLqJxAiJWZnXCf7AmZcYEaaWP19zALljc/YPY4TxcxXhjhCwD47JRxL7f+zg4xId7gW7juAiGkaVSg68C4JRl7WH5p8cJxfSAUOOCnOw9x6ZbhfyqGckL9xQ40O4EQazM74dbf/fbJbcR2LJmh8CQvVzGeH2ZYQP2zU9cBAOduPECFug6eLiKE9aJyAiHWZvalvYwxvPTSSxCJRAAAtVqNhQsXQiqVAoBRfZe0vTqdHleVTV/S25w/ju+D/Sl5uHSrFD/eKsFRrpwgo3ICIW3A7IQ7f/58o9dz585t0ObFF19sfUTELLeKq6Cu1UMq5MO/u7RFx5C5SfBsWE/868dcfJp4HVfvGxI4LcVISNswO+HGx8e3ZRzEQun3DOUEhdzN4kXFH/fH8X1w4Kc8pNw21HG9XEUI9e1qlRgJIcYsWi3s7t27OHHiBGprazF+/Hiad2tD9QNmihYMmD2uZ1cnzArtif0/5QEwlBNak8AJIaaZnXDPnDmDmJgYVFdXG94oEGD37t2YM2dOmwVHTHt8SlhrLYrsi/9cvgedntGNIglpQ2bPUnjvvfcQGRmJe/fuoaSkBC+//DL+8pe/tGVsxAS9nnGX9Ib0bH3C9enmhM1zhmDV0woMpXICIW3G7B5ueno6zpw5A7nccLnnp59+ih07duDhw4fo2pW+pO3pdkkVqrQ6iB0d0Nu9ZQNmT5pC824JaXNm93DLysrg6enJvZZKpXBycjJaT4G0j/pb6gyQudLiMoR0IBYNmmVlZaGgoIB7zRjD1atXUVFRwW0bOHCg9aIjjcp8NH3L3AVrCCH2waKEGxUV1WANhalTp4LH44ExBh6PB51OZ9UASUP1U8KsMWBGCGk/Zifc27dvt2UcxEyMMatNCSOEtC+zE66fn19bxkHMlFdagwp1HYR8B/TzdLF1OIQQC5g94hIVFYWDBw+a3F9cXIzevXtbJShiWv2CNf1lLhAKaMCMkI7E7G9sUlISnnvuOcTFxTW6X6fT4e7du1YLjDSOKydQ/ZaQDseiLtLWrVvx2WefYcaMGaisrGyrmEgTMvJbvgYuIcS2LEq4zzzzDC5evIisrCxERETg1q1brQ5gy5Yt8Pf3h1gsRmhoKM6ePWuybXJyMrcu7+OPa9eucW127NiBMWPGoGvXrujatSsmTJiAlJQUo+OsXLmywTG8vb1b/VnaGmOMS7g0JYyQjsfiIuCAAQOQkpICHx8fDBs2DCdPnmzxyQ8cOIClS5dixYoVSE1NxZgxYzBlyhTk5uY2+b7s7GwolUru0a9fP25fcnIy5syZg6SkJFy8eBG+vr6Ijo5Gfn6+0TEUCoXRMdLT01v8OdrL/XI1HlbXQuDAQ4AXDZgR0uEwM/F4PFZYWMi91uv17O2332aOjo5sw4YNrKCggDk4OJh7OMYYY8OHD2cLFy402ta/f3/2zjvvNNo+KSmJAWAPHz40+xx1dXXMxcWF7d69m9sWFxfHBg0aZFGsTyovL2cAWHl5eauOY4nvMpTM7+2jbPLGM+12TkJI88zNBxbdYufJ1+vWrcPevXvx3nvvYcGCBRYleq1Wi8uXLyM6Otpoe3R0NC5cuNDke4cMGQKZTIaoqCgkJSU12ba6uhq1tbXo1s34ljE5OTmQy+Xw9/fH888/b5XySFv7tZxA9VtCOqIW37W33uzZs3Hu3DmL/0teXFwMnU4HLy8vo+1eXl5Glw8/TiaTYfv27UhISMDBgwcRGBiIqKgonDlzxuR53nnnHfTo0QMTJkzgto0YMQJ79uzB999/jx07dqCgoAAjR45ESUmJyeNoNBqoVCqjR3vLaMU9zAghtmf2hQ9JSUkNeon1Bg8ejMuXL+PYsWMWB/Bkz5k9ukS4MYGBgQgMDOReR0REIC8vD5988gnGjh3boP369evx1VdfITk5GWKxmNs+ZcoU7nlISAgiIiLQp08f7N69G7GxsY2ee+3atVi1apVFn83aMh6toUBTwgjpmMzu4Y4bNw4Cgen83L17d4vuaebu7g4+n9+gN1tUVNSg19uU8PBw5OTkNNj+ySef4MMPP8SJEyeaXVBHKpUiJCSk0ePUW758OcrLy7lHXl6e2TFaQ6FKjQcVGjjwgCAZlRQI6YhsdqmSUChEaGgoEhMTjbYnJiZi5MiRZh8nNTUVMpnxWq4ff/wx/va3v+G7775DWFhYs8fQaDS4evVqg+M8TiQSwdXV1ejRnurLCX09nSER8tv13IQQ67BotTBri42Nxbx58xAWFoaIiAhs374dubm5WLhwIQBDrzI/Px979uwBAGzcuBG9evWCQqGAVqvFvn37kJCQgISEBO6Y69evx3vvvYd///vf6NWrF9eDdnZ2hrOzMwBg2bJlmDZtGnx9fVFUVIQ1a9ZApVI1uDOxPcmw4i11CCG2YdOEO3v2bJSUlGD16tVQKpUIDg7G8ePHuYVylEql0ZxcrVaLZcuWIT8/HxKJBAqFAseOHUNMTAzXZsuWLdBqtZg1a5bRueLi4rBy5UoAwL179zBnzhwUFxfDw8MD4eHhuHTpkl0v0JNOA2aEdHg8Zmr6AWmSSqWCm5sbysvL26W8ELH2FJTlanz9egSG+zc+eEkIsQ1z8wEtN9UBFFdqoCxXAwCC5DRgRkhHRQm3A6gfMOvtLoWzyKZVIEJIK1DC7QDq72FG9VtCOjZKuB0ALclISOdACbcD4GYo0JQwQjo0Srh2rqxai3sPawAACiopENKhUcK1c/X1W99uTnCTONo4GkJIa1DCtXNUvyWk86CEa+foCjNCOg9KuHaOmxJGA2aEdHiUcO1YhboWt4urAFAPl5DOgBKuHavv3fboIkE3qdDG0RBCWosSrh2rHzBT0PoJhHQKlHDtGF3SS0jnQgnXjqXTlDBCOhVKuHaqWluHmw8qAVAPl5DOghKunbqqVIExwNNFBE8XcfNvIITYPUq4dir9Hl3wQEhnQwnXTmXQgBkhnQ4lXDvFraFAU8II6TQo4dohda0OOUU0YEZIZ0MJ1w5dK6iATs/QXSqEzI0GzAjpLCjh2iFlWQ3Ejg5Q9HADj8ezdTiEECuhW8DaoSkhMkwM8oJKXWfrUAghVkQ9XDsl4DvQgjWEdDKUcAkhpJ1QwiWEkHZCNdwWYowBAFQqlY0jIYTYWn0eqM8LplDCbaGSkhIAgI+Pj40jIYTYi4qKCri5mZ47Twm3hbp16wYAyM3NbfIHbI9UKhV8fHyQl5cHV9eOdyVbR46fYreNto6dMYaKigrI5fIm21HCbSEHB0P5283NrcP946vn6uraYWMHOnb8FLtttGXs5nS8aNCMEELaCSVcQghpJ5RwW0gkEiEuLg4ikcjWoVisI8cOdOz4KXbbsJfYeay5eQyEEEKsgnq4hBDSTijhEkJIO6GESwgh7YQSbgtt2bIF/v7+EIvFCA0NxdmzZ20dUgNnzpzBtGnTIJfLwePxcPjwYaP9jDGsXLkScrkcEokE48ePR2Zmpm2CfcLatWsxbNgwuLi4wNPTE9OnT0d2drZRG3uNf+vWrRg4cCA35zMiIgLffvstt99e427M2rVrwePxsHTpUm6bvca/cuVK8Hg8o4e3tze33y7iZsRi+/fvZ46OjmzHjh0sKyuLLVmyhEmlUnb37l1bh2bk+PHjbMWKFSwhIYEBYIcOHTLav27dOubi4sISEhJYeno6mz17NpPJZEylUtkm4MdMmjSJxcfHs4yMDJaWlsaeeuop5uvryyorK7k29hr/kSNH2LFjx1h2djbLzs5mf/3rX5mjoyPLyMiw67iflJKSwnr16sUGDhzIlixZwm231/jj4uKYQqFgSqWSexQVFXH77SFuSrgtMHz4cLZw4UKjbf3792fvvPOOjSJq3pMJV6/XM29vb7Zu3Tpum1qtZm5ubuyLL76wQYRNKyoqYgDY6dOnGWMdL/6uXbuyL7/8ssPEXVFRwfr168cSExPZuHHjuIRrz/HHxcWxQYMGNbrPXuKmkoKFtFotLl++jOjoaKPt0dHRuHDhgo2istzt27dRUFBg9DlEIhHGjRtnl5+jvNxwF+P6NSw6Svw6nQ779+9HVVUVIiIiOkzcixYtwlNPPYUJEyYYbbf3+HNyciCXy+Hv74/nn38et27dAmA/cdNaChYqLi6GTqeDl5eX0XYvLy8UFBTYKCrL1cfa2Oe4e/euLUIyiTGG2NhYjB49GsHBwQDsP/709HRERERArVbD2dkZhw4dQlBQEPfltte4AWD//v345Zdf8NNPPzXYZ88/9xEjRmDPnj0ICAhAYWEh1qxZg5EjRyIzM9Nu4qaE20JP3tyRMdYhb/jYET7H4sWLceXKFZw7d67BPnuNPzAwEGlpaSgrK0NCQgLmz5+P06dPc/vtNe68vDwsWbIEJ06cgFhs+o7R9hj/lClTuOchISGIiIhAnz59sHv3boSHhwOwfdxUUrCQu7s7+Hx+g95sUVFRg9+e9qx+9NbeP8ebb76JI0eOICkpCT179uS223v8QqEQffv2RVhYGNauXYtBgwbhs88+s/u4L1++jKKiIoSGhkIgEEAgEOD06dPYtGkTBAIBF6O9xv84qVSKkJAQ5OTk2M3PnRKuhYRCIUJDQ5GYmGi0PTExESNHjrRRVJbz9/eHt7e30efQarU4ffq0XXwOxhgWL16MgwcP4ocffoC/v7/RfnuP/0mMMWg0GruPOyoqCunp6UhLS+MeYWFh+P3vf4+0tDT07t3bruN/nEajwdWrVyGTyezn595uw3OdSP20sJ07d7KsrCy2dOlSJpVK2Z07d2wdmpGKigqWmprKUlNTGQC2YcMGlpqayk1fW7duHXNzc2MHDx5k6enpbM6cOXYxvYcxxv74xz8yNzc3lpycbDTNp7q6mmtjr/EvX76cnTlzht2+fZtduXKF/fWvf2UODg7sxIkTdh23KY/PUmDMfuP/05/+xJKTk9mtW7fYpUuX2NSpU5mLiwv3vbSHuCnhttDnn3/O/Pz8mFAoZEOHDuWmK9mTpKQkBqDBY/78+Ywxw1SZuLg45u3tzUQiERs7dixLT0+3bdCPNBY3ABYfH8+1sdf4X375Ze7fhoeHB4uKiuKSLWP2G7cpTyZce42/fl6to6Mjk8vlbObMmSwzM5Pbbw9x02phhBDSTqiGSwgh7YQSLiGEtBNKuIQQ0k4o4RJCSDuhhEsIIe2EEi4hhLQTSriEENJOKOESQkg7oYRLSDvZtWsXunTpYuswiA1RwiWdyksvvdTgvlY8Hg+TJ0+2dWiYPXs2rl+/buswiA3Rerik05k8eTLi4+ONtolEojY7n1arhVAobLadRCKBRCJpsziI/aMeLul0RCIRvL29jR5du3YFYFiA+ssvv8SMGTPg5OSEfv364ciRI0bvz8rKQkxMDJydneHl5YV58+ahuLiY2z9+/HgsXrwYsbGxcHd3x8SJEwEAR44cQb9+/SCRSBAZGYndu3eDx+OhrKwMQOMlhW+++QahoaEQi8Xo3bs3Vq1ahbq6Om7/ypUr4evrC5FIBLlcjrfeeqsNfmKkvVDCJb85q1atwnPPPYcrV64gJiYGv//971FaWgoAUCqVGDduHAYPHoyff/4Z3333HQoLC/Hcc88ZHWP37t0QCAQ4f/48tm3bhjt37mDWrFmYPn060tLS8Prrr2PFihVNxvH9999j7ty5eOutt5CVlYVt27Zh165d+OCDDwAA//3vf/H3v/8d27ZtQ05ODg4fPoyQkJC2+aGQ9tGua5MR0sbmz5/P+Hw+k0qlRo/Vq1czxgzLPr777rtc+8rKSsbj8di3337LGGPsvffeY9HR0UbHzMvLYwBYdnY2Y8ywXOHgwYON2rz99tssODjYaNuKFSsYAPbw4UPGGGPx8fHMzc2N2z9mzBj24YcfGr1n7969TCaTMcYY+/TTT1lAQADTarUt/GkQe0M1XNLpREZGYuvWrUbb6u/2CwADBw7knkulUri4uKCoqAiA4RYzSUlJcHZ2bnDcmzdvIiAgAAAQFhZmtC87OxvDhg0z2jZ8+PAm47x8+TJ++uknrkcLGO7yq1arUV1djWeffRYbN25E7969MXnyZMTExGDatGkQCOhr21HR3xzpdKRSKfr27Wtyv6Ojo9FrHo8HvV4PANDr9Zg2bRo++uijBu+TyWRG53gca+RmhKyZpab1ej1WrVqFmTNnNtgnFovh4+OD7OxsJCYm4uTJk3jjjTfw8ccf4/Tp0w0+A+kYKOES8pihQ4ciISEBvXr1sqgn2b9/fxw/ftxo288//9zsubKzs5v85SCRSPD000/j6aefxqJFi9C/f3+kp6dj6NChZsdG7AclXNLpaDSaBndnFQgEcHd3b/a9ixYtwo4dOzBnzhz8+c9/hru7O27cuIH9+/djx44d4PP5jb7v9ddfx4YNG/D222/jlVdeQVpaGnbt2gWg4a25673//vuYOnUqfHx88Oyzz8LBwQFXrlxBeno61qxZg127dkGn02HEiBFwcnLC3r17IZFI4OfnZ9kPhNgNmqVAOp3vvvsOMpnM6DF69Giz3iuXy3H+/HnodDpMmjQJwcHBWLJkCdzc3ODgYPrr4u/vj//+9784ePAgBg4ciK1bt3KzFEzNAZ40aRKOHj2KxMREDBs2DOHh4diwYQOXULt06YIdO3Zg1KhRGDhwIE6dOoVvvvkG3bt3t/AnQuwF3dOMkDbywQcf4IsvvkBeXp6tQyF2gkoKhFjJli1bMGzYMHTv3h3nz5/Hxx9/jMWLF9s6LGJHKOESYiU5OTlYs2YNSktL4evriz/96U9Yvny5rcMidoRKCoQQ0k5o0IwQQtoJJVxCCGknlHAJIaSdUMIlhJB2QgmXEELaCSVcQghpJ5RwCSGknVDCJYSQdkIJlxBC2sn/A64Nypj2M0WwAAAAAElFTkSuQmCC", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAVMAAAC9CAYAAAAOcp4yAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAAhi0lEQVR4nO3deXRTZf4/8PfN1qZpGmhLaSsVWwERStkFHL6WDrJMRUc5qIOoqIxHR0A8qMMgo1WPTlGP6PhzRAZnCuOcGXRkObghRQsIsmihQ6FQGGmhLKUsbdI16/P7I82loQtpe9M04f06J6fJzSfJ5wnl0+e597nPlYQQAkRE1CmqQCdARBQKWEyJiBTAYkpEpAAWUyIiBbCYEhEpgMWUiEgBLKZERArQBDqBznC5XDhz5gyMRiMkSQp0OkQUhIQQqK6uRmJiIlSqjvcvg7qYnjlzBklJSYFOg4hCQFlZGfr06dPh1wd1MTUajQDcX0JUVFSAsyGiYGSxWJCUlCTXk44K6mLqGdpHRUWxmBJRp3R2VyEPQBERKSCoe6YetbW1UKvVzbar1WqEh4d7xbVGpVJBr9d3KLaurg6trRcjSRIiIiI6FFtfXw+Xy9VqHgaDoUOxDQ0NcDqdisRGRETIf9GtViscDocisXq9Xj4YYLPZYLfbFYkNDw+Xf1faE2u322Gz2VqNDQsLg0ajaXesw+GA1WptNVan00Gr1bY71ul0oqGhodVYrVYLnU7X7liXy4X6+npFYjUaDcLCwgC4DwLV1dUpEtue//ct1Y0OE0HMbDYLAK3eMjMzveIjIiJajU1PT/eKjY2NbTV21KhRXrF9+/ZtNXbQoEFesYMGDWo1tm/fvl6xo0aNajU2NjbWKzY9Pb3V2IiICK/YzMzMNr+3pmbMmNFmbE1NjRw7e/bsNmMrKirk2KeeeqrN2JKSEjn2ueeeazP24MGDcmxWVlabsXv37pVj33zzzTZj8/Ly5Nj333+/zdgvvvhCjs3JyWkz9tNPP5VjP/300zZjc3Jy5Ngvvviizdj3339fjs3Ly2sz9s0335Rj9+7d22ZsVlaWHHvw4ME2Y5977jk5tqSkpM3Yp556So6tqKhoM3b27NlybE1NTZuxM2bM8Podbis2MzNTriNms1l0Rkj0TImoe7hUa8X+k5WotztRWnq+zdj/llXh/317DCqVhFrzpTZjSy/UYm3+KahUgK2h9d4uAJyurMfa/FNwCYGrLTBadqkOq38oaTvIR1Jj5Q5KFosFJpMJZ86cafEAFIf5LcdymN81w3yXS8DmdMHqcMHmcMHmcMLmcAEaLVxQwepwoa7Bhtq6etgcAlan+3mb0wW7wwUhALVWC42mcejucMBht6G14ySaJrEupxN2m3uXgFMADqcLdqcLNqeA0+WCEyoIlRZ2pwtWmx3WBitsThccLgG70wW7S8DhcL/GKanhhBo2pwv1Vjvq6uvQYHOh3u6E64pfZUmthqR25yCEC8Le+ncmqdSQNJ5YAWFvfRdG+2JVkDQ6+bHL1vouDEmlgnA6UPbufTCbzZ06kB0SPVODweBVANqKa897+qppAVQytmnB7kiszeFC4ekq/Hy+FlHhFvSI0KFnhA49IrToERGOMM3V9xc1/WN0NWFhYfJ+LV9jrQ4nLPUOWBrssNTbYWlwNP680Mp278d2pwsatQpalQStRgWNSgWtWoJWrYJGLUGrUkGrka7Y7o7XND7WqlVQSRKcLncxcboE7I1Fx/PY4RJwOF3yfWdj0Wn62BMjF0Rn0PZTrkLtvmm18hFstUqCXqtGuFYNvU6FcI0ael3jY60a4VoV9Fr3Np1aBQHAJQScLsDlEnAKAZcQjffd29zPN/0J+b5nuyRJkACoJAkqFSBBgiS5H8s/4e6oqCTI21SSBDTet9fVYIUC30pIFFNyszqc+G+ZGXuOX8SekkvIP+EebrVGr1WjZ4QWpggdekZoG4ts4329zvtxhBZR4Vp378TmRF3jrd7uQL3NhTqbA/X2JtttjsafnrjG+3b39poGd6FssLfem/aVzeGCDQBsrbe1OwjTqKDTqNw/1SqEad2FRdd0e+NzOo27wDd1ZWluaYTTbIsAVCoJWrUEnecPTOMfEG0L9zVqFXRX3NeoVNBq3DEROs3lwqhVI1ynRrhGDa1aCtqzEC0WC1b8tvPvw2IaxBrsThSUVWH38YvYc/wS9p2shNXhXZyiDToMSohCnc2Bqno7qursqKqzwSXgLnBmJ86YWx8GdRVjuAZR4VpE6bWICtcgSq+FSa9t3Nb8Oc92nVrlHpI29gTtThccTgG7yz1U9gxbHc4mw9fGIa/defl1TiGgUUnQqCSo1Sr3z8YipFa1/lijkqBRez++XBzVcnEM5mJDvmExDSL1Nif2n6zE7uMXsbvkEgrKqtz74JqIjdRhTHIMxqREY0xyDPrHRUKl8v5P7HIJVFsdqKqzoarOjso6G8z1dlTW2rwKbmXjz6rG56qtDmjVKkTo1IhoHLLpdWpEaDXun57HOjUidBrotZe3ue9r5MdNi2dkmAZqFQsNBTcW026szuZA/olK7Dl+CbuPX8R/T1U12w8XZwzDmJQYjEmOxtiUGNzYy3DVHpBKJcHU2PPrG+N7PqJxHxURNcdi2o3UWB34qfQS9pRcwp7jF3HglBmOKw6XJpjCMSY5GmNSYjA2JQY3xER0WYFjISVqHYtpAFka7O7i2djzPHjGAucVxfO6HnqMSYnG2GR38UyK1rOoEXVDLKZdyFxnx95Sd69zd8lFFJ2xNJund310hNzzHJMcjaRo36dSEVHgsJj6UWWtDXtK3L3OPSWXcKTc0uyMjORYQ2PxdB8wSuzh+9xSIuo+WEwb2Z0uNNid8rQah1PI951Nptc4XJ4pNY3PNd3mck/YPnTGgj3HL6H4XHWzz0npZcDYJgeMekf5PimeiLovFlO4zxGeuXI36vww6bt/XKS7eKZE45bkaMQZWTyJQhGLKYDconNehVSSAK1K5Z6E7TkbRHX5FEW1yn2aokYteZ2aqGnc1jc6AmNTYnBLcjRiIn07vZKIghuLKYAj5e7h+IvTBuGRW2/gBHIiajeutA/gSLkFAJCaGMVCSkQdcs0XU0uDHacq3esjDozndaSIqGO6TTHNzs6GJEl45plnuvRzjzYO8RNM4TBFaLv0s4kodHSLYvrjjz/ir3/9K9LS0rr8sz37SwfGd+4yr0R0bQt4Ma2pqcGsWbOwcuVK9OzZs81Yq9UKi8Xidessz/7SgQkc4hNRxwW8mM6dOxd33HEHbr/99qvGZmdnw2QyybekpKROf/6Rs+yZElHnBbSYrlmzBvv27UN2drZP8YsXL4bZbJZvZWVlnfp8IQSK5WE+e6ZE1HEBm2daVlaGBQsWYPPmzT5fZ6g91xjyxemq+sYFjyWk9PL9mk9ERFcKWDHNz89HRUUFRo4cKW9zOp3Yvn073n//fVitVvnKkP7iGeLf2CsSWnXA93gQURALWDGdOHEiCgsLvbY9+uijGDhwIBYtWuT3QgpcPvh0Mw8+EVEnBayYGo1GpKamem0zGAyIiYlptt1fDnNaFBEp5Joe28oHn9gzJaJO6lYLnWzdurXLPqvB7sTx8zUA2DMlos67Znum/6uogUsAPSO0iDNymTwi6pxrtpgePtt45lN8FC9QR0Sd1q5hvtlsxvr16/H999+jtLQUdXV16NWrF4YPH44pU6bg1ltv9Veeiru8v5RDfCLqPJ96pmfPnsXjjz+OhIQEvPrqq6itrcWwYcMwceJE9OnTB3l5eZg0aRIGDRqETz75xN85K4ILnBCRknzqmQ4dOhQPP/ww9u7d2+q0pfr6emzYsAHLli1DWVkZnnvuOUUTVZq8wAlPIyUiBfhUTA8dOoRevXq1GaPX6zFz5kzMnDkT58+fVyQ5fzlfbcWFGhskCRjQmz1TIuo8n4b5vXr1woYNG+B0+nb1zqsV3kDz7C9NjjFAr/P/mVZEFPp8Ppo/Y8YMXHfddVi0aBGOHDniz5z8zjPEv4n7S4lIIT4X05MnT2L+/PlYv349Bg8ejPHjxyMnJwe1tbX+zM8vDp/lsntEpCyfi2liYiKWLFmCo0eP4rvvvsONN96Ip59+GgkJCfjtb3+LXbt2+TNPRV1eXZ89UyJSRocm7aenp2P16tU4e/Ysli1bhsOHD2P8+PEYPHiw0vkpzuF04ViF+zTSm9kzJSKFdOrc/MjISGRkZKC0tBRHjhzB0aNHlcrLb0ov1sLmcCFCp0afnvpAp0NEIaJDPdO6ujqsXr0a6enpGDBgAD755BMsXLgQpaWlCqenPM/+0pvijVCpeBopESmjXT3TnTt34u9//zv+85//wOFwYPr06diyZQsyMjL8lZ/iOFmfiPzB52I6YMAA/Pzzzxg+fDjeeOMNPPDAAzCZTP7MzS88c0xv5sEnIlKQz8V06tSpmDNnDoYOHerPfPxOHubzzCciUpDPxfS9996T7zscDmzduhU///wzHnjgARiNRpw5cwZRUVGIjIz0S6JKsDTYcbqqHgCH+USkrHYfzT9x4gSmTp2KkydPwmq1YtKkSTAajXjzzTfR0NCADz/80B95KuJo4xA/0RQOU4Q2wNkQUShp99H8BQsWYNSoUaisrIRef3lq0T333INvv/1W0eSUdpjXfCIiP2l3z3THjh3YuXMndDqd1/a+ffvi9OnTiiXmD0fO8px8IvKPdvdMXS5Xi6tHnTp1CkZj9y5SXBCaiPyl3cV00qRJePfdd+XHkiShpqYGWVlZyMzMVDI3RQkhmkyL4jCfiJTV7mH+O++8g4yMDAwaNAgNDQ144IEHcOzYMcTGxuLf//63P3JUxKnKetRYHdCqJSTHGgKdDhGFmHYX08TERBQUFGDNmjXIz8+Hy+XCnDlzMGvWLK8DUt2NZ4jfL84IrfqavSgrEflJhxY60ev1ePTRR/Hoo48qnY/feA4+3cz9pUTkBz510dqzVmltbS0OHTrU4YT85cg5XtqZiPzHp2L68MMPY9KkSfj0009RU1PTYkxRURFeeOEF9OvXD/v27VM0SSVcnhbFg09EpDyfhvlFRUVYsWIFXnrpJcyaNQsDBgxAYmIiwsPDUVlZiSNHjqC2thbTp09Hbm5uq5eDDpQGuxMlF9yXV+Ewn4j8QRJCiPa8YN++ffj+++9RWlqK+vp6xMbGYvjw4cjIyEB0dLS/8myRxWKByWSC2WxGVFTrPc7CU2bc+f4ORBt0yP/j7ZAkrmNKRG6+1pGrafcBqBEjRmDEiBEd/sBAuLyGqZGFlIj84pqYI+SZFsXTSInIX66RYuqZFsWDT0TkH9dGMT3LaVFE5F8hX0zPV1txsdYGlQT0j2MxJSL/CPli6hni3xBjgF6nDnA2RBSqfC6mmZmZMJvN8uPXX38dVVVV8uOLFy9i0KBBiianBA7xiagr+FxMv/nmG1itVvnxG2+8gUuXLsmPHQ4HiouLlc1OAYd5aWci6gI+F9Mr5/a3c65/wBRzQWgi6gIhvc/U4XTh2Dn3WgLsmRKRP/lcTCVJanb2UHc/m6jkQi1sThcMOjX69Oy+a60SUfDz+XRSIQQeeeQRhIWFAQAaGhrw5JNPwmBwr1rfdH9qd3G4yZlPKlX3LvxEFNx8LqazZ8/2evzggw82i3n44Yc7n5GCij0Hn3jNJyLyM5+LaU5Ojj/z8At5WhQPPhGRn7Vr1agTJ05g8+bNsNvtmDBhQrecV9rU5Us7s2dKRP7l8wGo7du3Y/DgwXjiiScwb948DBs2rNNXI83Ozsbo0aNhNBoRFxeHu+++W7G5quZ6O05X1QPgalFE5H8+F9MXX3wRGRkZOHXqFC5evIjHHnsMv//97zv14du2bcPcuXOxe/du5ObmwuFwYPLkyaitre3U+wLA0cZrPl3XQw+TXtvp9yMiaovPw/zCwkJs374diYmJAIC3334bK1euRGVlJXr27NmhD9+0aZPX45ycHMTFxSE/Px+33XZbs3ir1eo1a8BisbT63pev+cReKRH5n88906qqKsTFxcmPDQYDIiIivM7P7yzPuf+tXf4kOzsbJpNJviUlJbX6Xod55hMRdaF2HYAqKipCeXm5/FgIgcOHD6O6ulrelpaW1qFEhBBYuHAhxo8f3+oF+RYvXoyFCxfKjy0WS6sF1dMz5bQoIuoK7SqmEydObHZO/rRp0yBJEoQQkCQJTqezQ4nMmzcPBw4cwI4dO1qNCQsLk08aaIvLJXBUPo2UPVMi8j+fi2lJSYnfkpg/fz42btyI7du3o0+fPp1+v9NV9aixOqBTq5Aca1AgQyKitvlcTPv27av4hwshMH/+fKxfvx5bt25FcnKyIu97uHGI3y8uElp1SK/lQkTdhM+VZuLEiVi3bl2rz1+4cAEpKSnt+vC5c+fin//8J/71r3/BaDSivLwc5eXlqK+vb9f7XEledo8LQhNRF/G5mObl5eG+++5DVlZWi887nU6cOHGiXR++fPlymM1mTJgwAQkJCfLtk08+adf7XOkIj+QTURdr1wGo5cuX4/nnn8eBAwfw8ccfIzIyslMf7q8Fprm6PhF1tXbtUPz1r3+NXbt2oaioCOPGjcPx48f9lVeHNdidKL3gPoOKw3wi6irtPjpz8803Y+/evUhKSsLo0aOxZcsWf+TVYcfO1cAlgBiDDr0irz6NiohICR061G0ymfDll1/i8ccfR2ZmJt555x2l8+owzxD/pnhjt78SABGFDp/3mbZ0yZKlS5di+PDhmDNnDr777jvFk+uIy2uYcn8pEXWdDl+d1OP+++/Hjh07UFhYqFhSnXFEXl2f+0uJqOv43DPNy8trdQGSYcOGIT8/H19++aViiXWEEEKeFnUze6ZE1IV8Lqbp6eltPh8TExPwa0Cdr7HiUq0NKgno37tz07aIiNojpM619OwvvSHWgHCtOsDZENG1JLSKaeP+Ug7xiairhVgx5WmkRBQYoVVMG4f5vFQJEXW1kCmmdqcL/6twLwh9M1fXJ6IuFjLFtORCLWxOFyLDNLiuhz7Q6RDRNSZkiqlnf+lN8UaoVDyNlIi6VugUU17amYgCKHSKqXzmE4spEXW90CmmvLQzEQVQSBRTc70dZ8wNADjMJ6LACIlieuyce4h/XQ89osK1Ac6GiK5FoVFMK3jmExEFVkgU0+Jy92R9rmFKRIESEsXUM8y/iQucEFGAhEYxreC0KCIKrJAopnU2F3RqFZJjDYFOhYiuUSFRTAH3yvoadcg0h4iCTMhUH84vJaJACpliytX1iSiQQqaYcloUEQVS6BRT9kyJKIBCophGR2jRyxgW6DSI6BoWEsV0AA8+EVGAhUQx7d+bxZSIAiskiumAuMhAp0BE17iQKKZTUuMDnQIRXeNCophG6DSBToGIrnEhUUyJiAKNxZSISAFBPT4WQgAALBZLgDMhomDlqR+eetJRQV1ML168CABISkoKcCZEFOyqq6thMpk6/PqgLqbR0dEAgJMnT3bqS+hOLBYLkpKSUFZWhqio0DhFlm3q/kKtPYDvbRJCoLq6GomJiZ36vKAupiqVe5evyWQKmV8Aj6ioKLYpCIRam0KtPYBvbVKiM8YDUERECmAxJSJSQFAX07CwMGRlZSEsLHRWjGKbgkOotSnU2gN0fZsk0dn5AEREFNw9UyKi7oLFlIhIASymREQKYDElIlJAUBfTDz74AMnJyQgPD8fIkSPx/fffBzolAMD27dtx5513IjExEZIkYcOGDV7PCyHw8ssvIzExEXq9HhMmTMChQ4e8YqxWK+bPn4/Y2FgYDAbcddddOHXqlFdMZWUlHnroIZhMJphMJjz00EOoqqpSvD3Z2dkYPXo0jEYj4uLicPfdd6O4uDio27R8+XKkpaXJE7rHjRuHr7/+Omjbc6Xs7GxIkoRnnnkmaNv08ssvQ5Ikr1t8/OW1i7tde0SQWrNmjdBqtWLlypWiqKhILFiwQBgMBnHixIlApya++uorsWTJErF27VoBQKxfv97r+aVLlwqj0SjWrl0rCgsLxf333y8SEhKExWKRY5588klx3XXXidzcXLFv3z6RkZEhhg4dKhwOhxwzdepUkZqaKn744Qfxww8/iNTUVDFt2jTF2zNlyhSRk5MjDh48KAoKCsQdd9whrr/+elFTUxO0bdq4caP48ssvRXFxsSguLhYvvPCC0Gq14uDBg0HZnqb27t0rbrjhBpGWliYWLFggbw+2NmVlZYnBgweLs2fPyreKiopu256gLaa33HKLePLJJ722DRw4UPzhD38IUEYtu7KYulwuER8fL5YuXSpva2hoECaTSXz44YdCCCGqqqqEVqsVa9askWNOnz4tVCqV2LRpkxBCiKKiIgFA7N69W47ZtWuXACCOHDni1zZVVFQIAGLbtm0h0yYhhOjZs6f46KOPgro91dXVon///iI3N1ekp6fLxTQY25SVlSWGDh3a4nPdsT1BOcy32WzIz8/H5MmTvbZPnjwZP/zwQ4Cy8k1JSQnKy8u9cg8LC0N6erqce35+Pux2u1dMYmIiUlNT5Zhdu3bBZDJhzJgxcszYsWNhMpn8/h2YzWYAlxeaCfY2OZ1OrFmzBrW1tRg3blxQt2fu3Lm44447cPvtt3ttD9Y2HTt2DImJiUhOTsZvfvMbHD9+vNu2JygXOrlw4QKcTid69+7ttb13794oLy8PUFa+8eTXUu4nTpyQY3Q6HXr27NksxvP68vJyxMXFNXv/uLg4v34HQggsXLgQ48ePR2pqqpyLJ78r8+3ObSosLMS4cePQ0NCAyMhIrF+/HoMGDZL/EwVbe9asWYN9+/bhxx9/bPZcMP4bjRkzBv/4xz8wYMAAnDt3Dq+99hpuvfVWHDp0qFu2JyiLqYckSV6PhRDNtnVXHcn9ypiW4v39HcybNw8HDhzAjh07mj0XbG266aabUFBQgKqqKqxduxazZ8/Gtm3bWs2lO7enrKwMCxYswObNmxEeHt5qXDC16Ve/+pV8f8iQIRg3bhxuvPFGrF69GmPHjm0xl0C2JyiH+bGxsVCr1c3+clRUVDT7S9XdeI5GtpV7fHw8bDYbKisr24w5d+5cs/c/f/68376D+fPnY+PGjcjLy0OfPn3k7cHaJp1Oh379+mHUqFHIzs7G0KFD8ec//zko25Ofn4+KigqMHDkSGo0GGo0G27Ztw3vvvQeNRiN/XjC16UoGgwFDhgzBsWPHuuW/UVAWU51Oh5EjRyI3N9dre25uLm699dYAZeWb5ORkxMfHe+Vus9mwbds2OfeRI0dCq9V6xZw9exYHDx6UY8aNGwez2Yy9e/fKMXv27IHZbFb8OxBCYN68eVi3bh2+++47JCcnB32bWiKEgNVqDcr2TJw4EYWFhSgoKJBvo0aNwqxZs1BQUICUlJSga9OVrFYrDh8+jISEhO75b9Suw1XdiGdq1N/+9jdRVFQknnnmGWEwGERpaWmgUxPV1dVi//79Yv/+/QKAWLZsmdi/f788bWvp0qXCZDKJdevWicLCQjFz5swWp3T06dNHbNmyRezbt0/88pe/bHFKR1pamti1a5fYtWuXGDJkiF+mqPzud78TJpNJbN261WuaSl1dnRwTbG1avHix2L59uygpKREHDhwQL7zwglCpVGLz5s1B2Z6WND2aH4xtevbZZ8XWrVvF8ePHxe7du8W0adOE0WiU/493t/YEbTEVQoi//OUvom/fvkKn04kRI0bIU3UCLS8vTwBodps9e7YQwj2tIysrS8THx4uwsDBx2223icLCQq/3qK+vF/PmzRPR0dFCr9eLadOmiZMnT3rFXLx4UcyaNUsYjUZhNBrFrFmzRGVlpeLtaaktAEROTo4cE2xteuyxx+TfnV69eomJEyfKhTQY29OSK4tpsLXJM29Uq9WKxMREMX36dHHo0KFu2x4uwUdEpICg3GdKRNTdsJgSESmAxZSISAEspkRECmAxJSJSAIspEZECWEyJiBTAYkpEpAAWU6JGq1atQo8ePQKdBgUpFlPqUo888kiz6/pIkoSpU6cGOjXcf//9OHr0aKDToCAV1OuZUnCaOnUqcnJyvLaFhYX57fNsNht0Ot1V4/R6PfR6vd/yoNDGnil1ubCwMMTHx3vdPKuhS5KEjz76CPfccw8iIiLQv39/bNy40ev1RUVFyMzMRGRkJHr37o2HHnoIFy5ckJ+fMGEC5s2bh4ULFyI2NhaTJk0CAGzcuBH9+/eHXq9HRkYGVq9eDUmS5CtRtjTM//zzzzFy5EiEh4cjJSUFr7zyChwOh/z8yy+/jOuvvx5hYWFITEzE008/7YdvjIIBiyl1O6+88gruu+8+HDhwAJmZmZg1axYuXboEwL0eZXp6OoYNG4affvoJmzZtwrlz53Dfffd5vcfq1auh0Wiwc+dOrFixAqWlpZgxYwbuvvtuFBQU4IknnsCSJUvazOObb77Bgw8+iKeffhpFRUVYsWIFVq1ahddffx0A8Nlnn+Gdd97BihUrcOzYMWzYsAFDhgzxz5dC3V8HVsYi6rDZs2cLtVotDAaD1+3VV18VQriX+/vjH/8ox9fU1AhJksTXX38thBDixRdfFJMnT/Z6z7KyMgFAFBcXCyHcS88NGzbMK2bRokUiNTXVa9uSJUsEAHm5tZycHGEymeTn/+///k/86U9/8nrNxx9/LBISEoQQQrz99ttiwIABwmazdfDboFDCfabU5TIyMrB8+XKvbZ4rnQJAWlqafN9gMMBoNKKiogKA+/IceXl5iIyMbPa+P//8MwYMGAAAGDVqlNdzxcXFGD16tNe2W265pc088/Pz8eOPP8o9UcB9JdOGhgbU1dXh3nvvxbvvvouUlBRMnToVmZmZuPPOO6HR8L/VtYj/6tTlDAYD+vXr1+rzWq3W67EkSXC5XAAAl8uFO++8E2+88Uaz1yUkJHh9RlOihQukiass5etyufDKK69g+vTpzZ4LDw9HUlISiouLkZubiy1btuCpp57CW2+9hW3btjVrA4U+FlMKKiNGjMDatWtxww03tKsHOHDgQHz11Vde23766aerflZxcXGbhV+v1+Ouu+7CXXfdhblz52LgwIEoLCzEiBEjfM6NQgOLKXU5q9Xa7KqSGo0GsbGxV33t3LlzsXLlSsycORPPP/88YmNj8b///Q9r1qzBypUroVarW3zdE088gWXLlmHRokWYM2cOCgoKsGrVKgAtX+oXAF566SVMmzYNSUlJuPfee6FSqXDgwAEUFhbitddew6pVq+B0OjFmzBhERETg448/hl6vR9++fdv3hVBI4NF86nKbNm1CQkKC1238+PE+vTYxMRE7d+6E0+nElClTkJqaigULFsBkMkGlav3XOTk5GZ999hnWrVuHtLQ0LF++XD6a39oc1ylTpuCLL75Abm4uRo8ejbFjx2LZsmVysezRowdWrlyJX/ziF0hLS8O3336Lzz//HDExMe38RigU8BpQdM16/fXX8eGHH6KsrCzQqVAI4DCfrhkffPABRo8ejZiYGOzcuRNvvfUW5s2bF+i0KESwmNI149ixY3jttddw6dIlXH/99Xj22WexePHiQKdFIYLDfCIiBfAAFBGRAlhMiYgUwGJKRKQAFlMiIgWwmBIRKYDFlIhIASymREQKYDElIlLA/wfnmwzhEO/9jwAAAABJRU5ErkJggg==", "text/plain": [ "
" ] @@ -708,7 +717,10 @@ "id": "e8788f74-c208-4939-aa9b-3bbdfd8310ee", "metadata": {}, "source": [ - "We consistently underestimate the ZPE. Is it because we have too few oscillators?" + "We consistently underestimate the ZPE. Is it because we have too few oscillators?\n", + "\n", + "TODO: \n", + "- Save this data so we can compare it in a later notebook" ] }, {