solutions_tenpy.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "74eb0695",
   "metadata": {},
   "source": [
    "# Solutions for the `exercise_tenpy.ipynb` notebook"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "39301469",
   "metadata": {},
   "source": [
    "Uncomment and run the cells below (removing the `#`) when running this notebook on https://colab.research.google.com.\n",
    "\n",
    "Alternatively, you can run this notebook locally with jupyter, provided that you have the `toycodes` subfolder from \n",
    "\n",
    "https://github.com/jhauschild/WinterschoolBarcelona2021\n",
    "\n",
    "in the same folder as your notebook."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "d52f6ce1",
   "metadata": {},
   "outputs": [],
   "source": [
    "#!pip install git+git://github.com/tenpy/tenpy"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "8d42e939",
   "metadata": {},
   "outputs": [],
   "source": [
    "#!pip install git+git://github.com/jhauschild/WinterschoolBarcelona2021\n",
    "\n",
    "# use `pip uninstall tenpy-toycodes` to remove it a gain."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "e7d5c720",
   "metadata": {},
   "source": [
    "This tutorial focuses on a set of toy codes (using only python with numpy + scipy) that provide a simple implementation of the various MPS algorithms.\n",
    "\n",
    "You can add your code below by inserting additional cells as neccessary and running them (press Shift+Enter).\n",
    "\n",
    "**DISCLAIMER**: Like for the toy codes, we only use very small bond dimensions here. For state-of-the-art MPS calculations (especially for cylinders towards 2D), `chi` should be significantly larger, often on the order of several 1000s."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "b7d6aedd",
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import scipy\n",
    "import matplotlib.pyplot as plt\n",
    "from pprint import pprint\n",
    "\n",
    "np.set_printoptions(precision=5, suppress=True, linewidth=100)\n",
    "plt.rcParams['figure.dpi'] = 150"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "2932ad48",
   "metadata": {},
   "outputs": [],
   "source": [
    "import tenpy\n",
    "import tenpy.linalg.np_conserved as npc\n",
    "from tenpy.algorithms import tebd, dmrg, tdvp\n",
    "from tenpy.networks.site import SpinHalfSite, SpinSite, FermionSite\n",
    "from tenpy.networks.mps import MPS\n",
    "from tenpy.models.tf_ising import TFIChain\n",
    "\n",
    "tenpy.tools.misc.setup_logging(to_stdout=\"INFO\")"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "9d536bdf",
   "metadata": {},
   "source": [
    "# Overview\n",
    "\n",
    "The source code of TeNPy is at https://github.com/tenpy/tenpy/; \n",
    "you can find links to the documentation and the forum in the Readme there.\n",
    "\n",
    "The [**documentation**](https://tenpy.readthedocs.io) is roughly split into the \"User guide\" (upper part in the left side-bar) and the reference of all the functions and classes (lower part).\n",
    "\n"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "6a01aec3",
   "metadata": {},
   "source": [
    "## Exercise(s)\n",
    "\n",
    "Read the [overview](https://tenpy.readthedocs.io/en/latest/intro/overview.html) of the TeNPy documentation.\n",
    "\n",
    "Whenever you hit an example code, try to copy it here and run it.\n",
    "\n",
    "Try to modify it slightly and try to rerun it; for example try to calculate the overlap `<v|v>` in the first example."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "4f1b822f",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "id": "337e0452",
   "metadata": {},
   "source": [
    "# Initializing a Model\n",
    "\n",
    "In TeNPy, the model defines the Hilbert space and local operators, and ultimately fixes whether charge conservation is used. Therefore, you should usually start with the initialization of the model.\n",
    "There are many predefined models in `tenpy.models`, that you can often just use.\n",
    "\n",
    "We will first initialize the transverse field Ising model. Note, however, that TeNPy switches X and Z, which allows the conservation of the overall spin parity:\n",
    "$$H = - J \\sum_{i} \\sigma^z_i \\sigma^z_{i+1} - g \\sum_{i} \\sigma^x_i  \\textrm{ in toy codes} \\\\ \n",
    "\\rightarrow H = - J \\sum_{i} \\sigma^x_i \\sigma^x_{i+1} - g \\sum_{i} \\sigma^z_i  \\textrm{ in TeNPy}$$\n",
    "\n",
    "In TeNPy, allmost all parmaters can be changed dynamically through options. Default parameters are written back into the dictionaries."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "a4f0ecc0",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "INFO    : TFIChain: reading 'bc_MPS'='finite'\n",
      "INFO    : TFIChain: reading 'conserve'='best'\n",
      "INFO    : TFIChain: set conserve to parity\n",
      "INFO    : TFIChain: reading 'L'=20\n",
      "INFO    : TFIChain: reading 'g'=1.0\n",
      "used parameters, including default/not specified ones:\n",
      "{'J': 1.0,\n",
      " 'L': 20,\n",
      " 'bc_MPS': 'finite',\n",
      " 'bc_x': 'open',\n",
      " 'conserve': 'best',\n",
      " 'explicit_plus_hc': False,\n",
      " 'g': 1.0,\n",
      " 'lattice': <class 'tenpy.models.lattice.Chain'>,\n",
      " 'order': 'default',\n",
      " 'sort_mpo_legs': False}\n"
     ]
    }
   ],
   "source": [
    "model_params = {\n",
    "    'L': 20,\n",
    "    'g': 1.0,\n",
    "    'bc_MPS': 'finite',\n",
    "    'conserve': 'best'\n",
    "}\n",
    "model = TFIChain(model_params)\n",
    "# you can now print the default parameters used:\n",
    "print(\"used parameters, including default/not specified ones:\")\n",
    "pprint(model_params)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "9fe62909",
   "metadata": {},
   "source": [
    "Given the model, one can easily initialize a product state, e.g. for the Neel state:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "236024c3",
   "metadata": {},
   "outputs": [],
   "source": [
    "p_state = ['up', 'down'] * (model.lat.N_sites//2)\n",
    "psi = MPS.from_product_state(model.lat.mps_sites(), p_state, bc=model.lat.bc_MPS)\n"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "e3f4a8cf",
   "metadata": {},
   "source": [
    "Measuring expectation values is also similar to the toycode. However, we can even specify the local operators (defined in the sites) as strings:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "adf6d283",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "<sigmaz> =  [1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1.]\n",
      "S =  [-0. -0. -0. -0. -0. -0. -0. -0. -0. -0. -0. -0. -0. -0. -0. -0. -0. -0. -0.]\n"
     ]
    }
   ],
   "source": [
    "print(\"<sigmaz> = \", psi.expectation_value('Sigmaz'))\n",
    "print(\"S = \", psi.entanglement_entropy())"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "017d4d67",
   "metadata": {},
   "source": [
    "### Exercise\n",
    "\n",
    "Check the [`model.bond_energies`](https://tenpy.readthedocs.io/en/latest/reference/tenpy.models.model.NearestNeighborModel.html#tenpy.models.model.NearestNeighborModel.bond_energies) for the Neel state and make sure it matches what you expect.\n",
    "\n",
    "**Note** that you get different results than in the toycode notebook due to the switch of X <-> Z."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "d0ef22db",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "energy Neel: [-1.5 -1.  -1.  -1.  -1.  -1.  -1.  -1.  -1.  -1.  -1.  -1.  -1.  -1.  -1.  -1.  -1.  -1.  -1.5]\n"
     ]
    }
   ],
   "source": [
    "E = model.bond_energies(psi)\n",
    "print(\"energy Neel:\", E)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "99903e46",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "id": "bd218561",
   "metadata": {},
   "source": [
    "# Running DMRG\n",
    "\n",
    "Given the model and state, running DMRG isn't hard.\n",
    "Again, there are many (default) parameters for fine-tuning, see [this full option list](https://tenpy.readthedocs.io/en/latest/reference/tenpy.algorithms.dmrg.TwoSiteDMRGEngine.html#cfg-config-TwoSiteDMRGEngine) for details."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "ffc67843",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "INFO    : TwoSiteDMRGEngine: subconfig 'trunc_params'=Config(<2 options>, 'trunc_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_sweeps'=40\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : trunc_params: reading 'chi_max'=30\n",
      "INFO    : trunc_params: reading 'svd_min'=1e-07\n",
      "INFO    : checkpoint after sweep 1\n",
      "energy=-25.1067276087268887, max S=0.3801726459659367, age=20, norm_err=5.0e-02\n",
      "Current memory usage 130.2MB, wall time: 0.3s\n",
      "Delta E = nan, Delta S = 3.5048e-01 (per sweep)\n",
      "max trunc_err = 0.0000e+00, max E_trunc = 1.7764e-14\n",
      "chi: [2, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 2]\n",
      "================================================================================\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : checkpoint after sweep 2\n",
      "energy=-25.1077971109763851, max S=0.4437610868090598, age=20, norm_err=3.0e-04\n",
      "Current memory usage 130.6MB, wall time: 0.6s\n",
      "Delta E = -1.0695e-03, Delta S = 4.3710e-02 (per sweep)\n",
      "max trunc_err = 9.4252e-15, max E_trunc = 1.1013e-13\n",
      "chi: [2, 4, 8, 12, 12, 13, 13, 14, 14, 14, 14, 14, 13, 13, 12, 12, 8, 4, 2]\n",
      "================================================================================\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : checkpoint after sweep 3\n",
      "energy=-25.1077971116229968, max S=0.4437994305126334, age=20, norm_err=1.4e-09\n",
      "Current memory usage 130.6MB, wall time: 0.5s\n",
      "Delta E = -6.4661e-10, Delta S = 4.8294e-05 (per sweep)\n",
      "max trunc_err = 1.0328e-14, max E_trunc = 1.3500e-13\n",
      "chi: [2, 4, 8, 12, 14, 15, 16, 17, 17, 17, 17, 17, 16, 15, 14, 12, 8, 4, 2]\n",
      "================================================================================\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : checkpoint after sweep 4\n",
      "energy=-25.1077971116230394, max S=0.4437994305913012, age=20, norm_err=2.4e-11\n",
      "Current memory usage 130.6MB, wall time: 0.4s\n",
      "Delta E = -4.2633e-14, Delta S = 1.6891e-10 (per sweep)\n",
      "max trunc_err = 1.0277e-14, max E_trunc = 1.1369e-13\n",
      "chi: [2, 4, 8, 12, 14, 15, 16, 17, 17, 17, 17, 17, 16, 15, 14, 12, 8, 4, 2]\n",
      "================================================================================\n",
      "INFO    : DMRG finished after 4 sweeps, max chi=17\n"
     ]
    }
   ],
   "source": [
    "p_state = ['up'] * model.lat.N_sites\n",
    "psi = MPS.from_product_state(model.lat.mps_sites(), p_state, bc=model.lat.bc_MPS)\n",
    "algorithm_params = {\n",
    "    'trunc_params': {\n",
    "        'chi_max': 30,\n",
    "        'svd_min': 1.e-7,\n",
    "    },\n",
    "    'max_sweeps': 40,\n",
    "}\n",
    "eng = dmrg.TwoSiteDMRGEngine(psi, model, algorithm_params)\n",
    "E, psi = eng.run()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "3d9130aa",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "id": "b0eee6bf",
   "metadata": {},
   "source": [
    "### Exercise\n",
    "\n",
    "Run DMRG for `'infinite'` MPS. \n",
    "(You need to initialize a new model, state, and DMRG engine for this.)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "90538057",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "E_exact = -1.2732395447351628\n"
     ]
    }
   ],
   "source": [
    "import tfi_exact\n",
    "\n",
    "print(\"E_exact =\", tfi_exact.infinite_gs_energy(model_params['J'], model_params['g']))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "a9503ce9",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "INFO    : TFIChain: reading 'bc_MPS'='infinite'\n",
      "INFO    : TFIChain: reading 'conserve'='best'\n",
      "INFO    : TFIChain: set conserve to parity\n",
      "INFO    : TFIChain: reading 'L'=2\n",
      "INFO    : TFIChain: reading 'g'=1.2\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'trunc_params'=Config(<2 options>, 'trunc_params')\n",
      "INFO    : start environment_sweep\n",
      "INFO    : trunc_params: reading 'chi_max'=30\n",
      "INFO    : trunc_params: reading 'svd_min'=1e-07\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_sweeps'=40\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 10\n",
      "energy=-1.4196192748981424, max S=0.2533554575580327, age=128, norm_err=2.2e-12\n",
      "Current memory usage 131.5MB, wall time: 0.6s\n",
      "Delta E = nan, Delta S = 2.5336e-02 (per sweep)\n",
      "max trunc_err = 3.6232e-15, max E_trunc = 1.1369e-13\n",
      "chi: [18, 18]\n",
      "================================================================================\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 20\n",
      "energy=-1.4196192748981744, max S=0.2533554570965122, age=248, norm_err=8.1e-15\n",
      "Current memory usage 131.5MB, wall time: 0.6s\n",
      "Delta E = -3.1974e-15, Delta S = -4.5210e-11 (per sweep)\n",
      "max trunc_err = 3.7398e-15, max E_trunc = -5.6843e-14\n",
      "chi: [18, 18]\n",
      "================================================================================\n",
      "INFO    : DMRG finished after 20 sweeps, max chi=18\n"
     ]
    }
   ],
   "source": [
    "model_params = {\n",
    "    'L': 2,\n",
    "    'g': 1.2,\n",
    "    'bc_MPS': 'infinite',\n",
    "    'conserve': 'best'\n",
    "}\n",
    "model = TFIChain(model_params)\n",
    "\n",
    "p_state = ['up', 'up'] * (model.lat.N_sites // 2)\n",
    "psi = MPS.from_product_state(model.lat.mps_sites(), p_state, bc=model.lat.bc_MPS)\n",
    "algorithm_params = {\n",
    "    'trunc_params': {\n",
    "        'chi_max': 30,\n",
    "        'svd_min': 1.e-7,\n",
    "    },\n",
    "    'max_sweeps': 40,\n",
    "}\n",
    "eng = dmrg.TwoSiteDMRGEngine(psi, model, algorithm_params)\n",
    "E, psi = eng.run()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "e0b86afc",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "id": "fa17e5ea",
   "metadata": {},
   "source": [
    "### Exercise\n",
    "\n",
    "Reproduce the phase-diagram plot of the transverse field Ising model from the toy code noteboook with TeNPy.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "e9a35a7b",
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "INFO    : TFIChain: reading 'bc_MPS'='infinite'\n",
      "INFO    : TFIChain: reading 'L'=2\n",
      "INFO    : TFIChain: reading 'J'=1.0\n",
      "INFO    : TFIChain: reading 'g'=0.0\n",
      "INFO    : TwoSiteDMRGEngine: reading 'diag_method'='default'\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'trunc_params'=Config(<5 options>, 'trunc_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'combine'=False\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'lanczos_params'=Config(<9 options>, 'lanczos_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'start_env'=1\n",
      "INFO    : start environment_sweep\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list_reactivates_mixer'=True\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_to_trunc'=0.05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_min'=4.999999999999999e-16\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_max'=0.0001\n",
      "INFO    : TwoSiteDMRGEngine: reading 'E_tol_to_trunc'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'N_sweeps_check'=10\n",
      "INFO    : TwoSiteDMRGEngine: reading 'min_sweeps'=15\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_sweeps'=40\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_E_err'=1e-08\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_S_err'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_hours'=8760\n",
      "INFO    : TwoSiteDMRGEngine: reading 'update_env'=5\n",
      "INFO    : TwoSiteDMRGEngine: reading 'mixer'=None\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_N_for_ED'=400\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 10\n",
      "energy=-1.0000000000000160, max S=0.6931471805599453, age=128, norm_err=3.4e-16\n",
      "Current memory usage 131.5MB, wall time: 0.4s\n",
      "Delta E = nan, Delta S = 6.9315e-02 (per sweep)\n",
      "max trunc_err = 9.4504e-30, max E_trunc = 4.2633e-14\n",
      "chi: [2, 2]\n",
      "================================================================================\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 20\n",
      "energy=-1.0000000000000249, max S=0.6931471805599452, age=248, norm_err=4.4e-16\n",
      "Current memory usage 131.5MB, wall time: 0.3s\n",
      "Delta E = -8.8818e-16, Delta S = 0.0000e+00 (per sweep)\n",
      "max trunc_err = 4.8993e-28, max E_trunc = 8.5265e-14\n",
      "chi: [2, 2]\n",
      "================================================================================\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol_iter'=5\n",
      "INFO    : DMRG finished after 20 sweeps, max chi=2\n",
      "INFO    : TFIChain: reading 'bc_MPS'='infinite'\n",
      "INFO    : TFIChain: reading 'L'=2\n",
      "INFO    : TFIChain: reading 'J'=1.0\n",
      "INFO    : TFIChain: reading 'g'=0.1\n",
      "INFO    : TwoSiteDMRGEngine: reading 'diag_method'='default'\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'trunc_params'=Config(<5 options>, 'trunc_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'combine'=False\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'lanczos_params'=Config(<9 options>, 'lanczos_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'start_env'=1\n",
      "INFO    : start environment_sweep\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list_reactivates_mixer'=True\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_to_trunc'=0.05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_min'=4.999999999999999e-16\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_max'=0.0001\n",
      "INFO    : TwoSiteDMRGEngine: reading 'E_tol_to_trunc'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'N_sweeps_check'=10\n",
      "INFO    : TwoSiteDMRGEngine: reading 'min_sweeps'=15\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_sweeps'=40\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_E_err'=1e-08\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_S_err'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_hours'=8760\n",
      "INFO    : TwoSiteDMRGEngine: reading 'update_env'=5\n",
      "INFO    : TwoSiteDMRGEngine: reading 'mixer'=None\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_N_for_ED'=400\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 10\n",
      "energy=-1.0025015664216266, max S=0.6931533932074605, age=128, norm_err=4.6e-16\n",
      "Current memory usage 131.5MB, wall time: 0.4s\n",
      "Delta E = nan, Delta S = 6.9315e-02 (per sweep)\n",
      "max trunc_err = 1.2209e-19, max E_trunc = 1.4211e-14\n",
      "chi: [6, 6]\n",
      "================================================================================\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 20\n",
      "energy=-1.0025015664216532, max S=0.6931533932074607, age=248, norm_err=8.0e-16\n",
      "Current memory usage 131.5MB, wall time: 0.4s\n",
      "Delta E = -2.6645e-15, Delta S = 0.0000e+00 (per sweep)\n",
      "max trunc_err = 1.2209e-19, max E_trunc = -1.1369e-13\n",
      "chi: [6, 6]\n",
      "================================================================================\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol_iter'=5\n",
      "INFO    : DMRG finished after 20 sweeps, max chi=6\n",
      "INFO    : TFIChain: reading 'bc_MPS'='infinite'\n",
      "INFO    : TFIChain: reading 'L'=2\n",
      "INFO    : TFIChain: reading 'J'=1.0\n",
      "INFO    : TFIChain: reading 'g'=0.2\n",
      "INFO    : TwoSiteDMRGEngine: reading 'diag_method'='default'\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'trunc_params'=Config(<5 options>, 'trunc_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'combine'=False\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'lanczos_params'=Config(<9 options>, 'lanczos_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'start_env'=1\n",
      "INFO    : start environment_sweep\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list_reactivates_mixer'=True\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_to_trunc'=0.05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_min'=4.999999999999999e-16\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_max'=0.0001\n",
      "INFO    : TwoSiteDMRGEngine: reading 'E_tol_to_trunc'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'N_sweeps_check'=10\n",
      "INFO    : TwoSiteDMRGEngine: reading 'min_sweeps'=15\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_sweeps'=40\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_E_err'=1e-08\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_S_err'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_hours'=8760\n",
      "INFO    : TwoSiteDMRGEngine: reading 'update_env'=5\n",
      "INFO    : TwoSiteDMRGEngine: reading 'mixer'=None\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_N_for_ED'=400\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 10\n",
      "energy=-1.0100252539846153, max S=0.6932314274512770, age=128, norm_err=3.0e-16\n",
      "Current memory usage 131.5MB, wall time: 0.4s\n",
      "Delta E = nan, Delta S = 6.9323e-02 (per sweep)\n",
      "max trunc_err = 5.3784e-16, max E_trunc = -2.8422e-14\n",
      "chi: [6, 6]\n",
      "================================================================================\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 20\n",
      "energy=-1.0100252539846473, max S=0.6932314274512774, age=248, norm_err=4.5e-16\n",
      "Current memory usage 131.5MB, wall time: 0.4s\n",
      "Delta E = -3.1974e-15, Delta S = 2.2204e-17 (per sweep)\n",
      "max trunc_err = 5.3784e-16, max E_trunc = 8.5265e-14\n",
      "chi: [6, 6]\n",
      "================================================================================\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol_iter'=5\n",
      "INFO    : DMRG finished after 20 sweeps, max chi=6\n",
      "INFO    : TFIChain: reading 'bc_MPS'='infinite'\n",
      "INFO    : TFIChain: reading 'L'=2\n",
      "INFO    : TFIChain: reading 'J'=1.0\n",
      "INFO    : TFIChain: reading 'g'=0.30000000000000004\n",
      "INFO    : TwoSiteDMRGEngine: reading 'diag_method'='default'\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'trunc_params'=Config(<5 options>, 'trunc_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'combine'=False\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'lanczos_params'=Config(<9 options>, 'lanczos_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'start_env'=1\n",
      "INFO    : start environment_sweep\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list_reactivates_mixer'=True\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_to_trunc'=0.05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_min'=4.999999999999999e-16\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_max'=0.0001\n",
      "INFO    : TwoSiteDMRGEngine: reading 'E_tol_to_trunc'=None\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "INFO    : TwoSiteDMRGEngine: reading 'N_sweeps_check'=10\n",
      "INFO    : TwoSiteDMRGEngine: reading 'min_sweeps'=15\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_sweeps'=40\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_E_err'=1e-08\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_S_err'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_hours'=8760\n",
      "INFO    : TwoSiteDMRGEngine: reading 'update_env'=5\n",
      "INFO    : TwoSiteDMRGEngine: reading 'mixer'=None\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_N_for_ED'=400\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 10\n",
      "energy=-1.0226295149208635, max S=0.6935386402717904, age=128, norm_err=8.6e-16\n",
      "Current memory usage 131.5MB, wall time: 0.5s\n",
      "Delta E = nan, Delta S = 6.9354e-02 (per sweep)\n",
      "max trunc_err = 2.6584e-18, max E_trunc = 1.4211e-14\n",
      "chi: [10, 10]\n",
      "================================================================================\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 20\n",
      "energy=-1.0226295149208866, max S=0.6935386402717908, age=248, norm_err=3.5e-16\n",
      "Current memory usage 131.5MB, wall time: 0.5s\n",
      "Delta E = -2.3093e-15, Delta S = 3.3307e-17 (per sweep)\n",
      "max trunc_err = 2.6584e-18, max E_trunc = 2.8422e-14\n",
      "chi: [10, 10]\n",
      "================================================================================\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol_iter'=5\n",
      "INFO    : DMRG finished after 20 sweeps, max chi=10\n",
      "INFO    : TFIChain: reading 'bc_MPS'='infinite'\n",
      "INFO    : TFIChain: reading 'L'=2\n",
      "INFO    : TFIChain: reading 'J'=1.0\n",
      "INFO    : TFIChain: reading 'g'=0.4\n",
      "INFO    : TwoSiteDMRGEngine: reading 'diag_method'='default'\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'trunc_params'=Config(<5 options>, 'trunc_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'combine'=False\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'lanczos_params'=Config(<9 options>, 'lanczos_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'start_env'=1\n",
      "INFO    : start environment_sweep\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list_reactivates_mixer'=True\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_to_trunc'=0.05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_min'=4.999999999999999e-16\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_max'=0.0001\n",
      "INFO    : TwoSiteDMRGEngine: reading 'E_tol_to_trunc'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'N_sweeps_check'=10\n",
      "INFO    : TwoSiteDMRGEngine: reading 'min_sweeps'=15\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_sweeps'=40\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_E_err'=1e-08\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_S_err'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_hours'=8760\n",
      "INFO    : TwoSiteDMRGEngine: reading 'update_env'=5\n",
      "INFO    : TwoSiteDMRGEngine: reading 'mixer'=None\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_N_for_ED'=400\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 10\n",
      "energy=-1.0404170862358413, max S=0.6943390143742775, age=128, norm_err=5.7e-16\n",
      "Current memory usage 131.5MB, wall time: 0.5s\n",
      "Delta E = nan, Delta S = 6.9434e-02 (per sweep)\n",
      "max trunc_err = 3.4098e-16, max E_trunc = 2.8422e-14\n",
      "chi: [10, 10]\n",
      "================================================================================\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 20\n",
      "energy=-1.0404170862358342, max S=0.6943390143742776, age=248, norm_err=4.3e-16\n",
      "Current memory usage 131.5MB, wall time: 0.5s\n",
      "Delta E = 7.1054e-16, Delta S = -8.8818e-17 (per sweep)\n",
      "max trunc_err = 3.4098e-16, max E_trunc = 8.5265e-14\n",
      "chi: [10, 10]\n",
      "================================================================================\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol_iter'=5\n",
      "INFO    : DMRG finished after 20 sweeps, max chi=10\n",
      "INFO    : TFIChain: reading 'bc_MPS'='infinite'\n",
      "INFO    : TFIChain: reading 'L'=2\n",
      "INFO    : TFIChain: reading 'J'=1.0\n",
      "INFO    : TFIChain: reading 'g'=0.5\n",
      "INFO    : TwoSiteDMRGEngine: reading 'diag_method'='default'\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'trunc_params'=Config(<5 options>, 'trunc_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'combine'=False\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'lanczos_params'=Config(<9 options>, 'lanczos_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'start_env'=1\n",
      "INFO    : start environment_sweep\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list_reactivates_mixer'=True\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_to_trunc'=0.05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_min'=4.999999999999999e-16\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_max'=0.0001\n",
      "INFO    : TwoSiteDMRGEngine: reading 'E_tol_to_trunc'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'N_sweeps_check'=10\n",
      "INFO    : TwoSiteDMRGEngine: reading 'min_sweeps'=15\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_sweeps'=40\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_E_err'=1e-08\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_S_err'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_hours'=8760\n",
      "INFO    : TwoSiteDMRGEngine: reading 'update_env'=5\n",
      "INFO    : TwoSiteDMRGEngine: reading 'mixer'=None\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_N_for_ED'=400\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 10\n",
      "energy=-1.0635444099732894, max S=0.6960673544735364, age=128, norm_err=8.6e-16\n",
      "Current memory usage 131.5MB, wall time: 0.7s\n",
      "Delta E = nan, Delta S = 6.9607e-02 (per sweep)\n",
      "max trunc_err = 1.8068e-14, max E_trunc = 2.1316e-13\n",
      "chi: [10, 10]\n",
      "================================================================================\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 20\n",
      "energy=-1.0635444099733178, max S=0.6960673544735387, age=248, norm_err=8.0e-16\n",
      "Current memory usage 131.5MB, wall time: 0.5s\n",
      "Delta E = -2.8422e-15, Delta S = 3.5527e-16 (per sweep)\n",
      "max trunc_err = 1.8068e-14, max E_trunc = 2.5580e-13\n",
      "chi: [10, 10]\n",
      "================================================================================\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol_iter'=5\n",
      "INFO    : DMRG finished after 20 sweeps, max chi=10\n",
      "INFO    : TFIChain: reading 'bc_MPS'='infinite'\n",
      "INFO    : TFIChain: reading 'L'=2\n",
      "INFO    : TFIChain: reading 'J'=1.0\n",
      "INFO    : TFIChain: reading 'g'=0.6000000000000001\n",
      "INFO    : TwoSiteDMRGEngine: reading 'diag_method'='default'\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'trunc_params'=Config(<5 options>, 'trunc_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'combine'=False\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'lanczos_params'=Config(<9 options>, 'lanczos_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'start_env'=1\n",
      "INFO    : start environment_sweep\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list_reactivates_mixer'=True\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_to_trunc'=0.05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_min'=4.999999999999999e-16\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_max'=0.0001\n",
      "INFO    : TwoSiteDMRGEngine: reading 'E_tol_to_trunc'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'N_sweeps_check'=10\n",
      "INFO    : TwoSiteDMRGEngine: reading 'min_sweeps'=15\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_sweeps'=40\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_E_err'=1e-08\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_S_err'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_hours'=8760\n",
      "INFO    : TwoSiteDMRGEngine: reading 'update_env'=5\n",
      "INFO    : TwoSiteDMRGEngine: reading 'mixer'=None\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_N_for_ED'=400\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 10\n",
      "energy=-1.0922385835547068, max S=0.6994879401141262, age=128, norm_err=1.4e-13\n",
      "Current memory usage 131.5MB, wall time: 0.6s\n",
      "Delta E = nan, Delta S = 6.9949e-02 (per sweep)\n",
      "max trunc_err = 6.6755e-16, max E_trunc = 0.0000e+00\n",
      "chi: [14, 14]\n",
      "================================================================================\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "INFO    : Running sweep with optimization\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 20\n",
      "energy=-1.0922385835547885, max S=0.6994879401141427, age=248, norm_err=2.5e-15\n",
      "Current memory usage 131.5MB, wall time: 0.5s\n",
      "Delta E = -8.1712e-15, Delta S = 1.7097e-15 (per sweep)\n",
      "max trunc_err = 6.6754e-16, max E_trunc = 5.6843e-14\n",
      "chi: [14, 14]\n",
      "================================================================================\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol_iter'=5\n",
      "INFO    : DMRG finished after 20 sweeps, max chi=14\n",
      "INFO    : TFIChain: reading 'bc_MPS'='infinite'\n",
      "INFO    : TFIChain: reading 'L'=2\n",
      "INFO    : TFIChain: reading 'J'=1.0\n",
      "INFO    : TFIChain: reading 'g'=0.7000000000000001\n",
      "INFO    : TwoSiteDMRGEngine: reading 'diag_method'='default'\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'trunc_params'=Config(<5 options>, 'trunc_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'combine'=False\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'lanczos_params'=Config(<9 options>, 'lanczos_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'start_env'=1\n",
      "INFO    : start environment_sweep\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list_reactivates_mixer'=True\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_to_trunc'=0.05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_min'=4.999999999999999e-16\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_max'=0.0001\n",
      "INFO    : TwoSiteDMRGEngine: reading 'E_tol_to_trunc'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'N_sweeps_check'=10\n",
      "INFO    : TwoSiteDMRGEngine: reading 'min_sweeps'=15\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_sweeps'=40\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_E_err'=1e-08\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_S_err'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_hours'=8760\n",
      "INFO    : TwoSiteDMRGEngine: reading 'update_env'=5\n",
      "INFO    : TwoSiteDMRGEngine: reading 'mixer'=None\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_N_for_ED'=400\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 10\n",
      "energy=-1.1268286675867394, max S=0.7061087635509824, age=128, norm_err=3.0e-10\n",
      "Current memory usage 131.5MB, wall time: 0.9s\n",
      "Delta E = nan, Delta S = 7.0611e-02 (per sweep)\n",
      "max trunc_err = 3.7402e-14, max E_trunc = 2.9843e-13\n",
      "chi: [14, 14]\n",
      "================================================================================\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 20\n",
      "energy=-1.1268286675867429, max S=0.7061087635493897, age=248, norm_err=5.0e-15\n",
      "Current memory usage 131.5MB, wall time: 0.6s\n",
      "Delta E = -3.5527e-16, Delta S = -1.6247e-13 (per sweep)\n",
      "max trunc_err = 3.7402e-14, max E_trunc = 3.6948e-13\n",
      "chi: [14, 14]\n",
      "================================================================================\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol_iter'=5\n",
      "INFO    : DMRG finished after 20 sweeps, max chi=14\n",
      "INFO    : TFIChain: reading 'bc_MPS'='infinite'\n",
      "INFO    : TFIChain: reading 'L'=2\n",
      "INFO    : TFIChain: reading 'J'=1.0\n",
      "INFO    : TFIChain: reading 'g'=0.8\n",
      "INFO    : TwoSiteDMRGEngine: reading 'diag_method'='default'\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'trunc_params'=Config(<5 options>, 'trunc_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'combine'=False\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'lanczos_params'=Config(<9 options>, 'lanczos_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'start_env'=1\n",
      "INFO    : start environment_sweep\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list_reactivates_mixer'=True\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_to_trunc'=0.05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_min'=4.999999999999999e-16\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_max'=0.0001\n",
      "INFO    : TwoSiteDMRGEngine: reading 'E_tol_to_trunc'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'N_sweeps_check'=10\n",
      "INFO    : TwoSiteDMRGEngine: reading 'min_sweeps'=15\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_sweeps'=40\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_E_err'=1e-08\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_S_err'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_hours'=8760\n",
      "INFO    : TwoSiteDMRGEngine: reading 'update_env'=5\n",
      "INFO    : TwoSiteDMRGEngine: reading 'mixer'=None\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_N_for_ED'=400\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 10\n",
      "energy=-1.1678095085200795, max S=0.7195176070062523, age=128, norm_err=6.2e-07\n",
      "Current memory usage 131.5MB, wall time: 1.3s\n",
      "Delta E = nan, Delta S = 7.1952e-02 (per sweep)\n",
      "max trunc_err = 1.2023e-14, max E_trunc = 1.7053e-13\n",
      "chi: [20, 20]\n",
      "================================================================================\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 20\n",
      "energy=-1.1678095085204987, max S=0.7195176070786253, age=248, norm_err=1.6e-12\n",
      "Current memory usage 131.5MB, wall time: 0.7s\n",
      "Delta E = -4.1922e-14, Delta S = 1.3946e-11 (per sweep)\n",
      "max trunc_err = 1.1706e-14, max E_trunc = 3.1264e-13\n",
      "chi: [20, 20]\n",
      "================================================================================\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol_iter'=5\n",
      "INFO    : DMRG finished after 20 sweeps, max chi=20\n",
      "INFO    : TFIChain: reading 'bc_MPS'='infinite'\n",
      "INFO    : TFIChain: reading 'L'=2\n",
      "INFO    : TFIChain: reading 'J'=1.0\n",
      "INFO    : TFIChain: reading 'g'=0.9\n",
      "INFO    : TwoSiteDMRGEngine: reading 'diag_method'='default'\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'trunc_params'=Config(<5 options>, 'trunc_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'combine'=False\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'lanczos_params'=Config(<9 options>, 'lanczos_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'start_env'=1\n",
      "INFO    : start environment_sweep\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list_reactivates_mixer'=True\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_to_trunc'=0.05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_min'=4.999999999999999e-16\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_max'=0.0001\n",
      "INFO    : TwoSiteDMRGEngine: reading 'E_tol_to_trunc'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'N_sweeps_check'=10\n",
      "INFO    : TwoSiteDMRGEngine: reading 'min_sweeps'=15\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_sweeps'=40\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_E_err'=1e-08\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_S_err'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_hours'=8760\n",
      "INFO    : TwoSiteDMRGEngine: reading 'update_env'=5\n",
      "INFO    : TwoSiteDMRGEngine: reading 'mixer'=None\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_N_for_ED'=400\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 10\n",
      "energy=-1.2160008125853352, max S=0.7519468029822860, age=128, norm_err=3.4e-04\n",
      "Current memory usage 131.5MB, wall time: 1.7s\n",
      "Delta E = nan, Delta S = 7.5195e-02 (per sweep)\n",
      "max trunc_err = 2.9205e-14, max E_trunc = 2.2737e-13\n",
      "chi: [28, 28]\n",
      "================================================================================\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 20\n",
      "energy=-1.2160009141093227, max S=0.7519576065770814, age=248, norm_err=6.2e-07\n",
      "Current memory usage 131.5MB, wall time: 0.8s\n",
      "Delta E = -1.0152e-08, Delta S = 1.1780e-06 (per sweep)\n",
      "max trunc_err = 3.7089e-14, max E_trunc = 3.6948e-13\n",
      "chi: [28, 28]\n",
      "================================================================================\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol_iter'=5\n",
      "INFO    : DMRG finished after 20 sweeps, max chi=28\n",
      "INFO    : TFIChain: reading 'bc_MPS'='infinite'\n",
      "INFO    : TFIChain: reading 'L'=2\n",
      "INFO    : TFIChain: reading 'J'=1.0\n",
      "INFO    : TFIChain: reading 'g'=1.0\n",
      "INFO    : TwoSiteDMRGEngine: reading 'diag_method'='default'\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'trunc_params'=Config(<5 options>, 'trunc_params')\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "INFO    : TwoSiteDMRGEngine: reading 'combine'=False\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'lanczos_params'=Config(<9 options>, 'lanczos_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'start_env'=1\n",
      "INFO    : start environment_sweep\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list_reactivates_mixer'=True\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_to_trunc'=0.05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_min'=4.999999999999999e-16\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_max'=0.0001\n",
      "INFO    : TwoSiteDMRGEngine: reading 'E_tol_to_trunc'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'N_sweeps_check'=10\n",
      "INFO    : TwoSiteDMRGEngine: reading 'min_sweeps'=15\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_sweeps'=40\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_E_err'=1e-08\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_S_err'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_hours'=8760\n",
      "INFO    : TwoSiteDMRGEngine: reading 'update_env'=5\n",
      "INFO    : TwoSiteDMRGEngine: reading 'mixer'=None\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_N_for_ED'=400\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 10\n",
      "energy=-1.2732297890099140, max S=0.5971035592084446, age=128, norm_err=1.1e-03\n",
      "Current memory usage 132.0MB, wall time: 1.8s\n",
      "Delta E = nan, Delta S = 5.9652e-02 (per sweep)\n",
      "max trunc_err = 1.7558e-13, max E_trunc = 1.2363e-12\n",
      "chi: [30, 30]\n",
      "================================================================================\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 20\n",
      "energy=-1.2732371980491379, max S=0.6572441210794768, age=248, norm_err=6.3e-04\n",
      "Current memory usage 132.0MB, wall time: 2.0s\n",
      "Delta E = -7.4090e-07, Delta S = 6.0398e-03 (per sweep)\n",
      "max trunc_err = 1.5790e-11, max E_trunc = 1.0886e-10\n",
      "chi: [30, 30]\n",
      "================================================================================\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 30\n",
      "energy=-1.2732385090977729, max S=0.6913505799955206, age=368, norm_err=4.2e-04\n",
      "Current memory usage 132.0MB, wall time: 1.4s\n",
      "Delta E = -1.3110e-07, Delta S = 3.4213e-03 (per sweep)\n",
      "max trunc_err = 7.7968e-11, max E_trunc = 5.2478e-10\n",
      "chi: [30, 30]\n",
      "================================================================================\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 40\n",
      "energy=-1.2732389601518435, max S=0.7152926339808676, age=488, norm_err=3.1e-04\n",
      "Current memory usage 132.0MB, wall time: 1.4s\n",
      "Delta E = -4.5105e-08, Delta S = 2.3998e-03 (per sweep)\n",
      "max trunc_err = 1.8960e-10, max E_trunc = 1.2575e-09\n",
      "chi: [30, 30]\n",
      "================================================================================\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol_iter'=5\n",
      "WARNING : final DMRG state not in canonical form up to norm_tol=1.00e-05: norm_err=3.12e-04\n",
      "INFO    : start environment_sweep\n",
      "INFO    : start environment_sweep\n",
      "INFO    : start environment_sweep\n",
      "INFO    : start environment_sweep\n",
      "INFO    : start environment_sweep\n",
      "WARNING : norm_err=1.20e-04 still too high after environment_sweeps, call psi.canonical_form()\n",
      "INFO    : DMRG finished after 40 sweeps, max chi=30\n",
      "INFO    : TFIChain: reading 'bc_MPS'='infinite'\n",
      "INFO    : TFIChain: reading 'L'=2\n",
      "INFO    : TFIChain: reading 'J'=1.0\n",
      "INFO    : TFIChain: reading 'g'=1.1\n",
      "INFO    : TwoSiteDMRGEngine: reading 'diag_method'='default'\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'trunc_params'=Config(<5 options>, 'trunc_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'combine'=False\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'lanczos_params'=Config(<9 options>, 'lanczos_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'start_env'=1\n",
      "INFO    : start environment_sweep\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list_reactivates_mixer'=True\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_to_trunc'=0.05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_min'=4.999999999999999e-16\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_max'=0.0001\n",
      "INFO    : TwoSiteDMRGEngine: reading 'E_tol_to_trunc'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'N_sweeps_check'=10\n",
      "INFO    : TwoSiteDMRGEngine: reading 'min_sweeps'=15\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_sweeps'=40\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_E_err'=1e-08\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_S_err'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_hours'=8760\n",
      "INFO    : TwoSiteDMRGEngine: reading 'update_env'=5\n",
      "INFO    : TwoSiteDMRGEngine: reading 'mixer'=None\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_N_for_ED'=400\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 10\n",
      "energy=-1.3428640227248607, max S=0.3288585829175101, age=128, norm_err=1.2e-08\n",
      "Current memory usage 132.0MB, wall time: 0.6s\n",
      "Delta E = nan, Delta S = 3.2886e-02 (per sweep)\n",
      "max trunc_err = 1.3053e-14, max E_trunc = 1.7053e-13\n",
      "chi: [23, 23]\n",
      "================================================================================\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 20\n",
      "energy=-1.3428640227249247, max S=0.3288579907169442, age=248, norm_err=1.6e-12\n",
      "Current memory usage 132.0MB, wall time: 0.7s\n",
      "Delta E = -6.3949e-15, Delta S = -5.8017e-08 (per sweep)\n",
      "max trunc_err = 1.5356e-14, max E_trunc = 2.8422e-13\n",
      "chi: [23, 23]\n",
      "================================================================================\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol_iter'=5\n",
      "INFO    : DMRG finished after 20 sweeps, max chi=23\n",
      "INFO    : TFIChain: reading 'bc_MPS'='infinite'\n",
      "INFO    : TFIChain: reading 'L'=2\n",
      "INFO    : TFIChain: reading 'J'=1.0\n",
      "INFO    : TFIChain: reading 'g'=1.2000000000000002\n",
      "INFO    : TwoSiteDMRGEngine: reading 'diag_method'='default'\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'trunc_params'=Config(<5 options>, 'trunc_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'combine'=False\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'lanczos_params'=Config(<9 options>, 'lanczos_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'start_env'=1\n",
      "INFO    : start environment_sweep\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list_reactivates_mixer'=True\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_to_trunc'=0.05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_min'=4.999999999999999e-16\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_max'=0.0001\n",
      "INFO    : TwoSiteDMRGEngine: reading 'E_tol_to_trunc'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'N_sweeps_check'=10\n",
      "INFO    : TwoSiteDMRGEngine: reading 'min_sweeps'=15\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_sweeps'=40\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_E_err'=1e-08\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_S_err'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_hours'=8760\n",
      "INFO    : TwoSiteDMRGEngine: reading 'update_env'=5\n",
      "INFO    : TwoSiteDMRGEngine: reading 'mixer'=None\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_N_for_ED'=400\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 10\n",
      "energy=-1.4196192748981424, max S=0.2533554574789748, age=128, norm_err=4.5e-13\n",
      "Current memory usage 132.0MB, wall time: 0.9s\n",
      "Delta E = nan, Delta S = 2.5336e-02 (per sweep)\n",
      "max trunc_err = 3.6208e-15, max E_trunc = 4.2633e-14\n",
      "chi: [18, 18]\n",
      "================================================================================\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 20\n",
      "energy=-1.4196192748981886, max S=0.2533554571105219, age=248, norm_err=9.8e-15\n",
      "Current memory usage 132.0MB, wall time: 0.6s\n",
      "Delta E = -4.6185e-15, Delta S = -3.6706e-11 (per sweep)\n",
      "max trunc_err = 3.7398e-15, max E_trunc = 1.1369e-13\n",
      "chi: [18, 18]\n",
      "================================================================================\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol_iter'=5\n",
      "INFO    : DMRG finished after 20 sweeps, max chi=18\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "INFO    : TFIChain: reading 'bc_MPS'='infinite'\n",
      "INFO    : TFIChain: reading 'L'=2\n",
      "INFO    : TFIChain: reading 'J'=1.0\n",
      "INFO    : TFIChain: reading 'g'=1.3\n",
      "INFO    : TwoSiteDMRGEngine: reading 'diag_method'='default'\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'trunc_params'=Config(<5 options>, 'trunc_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'combine'=False\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'lanczos_params'=Config(<9 options>, 'lanczos_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'start_env'=1\n",
      "INFO    : start environment_sweep\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list_reactivates_mixer'=True\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_to_trunc'=0.05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_min'=4.999999999999999e-16\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_max'=0.0001\n",
      "INFO    : TwoSiteDMRGEngine: reading 'E_tol_to_trunc'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'N_sweeps_check'=10\n",
      "INFO    : TwoSiteDMRGEngine: reading 'min_sweeps'=15\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_sweeps'=40\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_E_err'=1e-08\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_S_err'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_hours'=8760\n",
      "INFO    : TwoSiteDMRGEngine: reading 'update_env'=5\n",
      "INFO    : TwoSiteDMRGEngine: reading 'mixer'=None\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_N_for_ED'=400\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 10\n",
      "energy=-1.5008232436750149, max S=0.2083921498315460, age=128, norm_err=3.9e-16\n",
      "Current memory usage 132.0MB, wall time: 0.6s\n",
      "Delta E = nan, Delta S = 2.0839e-02 (per sweep)\n",
      "max trunc_err = 7.1131e-15, max E_trunc = 1.7053e-13\n",
      "chi: [14, 14]\n",
      "================================================================================\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 20\n",
      "energy=-1.5008232436750006, max S=0.2083921498317719, age=248, norm_err=9.5e-16\n",
      "Current memory usage 132.0MB, wall time: 0.6s\n",
      "Delta E = 1.4211e-15, Delta S = 2.2982e-14 (per sweep)\n",
      "max trunc_err = 7.1131e-15, max E_trunc = 3.4106e-13\n",
      "chi: [14, 14]\n",
      "================================================================================\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol_iter'=5\n",
      "INFO    : DMRG finished after 20 sweeps, max chi=14\n",
      "INFO    : TFIChain: reading 'bc_MPS'='infinite'\n",
      "INFO    : TFIChain: reading 'L'=2\n",
      "INFO    : TFIChain: reading 'J'=1.0\n",
      "INFO    : TFIChain: reading 'g'=1.4000000000000001\n",
      "INFO    : TwoSiteDMRGEngine: reading 'diag_method'='default'\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'trunc_params'=Config(<5 options>, 'trunc_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'combine'=False\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'lanczos_params'=Config(<9 options>, 'lanczos_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'start_env'=1\n",
      "INFO    : start environment_sweep\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list_reactivates_mixer'=True\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_to_trunc'=0.05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_min'=4.999999999999999e-16\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_max'=0.0001\n",
      "INFO    : TwoSiteDMRGEngine: reading 'E_tol_to_trunc'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'N_sweeps_check'=10\n",
      "INFO    : TwoSiteDMRGEngine: reading 'min_sweeps'=15\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_sweeps'=40\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_E_err'=1e-08\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_S_err'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_hours'=8760\n",
      "INFO    : TwoSiteDMRGEngine: reading 'update_env'=5\n",
      "INFO    : TwoSiteDMRGEngine: reading 'mixer'=None\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_N_for_ED'=400\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 10\n",
      "energy=-1.5851883000539750, max S=0.1770262953951441, age=128, norm_err=7.0e-16\n",
      "Current memory usage 132.2MB, wall time: 0.6s\n",
      "Delta E = nan, Delta S = 1.7703e-02 (per sweep)\n",
      "max trunc_err = 1.7903e-15, max E_trunc = 1.1369e-13\n",
      "chi: [13, 13]\n",
      "================================================================================\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 20\n",
      "energy=-1.5851883000539573, max S=0.1770262953895704, age=248, norm_err=3.4e-16\n",
      "Current memory usage 132.2MB, wall time: 0.5s\n",
      "Delta E = 1.7764e-15, Delta S = -5.5357e-13 (per sweep)\n",
      "max trunc_err = 1.7903e-15, max E_trunc = 3.4106e-13\n",
      "chi: [13, 13]\n",
      "================================================================================\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol_iter'=5\n",
      "INFO    : DMRG finished after 20 sweeps, max chi=13\n",
      "INFO    : TFIChain: reading 'bc_MPS'='infinite'\n",
      "INFO    : TFIChain: reading 'L'=2\n",
      "INFO    : TFIChain: reading 'J'=1.0\n",
      "INFO    : TFIChain: reading 'g'=1.5\n",
      "INFO    : TwoSiteDMRGEngine: reading 'diag_method'='default'\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'trunc_params'=Config(<5 options>, 'trunc_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'combine'=False\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'lanczos_params'=Config(<9 options>, 'lanczos_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'start_env'=1\n",
      "INFO    : start environment_sweep\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list_reactivates_mixer'=True\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_to_trunc'=0.05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_min'=4.999999999999999e-16\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_max'=0.0001\n",
      "INFO    : TwoSiteDMRGEngine: reading 'E_tol_to_trunc'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'N_sweeps_check'=10\n",
      "INFO    : TwoSiteDMRGEngine: reading 'min_sweeps'=15\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_sweeps'=40\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_E_err'=1e-08\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_S_err'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_hours'=8760\n",
      "INFO    : TwoSiteDMRGEngine: reading 'update_env'=5\n",
      "INFO    : TwoSiteDMRGEngine: reading 'mixer'=None\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_N_for_ED'=400\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 10\n",
      "energy=-1.6719262215361468, max S=0.1534912555174990, age=128, norm_err=5.0e-16\n",
      "Current memory usage 132.5MB, wall time: 0.6s\n",
      "Delta E = nan, Delta S = 1.5349e-02 (per sweep)\n",
      "max trunc_err = 7.5713e-15, max E_trunc = 1.9895e-13\n",
      "chi: [11, 11]\n",
      "================================================================================\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 20\n",
      "energy=-1.6719262215361752, max S=0.1534912555174835, age=248, norm_err=8.3e-16\n",
      "Current memory usage 132.5MB, wall time: 0.6s\n",
      "Delta E = -2.8422e-15, Delta S = -2.1955e-15 (per sweep)\n",
      "max trunc_err = 7.5713e-15, max E_trunc = 2.8422e-13\n",
      "chi: [11, 11]\n",
      "================================================================================\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol_iter'=5\n",
      "INFO    : DMRG finished after 20 sweeps, max chi=11\n",
      "INFO    : TFIChain: reading 'bc_MPS'='infinite'\n",
      "INFO    : TFIChain: reading 'L'=2\n",
      "INFO    : TFIChain: reading 'J'=1.0\n",
      "INFO    : TFIChain: reading 'g'=1.6\n",
      "INFO    : TwoSiteDMRGEngine: reading 'diag_method'='default'\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'trunc_params'=Config(<5 options>, 'trunc_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'combine'=False\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'lanczos_params'=Config(<9 options>, 'lanczos_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'start_env'=1\n",
      "INFO    : start environment_sweep\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list_reactivates_mixer'=True\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_to_trunc'=0.05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_min'=4.999999999999999e-16\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_max'=0.0001\n",
      "INFO    : TwoSiteDMRGEngine: reading 'E_tol_to_trunc'=None\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "INFO    : TwoSiteDMRGEngine: reading 'N_sweeps_check'=10\n",
      "INFO    : TwoSiteDMRGEngine: reading 'min_sweeps'=15\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_sweeps'=40\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_E_err'=1e-08\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_S_err'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_hours'=8760\n",
      "INFO    : TwoSiteDMRGEngine: reading 'update_env'=5\n",
      "INFO    : TwoSiteDMRGEngine: reading 'mixer'=None\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_N_for_ED'=400\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 10\n",
      "energy=-1.7605081222046017, max S=0.1350487935034509, age=128, norm_err=5.1e-16\n",
      "Current memory usage 132.5MB, wall time: 0.5s\n",
      "Delta E = nan, Delta S = 1.3505e-02 (per sweep)\n",
      "max trunc_err = 1.4088e-15, max E_trunc = 1.7053e-13\n",
      "chi: [11, 11]\n",
      "================================================================================\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 20\n",
      "energy=-1.7605081222045982, max S=0.1350487935032368, age=248, norm_err=2.6e-16\n",
      "Current memory usage 132.5MB, wall time: 0.5s\n",
      "Delta E = 3.5527e-16, Delta S = -2.1552e-14 (per sweep)\n",
      "max trunc_err = 1.4088e-15, max E_trunc = 1.7053e-13\n",
      "chi: [11, 11]\n",
      "================================================================================\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol_iter'=5\n",
      "INFO    : DMRG finished after 20 sweeps, max chi=11\n",
      "INFO    : TFIChain: reading 'bc_MPS'='infinite'\n",
      "INFO    : TFIChain: reading 'L'=2\n",
      "INFO    : TFIChain: reading 'J'=1.0\n",
      "INFO    : TFIChain: reading 'g'=1.7000000000000002\n",
      "INFO    : TwoSiteDMRGEngine: reading 'diag_method'='default'\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'trunc_params'=Config(<5 options>, 'trunc_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'combine'=False\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'lanczos_params'=Config(<9 options>, 'lanczos_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'start_env'=1\n",
      "INFO    : start environment_sweep\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list_reactivates_mixer'=True\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_to_trunc'=0.05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_min'=4.999999999999999e-16\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_max'=0.0001\n",
      "INFO    : TwoSiteDMRGEngine: reading 'E_tol_to_trunc'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'N_sweeps_check'=10\n",
      "INFO    : TwoSiteDMRGEngine: reading 'min_sweeps'=15\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_sweeps'=40\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_E_err'=1e-08\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_S_err'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_hours'=8760\n",
      "INFO    : TwoSiteDMRGEngine: reading 'update_env'=5\n",
      "INFO    : TwoSiteDMRGEngine: reading 'mixer'=None\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_N_for_ED'=400\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 10\n",
      "energy=-1.8505593362793462, max S=0.1201649555165585, age=128, norm_err=3.2e-16\n",
      "Current memory usage 132.5MB, wall time: 0.4s\n",
      "Delta E = nan, Delta S = 1.2016e-02 (per sweep)\n",
      "max trunc_err = 1.0356e-14, max E_trunc = 2.5580e-13\n",
      "chi: [9, 9]\n",
      "================================================================================\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 20\n",
      "energy=-1.8505593362793249, max S=0.1201649555165666, age=248, norm_err=2.2e-16\n",
      "Current memory usage 132.5MB, wall time: 0.4s\n",
      "Delta E = 2.1316e-15, Delta S = 8.0491e-16 (per sweep)\n",
      "max trunc_err = 1.0356e-14, max E_trunc = 4.5475e-13\n",
      "chi: [9, 9]\n",
      "================================================================================\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol_iter'=5\n",
      "INFO    : DMRG finished after 20 sweeps, max chi=9\n",
      "INFO    : TFIChain: reading 'bc_MPS'='infinite'\n",
      "INFO    : TFIChain: reading 'L'=2\n",
      "INFO    : TFIChain: reading 'J'=1.0\n",
      "INFO    : TFIChain: reading 'g'=1.8\n",
      "INFO    : TwoSiteDMRGEngine: reading 'diag_method'='default'\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'trunc_params'=Config(<5 options>, 'trunc_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'combine'=False\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'lanczos_params'=Config(<9 options>, 'lanczos_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'start_env'=1\n",
      "INFO    : start environment_sweep\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list_reactivates_mixer'=True\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_to_trunc'=0.05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_min'=4.999999999999999e-16\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_max'=0.0001\n",
      "INFO    : TwoSiteDMRGEngine: reading 'E_tol_to_trunc'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'N_sweeps_check'=10\n",
      "INFO    : TwoSiteDMRGEngine: reading 'min_sweeps'=15\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_sweeps'=40\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_E_err'=1e-08\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_S_err'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_hours'=8760\n",
      "INFO    : TwoSiteDMRGEngine: reading 'update_env'=5\n",
      "INFO    : TwoSiteDMRGEngine: reading 'mixer'=None\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_N_for_ED'=400\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 10\n",
      "energy=-1.9418043043913684, max S=0.1078900285324606, age=128, norm_err=2.7e-16\n",
      "Current memory usage 132.5MB, wall time: 0.4s\n",
      "Delta E = nan, Delta S = 1.0789e-02 (per sweep)\n",
      "max trunc_err = 3.0124e-15, max E_trunc = 1.1369e-13\n",
      "chi: [9, 9]\n",
      "================================================================================\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 20\n",
      "energy=-1.9418043043913968, max S=0.1078900285324698, age=248, norm_err=2.7e-16\n",
      "Current memory usage 132.5MB, wall time: 0.4s\n",
      "Delta E = -2.8422e-15, Delta S = 4.5797e-16 (per sweep)\n",
      "max trunc_err = 3.0124e-15, max E_trunc = 2.8422e-13\n",
      "chi: [9, 9]\n",
      "================================================================================\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol_iter'=5\n",
      "INFO    : DMRG finished after 20 sweeps, max chi=9\n",
      "INFO    : TFIChain: reading 'bc_MPS'='infinite'\n",
      "INFO    : TFIChain: reading 'L'=2\n",
      "INFO    : TFIChain: reading 'J'=1.0\n",
      "INFO    : TFIChain: reading 'g'=1.9000000000000001\n",
      "INFO    : TwoSiteDMRGEngine: reading 'diag_method'='default'\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'trunc_params'=Config(<5 options>, 'trunc_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'combine'=False\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'lanczos_params'=Config(<9 options>, 'lanczos_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'start_env'=1\n",
      "INFO    : start environment_sweep\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list_reactivates_mixer'=True\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_to_trunc'=0.05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_min'=4.999999999999999e-16\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_max'=0.0001\n",
      "INFO    : TwoSiteDMRGEngine: reading 'E_tol_to_trunc'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'N_sweeps_check'=10\n",
      "INFO    : TwoSiteDMRGEngine: reading 'min_sweeps'=15\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_sweeps'=40\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_E_err'=1e-08\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_S_err'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_hours'=8760\n",
      "INFO    : TwoSiteDMRGEngine: reading 'update_env'=5\n",
      "INFO    : TwoSiteDMRGEngine: reading 'mixer'=None\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_N_for_ED'=400\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 10\n",
      "energy=-2.0340345587409665, max S=0.0975950958265605, age=128, norm_err=1.8e-15\n",
      "Current memory usage 132.5MB, wall time: 0.4s\n",
      "Delta E = nan, Delta S = 9.7595e-03 (per sweep)\n",
      "max trunc_err = 9.6326e-16, max E_trunc = 8.5265e-14\n",
      "chi: [9, 9]\n",
      "================================================================================\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "INFO    : Running sweep with optimization\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 20\n",
      "energy=-2.0340345587409985, max S=0.0975950958265605, age=248, norm_err=1.5e-15\n",
      "Current memory usage 132.5MB, wall time: 0.4s\n",
      "Delta E = -3.1974e-15, Delta S = -4.3437e-16 (per sweep)\n",
      "max trunc_err = 9.6326e-16, max E_trunc = -1.1369e-13\n",
      "chi: [9, 9]\n",
      "================================================================================\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol_iter'=5\n",
      "INFO    : DMRG finished after 20 sweeps, max chi=9\n",
      "INFO    : TFIChain: reading 'bc_MPS'='infinite'\n",
      "INFO    : TFIChain: reading 'L'=2\n",
      "INFO    : TFIChain: reading 'J'=1.0\n",
      "INFO    : TFIChain: reading 'g'=2.0\n",
      "INFO    : TwoSiteDMRGEngine: reading 'diag_method'='default'\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'trunc_params'=Config(<5 options>, 'trunc_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'combine'=False\n",
      "INFO    : TwoSiteDMRGEngine: subconfig 'lanczos_params'=Config(<9 options>, 'lanczos_params')\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'start_env'=1\n",
      "INFO    : start environment_sweep\n",
      "INFO    : TwoSiteDMRGEngine: reading 'chi_list_reactivates_mixer'=True\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_to_trunc'=0.05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_min'=4.999999999999999e-16\n",
      "INFO    : TwoSiteDMRGEngine: reading 'P_tol_max'=0.0001\n",
      "INFO    : TwoSiteDMRGEngine: reading 'E_tol_to_trunc'=None\n",
      "INFO    : TwoSiteDMRGEngine: reading 'N_sweeps_check'=10\n",
      "INFO    : TwoSiteDMRGEngine: reading 'min_sweeps'=15\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_sweeps'=40\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_E_err'=1e-08\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_S_err'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_hours'=8760\n",
      "INFO    : TwoSiteDMRGEngine: reading 'update_env'=5\n",
      "INFO    : TwoSiteDMRGEngine: reading 'mixer'=None\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : TwoSiteDMRGEngine: reading 'max_N_for_ED'=400\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 10\n",
      "energy=-2.1270888199466818, max S=0.0888429735319253, age=128, norm_err=5.7e-16\n",
      "Current memory usage 132.5MB, wall time: 0.4s\n",
      "Delta E = nan, Delta S = 8.8843e-03 (per sweep)\n",
      "max trunc_err = 3.3343e-16, max E_trunc = 1.1369e-13\n",
      "chi: [9, 9]\n",
      "================================================================================\n",
      "INFO    : Running sweep with optimization\n",
      "INFO    : start environment_sweep\n",
      "INFO    : checkpoint after sweep 20\n",
      "energy=-2.1270888199466711, max S=0.0888429735319024, age=248, norm_err=8.8e-16\n",
      "Current memory usage 132.5MB, wall time: 0.4s\n",
      "Delta E = 1.0658e-15, Delta S = -1.3864e-15 (per sweep)\n",
      "max trunc_err = 3.3343e-16, max E_trunc = 0.0000e+00\n",
      "chi: [9, 9]\n",
      "================================================================================\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol'=1e-05\n",
      "INFO    : TwoSiteDMRGEngine: reading 'norm_tol_iter'=5\n",
      "INFO    : DMRG finished after 20 sweeps, max chi=9\n"
     ]
    }
   ],
   "source": [
    "#gs = [0.1, 0.5, 1.0, 1.5, 2.0]\n",
    "gs = np.linspace(0., 2., 21)\n",
    "entropies = []\n",
    "vals_X = []\n",
    "vals_Z = []\n",
    "corrs_XX = []\n",
    "\n",
    "\n",
    "for g in gs:\n",
    "    model = TFIChain({'L': 2, 'J': 1., 'g': g, 'bc_MPS': 'infinite'})\n",
    "    psi = MPS.from_product_state(model.lat.mps_sites(), ['up', 'up'], bc=model.lat.bc_MPS)\n",
    "    eng = dmrg.TwoSiteDMRGEngine(psi, model, algorithm_params)\n",
    "    E, psi = eng.run()\n",
    "    entropies.append(np.max(psi.entanglement_entropy()))\n",
    "    vals_X.append(np.mean(psi.expectation_value(\"Sigmax\")))\n",
    "    vals_Z.append(np.mean(psi.expectation_value(\"Sigmaz\")))\n",
    "    corrs_XX.append(psi.correlation_function(\"Sigmax\", \"Sigmax\", [0], [20]).item())\n",
    "\n",
    "    "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "id": "4ce49bc6",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.legend.Legend at 0x7f15e7fb2250>"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAyMAAAIlCAYAAADG/SdsAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAABcSAAAXEgFnn9JSAADd1ElEQVR4nOzdd3iUxfbA8e/sZtMTSihJSELoHUQEBKQGARHQi6hIF1HsFQviBVGviKI/7PeqIFWwYEcBwQ5IERQk9JaQBEgCpJfN7vz+2GTT+yabhPN5nn123/edfecsCtmTmTOjtNYIIYQQQgghRHUzODsAIYQQQgghxOVJkhEhhBBCCCGEU0gyIoQQQgghhHAKSUaEEEIIIYQQTiHJiBBCCCGEEMIpJBkRQgghhBBCOIUkI0IIIYQQQginkGRECCGEEEII4RSSjAghhBBCCCGcQpIRIYQQQgghhFNIMiKEEEIIIYRwCklGhBBCCCGEEE7h4uwALhdKqbOAJxDp7FiEEEIIIYRwoGAgVWvtX943Kq11FcQjClJKJbq5ufm0atXK2aEIIYQQQgjhMMePHycjIyNJa+1b3vfKyEj1iWzVqlXHAwcOODsOIYQQQgghHKZTp06Eh4dXaPaP1IwIIYQQQgghnEKSESGEEEIIIYRTSDIihBBCCCGEcApJRoQQQgghhBBOIcmIEEIIIYQQwikkGRFCCCGEEEI4hSQjQgghhBBCCKeQZEQIIYQQQgjhFJKMCCGEEEIIIZxCkhEhhBBCCCGEU0gyIoQQQgghhHAKSUaEEEIIIYQQTlFrkxGlVA+l1FNKqc+VUlFKKa2USq/E/eorpRYrpU4rpTKyn19XStV3YNjV4sKKFZijo0ttZ46O5sKKFdUQkRBCCCGEEIXV2mQE+DewAPgXEFiZGyml/ICdwENAFvAlkAQ8COzKvl4rXFixgnMvLuD01GklJiTm6GhOT53GuRcXODQhkURICCGEEEKUVW1ORrYDzwGjAf9K3uv/gDbA50A7rfWtWuvOwJtAa+C1St6/2vgMHYopOBhzZGSxCUlOImKOjMQUHIzP0KEO6dvZiZAQQgghhKhdam0yorVeqLWep7X+Vmt9rqL3UUr5AxMBM3Cv1jorz+XHgVhgolKqaeUirh6mwECaL19WbEJSMBFpvnwZpsBKDSzZOTMREkIIIYQQtU+tTUYc6Dpsfw6/FkxqtNYZwDeAMbtdrVBcQlKViUhJ/eao6v5BpokJIYQQQtQmLs4OoAbolv28p5jre4DpedrVCjmJQc6X/2NDwuzXlMmEwdeX6DlzMLi6odzcUO5uGNzcUK5leJ39KOq1wbceIUs+IOKOGfaEpPnyZQDVkoice3EBF1auKvH+eZMigIZTpjg0DiGEEEIIUTaSjEBI9vOZYq6fKdCuREqpA8VcalWeoBwhJyHJm4gAaLOZjAPFhekgBgMoVTgR8vDAs2dPEtavx9SkCS6NG+OS/Wzw9UUpVeEufYYO5cLKVfmSoIIJiUwTE0IIIYSoOSQZAe/s59RirqcUaCfKwmot8rROSyPh88+LvKbc3GzJSZ4Exf7cuDEuTWzPxvr1i0xaCo4GFUxIqmOamBBCCCGEKDtJRiDnW60u5XqZaK07FXkT24hJx/Lcq7JyvnwDYDTani0WjH5++N11JwZPT3RGJjojHWtGhu11ejrWzAx0egY6I6P41xm255zXZGUVG0dZ6YwMzGfOYD5T3CCVjTKZik9amjTGf95czs6dV+3TxIQQQgghRPlIMmLbTwTAq5jrntnPydUQi8MUNQoAuV/IL65a7dAv5DorC52ZaU9SMiMjiX78CbLOns2XCBnq+eIzaBDW1FTM58+TFRtLVmwcmM1l78tsthfkl8YcGcmxsOypWFpj9PMj4Ln5uDRuXJGPWaoLK1bYVhUr5c/VHB1N0ubNUq8ihBBCiMuaJCMQkf0cVMz1oALtarySpiOVNI2pMpSLC8rFBYOnJ+boaGKenkPW2bNFJkKpe/bm61dbrVgSEsg6H0tWToJS8Dn7tS5H0mKncwe9LPHxRNw+HUwm3Fq0wK1Nm+xHa9zatMEUFIQyVGyROSmgF0IIIYQoH0lG4O/s5yuLuZ5zfl81xFJppdVFlFZXUZX9F9evMhhwadAAlwYNoF3bYu+ttcaakJBnRCXWlsAUlbRkZJQSqJmMI0fIOHIk32nl4YFbq1a4tbYlJ25tbcmKS9OmpRbXSwG9EEIIIUT5KK2LK5WoXZRSGsjQWruX830B2FbMygKCtdbn81xzAyIBP6CZ1vpsJeI70LFjx44HqngVq5zfzpdWF5H3S3HTp2c75LfzZSkQr44icq01GUePEnnXzNxpYlrbiupdXGzPxRTYF8fg45OboOQZSXHx88vXrqTPJwX0QgghhKiLOnXqRHh4eHhxtdMluWxGRpRS9wP3A19orWfnnNdaxyil1mDbhf0dpdT4PLuwvww0BlZVJhGpTjlJRWl1CzkjJI6sW0javLnUL9oFR2aqom4iKyaGM/fdX+w0MVNQEP7z5mK5lEDG0aP2R0mF89akJNL27iVt7958540NG9qSkzwjKUFvv8WZ++6XAnohhBBCiFLU2pERpdT1wL/znOqNbUWsnXnOPa+1Xp/d/llgHrBcaz2twL0aAX9g2wvkOLAb6AR0zj6+WmsdV8l4q2VkxNmcXcBdmZEJa0oKGSdOkHEkO0E5doyMo0fJOneu3HEYGzXCmpyMTk+3jcwoBVlZkogIIYQQos65XEdGGmNLQPJSBc6VackkrXWcUqonMB+4EfgXcA54C5intb5Q6WgvE2VNLkyBgdWaiOT0WVK9jMHLC48uXfDo0iXffS0JCWQcP54/STlyBMvFi8XGYonLk7taLPaXrs2bk7TlRzx79cKtTesKF8sLIYQQQtQFtXZkpLa5XEZGnKm662Wy4uOzp3gdy53udewY1qSk0t8MGBs0wLNnTzx79cKrdy9cW7eu1A70QgghhBDOUJmREUlGqokkI9XD2dPEtNZknTtHyh9/cG7BS1gTEsr8XmPDhrbkpHcvvHr1wrVVq3InJ87+/EIIIYS4/Fyu07SEKMSZ08QAW/JgtRL39jtYExLsozSWlBQibp+OJS4OZTKhs7Ly7X8CYLlwgaSNG0nauBEAo58fnr164tWrF569euHasmWJyYnscyKEEEKI2kaSESEcqLi6FRPQ4pOPc681a4bfvfeQefQYqTt3kn7wYOHkJD6epO83kPT9BsBWFO/Vyzaty7NXb1xbhOZLTmSfEyGEEELUNjJNq5rINK26rzL7rFgSEkj9809Sd+wkZddOMg4eKpScFGRs3Aivnr3w7N0bz149cQ0NJSsmRvY5EUIIIUS1kpqRWkCSkbrPkQX0lkuXSP3zT1J27CB15y4yDh0qtX+XJk1sq3S1bcPFNWvIiilmnxVJRIQQQgjhQJKM1AKSjFweqqqAPOviRVJ37yZ15y5Sd+4k4/Dh0t9kNNqWFZZ9ToQQQghRhaSAXYgaoqoK6F0aNMD32mvxvfZaIDs52bXLlpzs2EHG0aOF35Szv0nOs1LUHzcO5eFR5n6FEEIIIaqSjIxUExkZEVUp68IFUnftJnXHDlJ37STj6LHiG7u44N2vH76jrsdnyBAMXl7VF6gQQggh6hwZGRHiMufSsCG+w4fhO3wYAGkHDhB510ws8fGFG2dlkfzLLyT/8gvK3R3vwYOod/31eA0YgMHVtXoDF0IIIcRlTZIRIeoYc3Q0UQ8/giU+HlNwMCHLPiTj4EGin5qNNTk5X1udnm5fPtjg44PPtdfie/1IvHr3RrnIPw9CCCGEqFoGZwcghHCcopbvdW3WDJ+hQ2n59VeYgoMAMHh5oby9873XmpREwuefE3nHDI4OHMTZ554ndc8etNXqjI8ihBBCiCKsCl9FTHJMqe1ikmNYFb6qGiKqHElGhKgjSttHxBQYSPPlyzEFB2NNScGlfn38X3ge35EjCxW1W+LjufjRR5yeMJFjQ4dyftEi0g8epKQaswsrVmCOji5TnBdWrKj4BxVCCCEuU6vCV7Fw10Kmb5xeYkISkxzD9I3TWbhrYY1PSCQZEaKOSNq8udR9RGwJyTJMwcGYz5xBp6bS7LVXabv1dwIXLcJ7yBAwmfK9Jys6hvgPlnDyX2M5cf0oYt9+m8xTp/K1ydlj5fTUaSUmJDkJ07kXF0hCIoQQQpRTWEgYQd5BnEk+U2xCkpOInEk+Q5B3EGEhYU6ItOxkNa1qIqtpiergiH1OLAkJJP3wAwnr15O6YycUM03LvVMnfK+/Ht+R14HVWuHd54UQQojaZFX4KsJCwgjwDiixXUxyDFsitjCp4ySH9l8w2Vg6fKk9lpKuVSXZ9LAWkGRE1EZZsbEkfr+BxPXrSfv776IbKYVnjx54XXMNlz79BHNUdKFkQxIRIYQQdUHONKnSvujnTQqe7PlktSQkgFMSEZBkpFaQZETUdplnzpC4/jsS168n48iRohsZDChXV3R6OqZmzWi+0jYVSxIRIYQQdUFZRh6qYnTCqq1cTL9IbFossamxxKXFcSLhBJ8e+ZQUcwoKhcb2nb66ExGQZKRWkGRE1CUZR4+SsH49ieu/wxwZWXxDpcBgAItFEhEhhBB1giOnSWVZs4hPiycuLY7YtFjOp563v45LtT3HpsVyIe0CWTqr1NiMysj3Y7+v1kQEZNNDIUQ1c2vThiYPP0zjhx4iff9+EtevJ/G778mKjc3fUGuwWACod8MNGBs0cEK0Qggh6hJn12wEeAewdPhSe9IxfeP0QtOkmnk348X+L3Iu9RwH4g/YRzRykoucRONi+kX7iIYjOPJe1UVGRqqJjIyIuk5bLKTu2k3id9+R+P33WJOSCrUxNmhAwymTaTBhAsZ69ZwQpRBCiNqsJtRspJpTOZt6lvC4cF7e9TIXMy6iUIAtGcg7ZcrRPFw8aOTRiMYejfEyebH3/F6SzckYlREAi7bINC1RNElGxOXCHB3N6SlTMZ85Y5umVcS/MQZPT+rfeisNp03D1LSJE6IUQghRG1V1zYbFaiEuLY6YlBjOppwlJiXG/jibcpazKWe5lHHJ4Z/Lx+RDI09bkpGTbDT2bGx/zpuAKKWkgF2UnyQj4nJQ1KpZ5vPniZxxJ9bk5ELtlclEvRtvoOH06bi1aOGEiIUQQtQ2lanZSMpMsicZ+ZKNZNvx+dTzZarNKKsGbg2KTDIKvvZw8Sj9Zg74/FVFkpFaQJIRUdeVtHyvOTqa05MnY46KthW0F9y7RCl8hg3D78478ehc7n/HhBBCVDNn120U9aXbrM3M2DiDmJQYGro35IbWN5CSmZJvZCPZXPgXYxXholzw8/DjYsZFMi2ZGJRtH3GrttLEswmvDXyNjn4dMRlNpdypfJy1mldpJBmpBSQZEXVZWfYRydvGWL8+mExYCha8A159++J315149u6NUqqaPoEQQoiycmbdRlpWGmeSzhCRFMGBuAN8dOgjUswplb5vQfXd6hPgFYC/lz8BXgG2197++Hvajs1WM3duurPap0nVhJqZokgyUgtIMiLqsgsrVnDuxQWlLt+bNyFp8sQTGOvVI/6DD8g8ebJQW/euXfG7cwY+YWEog6GqP4IQQogyqurfzqeYU4hIjCAyKZKIpOznxAgikiI4n3q+0vG7GlztSYa/lz8B3gH21/5etoTD0+RZ7PudPU3K2aNSRZFkpBaQZETUdRdWrMBn6NBS9xExR0eTtHkzDadMAUBbrSRt2UL8e++Tvn9/ofauLVviN2MG9UZdj3J1rZLYhRBClE9lv5AnZCTkSzLyvr6QfqFSsdV3q0+QdxAB3vlHNgK8Amjq1ZSG7g3t06oc+bnL06aukWSkFpBkRIiSaa1J3bGD+PfeJ2XbtkLXXfz98bt9GvXHjcPg5eWECIUQQuRV2opOgV6BPN7zcVKzUu2JRs4Uq4SMhAr362JwIcg7iCaeTQiPD7cvbatQZOmsKk0Aauo0KWeTZKQWkGREiLJL2/8P8R98QNKmTYWWBjbWq0eDyZNpMHECLnk2UazoyIwQQtRGNWWqzunE09yx8Q7OpZ7DoAxorR2y14ab0Y1gn2BCfEII8Q0h2CfYduwbgr+nP+dTzzttadua8mdfk0gyUgtIMiJE+WWcPEn8kiUkfPU1mM35rilPTxrcfDMNb59G0qZN5a5Zafr0bElIhBC1kjN+O59lzSIiKYJjF49x/NJxjl46yrFLx4hIjMCiLRW6p6eLpz3RyJt0hPiE0NizcbFTqZxdsyEKk2SkFpBkRIiKM587x4Vly7n48cfo1NT8F00mfMLCSPvrL7LOni3Tal6lJS1CCFGTVWXdglVbiU6O5tilY7mPi8c4kXACs9Vc6vsL8nH1oblPc9vIhm/+pMPP3a/cqyZKzUbNJMlILSDJiBCVZ7l0iQsffcTFFSuxXLpU6Lry8ECnpRW9z4kkIkKIOqSyowNaa2LTYjl28Vj+xOPSMdKy0ioUk0KhlLLvtfFO2Du0a9iu0p81L6nZqJkkGakFJBkRwnGsqalc+mwd8R9+SFZMTJFtXBo3pvnaNSilJBERQtRJpRWQ55zzcPGwT6s6fuk4Ry/aXidmJparP4MyEOITQuv6rWndoDV+7n4s2b+Es6lnq32vDanZqFkkGakFJBkRwvG02UzC+vW2vUqOHS/cQCnbju8WiyQiQog6KW9C4mJwQWuNRVvwNnnTpkEbIpMiiUuLK/d9A70Cad2gtS3xqN+aNg3a0KJeC9yMboX6lboNIclILSDJiBBVR1utJP/0E3HvvUf63/sKN1CKkFWr8OpxZfUHJ4QQVSDLmsWxS8fYH7efHdE72Hh6Y4Xu09ijMa3qt7InHK3rt6ZV/VZ4mYpfQl3qNkRBlUlGXKoiICGEqE7KYMAnLAzvIUNI3bWL2DfeJG337twGWhMxfTpNHn6YhpMnoVzknz4hROVU51QhrTXRKdHsj9vP/tj9/BP3D+Hx4aRb0st8D19X33wJR86jvnv9csezJWJLqUlGgHcAS4cvtSckMl1KFEd+Igsh6gylFK5BQWSdO5dzInefkowMzi9cSOL69QS88Dzu7ds7L1AhRK2WU0S9+uDqMhdRA2X+Mp6QkcCBuAPsi9vHP3H/sD9uf7l3Jfdy8WJSh0lc6X8lbeq3oZFHo3KvXFWcnM9RWjKWk5BIIiJKItO0qolM0xKi6hVcNStk2Yck/fgT5xcsAKs1t6GLC3533EGje+/B4ObmvICFELWSI6cpZVoyOXzhsD3x+CfuH04lnipTHEaMGAwGzFYzDdwb8Er/VwjyCWLGphkyPUpUK6kZqQUkGRGiapW0fG/awYNETJyEtcAeJa6hofg/Nx+vXr2cEbIQoharSAG3VVuJSIywTbeKs023OnThUJn37wj2CaZzo850bdQVfy9/Fu1eRFRylBSQC6eTZKQWkGREiKpTln1EzNHRnLx1PJbY2ELvr3/LLTSZ9RhGX9/qClkIUQeUtrRuoFcgd3e7mzPJZ+zTrZIyk8p07/pu9e2JR+dGnenSqIu9vkMKyEVNI8lILSDJiBBV58KKFZx7cUGpy/eao6M5PWUK5jNR+etJsO1L0nTuv/G99trqClsIUQfk/dJvVEY0Gqu2YlRGLNpSpnu4Glzp4NeBLo262B9BPkHF1njIxn+ippHVtIQQl7WGU6YA4DN0aIn7iJgCA2m+YgVJmzfj0b07Mc/8m4zDhwHIio0l6oEHSbz2Wpo+8wympk2qJXYhRO0VmRTJ9ujthNYL5UzymXzJR3GJiELRol6L3MSjcRfaNGiDyWAqc79SQC7qEhkZqSYyMiJEzaPNZuKXLCXunXfQmZn28wYfH5o8Pov6N9/ssNVnhBC1X0JGAjvP7mR79Ha2R2+3r5JVksYeje1JR5dGXejo1xEfV59qiFaI6iMjI0IIUQHKZKLR3TPxGTaMs3Pnkpq9N4k1KYmzc+eR+O16Ap6bj2toqHMDFUI4hdli5u/Yv9kWvY0/Yv7gQPwBrNpa4nsUCo3Gz92PN4a8QZdGXeSXGkKUwODsAIQQwtncWrYgZMVy/OfPx+DtbT+funMnJ8bcQNx776PNZVvtRghRfVaFryImOabUdjHJMawKX1VqO601xy8dZ1X4Ku7bch/91vbj9o238/7+99kft7/IRMTFYPu9biP3Rnx5w5dsvGkjQd5BxKfH8+SvT3I25Wz5P5gQlxGZplVNZJqWELWD+dx5zj7/HMmbt+Q779a+PQHPP49Hl85OikwIkZejirjj0uL4I+YPtkdv54/oPzifdr7Efhu6N6R3QG86NOzAmkNriEmJkaV1xWVPVtOqBSQZEaJ2Sdy0ibPPP48lNi73pMFAw6lTafzA/Rg8PZ0XnBCiwsvb1nevz55ze2x1HzHbOXLxSIn9uBnduLLJlfQJ7EOfwD60bdCWcynnZGldIfKQZKQWkGREiNrHkpjI+VcWcenTT/OdNwUF4T//Wbz79XNSZEIIKPvGg008m3B9y+sJjwtnz/k9pW4y2KFhB64OvJo+AX3o3qQ77i7u+a7L0rpC5CfJSC0gyYgQtVfKjp2cnTuXzNOn852vd+ONNHnyCVwaNODCihWlLi0Mtr1OkjZvti9HLISonKISknOp53jwxwe5mHERAwaslFx03tSzKX0C+9A3sC+9A3rT0L1hqf2uCl9V6tK6OfHJ0rqirpNkpBaQZESI2s2ank7cO+8Sv2QJWHL3DzD6+eHd/xoSvvyqbJsuZu8U3/Tp2ZKQCOEgMckxTP5+MudSz5WpvaeLJ738e9lGPwL70MK3hax4JUQlyNK+QghRxQzu7jR59BF8rxtBzDP/Jj37FwuW+HgSvvwK5e6OOTKS01OnFZmQ5E1ETMHB+Awd6oyPIUSdEpUcxfcnv+e7k9+VmIgYlZHOjTrb6j4C+tClcZdybTIohKg6kowIIUQ5uHfoQOjHa7mwYiWxb7yBTk8HsD0rVWRCUjARKWn0RAhRsgvpF9h0ahPfnfyOvef3lth2VMtRDG0+lF7+vWSjQSFqKElGhBCinJSLC37Tb8fn2qGcnTePlG3bbReyp72aIyM5NWEioR+tBpBERIhKSjWn8mPkj6w/sZ7t0duxaEuxbQ3KgAEDWTqLv87/xYPdH5RERIgaTGpGqonUjAhRN2mtSfjyK8699BLWhIT8Fw0G2yMrSxIRIcrJbDGzNXor3534jp8ifyLdkl5kuyDvIBIzE0nMTLQXsAOyrK4Q1agyNSOyA7sQQlSCUor6/7qRVuu/xXfkyPwXrVbIygIXF0lEhCgDq7ay6+wu5m+fz6BPBvHAjw/w/anvCyUiTTyaMKXjFN4KewsgXyIS4B1AgHcAS4cvJcg7iDPJZ5i+cXqZdmoXQlQ/SUaEEMIBXBo1otlrrxL07jsYGzfOfzEri8zISOcEJkQVWhW+qkxf8mOSY1gVvqrIa1prDl04xGu7X2PYZ8OYvnE6nx35jMTMxHztfFx9uKnNTSwZtoRN4zYxqcMkXtrxUrGjH5KQCFE7SM2IEEI4kHu7dhhcXSk4oz3i9ukEvf02PoMHOSEqIRwvZ+O/1QdXl3njP8C+30ZkUiTfn/ye9SfWcyLhRJHvdTO6MTBoICNbjqR/s/64Gl3t17ZEbCl1GlZOQpLTv+z3IUTNIzUj1URqRoSo+wqumlVvzGji3n4nt4FSBC56hXrXX++8IIVwkJJ2Py+uzaKBi/gr9i++O/kd+2L3FXlfozJydcDVjGw5kiHBQ/B29S42Btl4UIiaQTY9rAUkGRGibitu+d6Eb74h+qnZ+TZKbDpvHg1vG+/EaIVwjJISkrzXGrg3oKVvS/bG7sWqi94NvVvjboxsMZJhocNo5NGoOj+GEKKSJBmpBSQZEaLuKm0fkeRffyXygQchI8N+rvHDD9Po7pnOCFcIhyoqIcm0ZjL1+6nEp8eX+N5W9VpxfcvrGdFiBME+wdUUsRDC0WQHdiGEcKKkzZtL3EfEe8AAmi/7kIi7ZqKTkgCIXbwYtBW/u+9GKeWEqIVwjIJ1GcPXDUdT/C86/b38ua7FdVzf4nraNmgr//8LcZmTZEQIISqp4ZQpAPgMHVrs8r2e3bvT4qPVnL59Opa4OABiX3+DrAsXaTr7KZRBFjcUtdeFjAs0923OmeQzRSYi9dzqMbz5cEa2HEn3Jt0xKPn/XQhhI8mIEEI4QE5CUhK3Nm1o8fFaIu6YQeapUwBcXLkSy4ULBC54EeXqWvINhKhhdp/dzQf7P2Br9NYirw8OHsy4tuPoE9AHk9FUzdEJIWoDSUaEEKIamZo1o/lHq4m8aybp//wDQOL69VgSEgh643UMnp5OjlCIkmmt2Rq9lff3vc+e83uKbGNURizawtGLR2lTv40kIkKIYsk4qRBCVDOXhg0JWbYMzz5X28+l/P47p2+/nayLF50YmRDFs2orP5z+gVu/vZV7Nt9TZCLi5+7HhrEb+H7s97LZoBCiTCQZEUIIJzB6exH8v//hM2KE/Vz63/s4PWky5rNnnRiZEPmZrWa+Pv41N351I4/+/CgHLxzMd93N6AZAkHcQa65fQzOfZrL7uRCizCQZEUIIJzG4utLs1UXUz7PnSObx45y6bQIZJ4rekVqI6pJhyeCTw58w+ovRzPl9DicTTua73rNpT5p4NCHDklHkpoeSkAghyqJWJyNKKXel1Hyl1BGlVLpSKloptVQpFVSBe41QSn2vlIpTSpmVUueVUt8qpcKqInYhhABQRiP+c+fS6L777OeyYmI4PWEiafuK3qFaiKqUak5l+YHlXLfuOp7/43mikqPs1xSKYc2H8enoTxkSMoTzaeeL3X0dCickWyK2VOdHEULUArW2gF0p5Q5sAfoCMcBXQChwOzBKKdVHa328jPd6FHgV0MBWIApoCVwPXK+Uukdr/V+HfwghhACUUjR+4H6MDRpw7j//Aa2xXLrE6Wm3E/TmG3j36+fsEMVlICEjgY8OfcTqg6tJyEjId81FuXB9y+uZ3mU6Leu1BKB9w/YAhIWEFZmI5MhJSLZEbGFSx0lV9wGEELVSrd2BXSn1HPBvYDswTGudnH0+J7H4VWs9sAz3aQycyT4M01r/nufaTcCnQBrQNKePCsYrO7ALIUqVsH490U/NBrPZdsJkotnLC/G97jrnBibqrLi0OFaEr+DjQx+TmpWa75qrwZWxbcZye+fbCfQueg8dIYS47HZgV0qZgAeyD+/LmyRorV9TSk0FBiilemit/yzldr0BV2BD3kQk+17rlFL7gG5AR2Cnwz6EEEIUod7112OsV58zDz6ITk0Fs5moRx8j6+JFGk6Y4OzwRA2zKnxVqSMTADHJMYVGJmKSY/jwwId8fvRzMiwZ+dp7unhya7tbmdJpCo08GlVJ7EIIAbU0GQGuAeoDx7XWe4u4/hnQFRgNlJaMZJRyPceFMkcnhBCV4H1NP5p/uJTImXdjuXQJtObcc89jib9Ao/vvQynl7BBFDbAqfBULdy1k9cHVxdZsgC3pmL5xOmeSbZMArml2DUv+WcK3x78lS2fla+vr6sukDpOY0GEC9dzqVflnEEKI2lrA3i37uejdlnLPdyvmel67gARgiFLqmrwXlFJjsSU127TWxyoSqBBCVIRHt240X70KF39/+7m4t9/m3PPPoy0WJ0YmaoqwkLBSV6rKm4g09WzKHzF/cMNXN/DlsS/zJSKNPBrxWI/H2DRuE/dccY8kIkKIalNbk5GQ7OczxVw/U6BdsbTWl4AZ2Ye/KqV+U0qtVUrtwDbCsgH4V1kDU0odKOoBtCrrPYQQAsCtVStC13yEa8uW9nMXP1pD1KxZWDMznRiZqAlKWzo3byLibnTnXOo5fjnzC1ZttbcJ9ApkTu85bLhpA9M6T8PL5OWMjyKEuIzV1mTEO/s5tZjrKQXalUhr/RlwHRCPbQrYrUAv4DzwY/Z5IYSodqaAAJqvXoV7t672c0nfb+DM3XdjTUkp4Z3iclBcQhKdFM3E7ybap2alW9LzvS/UN5QX+r3At2O/ZXz78faNC4UQorrV1mQkZ8J0cUuBlWtCtVLqMeAH4Fds07K8s5+3A68AH5f1XlrrTkU9gDItMyyEEAW5NGhA86VL8cqzxG/Ktu2cnnY7WRcvOjEyURMUTEiu+/w6Rnw+gti02EJt2zdsz6sDX+XLG77khtY3YDKYnBCxEELkqq3JSFL2c3HjyZ7Zz6UuxauUGggsAv4CbtZa79dap2it9wPjgL3ATUqpYZULWQghKs7g5UXwu+/gO3Kk/Vz6/v2cnjARc3S0EyMTNUGAdwBvDnkThcKiLegCv6u7ovEVvB32Np+M+oRhocMwGoxOilQIIfKrrclIRPZzcTutBxVoV5Ip2c+fa51nIi2gtbYAn2cfDipPgEII4WjK1ZXARa/QYOJE+7nMkyc5ccONpGzbVur7zdHRXFixoipDFE7y57k/uXvz3YWSkCubXMnS4UtZcd0KBgQNkJXYhBA1Tm1NRv7Ofr6ymOs55/eV4V45iUtiMddzzjcsw72EEKJKKYOBps/ModGDD9jPWZOSiJhxJ0mbNxf7PnN0NKenTuPciwskIalDzBYzr+95nWkbpnEu9Zz9vMqerXw+9TxB3kGShAghaqzamoxsxbYcbyulVPciro/Lfv62DPc6m/18VTHXe2Y/nypzdEIIUYWUUjS+9178n52Xe9Jq5cwDD3Lpyy8Ltc9JRMyRkZiCg/EZOrT6ghVV5sSlE0z8biIf7P8g3/nJHSezfuz6Upf9FUKImqBWJiNa60zgrezDt5RS9toRpdSj2IrPf9da78pz/n6l1CGl1IICt/sy+3miUmp03gtKqRuACYAV+MKxn0IIISqnwfjxNFv8f+CSvX+t1sQ8NZv4lavsbQomIs2XL8MUGOikiIUjaK1ZfXA1t3x7CwcvHLSfd1EuvD74dZ7o+QTBPsElLvsrhBA1Ra1MRrK9AOwA+gJHlVIfK6X+AF7FthTv7QXaNwLaAQW3qP0S+BQwAl8rpXYppT5RSu3KvmYA/q21PlxVH0QIISrKd8QIQt77H8rd3X7u/H/+Q+ybb0kiUgedTz3PPZvv4aWdL5FhybCf9zJ5sW7MOoaEDLGfK20fEiGEqAlqbTKitU4HBgPPY9tv5EYgFFgOdC/rjulaa41tX5E7sC3t2xrbJoehwHfAdVrrFx0bvRBCOI5X3740X7kSQ73cXbPj3n6bY0OvlUSkDvnh9A+M/XosW6O35jvfyL0RX4z5gpb1WxZ6T8GEZEvEluoKVwghykTZvouLqqaUOtCxY8eOBw4ccHYoQog6KuPkSU5PnYbl/PnckwYDrTf/IIlILZacmcxLO1/iq+Nf5Ts/KGgQnRt1ZkyrMQR4Fxz0zy8mOYYtEVuY1HFSVYYqhLhMderUifDw8PDsvfXKxaUqAhJCCFH93Fq0IOitNzl9y625J61WUnbtpv4NY5wXmKiwPef28PTvTxOVHGU/5+HiwRM9n+CmNjeVeZWsAO8ASUSEEDWSJCNCCFFHmKOjiX5sVqHzMU8/jXvHDri3aeOEqERFmC1m3vn7HZb+sxRrni2wujbqyov9X6S5b3MnRieEEI5Ta2tGhBBC5CpYrB78/nu5q2xZLJwef5vs1F5L5F2yNycRMSoj93a7l+XXLZdERAhRp8jIiBBC1HLFrZrV+L57iX39DQCsKSmcvOVWWnzycZXWj8QmZfDxrgh2nLxAckYW3m4uXN3Sj1uuCqaxj1uV9VsXaK1Zc2gNr/35Wr6VskJ8QljQfwFdG3d1YnRCCFE1JBkRQoharKTle/3uvJOkzVtIz144wxIXx+nJU2i+coXDE5J0s4X53xzgsz/PYLbkXxjlt6NxLN58hHE9gpk3uiPuJqND+86rtiZDsamx/Hvbv9kalX+lrHFtx/H4VY/jafJ0UmRCCFG1JBkRQohaLGnz5mKX71UuLgS+tICTY29Cm80AmKOiSNq8mYZTpjgshnSzhalLd7Lj5IVi25gtmjU7IzgRm8zy6b0cnpDUlGSoIjaf3sz87fO5lHHJfq6he0Pm953PoOBBTotLCCGqgyQjQghRi+UkFT5DhxY52uHWpg2NHnyA2Fdfs58zNWvm0Bjmf3OgxEQkrx0nLzD/m3AWjO3isP5rQjJUESUt2fts32fx8/BzUmRCCFF9pIBdCCFquYZTppQ47cpv+nQ8unWzH8fMnUfWxYsO6ft8Ujqf/XmmXO/57M9IYpMySm9YRhVJhpxtz7k9jPtmXL5ExMPFg3l95vHGkDckERFCXDYkGRFCiDpOGY0ELFiAcrPVTFji4zn3/AsOufcnuyILTYsqjdmieeKzv1mzM4Kv/opi44Gz/HY0lj9PX+BAdAIn41I4l5hOQpoZs8Va4r2cnQytCl9FTHJMqe1ikmNYFb4Ks8XMG3ve4PaNt+fbO6Rro658OvpTxrUdV+a9Q4QQoi6QaVpCCHEZcGvZgsaPPMz5lxYCkPjdd/gMG4bviOGVum9ZRyQK+ulwLD8dji1TW5NR4W4y4ulqxNPVJc9rI9GX0iqUDH2yO5L7BreuSOh2q8JXsXDXQlYfXM3S4UuL3QU9JjmG6Runcyb5DMsOLONc6jn7NaMyMrPrTO7seicuBvmRLIS4/MjIiBBCXCYaTp6MR48e9uOz8+eTFR9fqXteSMmsbFilMls0SelZnEvM4GRcCgdjEvnz9EV+OxrH8diUCt3zjxOV+9wAYSFhBHkHcSb5DNM3Ti9yhCRvIqJQ+RKREJ8QVly3gnuuuEcSESHEZUuSESGEuEwoo5HAF/+D8vAAwHLxImefnY/W5RtZ0Fqz7VgcM1fu5kB0YoViaehlom8rP7qH1Ke9vw8hDT1p7OOGt5sLhmqYpZSckVXpewR4B7B0+NJiE5KY5BimbpjKmWTbNDJN7p/zuLbj+HT0p7J3iBDisie/ihFCiMuIa/PmNHnsMc69YKsZSfrhBxLXf0e9UdeX+t7kjCy+2HOG5dtPc+x8cqXiuOOalsVOk9Jak2mxkpZpIc1sITXTUuB1lv31B7+d5GRc+UdHvN0c8+MvJyHJGf2YvnE6S4cvBeC29bcRn55/BKahe0Oe7fMsg0MGO6R/IYSo7SQZEUKIy0yDCbeR9MMPpO7YAcDZ55/Hs1dPTE2aFNn+2PlkVv1xms/+POOQEQWTUXHLVcHFXldK4eZixM3FSP1S7nUxJZNFm46UO4arWzputaqCCcnIL0ZisVryjYQADAwayLN9n6WRRyOH9S2EELWdTNMSQojLjDIYCPjPfzB42nb1tiYkcHbes/mma1msmk0HzjLpgx0Mfe0Xlm07VSgRGdSuMUunXcWtPYtPLIoyrofjdkO/pWcwJmP55nWVlgxVRE5C4mJwIcualS8R8XDx4N9X/5s3h7wpiYgQQhQgIyNCCHEZcg1qRpMnnuDss88CkPzTTyR89RXWa0fy8a5IVv1xmqhLaYXe5+Puwi1XBTPp6ua0aOQFQN9WjTgVl1KmlbV6t2jIvNEdHfY5mvi4M65HEGt2Rpb5PY5MhvLSWmOxWvKda9egHYsGLiK0XqjD+xNCiLpAkhEhhKhDYpMy+HhXBDtOXiA5IwtvNxeubunHLVcV/gJe/9ZbSNq0iZRt2wA49ezz3P1bGmdNvoXu297fhyl9QrmxeyCervl/dLibjCyf3ov534Tz2Z9F7ztiMirG9Qhm3uiODt/9fN7oTpyIdU4ylCMmOYbx68cXmpqVYk7Bzej4xEcIIeoKVd5VVCrUiVK+QFvgjNb6bJV3WAMppQ507Nix44EDB5wdihCiDko3W5j/zQE++/NMmZOBjCwLm37+m6BHZ+CeaRsF2dWkHXP7zAClMBoUIzr7M+Xq5vRq0bBMm/HFJmXwye5I/jgRX2oy5Ei2z++cZCgmOYbJ30/Ot2xvx4YdScxM5EzyGYK8g0rch0QIIWq7Tp06ER4eHq617lTe9zosGVFKDQPGA29qrffmOX8P8BrgCmhgsdZ6lkM6rUUkGRFCVJV0s4WpS3eWeWTgpbFdWLcnijU7I4hPyWTY6Z08svcTe5sPrp5As4m3MqFXCP713KsydIfLSYY++/OMfZWt5n6efHZ33ypJhmKSYwrtpt7ArQGfjP4ErbW9qF0SEiFEXVaZZMSRBewzgFuBYzknlFIdgTcBI/AHkAg8opQa7cB+hRDisjb/mwNl3gl9x8kLDHn1F9766Rjx2RsWbgrpyc6m7e1t7tz/FQ908q51iQhAYx837hvcmuduyP15mJKRVWWJyPSN0/MlIgrFS/1fwt/Lv9R9SIQQQjg2GbkS2Ku1Tspz7nZAAdO01v2A7kAGcK8D+xVCiMvW+aR0PvvzTLnek3c83M3FwK09Q+jx+ssYfG21IjolhZhnnin3Zog1SYeA3LqXuORMzielO7yPLRFb7Bsa5rir6130bdbXflwwIdkSscXhcQghRG3myAL2psDuAueGApeAtQBa69NKqV+Bcg/hCCGEKOyTXUXXSJSmnocL9w1uzc09gmng5QpAwpyniX7yKQBStm3n0scf02D8eIfGW10aebvR2MeN2KQMAA7GJNHEx7EjPVc0uQKDMmDVVgB6+/fmnm73FGqXk5BsidjCpI6THBqDEELUdo5MRrKw1YUAoJTyBjoD67XO/pfaJhZo7MB+hRCiRinPilaVkZqZxQ8Hz1fovV2a1eeuAa3ynfMdM4bETT+QvMX22/tzL7+CV79+uAY7dk+O6tIhwJfYpFgADsYkMrCt4370JGQk8NjPj9kTkUYejXhpwEsYDUUXxwd4B0giIoQQRXBkMnIK6JHn+HpstSI/FGjnB8Q7sF8hhKgRSlrR6rejcSzefKRCKzqZLVZOxaVw6GwSR84lcehsEofPJhF5MZWKzqRKySy8k7pSioBn53Fi924sCQno1FRinp5DyPJlKEPt2yO3Q4APvx7JTUYcxaqtzPl9DtEp0QAYlIGXB7wsGxoKIUQFODIZWQssUEqtA34BngAygS9zGijbupA9gMMO7FcIIfKprpGJvMqyopXZolmzM4ITscksn96rUEKitSbqUpo94Thy1vZ8IjaFTIu1mLtWjLdb0f/8uzRuTNO5/yb6Mduih6m7dnFx9Uc0nFz7fqvfMU/diCOTkWUHlvHLmV/sxw90f4Ce/j0ddn8hhLicODIZeRMYA/wr+2EFZmmto/K0GYJtitZ/HdivEEIAVTcyURblXdFqzhf7GdcjODfxOGdLPpIyCo9YlMTFoMiyln945OqWfsVe8x05kqRNP5C0cSMA5199Fe/+1+AaGlrufpwpbxH78dgU0s2WSv93//Pcn7yx5w37cf9m/ZneeXql7imEEJczhyUjWutUpVR/oD/QBPhLa320QDML8AjwjaP6FULUPLV1ZKKiKrKi1bo9UazbE1V6w2xGg6JFIy/aNfWhnb8PbZv60N7fB3eTkf4v/1iuInaTUXHLVcXXgSil8J83l9Rdu7BcuIBOTyf66Tk0X7kCZXRsEleVWjbywtXFQGaWFYtVc+x8Mp2b1avw/eLT4nnilyewaAsAAV4BvHjNixhU7ZvCJoQQNYUjR0bILlT/pYTrPwM/O7JPUTJnfCmsCX1L/85LCGrLyMT8b8JZMLZLhfqyWjWX0szEJ2cQl5zJ6h2nK7SiVXEC67nbEg5/W8LRrqkvLRt7FftnNq5HEGt2Rpb5/uN6lP7/gEvDhvjPm0fUQw8BkLZnDxdWrMTv9mll7sfZXIwG2jb15p8o2xSt8JjECicjFquFp357ivNptgUDXAwuLBq4iPru9R0VrhBCXJYctgN7oRsr1QZoBMRrrY9USSe1SHXvwF7Sl0Kw/Wa0qr4UOrNv6d95/Zd3F3BHj0z0e6n8owPbngqjsY8bWmtSMy3EJ2cSl5JBfHIm8ckZxKdkEpecfZx9Pi45k4upmVgqMDWqIBeDokfzBrTzt412tGtqS0B83U3luk9V/tlHPTaLxPXrAVCurrT48gvcWrYsV3zO9Pinf/Np9qjV7f1CmTe6YivLv/vXu7zz9zv246d6PcXEDhMdEqMQQtR2ldmB3aEjI0opD+BZbLux188+vRyYnn39duABYLrW+i9H9i1yOXO6ijP7lv6d2391jkwUVJG9NswWzeg3f8NoMBCfkkG62bEF4mXRpVk9Pp7Zp9L3cTcZWT69F/O/CeezP4v+s6hoEtr0mTmk7NyBJTYOnZlJ9FOzCf1oNcrFoT8+qkwHBxSxb4vexrt/v2s/vrb5tUxoP6HSsQkhhHBgMqKU8gJ+wrZa1jlgPTCqQLNfgSXArcBfjupb5OfML4XO7Fv6d17/FamZ+OzPSGYObImnq5G0TAspGRbSzFmkZFhIzcx9nZZpO07NzCI100JKZlahcydikysU99nEjAq9L4er0YCftyuJ6WZSMizlfr+3u+O+0LubjCwY24VHr23LJ7sj+eNEvEOm57k0aEDA/Pmcufc+ANL37SP+ww9pdOedDou9KuVPRpLQWmNb2LFszqWcY/Zvs9HZ+9aH+IQwv+/8ct1DCCFE8Rz5q60ngauA94CHtNYZSql8v2rUWh9XSh3CtjP7bAf2LbJV9Evho9e2rXQdgTP7lv4r1v+nuyMZ2cUfbzcXzBaN2WLNfhT9OjPLSpZVY87KPrZosixW/jx9sUIjE4Ne+blc76lqSkEDT1f8vFzx83bFz9uNRl62Zz9vV/y83GjknXvs4+aCUoq3fjzKok3ln41a0opWFdXYx437BrfmvsGtHXZPnyFDqHfDDSR89RUAcW+8iffAgbi3beuwPqpK3uV9E9LMxCSkE1jfo0zvzbJm8cSvT3Ah3ZbguxpceXXQq/i4+lRJrEIIcTlyZDJyK7aND+/TWpf0K8LTwBUO7FfkUdHpKvd/tIcezRtUqu/dFfxCet/qPXRvXj//BV3iIQVrnbSGvZGXKtT/jOW76NSsXvbmcRqtbffTOa/JPcZ+rPOctx0fOZdUof5vfPt3mvi6Y9W2wmiLVWPVtofFaovBkve1VWPRGp19Lud8epal3P1nWTWTl+ws13vqknZNfZg7uqM90WjgacLFWP6VkW7pGczrW446dEWrmqbp07NJ2b6drPPn0WYzMbOfJnTtGpSpfPUt1a2ep4nAeu5EJ6QDtqlaZU1G3tj7BnvO77Efz+49m/YN21dJnEIIcblyZDISAnxbSiICkAhU7luvKFZZp+gU9b6Kvreydp66wM5Tzukb4O8zCfx9JsFp/UddSifqUrrT+q8JlAJPkxFPNxc8XY14mIx4Zb+2PVzwcDXi5WrEw9UFr+zzOa9/OHiOz8uxTG6OMVcE0q915XfNbuLjXiUrWtUkxnr1CHjheSLvmglA+oEDxL3/Po3vvdfJkZWuQ4BvvmQkrEPTUt/zc+TPfPjPh/bj0S1Hc1Obm6oqRCGEuGw5MhlJwbZ6VmlaAPEO7FfkkVzODdOEMCqo7+mKyWjA5KJsz4Y8r40GTMbc165GAy75jhU/HDxHdAUSqt4tGrLs9l64mwyVmoPfI7QB3/wd7dSRiXmjO3EiNqXMK1rNG93RYX1XF+8BA6g37iYSPlsHQNw77+IzeDDuHTo4ObKSdQjwZcsh25K8B88mldo+KjmKOb/PsR+3qteKZ65+RupEhBCiCjgyGfkTuEYpFay1LvLXg0qpTkB34CsH9ivy8Har2H9Sf183BrdvUqm+fzp0vkIFwf6+7gzp0ISCP+YL/txXBVoUvP5D+DliEsr/hTiovgfXdwtAoVAKVPa98x6jVJHnlcL+BeXT3ZGcik8td/9tm3pz3+DWGJTCoBRGA3leKwwGhUGBUeW8trVRSmHMbqMU/PvLf9gTcanc/fdt3YiVd/Qu9/vyauzjVqGaiQFtG+PhWvnVvGrCyERVrmhVkzR96ilStm4jKyYGsrKInv00LT75GOXq6uzQitU+ILfGo7QVtTItmcz6eRaJmbZ2Hi4evDboNTxNnlUaoxBCXK4cmYy8BVwLfK6Uuk1rfSzvRaVUc2AFYMhuK6pA7xYN+e1oXLnfN7lPaKULXitaxDu5T3OHFNs2qeAX4tt6hzikf611hfq/4Ypm3HBFs0r3P6R9kwolI44ooq4JNRM1YWSiqla0qkmM3t4E/ucFIqbfAUDGoUPE/fe/NH7wQSdHVry8K2qdikshLdNSbBL86u5X+Sf+H/vx3D5zaVm/9uyrIoQQtU35KzWLobX+Bvg/bEv7HlZK/YOttneYUmo3cBTbqMjL2TuxiypwS89gTMbyTSVw1JdCZ/Yt/Tu3/5yRifKoqpGJ23qFFPvnYDIqbusV4vD9XQrKWdFq5R29+eLefqy8ozf3DW5d6xORHF59+1J//K3247j/vUfaP9WzoWtFhPp54W6y/bizajh8ruipWhtPbeSjQx/Zj29uezOjWhZcoV4IIYQjOSwZAdBaPwaMB/YDHbHNZAkErgSOA5O11rKkbxVy5pdCZ38hlf6d2/+80Z3o3aJhmdpW9cjEtqfCeHx4O/q3aUT3kPr0b9OIx4e3Y9tTYSwY26XWTpGqSZo+/jimZtkjehYLMbOfwpqZ6dygimE0KNr5l7z54amEU8zbNs9+3L5he57s9WS1xCeEEJczhyYjAFrrT7TWVwBNgV5AHyBYa91Ba73a0f2Jwpz5pdDZX0ilf+f1LyMTlxeDlxcBL75oP844eoy4N2vuDNyOJdSNpGel89gvj5FiTgHA2+TNawNfw80o/78IIURVUwX3axBVQyl1oGPHjh0PHKieqQzpZovTCmmd2bf07/z+AWKTMupszYTI7+wL/+HiqlW2A4OB0DUf4dGtm3ODKsKK7aeY+5Xt39+eoQ349O6+9mvzts3j86Of24//b9D/MbT50GqPUeSns/dTEkJUH6VUhVYO7NSpE+Hh4eFa607l7lP+oleP6k5GcjjzS6Gzv5BK/5IQiKpnTU3laFgY1ouXAHBt2ZIWn6/D4O5eqK05OpqkzZtpOGVKNUcJu05d4Ob/bgfAx82Ffc8OQynFV8e+4pmtz9jbTeowSaZnOZHFYiE+Pp6kpCQya+i0PyHqOldXV3x8fPDz88NoLNsvLWtEMqKU+rEczbXWOswhHdcSzkpGhBCiKl1YsYJzLy7Id67h7bfT9Mkn8p0zR0dzeuo0zJGRNH16drUnJEnpZro8u8l+/NsTg0lXUUxYP4F0i21J8K6Nu7Js+DJMxpq9q3xdZbFYiIiIID398t4EVoiawt3dnZCQkDIlJJVJRhy5tO+gMrTR2IraZThGCCHqAJ+hQ7mwchXmyNw9Xi4sW4bPtUPxvPJKIH8iYgoOxmdo9U+B8nE3EdzQg8gLaQDsPXOO948/ak9E6rnVY9GARZKIOFF8fDzp6ekYjUaaNm2Kl5cXBoPDS1uFECWwWq2kpKRw7tw50tPTiY+Pp0mTyu1DVxpHJiMtijlvAIKB4cBDwNvAOw7sVwghhJOYAgNpvnwZp6dMxXzmjO2k1kTPnk3LL77AculSvkSk+fJlmAIDnRJrB3/f7GREs+TQy5xKO2W/tuCaBQR4BzglLmGTlGRbcrlp06bUq1fPydEIcXkyGAz2v3/R0dEkJSXVnmREa326hMsngV+VUj8B3wN/ACW1F0IIUUuYAgNpvmI5p267jaxz5wEwn47g7Av/IXXXrhqRiIBt88NN4ecw1d/BibTf7efv7HIn/YP6Oy0uYStWz6kR8fLycnI0Qoicv4eZmZlorStU1F5W1Tr+qbXeDPwJPFWd/QohhKhapsBAQtesweCTu4Ruwuef15hEBGzJiME9Crem39jPXdX0Ku694l4nRiWAfKtmydQsIZwv79/Dql7syhl/4yOBche3CCGEqNlMgYGEfvZp/pNK1YhEBCCkkcKj2WqUwQJAQ3c/Xh7wMi4GR85YFkIIUR7VmowopTyAnoAslSGEEHWQwWSCvCuv1JDl47XW/Df8RQyuF7KPFTPaPkNjz8ZOjkwIIS5vDktGlFIhJTw6KqVuBDZhK2bf4Kh+hRBC1Aw5q2ZhseQ7f3rSZMzR0c4JKtuK8BX8FJm7An1m7LWQ3tqJEQkhhADHjoycwlaoXtRjP7AO6AccAR53YL9CCCGcrODyvUY/v8LXqighWRW+ipjkmGKv/3X+Lxb/udh+nJXclsz4QYTHJFVJPEIIIcrOkcnIryU8NgMrgTuA7lrr4n9qCCGEqFUKJiLNly/D88ru9usGX1/MkZFVkpCsCl/Fwl0Lmb5xepEJycX0i8z6ZRZZOst+zpIaChg4dDbRobEIIYQoP4clI1rrQVrrwcU8hmutp2mtP9RaS72IEELUEUUlIqbAQNy7dLW3cWvTGlNwcJUkJGEhYQR5B3Em+UyhhMSqrcz+bTbnUs/Zz9V3bYQ5wbYZ4+GzSVitNaOmRQghLleyfp4QQogKS9q8ucjlez265iYjGYePEPLhUntCkrR5s8P6D/AOYOnwpUUmJB/s/4Ct0Vvtbeu51WP5iBVgqQ9AaqaF0xdSHRaLEELUFDExMcyZM4crr7wSX19f3NzcCAoKom/fvsyZM4cDBw44O0Q7Wc9QCCFEhTWcMgUAn6FD8y3f6965EygFWmNNTkZnZtJ8+TKSNm+2v8dRchKS6Run2xOSB7s/yFt738qNx+jOJ9d/QqBPIKF+xzkZlwLAwZhEWjSSTfYuN7FJGXy8K4IdJy+QnJGFt5sLV7f045argmns4+bs8IQDnDp1ihYtWjBw4EB+/vlnZ4dTrb799lsmTZpEQkICoaGh9O/fH19fX06fPs2ff/7J9u3bCQwMpFOnmrHTRoWTEaXU3Er0q7XWz1fi/UIIIWqIopILo7c3bq1bkXH0GABpf++j/r9udHgikqNgQvLEb0/kxqKMrL5+NYE+tmSpQ4BPvmRkZJeAKolJ1DzpZgvzvznAZ3+ewWzJP0Xvt6NxLN58hHE9gpk3uiPuJmMxdxGi5vrrr7+46aab8PDw4LPPPmPs2LH5dk+/dOkSy5Yto3///k6MMr/KjIw8C2igIvvDa0CSESGEqMPcu3S1JyPp+/fBv26s0v4CvANYMmwJIz4fgSb3i+YbQ96gbYO29uMO/r58t/8sYEtGxOUh3Wxh6tKd7Dh5odg2Zotmzc4ITsQms3x6L0lIRK0zf/58MjMzWbx4MTfddFOh6/Xr1+fhhx+u/sBKUJmakduB6dnP5X1Mr0S/QgghaoG8dSNpf++rlj6/P/V9vkTEoAy0qd8mX5sOAb721wdled/LxvxvDpSYiOS14+QF5n8TXsURlc2pU6dQSjFo0CBSUlJ49NFHCQ4OxsPDgyuvvJJvvvnG3vbTTz+lV69eeHl50bRpUx588EHS0tIK3XP9+vVMnz6dDh064Ovri5eXF926dePFF18kIyOjUPtbb70VpRRPPvlkoWuHDh3C09MTX19fTpw4Ua7PNXPmTEJDQ3Fzc6Nx48aMGzeOffv2FWqX8/nT0tJ46qmnaN68OW5ubrRu3ZqFCxeiC2yu+uyzz9KiRQsAfvnlF5RS9se0adMK3TcxMZHHHnuMFi1aYDKZ8n1Z3759OzfccAONGzfGzc2N0NBQ7r33XqKLWIij4D0feughgoODcXd3p0OHDvzf//0fVqvV3n7Xrl0opejXr1+xf07z589HKcULL7xQpj/Xn376CYC+ffuWqX1NUOGREa31ckcGIoQQom7x6NrF/jr98GGsGRkY3KpuPv5f5//i9T2v248VCqu2Mn3jdJYOX0qAt206VofA3GQk6lIaCalm6nmaqiwu4Xznk9L57M8z5XrPZ39G8ui1bWtMDUlmZiZhYWEcP36cq6++muTkZH799Vf+9a9/sWHDBvbv388TTzxBz549GTZsGL/99htvvvkm8fHxrF69Ot+97rjjDlJSUujUqRNdunQhMTGRnTt3MmfOHLZs2cKmTZswGnNHhf773/+ybds2Fi1axHXXXcegQYMAMJvNTJgwgbS0ND788ENatmxZps/y+++/c/3115OYmEinTp0YM2YMUVFRfP7553z33XesX7+ewYMHF/r8w4YN48CBA/Tq1YsOHTrwyy+/8NRTT5GUlJTvy/oVV1zBTTfdxLp162jatCkjRoywX7vmmmvy3TctLY2BAwdy+vRpBg4cyJVXXkmDBg0AWLVqFdOmTcNqtdK3b1+Cg4PZs2cP7777Lp9//jk///wz7du3L/T5MjIyGDJkCMePH2fIkCFkZmayZcsWHn30Ufbt28eHH34IQM+ePenRowfbtm3jwIEDhWo4rFYrH374IUajkdtvv71Mf7Y5yeRvv/1Gt27dyvQep9Nay6MaHsCBjh07aiGEuFxYMzP1wW5X6PB27XV4u/Y6Zc+eKusrKjFK91jZQ3de1ll3XtZZX736ar0/dr8e8dkI3XlZZz3isxE6OinaFpfVqrvM26CbP/mtbv7kt3r78bgqi0uUjcVi0eHh4To8PFxbLJZC161Wq76UmlnhxysbDtn/e5fnsWjjoUr1a7VaK/1nc/LkSY1terseNGiQvnDhgv3ahx9+qAHdunVr3bBhQ/3rr7/ar0VFRekmTZpoQB8/fjzfPb/44gudnJyc71xiYqIeNWqUBvTy5csLxfHjjz9qg8Ggg4OD9cWLF7XWWj/xxBMa0OPGjSvz50lISND+/v7aZDLpTz/9NN+1H374Qbu6uupmzZrpjIyMQp+/f//+OjY21t5+165d2sXFRXt6euqkpKQi/9wGDhxYZBx579unTx/7Z8oRERGhPTw8tIuLi/7mm2/s5y0Wi3744Yc1oHv27FnsPbt27Zov1mPHjunAwEAN6K+++sp+/r333tOAfvjhhwvF+P3332tAjx49usjPUJQhQ4bYYwgLC9Nvv/22Pnz4cJnfn/dzlvR3sqCOHTtq4ICuwHdkpbWssV4dlFIHOnbs2LEmLaUmhBBV7dSEiaTt2QNA06dnV0kBe0xyDLd+eysXMy7azy3ov4BRLUcRkxxjL2oP8g6yj5Dc+r/t9ik780Z35PZ+LRwelyg7q9XK4cOHAWjXrh0GQ/5Z5AlpZrrN3+SM0Crl73nDqOdRuVG3nFWhjEYjhw4donXr1vZrVqsVf39/YmNjmTt3LvPnz8/33kcffZT/+7//48MPP7RPTyrJsWPHaNOmDWPHjmXdunWFrj/xxBO88sor3Hbbbdx1112EhYXh7+/P/v37adiwYZk+z+LFi3nkkUeYPXs2L774YqHrjzzyCIsXL2bdunWMHTvW/vkNBgMHDx6kbdu2+dqPGTOGb775hp9++sk+YgOlr6aVcx1s06WuuuqqfNfnzZvHc889x+TJk1mxYkW+axkZGbRs2ZLo6Gi2b9/O1VdfXeiemzZt4tprr833vv/+97/cc889DBs2jI0bNwKQkpJCYGAgJpOJqKgo3PKMHo8bN45169bx9ddfM3r06JL+WO2OHTvGuHHj+Pvvv/Odb926NTNmzODBBx/Ew8Oj1PuU9neyoE6dOhEeHh6utS73El0O32dEKXWNUuoVpdSXSqktSqkfi3hscXS/Qgghap6qrhuJSY5h6oap+RKRwcGDub7F9UDx+5DkrxuRInZR84WGhuZLRAAMBgPNmzcHKPTFF6BVq1aAbc+Jgo4ePcrrr7/OAw88wPTp05k2bRrPP/+8/VpRXnjhBbp3786aNWu44YYb0FqzfPnyMiciAD/88AMAN954Y5HXc6ZR7dq1K9/50NDQQokIYD9X1Gcsi4CAgEKJCNimOQFMnDix0DU3NzduvvnmfO3yatiwYZH/PSZMmADAtm3b7HUuXl5eTJw4kfj4eL744gt72/Pnz/P1118TGBjIyJEjy/x5Wrduzd69e/nxxx954IEH7NPIjh07xlNPPcXAgQPJysoq8/2qg8P2GVG2dcOWAFPJXWGr4GpbOccyHCOEEJeBvHUjafv3O/z+myM2E5OS+yWknls95vaZm28py4LL/m6J2ELHgIH261LELmqDZs2aFXney8ur2Os51/IWpWutmTVrFv/3f/9XqPA7R1JS0X8nXF1dWbZsGd26dSMxMZEHHniAoUOHlutznDp1CoDevXuX2C4uLi7fcVBQUJHtvL29AYosvC+LkJCQIs/nFKiHhoYWeT3nfFGF7DkJYkG+vr7Ur1+fS5cukZiYSL169QC4++67effdd3n//fcZP348AMuWLcNsNjN9+vR89TtloZRi8ODB9rqbU6dO8corr/DOO++wa9cuPv30U2677bZy3bMqOXLTw7uBacBuYDZwD/AvoB3QErgVmAy8BrzjwH6FEELUUO5dcwsozRERZF28iEt2cagjmAz5p8A83etpGnk0KtQuJyHZErGFSR0nsf9Mgv3a4XNJZFmsuBgdPllAOIivuwt/zxtW4fe//+sJ3vrpWLnf98CQ1szoX7ai7KL4ujvua1beBLsi13N8/PHHvPbaawQFBbF48WL69OlD48aNMZlMZGZm4ubmVmySkvP+HHv37sVqtZY6hScvi8UCwM0334ynp2ex7QomK2X9fOXl7u5e4nVH/bnnKOrPtmvXrlx99dX89NNPHD9+nFatWrFkyRKUUtxxxx3lun9RQkNDefvtt/nll184cOBAjdp9HRybjEwDUoDrtNbxSqlJAFrro8BRYKNS6jvgY2AbcNqBfQshhKiBTM0CMTZsiOWCrT4jff9+vAcMcMi9I5Miee3P1+zHYSFhXNfiumLbB3gHMKnjJADaNPXGaFBYrJrMLCsn41Jo09THIXEJx1NKVar2Ykrf5vzv1+OFNjosicmomNIntNI1HzVNzlSgd999l1GjRuW7VtrSvL/99hsvvfQSgYGBtG/fnh9//JGXXnqJp59+usz9BwUFcfjwYZ555hm65pnGWdMEBgZy+PBhTp48WeT0sNOnbV9jAwIKb5oaERFR5D0TExNJSEjAy8sLX1/ffNfuvvtu/vjjD5YsWcLw4cM5cuQIw4YNK3ZkpiKaNWvGgQMHSkwCncGRvwbqAGzXWsdnH2sApZR9bElr/RnwJzDLER0qpdyVUvOVUkeUUulKqWil1FKlVNFjeaXfr7VS6n2l1Kns+8UqpbYppR53RLxCCHG5UUpVSd2IVVuZu3UuaVm2PRTqu9XnmaufKfNvKd1NRlo28rIfh0vdSJ3WxMedcT3K99VgXI/gGrOsryNdvGirrwoODi507ZNPPin2fYmJiUyePBmtNR9++CGrVq3Cz8+PZ599lt27d5e5/5xpXV9++WX5Ai8nV1dXgArXR+TsUF5wWWSwLTP86aef5muXV3x8PJs3by50fs2aNYBtD5CC/1bdcsstNGjQgGXLlvHuu+8CcOedd5Y53tIWpIqPj2f79u0AhIWFlfm+1cGRyYgByDvBLzX7ueB4/FGgC5WklHIHtgBzAW/gKyAS26aKe5RSrcp5v38B+4E7gHjgC2Av0AKYWdl4hRDicuWer27EMcnI2kNr2X0u9wvQnN5zipyeVRLZ/PDyMm90J3q3KFuhde8WDZk3umMVR+QcOb/lf++99/J9gf3tt9945ZVXin3fvffey+nTp3nggQcYNmwYAQEBvP/++5jNZiZOnEhqamqx781r5syZNG7cmBdffJEPP/yw0JfolJQUVqxYwZkz5dsXpqBGjRphMpk4fvy4fWpYedxxxx14eHiwZs0a1q9fbz9vtVp5+umniYqKomfPnvaVtAp6/PHHiY+Ptx+fPHnSvkDAvffeW6i9h4cHU6ZMISYmho8//pjGjRtzww03lDneJ598klmzZhEZGVno2l9//cXIkSNJSkpi7NixpdbrVDdHJiNRQN5fO+RMw+peoF1bwBFl/E8DfYHtQFut9a1a697AY0BjYGlZb6SU6gasxTbNbIDWuofW+jat9TCgGTDeAfEKIcRlySNP3Uj6vv2l/gavNJGJkSzes9h+fG3zaxkeOrzc95EVtS4v7iYjy6f34rZeIZiMRY+gmYyK23qFsHx6L9xN5Ssari0efPBBvLy8eOedd+jcuTO33XYbAwYMYODAgdx9991Fvmft2rWsXr2aTp06sXDhQvv5f/3rX9x+++0cOXKERx99tEz9N2jQgC+++AIvLy+mT59OixYtGDVqFDfddBM9e/akadOmTJ06tVABe3m5uroyYsQIzp49S7du3ZgyZQozZsywbzhYmpCQEHvCNnr0aPr378+ECRPo2LEjr776Kk2bNi205G+Oq6++GoPBQJs2bRg3bhxjxoyhc+fOREVFMWnSpGJXEps5M/d339OmTcNkKvsUwW3btvHqq68SGhpKt27duPnmmxk7dixXXHEF3bt3Z+fOndx4442sXLmyzPesLo5MRvYAHZVSOXUom7CtnPWKUqqDUsone7pTD2wjDhWmlDIBD2Qf3qe1Ts65prV+DdgHDFBK9SjjLd8EXIFpWuvf817QWlu11mUffxRCCJGPR5fO9teWS5cwF/Gbu7Kyaiv/3vZv+/SsBm4NmNN7ToWKW9sH5NaIHDorycjlwN1kZMHYLmx7KozHh7ejf5tGdA+pT/82jXh8eDu2PRXGgrFd6mwiAraRkV27djF69Gji4uL4+uuvSU5O5n//+1+RIyORkZHcc889uLq6smrVqkIF32+88QYtW7bkf//7H19//XWZYujXrx/79+/nsccew8PDgx9//JFNmzaRmJjIqFGj+Pjjj+nYsfIjUx988AGTJ08mPj6ejz76iCVLlvDLL7+U+f2TJk3i119/ZdSoURw8eJDPPvuMtLQ07rnnHv78888id18H29K/P/74I7fddhvbt29n48aNBAcHs2jRIpYtW1Zsfx06dCAwMBCAGTNmlOuzvvfee7z88suEhYWRmprK999/z3fffcelS5e49dZb+f777/niiy9qXL0I4LhND5VStwGrgbFa6y+zz32EbVQhbycW4Bqt9c5K9DUY+BE4rrVuXcT1fwPPAfO11s+Wcq8OQDhwRGvdrqIxlUY2PRRCXM6Oj7iOzOwlPQMXLaLeqOsrdJ/VB1fz0s6X7MeLBi6q0KgIwLnEdHq/mLvt1Z5/X0tDL9cK3UtUTnk3WBOiJipto8XSbNu2jX79+lX4/Y5UKzc91FqvATyAb/KcnoptOtUu4BjwHRBWmUQkW86Y/55iru8p0K4kOVU8P2QXxE9VSr2plHpDKTVDKeVb4ruFEEKUKl/dyL6/S2hZvNOJp1n852L78fDQ4RVORACa+LjlSz5kqpYQwplydqS///77nRxJ9XLk0r5orTMKHJuBl7IfjpSzQ01x1U1nCrQrSU4Glwb8hW1flLwWKKVu0lr/WpbAlFLFDX2Uq6BeCCHqEo+u3Uj82va7qvR95d/8MGf1rHRLOgAN3RvydO+yLydaFKUUHQJ82HrMVmR6MCaRfq3LVwQvhBCVsW3bNpYsWcI///zDzp076dGjB2PHjnV2WNXKYSMjSilvR92rDHL6Km7phpQC7UqSs9rXw0BDYCxQH1tS8hHQCPhSKVV4IWkhhBBlkncn9vTwcHRmZrnev/rgavaczx0Mf+bqZ2joXraVkUrSwT938FuW9xVCVLcjR46wdOlSDh48yOjRo/n8888vu2mKjhwZOaeU+hpb3cgGrbUjVswqTk6lYnEFL+WpZMypUnMBJmmtN2UfJwATlVJtgJ7AfcAzpd2suLly2SMmdXOdQCGEKIVb+/YokwltNqMzM0k/chSPzmWbWnwq4RSv73ndfnxd6HVc2/xah8Qly/sKIRwlNDS03KsFTps2jWnTplVNQLWEI1OvDOBWbPt9nFVKvaOUusaB988r5yeGVzHXc5YKSC7melH3isqTiOSVswbcoLKFJoQQoiCDqytuHTrYj8taN2KxWvj31n+TYbHNAm7o3pDZvWc7LK68ycix80lkZlkddm8hhBClc2Qy0hTbFKd12ArZ7wZ+yd7N/D9KqXJX15cgIvu5uO1Ugwq0K8mp7OfTpVxvUoZ7CSGEKEbendjLWjey6uAq/or9y3489+q5NHAvuJduxbVu4m3fc8Js0RyPLcvvsIQQQjiKI1fTMmutv9Ra34ItMZkG/AAEArOBfUqpv5VSs5RSxSURZZXzK7Uri7mec74sW/3m7HlS3ORjv+xn+QklhBCV4JFvRa3S/3k+mXCSN/e+aT8e2WIkYc3DSnhH+bm6GGjVOLe8UFbUEkKI6lUlFTJa62St9Qqt9QhsyciDwA6gC7AQOFnJLrZiq+lopZQquMM7wLjs52/LcK8t2AreWymlgou4Pij7ubhlhIUQQpSBe5fcZCTzxAksScXXaFisFp7Z+ox9epafux+zezluelZeHWUndiGEcJoqL9fXWsdprd/CVk/yX2zF5ZXqV2udCbyVffiWUspeO6KUehToCvyutd6V5/z9SqlDSqkFBe6Vim0HdhPwboF7jcC2V4oG3qtMzEIIcblzDQ3F4Jv7xT99f/FTtVaEr2BfbO7oydw+c6nvXr9K4pIidiGEcJ4qTUaUUvWVUncqpX7GNhoyM/vSruLfVWYvYBtt6QscVUp9rJT6A3gViAduL9C+Ebbleotaonc+ttGW67Pv9YVSaiuwHtsqW884YKNGIYS4rCml8OiSd6pW0cnIiUsneGvvW/bjUS1HMSRkSJXF1aHAyEh5V8MRQghRcQ5PRrJ3Mb9FKfUlcBbbaMgA4DjwHNBWa311ZfvRWqcDg4Hnse03ciMQCiwHumutj5XzXkOAOcAl4DpsmyH+BIzSWr9Y2XiFEEKAR7fcIvai6kayrFk8s/UZMq22fUgaezTmqV5PVWlMHQJ87K/jUzKJTcooobUQQghHctg+I0qpYcBEbEmBN7bpWOeAj4HVeadMOYrWOg2Ym/0ore2zwLMlXM8EXsx+CCGEqALuXfIXsWutUSp3a6jlB5azPy53xGRen3nUc6tXpTH5ebvRxMeN89lJSHhMIk183au0TyGEEDaOHBnZAEzGloSsAkYAzbTWD1dFIiKEEKL2ybu8ryUujqyYGPvx8UvHefuvt+3HY1qNYWDwwGqJS+pGhBDCORyZjHwL3AY00VpP1Vpv0lrL7lFCCCHsXPz8MDVrZj/OqRvJsmYx5/c5mK1mAJp4NOGJnk9UW1wF60aEEEJUD0fuMzJGa/1xdv2FEEIIUaSi6kaWHVjGgfgD9vPz+lb99Ky88taNSDIihBDVp8qX9hVCCCHycu+Sdyf2fRy9eDTf9KwbW9/IgKAB1RpT3pGRE3EppJst1dq/ELXRDz/8gFKKrVu32s99/fXXeHh48O23ZdnqTQhJRoQQQlSzfCMjB/5h7q9zyLJmAdDEswmP93y82mNq2cgLVxfbj0SLVXPsfHK1xyBEbXPttdcyfPhwZs+2bUj6yy+/cNttt/HBBx8watQoJ0cnagtJRoQQQlQr9w4dwGgEQKelk3gk3H7t2T7P4uvqW9xbq4yL0UDbpt7243CZqiVEmbzyyits3bqVF198kTFjxvDyyy8zceJEZ4clahFJRoQQQlQrg4cHbu3a2o/bRNs2GRzbZiz9g/o7Kyw6+EsR+2UlMwWerWd7ZKY4O5oqdfz4cSIiIqrk3l26dGH8+PHMmTOHWbNmcd9995Xpfdu3byc9XcqMhSQjQgghnMCtc2f769bRmqaeTZl11SwnRiQraom6JTExkSVLltC/f39at27Nnj17Sn2P1ppVq1YxdOhQ/Pz8MBqNKKUKPSZNmmR/z99//813332HyWSicePGZY5v9uzZ+Pv7M3PmTLZt21ahzyjqBodteiiEEEKU1c6GF2mf/bp1tGZ+3/n4uPqU+J6qVnCvkYIbMgpR01mtVjZv3szy5cv54osvSEtLQynFwIEDad++fanvve222/jkk09wd3dn4MCB1K9fn99//52oqKh8bbtkb1567Ngxhg8fzr333ouLiwvPPvsskyZNwtvbu6gu8rnhhhs4cuQI7733Hu+99x5t2rRh6tSpTJ48mZCQkIr/IYhaR0ZGhBBCVKvDFw6zxPKr/Tg4TtGn/hXOCyhbxzzJSEKamZgEmUIiaodDhw4xe/ZsQkJCGD58OB999BEhISG88MILnDx5kp9//rnUZOT555/nk08+oVOnThw8eJANGzawdu1ajhw5Qvfu3QFYuHAhaWlpzJo1i6ioKK699lpuvPFG/vOf/zBr1iwsFguvvPJKmWJ+5JFHiIyMZMOGDUyYMIGoqCieeeYZQkNDCQsLY+XKlaSk1O3pc8KmSkZGlFIugB/gVlwbrXXVTF4UQghRY5ktZp7Z+gwRDS2kuoJnJhi0Ju3AAbx69XJqbPU8TQTWcyc6Owk5GJNIYH0Pp8YksmkN6QmVu0dmav5jc57jxBgweea/7lrguCLc60EVja5dvHiRtWvXsmzZMnbu3AlAo0aNuO+++5g8eTK9e/cu873OnTvHggULcHNz45NPPiE0NNR+zdPTk+nTp/PAAw/w+++/88QTTxAfH8+wYcO46qqreOeddwDw8fFhzpw5PPPMM9x9990EBASU2q/RaGT48OEMHz6cpKQk1q1bx4oVK/jpp5/48ccfuffeexk3bhxTp05l4MCBMlJZRzk0GVFKDQWeAa4GTCU01Y7uWwghRM33/v73OXThEBgUJwIUnU/bitfT9+93ejICtqlaeZORsA5NnRyRAGyJyMLmVXf/t3pUzX2fPA0e9R1/2yef5PXXXycjIwM3NzfGjh3LlClTGDlyJCZTSV+/irZy5UoyMjKYOXMmHTt2LHQ959zFixcB8PPz48CBA4XaPfzwwzz88MPl7h9sycy0adOYNm0akZGRrFy5kpUrV7Js2TKWLVtGaGgo69at48orr6zQ/UXN5bCEQCk1CvgCMAIXgROALNQuhBACgIPxB3l/3/v2Y2On9nD6IABp+/Y7K6x8OgT4suXQecBWNyJETbRjxw4yMjIwGo08/vjj3Hffffj7+1f4fps3bwZsdRxFSUiwjUr5+flVuI/yCA4OZvbs2fTr149Zs2axe/duTp06RUREhCQjdZAja0bmZd/vYaCJ1rqn1npwcQ8H9iuEEKKGy5melaVtmxsGegXSf9h0+/W0ffucFVo+sqKWqA0WLFjAtGnT8PDw4IUXXiAoKIgRI0awevXqCtVZ7Mv++9ermNHJXbt2AdCjRxWNIOVx6NAh5syZQ2hoKIMGDWL37t106tSJl19+mYEDB1Z5/6L6OXKqVCdgu9b6DQfeUwghRB3wv33/48jFI/bj+f3mU9/Qkrjs46yYGMznz2Nq0sQ5AWbrEJC7otfJ+BRSM7PwdJVZxU7nXs825cmRMlPh/zrYXj9y0DE1IgW513P8PYE+ffrQp08f3nrrLdatW8fy5cvZtGkTGzduxMvLi7FjxzJ58mTCwsIwGEr/vXNsbCwAvr5Fbzj65ZdfAlTZrurnz59n7dq1rFy5kt27dwO2UZj777+fqVOnctVVV1VJv6JmcOS/sMnAOQfeTwghRC2wKnwVYSFhBHgXXbAaHh/OB/s/sB/f2u5Wrg64GgCXpk3JOmf70ZG+fz+msLCqD7gEzf288DAZSTNb0BoOn02ie0gDp8YksBWBO7r2wpintsKjHrh6Ofb+1cDLy4spU6YwZcoUIiIiWLlyJcuXL7fXWwQGBjJhwgQmTZpEt27dir2Pj48PFy9e5NixY3To0CHftffee4+DBw8yYMAA+6pajpCWlsbXX3/NypUr2bhxI1lZWZhMJsaMGcPUqVMZNWoUrq6uDutP1FyOnKa1GbhSKSXLBQshxGViVfgqFu5ayPSN04lJjil0PdOSyZzf52DRFvu5QK9A+2uPrl3tr2tC3YjRoGjnnzs6InUjorYICQlhzpw5HDlyhK1bt3LXXXeRmprKokWLuOKKK1i/fn2x7x0wYABgW97XarXaz2/YsIHHHnsMLy8v3n77bYfGO2bMGMaPH8/69evp3LkzixcvJioqiq+++oqxY8dKInIZcWTi8CTgAbyqlDI68L5CCCFqqLCQMIK8gziTfKbIhOS/f/+XY5eO2Y8bezTmuhbX2Y/du3axv07b93fVB1wGUjciaru+ffvyv//9j5iYGNasWcOIESPQWhfb/tlnn8XDw4M1a9bQqVMnJkyYwNVXX811112Hi4sLX375JZ07d3ZojH5+fjz66KPs27ePvXv38tBDD5VrB3dRdzhymtbtwPfAg8AopdTPwBlsy/gWpLXWzzuwbyGEEE4Q4B3A0uFLmb5xuj0hWTp8KQHeARyIO8CS/Uvsbb1MXqweuTrfdC6PrrlTR9L3/4O2WlFlmONelfLWjUgyUoe5esGzldy7pIZzd3dn/PjxjB8/vsR2V1xxBb/99htz585l27ZtrFu3jpCQEB566CEef/xxmjVr5vDY1q5d6/B7itrJkcnIs9gSDwW0yn4URwOSjAghRB1QVELy32v/y5O/PYkV25QPozLy0ciPCtWVuHfqZKsH0BprcjKZJ0/i1qqkHx9VL+/IyKGzSVitGoNBNlsTdVuPHj1KnMolRFVx9MiIEEKIy1DBhGTMl2Ow6ty55y/1f4mW9VsWep/R2wu31q3JOHoUsNWNODsZaZ+nZiQ5I4uoS2kEN6yClZaEEEI4LhnRWi931L2EEELUPjkJyXWfX5evYP3G1jcyosWIYt/n3rVLnmTkb+r/68aqDrVEPu4mght6EHkhDYDwmERJRoQQoorIyldCCCEc5u/Yv/MlIgC3dyp54Dxf3UgNWFELoIO/FLELIUR1qJJkRCnVSyn1mFLqVaXUouzXRW/rKYQQok7YdGoTT/z6RKHz9225r8hlf3N45FlRK/3wYazp6VUSX3nIilpCCFE9HJqMKKXaKqX+ALYDLwOPAI9mv96ulNqulGrjyD6FEEI43w+nf2DWL7PQ2QsouhnceGXAKyUu+5vDrU0blLu77SAri/SDB6sr7GLlT0ZkrxEhhKgqDktGlFIBwC9ALyAGeANbMvII8DoQDfQGfs5uK4QQog7YErGFx35+zJ6IuBpceWvoW4xoMYKlw5eWmpAoFxfbqlrZ0vftq7bYi9MxTzIScSGVpHSzE6MRQoi6y5EjI88ATYH/A1pqrR/RWr+e/XgUaAm8BgQATzuwXyGEEE7yc+TPPPLTI/ZExKRMvDnkTa4OuBrILWovLSHx6JJ380Pn140ENfDA2y13jZfDZ2V0RAghqoIjk5GRwGGt9WNa68yCF7XWZuBx4DAwyoH9CiGEcIID8QfyTc1yUS68PuR1+jbrm69dwYRkS8SWQvfy6NbV/jqtBoyMGAwq3xK/UjcihBBVw5HJSACwp6QGWmud3UamaQkhRC126MIh7tp0FxmWDMC2qeHiwYvpH9S/yPY5CcmTPZ9kUsdJha67d8lNRsyRkWRdvFg1gZdD3rqRcKkbEUKIKuHIZCQRCC5Du+DstkIIIWqhwxcOc+emO0nMtP1T7qJceG3QawwMHlji+wK8A4pMRABMzQIx+vnZj2tC3YisqCWEEFXPkcnIdqCvUuq64hoopUYC/YBtDuxXCCFENTl68Sh3brqTSxmXANuIyCsDX2FIyJBK3VcpVePqRjoE5E7TOnw2CYtVOzEaIYSomxyZjLwEaOBLpdSHSqlrlVJtlFKts18vA74ALNlthRBC1CInLp1gxqYZXMywTaEyKiMLByxkaPOhDrl/Tasbaefvg1K212lmC6fjU5wbkBBC1EEOS0a01tuB24EsYCqwATiErWB9AzAl+9rtWus/HNWvEEKIqncy4SR3bLqDC+kXADAoAy9e8yLDQ4c7rI+8dSPp+/ZhKzN0Hk9XF1r4edmPZb8RIYRwPIdueqi1XgW0A54HfgaOZD9+Bp4D2mutVzuyTyGEEFUrIjGCGRtnEJcWB4BC8UK/FxjZcqRD+/Ho0tn+2pKQgDkiwqH3rwipGxFCiKrlUnqT8tFanwHmOfq+Qgghql9kUiTTN07nfNp5wJaIPNfvOUa3Gu3wvoz16uEaGkrmqVOArW7EtXlzh/dTHh0CfFi/37YviiQjQgjheA4dGRFCCFF3RCVHccfGOziXes5+bl6fedzY+sYq67Om1Y2095eRESGK895776GUKvTo2rVr6W8WIpskI0IIIQqJSY7hjo13EJOSu1v6v6/+Nze1valK+y1YN+JsHQJzk5HohHQupRba01eIy9aECROIiYmxP1auXIm7uztPPvmks0MTtUiFkxGllFUplaWUapt9bCnHI8txH0EIIYQjnU05y/SN04lKjrKfe7r309zS7pYq7zvvyEj6wYPoTOd++Q+s546ve+6M5kNnpYi9tloVvoqY5JhS28Ukx7AqfFU1RFT7eXt74+/vj7+/Pzt37uTee+9l5cqVTJw40dmhiVqkMiMjEUAkYM4+jsw+V5ZHZCX6FUIIUUXOp57njo13cCb5jP3ckz2f5Lb2t1VL/27t2qFMJgB0Zibph49US7/FUUpJEXsdsCp8FQt3LWT6xuklJiQxyTFM3zidhbsW1qmE5Pjx40RU4YIQa9euZdKkSaxZs4Zx48aV2n779u2kp6dXWTyidqlwMqK1DtVat9BanyxwXKaH4z6CEEIIR4hNjeWOjXcQkZT7pWXWVbOK3TW9KhhcXXHr2MF+nLa/BkzVkmSk1gsLCSPIO4gzyWeKTUhyEpEzyWcI8g4iLCTMCZE6TmJiIkuWLKF///60bt2aPXv2lPoerTWrVq1i6NCh+Pn5YTQai6wJmTQp99+EJUuWMHPmTL788kuuv/76MsU2e/Zs/P39mTlzJtu2yT7YlzupGRFCCEFcWhwzNs3gVOIp+7lHejzC1E5Tqz0Wj7x1I387PxnpmC8ZkWlatVGAdwBLhy8tNiEpmIgsHb6UAO8AJ0ZcMVarlU2bNjFx4kT8/f2ZMWMGW7duZeDAgbRv377U944fP57JkyezdetWevbsyc0330yzZs0Kte3SpQsAr7/+Oo8//jgbNmxgyJAhZY7zhhtuwNPTk/fee49+/frRtm1b/vOf/1Tp6I2ouRyWjCilliqlppeh3TSl1FJH9SuEEKJyLqRf4M5Nd3Ii4YT93IPdH2R651L/Sa8S+VbU2r/fKTHklXdk5PC5JLIsVidGIyqquISkLiQihw4dYvbs2YSEhDB8+HA++ugjQkJCeOGFFzh58iQ///xzqcnI888/zyeffEKnTp04ePAgGzZsYO3atRw5coTu3bsDsHDhQtLS0pg1axYLFy7kmWeeYfXq1bRo0YKzZ89y9uxZEhISSo33kUceITIykg0bNjBhwgSioqJ45plnCA0NJSwsjJUrV5KSkuKQPxtR8ylH7XCrlLICy7TWJf70Ukq9D0zXWhsd0nEtoZQ60LFjx44HDhxwdihCCGF3Mf0id2y6g6MXj9rP3XvFvdzT7R6nxZR56hTHR1xnP267cwdGX98S3lG10s0WOs3biMVq+3n5wyMDaNPUx2nx1EVWq5XDhw8D0K5dOwyG/L8r1VqTZHbMqNTZlLM88OMDRCdH42KwLU6QZc0i0DuQN4e8ib+Xv0P6AfAx+aCUctj98rp48SJr165l2bJl7Ny5E4BGjRpx6623MnnyZHr37l3me507d47m2Xv67Nmzh44dO+a7/tZbb/HAAw8wevRovv76a7TW1K9fn8TEwtMWn3rqKRYsWFCuz5KUlMS6detYsWIFP//8M1prvL29GTduHFOnTmXgwIFV9ucoilba38mCOnXqRHh4eLjWulN5+3L4podl4ApYnNCvEEKIPBIyErjrh7vyJSIzu850aiICYGreHEO9elizf8Oa/s8/ePXt67R43E1GWjby4uj5ZADCYxIlGalmSeYk+q3p5/D7ZllzF/eMTo7mpq8du3T11tu24uvq+ET6ySef5PXXXycjIwM3NzfGjh3LlClTGDlyJKbsBSDKY+XKlWRkZDBz5sxCiQhgP3fx4kXAtrBDWUZAysrHx4dp06Yxbdo0IiMjWblyJStXrmTZsmUsW7aM0NBQ1q1bx5VXXumwPkXNUa01I8qW1l4JxFZnv0IIIfJLyEjgzk13cujCIfu5OzrfwX1X3OfEqGyUUnhkz0mHmrH5YQepGxE1yI4dO8jIyMBoNPL444/z9ttvc8MNN1QoEQHYvHkzYKvlKEpO4uHn51exgMshODiY2bNn89///perrroKgFOnTkk9SR1WqWREKfVjziP71Ii85wo8fgXOAB2BHyobuBBCiIpJykzi7h/u5uCFg/Zz0zpN46ErH6oxUyE8uubdib1m1Y3IilrC2RYsWMC0adPw8PDghRdeICgoiBEjRrB69eoK1Vrsy074e/XqVeT1Xbt2AdCjR4+KB10Ghw4dYs6cOYSGhjJo0CB2795Np06dePnllxk4cGCV9i2cp7LTtAblea0B/+xHcczAt8CsSvYrhBCiApIzk7l78938E/+P/dykDpN4tMejNSYRAXDvmn9kRGvt1Pg6BOROy5JkpPr5mHzYetvWSt8nb71ITo0IUOico+pGfExVM52vT58+9OnTh7feeot169axfPlyNm3axMaNG/Hy8mLs2LFMnjyZsLCwUuf6A8TG2ias+BZTm/Xll18CMGrUKId9hhznz59n7dq1rFy5kt27dwO2EZj777+fqVOn2kdHRN1V2WQkZ78QBZwAPgMeL6ZtJhCntTYXc10IIUQVSjGncM/me9gXmzvt6bb2t/FEzydqVCIC+UdGLHFxZMXEYAoMdFo8eZf3PZ+UQXxyBn7ebk6L53KjlKp07UVMcgwP/fgQ0cnRhVbNWjZ8mX1FrYd+fKjWrKjl5eXFlClTmDJlChEREaxcuZLly5fbay4CAwOZMGECkyZNolu3bsXex8fHh4sXL3Ls2DE6dOiQ79p7773HwYMHGTBggH1VrcpKS0vj66+/ZuXKlWzcuJGsrCxMJhNjxoxh6tSpjBo1CldXV4f0JWq+Sk3T0lqfzn6cAuYDq/KcK/iIkURECCEcb1X4qhJ3lQZINacyY+MM/or9y37u1na3MrvX7BqXiAC4NGyIKSjIfuzsupHGPm74eeV+OZK6kdqltOV7S9uHpDYICQlhzpw5HDlyhK1bt3LXXXeRmprKokWLuOKKK1i/fn2x7x0wYABgW97Xas1dunrDhg089thjeHl58fbbbzss1jFjxjB+/HjWr19P586dWbx4MVFRUXz11VeMHTtWEpHLjMMK2LXW87XWXzvqfkIIIUq3KnwVC3ctLPHLU1pWGnf+cGe+qVk3tbmJp3s/XSMTkRw1qW5EKSV1I7XYlogtpe4jUjAh2RKxxQmROkbfvn353//+R0xMDGvWrGHEiBGUtJXDs88+i4eHB2vWrKFTp05MmDCBq6++muuuuw4XFxe+/PJLOnfu7LD4/Pz8ePTRR9m3bx979+7loYceonHjxg67v6hdqmxpX6VUfcAH2xSuQrTWsiyCEEJUUlhIGKsPrrb/NrfgF630rHTu+uGufFOzhjUfxtw+czGoal1Qsdzcu3Yh8bvvAEjb97eTo7HVjfx+LA6QZKS2mdRxEmD7+1LS9KuchGRLxBb7e2ozd3d3xo8fz/jx40tsd8UVV/Dbb78xd+5ctm3bxrp16wgJCeGhhx7i8ccfL3IX9spYu3atQ+8najeHJiNKKX/gBeAGoGEJTbWj+xZCiMtRzpennCkoeROSDEsGd/9wN3+d/8vefmjIUF4e8HKNT0QAPLrmznFPPxCOzspCuTjvR0d7/9yRkXBJRmqdsiYXAd4BdSIRKa8ePXqUOJVLiKrisJ9GSqkAYDcwHUjHtpeIAv4AzpM7QrId+M1R/QohxOWuqPnupxNPc88P9/Dn+T/t7QYHD2bRwEUYDUYnRlt27h07QHbyodPSyDh2zKnx5J2mdTw2mcwsawmthRBClIUjfzX2DBAIzNVaBwPfA1pr3U9rHYBtGeBD2EZFrnNgv0IIcdkrmJCM/mI0u87tsl8fGDSQ1wa9VmsSEQCDuzvubdvaj51dxN66iTcmo+33amaL5lj2juxCCCEqzpHJyAjgpNb6haIuaq1/BYYB3YF/O7BfIYQQ2BKSVwe+ikKhyS1WvabZNSwevBgXQ+2bHVtwvxFncnUx0Kqxt/340FmZqiWEEJXlyGSkGfBXnmMLgFLKvhC71joK+Am4xYH9CiGEAE4mnOShnx7Kl4goFLN7z66ViQgUqBupATuxd5QVtYQQwqEcmYwkkn/lrEvZzwWXYEgv4pwQQohK2HNuDxO/m8jZ1LP2cwYMaDQzN82sdXsm5PDIMzKScewY1pQUJ0ZDgeV9Za8RIYSoLEcmIxFAaJ7jnAXtR+acUEp5Av2A2vlTUQghaqCNpzYyY+MMkjJzvxzP6DKD78d+X6s3cQNwbdECg5eX7cBqJe3AAafGU3CvkZL2bhBCCFE6RyYjPwKdlVJNs4+/BlKARUqphUqpB7BN0WqKrbhdCCFEJWitWfbPMmb9MguzNtvPP93raR668iECfQJr/a7SymjEvUvu6Ei6k+tGOgT42F/Hp2QSm5ThxGiEEKL2c2QyshpYB3QA0FpfAGZmX3scWAz0BMKBOQ7sVwghLjtZ1iz+s+M/vPrnq/ZzCsWrA1/ltg632c8VtexvbUtIPLrkLWJ3bt2In7cbTXzspZCy34gQQlSSw5IRrfXfWuvbtNY/5zm3BmgL3Itt6d+bgSu11gmO6lcIIS43qeZUHv7pYT4+/LH9nFEZ+WDYBwwLHVaofcGEZEvEluoMt9I8unW1v3b2ilogdSNCCOFIVb68itY6AvhvVfcjhBCXg7i0OO7fcj8H4nNrJ/w9/Vk8eDGdGnUq9n05CcmWiC21bndp9y65yUjW2bOYz53H1LSJ0+LpEODLL0diAVlRSwghKsuRO7D/qJR6ogztZimlfnRUv0IIcbk4cekEk76blC8RGRQ8iK9u/KrERCRHgHdArUtEAExNm+Di728/Tv/HuVO18taNSDIihBCV48iakUFA+zK0awcMdGC/QghR5+0+u5tJ308iKjnKfm5C+wksHrQYT5OnEyOrHvnqRv527lStvHuNnIhLId1scWI0QghRuzkyGSkrdyDLCf0KIUSt9N2J77jrh7vsS/cqFI9f9ThP9XoKo8Ho5OiqR766kf3OTUZaNPLC1cX249Ni1Rw9l+zUeIQQojar1mREKeUL9EX2GRFCiFJprVmyfwlP/vYkZqtt6V5XgyuLBi5iSqcpKKVKuUPdkbduJH3/P2ir1WmxuBgNtGsqU7WEEMIRKpWMKKVO5DyyT43Le67AIwKIxbYx4leVjDunf3el1Hyl1BGlVLpSKloptVQpFVTJ+7ZRSqUppbRSaoMjYhVCiPLIsmbx/B/Ps3jPYvu5+m71WTJ8SZErZtV1Hp07gcH2I8uanEzmyZNOjSdv3Ygs7yuEEBVX2dW0QvO81oB39qMoZiAa22aIsyvZL0opd2ALuSMtX2XHczswSinVR2t9vIK3/x/gVmorIYSoAqnmVGb9Movfon6znwv2Cebdoe/S3Le5EyNzHoOXF26tWpFx9Chgqxtxa9XKafEU3IldCCFExVRqZERrbch5AApYlvdcgYeb1rqF1vohrXWqA2J/Glsish1oq7W+VWvdG3gMaAwsrchNlVJ3AIOB9x0QoxBClEtsaizTNkzLl4h0bdyVVSNXXbaJSA73GlQ30t4/fzKitXZiNEIIUXs5smbkdmCJA+9XLKWUCXgg+/A+rbW9elBr/RqwDxiglOpRzvs2AV4BNgNrHBSuEEKUybGLx5j43UQOXjhoPxcWEsaSYUto6N7QiZHVDB5560Zq0IpaielZRCekOzEaIZxn8ODBjBo1Kt+5hIQEunXrxsCBA0lPr/zfjeroQziPI3dgX6613uqo+5XiGqA+cFxrvbeI659lP48u533fADyAeyoemhBClN/OmJ1M+X4KMSm563tM6jCJVwe+iruLuxMjqznyrqiVfuQIVid+AannaSKwXu5/l4PRMlVLXJ7mz5/P+vXr2blzJwDp6emMGTMGpRRff/017u6V//erOvoQzuPITQ8nZheqX1tCm2HZbW6tZHfdsp/3FHN9T4F2pVJKjQRuBV7UWh+rRGxCCFEu3xz/hpmbZ5Jkzl2698meT/Jkrycvm6V7y8KtdWuUh4ftICuL9PCDJb+hiuWtGzl0VpIRcXkaMGAAYWFhzJs3D4vFwvjx4zlz5gwbNmygXr16taYP4TyVLWDPazLgBfxUQpsfsRW4TwU+rkRfIdnPZ4q5fqZAuxIppbyAd4DDwMJKxIVS6kAxl5xXaSmEqJG01ry//33e3Pum/Zyb0Y2F/RcS1jzMiZHVTMrFBfdOHUnb/ScA6fv34Xlld6fF0yHAly2HzgNwMCbJaXEI4WzPPfcc/fr1Y/jw4fzzzz9s3boVf3//WteHcA5H1ox0BvZprYvd0DD72t/ZbSsjZ8Wu4grhUwq0K80LQHPgHq11ZmUCE0KIsjBbzczfPj9fItLArQFLhi+RRKQEeetGnL0Tu6yoVbtcWLECc3R0qe3M0dFcWLGiGiKqO/r27UuXLl347bff2LBhA62qYKW76uhDOIcjR0YaAefL0O480L+SfeXs9FXc8iVl3glMKXUVtmL4FVrrkkZ1ykRr3amYfg4AHSt7fyFE7ZdiTuGxnx9ja3RumV1z3+a8G/Yuwb7BToys5su/E/t+J0aSf6+Rk/EppGZm4enqyB+rwlEurFjBuRcXcGHlKpovX4YpMLDIduboaE5PnYY5MhKAhlOmVGeYVeb48eOYTCZCQso0YaTcFi5cyNGjR8nMzCQ5Obn0N9SQPrZv30737t2l5sTJHDkyEk/ZpiK1Ai5Vsq+c8XCvYq57Zj+X+H+rUsoF2xK+CcCsSsYkhLiMrQpfRUxyTKnt9sftZ8yXY/IlIlc0voKV162URKQMPLp0sb82R0aSdeGC02Jp7ueFh8lW06M1HD4rU7VqKp+hQzEFB2OOjLQlG0WMkORNREzBwfgMHeqESB0nMTGRJUuW0L9/f1q3bs2ePcWV2ebSWrNq1SqGDh2Kn58fRqMRpVShx6RJk+zv+eCDD5g7dy5ffPEF1113HXPnznV6H2U1e/Zs/P39mTlzJtu2bXPIPUX5OTIZ2Qr0VEoVO+qhlLoG6AVU9r94RPZzcTutBxVoV5wg4AogE/hUKfVzzgNYnN2mV/a5byserhCiLlsVvoqFuxYyfeP0EhOSbVHbmPzdZM6n5g4iX9v8Wt4f9j4N3BtUR6i1nktgIMZGjezH6U4cHTEaFO38c0dHpG6k5jIFBtpGRIpJSAomIiWNntRkVquVTZs2MXHiRPz9/ZkxYwZbt25l4MCBtG/fvtT3jh8/nsmTJ7N161Z69uzJzTffTLNmzQq17ZL9S4F169Zxzz33sGzZMkaMGMFzzz3HTz/9xC+//OK0PsrjhhtuwNPTk/fee49+/frRtm1b/vOf/xARUdrXR+FQWmuHPIA+gAW4CDwMeOW55pV9Lj67Tf9K9jUY2xStY8Vc/3f29fml3Cc0u11ZHpcqGfOBjh07aiFE3ROdFK1HfDZCd17WWY/4bISOToou1Obb49/qLsu66M7LOtsfi3Yt0harxQkR124Rd9+jw9u11+Ht2uvzr7/h1FieWrdPN3/yW938yW/1M1/sd2ostZnFYtHh4eE6PDxcWyxV93ciMypKHx16rQ5v114fHXqtzoyKKvJcbXPw4EH91FNP6WbNmtm/t7Rr106/8MIL+tSpU2W6x7PPPqsB3alTJ33y5En7+ZSUFN29e3cN6IULF+q0tDSdlZWlf/jhB+3m5qbfeCP/38EbbrhBDxgwwGl9lFdWVpbesGGDnjBhgvb09NSAVkrpIUOG6BUrVujk5GSH9FPblPfvZMeOHTVwQFfkO3JF3lTszWwJhyXPIwaILnDuUQf044ptqpcGuhdx/e/saz0r0ceg7HtscNCfjSQjQtRhJSUky/5Zli8J6bKsi14dvtqJ0dZuse++a09GTs+406mxrNh20p6M3PTOVqfGUpuV9sXHarXqrIQEhzzSDh/WR4eE2f4f6tTZ9mjXXh8dEqbTDh92WD9ZCQnaarVW2Z/ZhQsX9DvvvKN79eplT0AaNWqk77vvPv3HH3+U615nz57Vbm5u2s3NTR84cKDQ9TfffFMDevTo0VprrXfs2KG9vb313LlzC7X966+/tFJKb968udr7qKzExET94Ycf6sGDB2ullAa0t7e3njZtmv7pp5+q9L9nTVOdyYhDK+201ouVUnuAp4CBQNPsS2nAz8BCrfWvDugnUyn1FjAHeEspNUxrnQKglHoU6Ar8rrXelfMepdT9wP3AF1rr2ZWNQQgh8grwDmDp8KVM3zidM8lnmL5xOkuGLWHVoVWsDF9pb+dmcOOVga8wOGSwE6Ot3dzz1I2k79uH1hqlyrxuiUPl32skCatVYzA4J5a6zJqUxJFevR1/46zcBUDNUVGcHHODQ2/fducOjL6+pTcspyeffJLXX3+djIwM3NzcGDt2LFOmTGHkyJGYTKZy32/lypVkZGQwc+ZMOnYsvNZOzrmLFy8C0KtXL5KSip6W2K1bN6xWq1P6qCwfHx+mTZvGtGnTiIyMZOXKlaxcuZJly5axbNkyQkNDWbduHVdeeaXD+76cOXzZj+xk41ellAHbClsAcVprR/9f8wIwFOgLHFVK/YZted7e2KaD3V6gfSOgHRDg4DiEEAIonJCM+HwEOs+if/Xc6vFu2Lt0adylhLuI0uQtYrckJGCOiMC1eXOnxNI+TzKSnJHFmYtphPh5lvAOISpvx44dZGRkYDQaefzxx7nvvvsqtefG5s2bAVsNRVESEhIA8PPzq9F9OFJwcDCzZ8+mX79+zJo1i927d3Pq1CkiIiIkGXEwRxaw56O1tmqtz2c/HJ6+aq3TsdWOPI9tv5EbsdWALMc2dUt2URdCVLsA7wBeG/QaCpUvEQnyDmLt9WslEXEAo68vri1a2I/T9jlvvxFvNxdCGuYmH+Gy34ioBgsWLGDatGl4eHjwwgsvEBQUxIgRI1i9ejUpKSml36CAfdl/h3r16lXk9V27bBNNevToUeGYq6MPRzl06BBz5swhNDSUQYMGsXv3bjp16sTLL7/MwIEDnR1enePwkZHs5XJHAT2xjUbs0Fovzb4WmH0uXJewOWJZaa3TgLnZj9LaPgs8W457/0w59isRQgiAf+L+4cEfH8yXiCgUiwcvJsinuAUARXl5dO1K5smTAKTt20+90aOdFkuHAB8iLtj24D0Yk8iIzrIrtKMZfHxou3NHpe9jPnuWM/fcizkqClOzZgS9+w5AoXMmB+3sbfDxKb1RBfTp04c+ffrw1ltvsW7dOpYvX86mTZvYuHEjXl5ejB07lsmTJxMWFobBUPrvnWNjYwHwLWZK2ZdffgnAqFGjKhxzdfRRGefPn2ft2rWsXLmS3bt3A7ZRmvvvv5+pU6dy1VVXOSWuy4FDR0aUUgOBE8A6YDYwA7gmT5MwYC/g2EmZQghRA3xx9AumfDeF2LRY+7mcEZKHf3q4TPuQiLJx75o7wpS2728nRiI7sVcHpRRGX99KPazJyZy5735b0hEcTPOVK3Bv2xb3tm1pvnKFbdnfqCjO3Hc/1uTkSvdn9PWt8lomLy8vpkyZwpYtWzh16hQvvPACgYGBrFy5kmHDhhEcHMzjjz/O33+X/HfEJztpOnas8KSS9957j4MHDzJgwAC6d++e79qmTZvo0aMH7u7uNGvWjGeeeQaLxeKwPj788EMGDx5M48aN8fHxoUePHqxevbrQ+8sTR15paWl8/PHHjBo1imbNmvHQQw/x999/M2bMGNatW0d0dDRvvvmmJCJVzGHJiFKqC/Ad0AR4HbiZwiML67BNqbrJUf0KIYSzmS1mnt/+PHO3zcWszfbzd3e9mw1jNxDkHWQvapeExDE8unazv84IP4jOzHRaLO398yQjZyUZqYlK20ektH1IaoOQkBDmzJnDkSNH2Lp1K3fddRepqaksWrSIK664gvXr1xf73gEDBgDw/PPP5ysM37BhA4899hheXl68/fbb+d6zd+9eRo0axYABA9i7dy9vvPEGb731FnPmzHFYH1u2bGHMmDF899137N27175Hyccff1zhOPIaM2YM48ePZ/369XTu3JnFixcTFRXFV199xdixY3F1dS31HsIBKrIEV1EP4FMgCxiS55wVWFqg3U7goKP6rS0PZGlfIeqkcynn9KT1k/It3dt1eVf9/Ynv7W3Ksg+JKB9rRoY+2KWrfYnf1H3O2+PjdFyKfXnf5k9+qxPTMp0WS21V1fuMxC9fXqZ9RPLuORK/fLnD46huaWlpes2aNXrEiBH6m2++Kbbd3r17tYeHhwZ0+/bt9W233aZ79+6tAV2/fn39ww8/FHrPhAkT9JVXXpnv3Ouvv649PDx0UlKSQ/ooynXXXafHjh1b4TjyuvXWW/Wjjz6q9+3bV6a+LyfVubSvI6dpDQT+0Fr/WEq7CKD2bWsqhBAF7D2/l1u/vZW/Yv+ynzMZTKy8biUjWoywn8tZZUtGSBxHubri3qGD/diZU7WCGnjg7ZZbgnnorOzEXtM0nDKFpk/PLnVn9ZwRkqZPz6bhlCnVGGHVcHd3Z/z48Xz//fcl1mJcccUV/Pbbb4wcOZKzZ8+ybt064uPjeeihh/jnn38YOnRoofds3bqVkSNH5jt3/fXXk5aWxp9//umQPopy6dIlGjVqZD8ubxx5rV27lldffdW+27twDkcWsPsCUWVo5wYYHdivEEJUK601Hx/+mIU7F5KVZy0OLxcvPh79Mc19Cy8zW3DZ3y0RW5jUcVJ1hl3nuHftSlr2XPj0ffthonPiMBgU7f192H3atj/CwZhEeoY2dE4wolhlTS5MgYF1IhEprx49epQ4laugmJgYAgLy75aQs7xwdDFT3MrbR0HLly9n9+7dvPX/7d13eJRV+v/x90ky6aG3BAIiiBhEEKQqKIJiw766WBG77lrXta0FC6w/Xdu67asiKIqFZRdZCwq6K0UBsYAEpUNCQg0lhfTz+2Mmk07azDyZzOd1XXNN5jzl3PMwPMk9p73ySpPikObFl8lIJnBcnXvB8cA2H9YrIhIwBSUFPPn1k8zbNM9bFhEWwenJp3Pv4HtJSqj9W9eyhESJiG/EnNCf/Z6fnZzeF9yD2MuTEbWMSGgqG7Dvj4H78+bN4+abb+Yf//hHnet8+DMO8T1fJiOfATcYYy6y1v6rph2MMZNwL0z4nA/rFREJiMycTO76712k7kv1lnWM6cjzpz3PwE4D63WOxPhEJSI+EnPCCd6fC7dsoeTQIb+sdl0fmlFLQk1iYiKZmZW7m5a9rtpS0VTvvvsukyZN4u9//zuTJk1yLA7xD1+OGZkK5ACzjTFPGmPK5kGLNcYcb4z5A/BX3KujP+/DekVE/G555nIu/8/llRKREzudyHvnvVfvRER8y9W9O+GtW3tfH16zxrFYjkssX0/il53ZlJTaI+wtEvxOPvlkPvnkk0plH3/8MTExMT5duPDVV19l0qRJzJgxo1oiEsg4xH98loxYa7cB5wL7gYeB5YDFPcXvj8ATQDZwgbV2p6/qFRHxJ2stM9fO5KbPb2J/wX5v+a+P/TWvn/k6HWM7OhhdaDPGEF2hdSTfwWTk2C4JlPUIOVxUwrZ9DV8FWySY3HvvvaxevZp7772XdevWMXfuXB599FHuuOMO4uPjfVLH888/z6233spLL73Eaaedxs6dO9m5c6d3AcVAxSH+5dNFD621S4A+wD3AJ8A6YD3wBe5FEI+11i7zZZ0iIv6SV5TH/V/dz3PfPkepdc+LHxkWyVMnP8XDwx/GFe5yOEKJqTALzuEfnRs3EhsZQc/2cd7XGjciLd2gQYOYP38+X375JQMHDuQ3v/kNt912G08//bTP6nj55ZcpKSnhlltuITEx0fsYMmRIQOMQ//LlmBEArLXZwIueh4hIUEo7lMad/72TDfs3eMsS4xJ5YcwL9Gvfz8HIpKKYAeUtI4fXrMFa69ig1eMSW7F5r7tFZF3mIc49Qf3VpWUbP34848eP99v5t27d2iziEP/yacuIiEhLsGTHEi7/6PJKiciwLsN497x3lYg0M9EVWkZK9u6l2MGpPCuOG9EgdhGR+vF5y4gxpjMwCTiF8sUNM4ClwExrrVb6EpFmqdSW8tqa13jl+1ewlA9AntRvEncOupOIMJ/fMqWJItq1w5WcTFFaGuBuHXF17epILEltYrw/f7VhD1e/vpzhR7fnspOS6ZgQ5UhMIiLNnU9/sxpjrgFeAeKAiu3kJwLnAA8bY+6w1r7hy3pFRJoqpzCHh5c8zBdpX3jLYiJieGLkE5VWU5fmJ6Z///Jk5MfVtDorsP9e+UUlTJm/lg++TfeWFZVYFm/Yy+INe3lx4XouHZzMYxNSiHZpzV8RkYp81k3LGHM28AYQA7wPXAgMxJ2IXOApiwZe8+wrItIsbD64mSs+vqJSItItvhuzzpmlRCQIVB43EthB7PlFJVw7fQWzV6RRXMt0vkUlltkrtnPt9BXkF5UEND4RkebOl2NGHsY9le8F1tqJ1toPrbWrrbU/WmvnW2sn4k5KAB7yYb0iIo22aPsirvjoCrYc3OItO6XrKbx73rv0advHwcikvqL7V5jed20qtrg4YHVPmb+W5Vuy6rXv8i1ZTJmfWveOIiIhxJfJyEBgsbX249p28Gz7CndriYiIY0pKS3j5u5e568u7yC0qXxPiphNu4pXTX6F1VOsjHC3NSXTKcRDh7nVsDx+mYOPGgNS7OzufOavS696xgjmr0tiTXeCniIJXxRnQSktLHYxERKDy/0N/z1Doy2QkH/dA9bpkevYVEfGZWamzyMype36MzJxMXl39Krd/cTuvrnnVWx7niuPFMS/y2xN/S3iY+vUHk7DoaKL7lLdiBWq9kfdXplFU0rCV1otKLO9/m+aniIKXMYbIyEgAcnO1YKSI08r+H0ZGRvo9GfHlAPb/AcOMMWHW2hq/1jDGhAPDcLeOiIj4xKzUWTyz8hneXvc208dPJzG+5vUdMnMyufqTq9mVt6tSec/WPXlxzIsc3froQIQrfhA94ATyU91doA6vWU3byy/ze5317Z5V1Teb93H7mN4+jib4JSQksG/fPnbtcv//jIuLIyxMKxCIBFJpaSm5ubne/4cJCQl1HNF0vkxGHgBW4B6gfre19mDFjcaYVsALQFvcq7GLiPjE2O5jeXvd26TnpDN5weQaE5LMnEx+/dGvycrPqnbsUyc/RXxkfCBDFh+L6X8CB2a/C0B+gFpGcgoaNzalsce1dO3btyc3N5f8/HwyHFwvRkTcoqOjad++vd/r8WUyMhGYj3uNkUuMMZ8B2zzbegBnAvHA28DlVZp8rLX2SR/GIiIhJDE+kenjpzN5weQaE5K07DQum38ZOUU53mMMht+e+Fuu7389YUbfvga7ijNqFWzcSElOLuHxcX6tMz6qcb9CG3tcSxceHk737t3Zt28f2dnZFBYWOh2SSEiKjIwkISGB9u3bEx7u/27LvrwjPg7eVcISgEtq2e+qGsosoGRERBqttoTkQMEBrv7kagpKygcNJ0Qm8MyoZxjVbZSDEYsvRfbsSVh8PKU5OWAt+WvXEjdsqF/rHNazHYs37G3wccOP9v83jcEqPDycTp060alTJ6y1WNuwMTki0jTGGL+PEanKl8nIdT48l4hIg1VNSM6eezYltvK6Dse0PYaXTnuJ5FbJDkUp/mDCwojufzx5X38DQP6a1X5PRi4bksxLizY0aBB7eJjhspP02asPJ/4oEpHA81kyYq2d6atziYg0VmJ8Iq+f+Tpn/6t6InL2UWfz+MjHiXXFOhSd+FNM/xO8yUggZtTqlBDNpYO7MXtF/WfHstby046DjOnbyY+RiYgED3WUFpEWJbswm6dXPE1plUn9bup/E8+MfkaJSAtWeSX2NQGp87EJ/RjWs1299y+1cPNbq1iYuqvunUVEQoDPkhFjzL313K+dMeYDX9UrIlJm7d61XDzvYr5Krz57+MdbPmZn7k4HohJ/y3rzTYoyMoju399bVrxzJ0W7dlfarygjg6w33/Rp3dGucGZOHsrEod1xhdfcpcgVbjjjuE7ERboHghaWlHLLrFV8+pM+jyIivmwZedYY87kxpmttOxhjxgI/Ahf7sF4RCXHWWmalzuLKj69kZ175H3gDOg5g9rmz6RbfzTuovT4LI0rwyHrzTXZNnca2aydBcTERieVTOuevKe+qVZSRwbZrJ7Fr6jS/JCTTLu7PsgfGct/4Yxl1TAdO7N6GUcd04L7xx7LsgbG8eu0Q3rphGAmembSKSy23v/MdH63W51FEQpsvk5EvgbHAj8aYSytuMMa4jDF/AhYAicBzPqxXRELYwYKD3PHlHTyz8plKY0SuTrmamWfN5PgOxzN9/HQlJC1UwrhxuJKTKUpLY9u1k4jqXb6Y4OHV7q5aZYlIUVoaruRkEsaN80ssHROiuH1Mb966fhj/uu1k3rp+GLeP6U3HhCgABnVvy9s3DqNVtDshKSm13PHu98z7YYdf4hERCQY+S0astWOB3+NeS+Q9Y8wbxph4Y0w/YCVwN5AOjLXW3u+rekUkdP2w+wcunX8p/037r7cszITx7Ohn+f2Q3xMe5u4WUzbLlhKSlseVlESPmTO8Ccnh1eWtIYdXr66WiPSYOQNXUpJj8Z7QrQ3v3DicNrEuwJ2Q3P3eD8z9Lt2xmEREnOTTAezW2ueAYcDPwDVAKu5E5ATgHWCAtfZ/vqxTREJPqS3l9TWvM+nTSZXGgUSHR/Puue9yVs+zqh1TNSFZtH1RIEMWP6qYkJQePOgtP/zjj80qESlzfNfWvHPDcNrFRQLuQe33fvAj76+s/6xcIiIthc9n07LW/ghMALKBbkAU8Ja19ipr7cEjHiwiUod9h/dx28LbePG7F73dssJNOKd2O5V/X/Bvjmt/XK3HliUk9w+5n6tSalp/VYKVNyHpWj5s0R4+3OwSkTIpSa2YfeNwOsS7ExJr4ff/XM07y7c7HJmISGD5PBkxxlwIfAO0AlYBhcBVxpiZxpgEX9cnIqFjReYKfjX/VyzNWOotS4xLZMZZM3hl7Ct0Tah1/ozy/eMTlYi0UK6kJHq8VWVwelhYs0tEyhzbJYF3bxruHVMC8NC/1vDm11udC0pEJMB8ObVvrDHmNeCfuBOR3wNDPY+1wNW4B7ef4qs6RSQ0lJSW8Ncf/soNn93AnsN7vOWnJ5/OBxM+YGCngc4FJ81PWIVfbbb+q6M7oXenBN67aThdWkV7yx6dt5bpS7Y4GJWISOD4smXkB2Ay7vEiw6y1z1m3NcBJwEtAD+BLY8xUH9YrIi3Y7rzd3Pj5jfztx79hcf9h6Qpz8cDQB3hxzIu0jmrtcITSXJQNVqe0woKX1rLtmmspyshwLK66HN0xnvduHk5S6/KE5In/pPJ/X21yMCoRkcDwZTLSG3gFGOwZN+JlrS201t4NnAXsBjSblojUacmOJVz64aWs3LnSW5ackMxb57zFlcddiTE1LzInoafSrFndukFkZPm29HT3tmackPRoH8d7N4+gW9sYb9nUj3/mL19udDAqERH/82Uycq619g5rbX5tO1hrP8c9s9a/fFiviLQwRaVFvLDqBW5deCv7C/Z7y88+6mzeP+99+rXv52B00txUm773zZnEDRvm3R6WkOBdh6Q5JyTJ7WJ57+YRdG8X6y17dsEvvLRwg4NRiYj4ly/XGfmknvvts9ZeWveeIhKKMnIyuO7T65j+03RvWVR4FI+PeJxnRj9DfGS8g9FJc1PbOiLxo0Z594no0KHSwojNOSHp2iaG924eTs8Ocd6yFxau5/nPfsE28/EvIiKN4Y/ZtCKMMRcaY542xvzDGDO5wrYkY8wJxpgIX9crIsHvi+1f8Kv5v+LHPeU9PY9ufTTvnPsOl/S5RN2ypJrshQtrnL43fnR5MlK4ZQtJzz3rTUiyFy50Ktx6SWwdw7s3DadXx/KE5OUvNvL/FighEZGWx6dJgTHmVOAtoCtgAAu4gLKvOMcCM4DLcM+6JSJCYUkhz696nrfXvV2p/MLeF/Lg0AeJdcXWcqSEunbXXANAwrhxlabvdfXo4U0+AArWr6fHzBlkL1zoPaY569wqmndvGsEVr37Dht05APztv5soLinloXOOU2IuIi2GL6f27Q98DHTCPXPWr3AnJBX9E8gDLvFVvSIS3LYf2s5VH19VKRGJiYhh6ilTefLkJ5WISJ3aXXNNtXVEjDGVumrlLl6CKykpKBKRMh0Tonj3puH07VK+RNeri7fwxH9S1UIiIi2GL7tpPYp7tfVzrLX3WGurtXxYa/OAdcCJPqxXRILUJ1s+4bL/XMa6rHXesmPbHst7573HhF4THIxMWoK4Cl21cpctwxYVORhN47SPj+KdG4eTktjKW/bG0q08Om8tpaVKSEQk+PkyGTkV+MZa+0Ud+20Hmt9SuCLSJLNSZ5GZk1nnfpk5mbzx0xtM+XoKv//q9+QW5Xq3XX7s5bx97tv0bN3Tn6FKiIgbOhTjcgFQmpPD4R9+cDagRmoXF8k7Nw6jf9fyNXXe+mYbD/97jRISEQl6vkxGWgE76rFfFBDuw3pFxGGzUmfxzMpnmLxg8hETksycTK755BqeX/U8c9bP8ZbHu+L506l/4g/D/0BUeFQgQpYQEBYbS+yQId7XOYuXOBhN07SJjWTWDcMYmNzGWzZ7RRr3/3M1JUpIRCSI+TIZyQSOq8d+xwPbfFiviDhsbPexdIvvRnpOeq0JSWZOJpd/dDk783ZWKj++/fG8P+F9zjzqzECFKyGkYletnMWLHYyk6VrHuHjr+qEM7tHWW/bBqnTu++BHSkote7ILeOWLDVz9+nIu+utSrn59OX/5ciN7sgscjFpE5Mh8mYx8BvQzxlxU2w7GmElAD+AjH9YrIg5LjE9k+vjptSYkmw9s5sJ5F7I/f3+l465JuYY3z36T5ITkQIcsISJ+9GjvzwXr1lG0e7eD0TRdQrSLmZOHMvSodt6yud/vYMxz/2XEHxfx3GfrWbxhL99vP8DiDXt5dsEvjPzjIh6cu4b8ohIHIxcRqZkvk5GpQA4w2xjzpDHmJE95rDHmeGPMH4C/AvuA531Yr4g0A7UlJEvSl3DJh5eQV5zn3bd1VGteOf0V7htyH65wl4NRS0sX2bNnpZm2coO4q1aZ+KgIZkwewoij23vLtmflUVxSc3etohLL7BXbuXb6CiUkItLs+HIF9m3AucB+4GFgOe51Rn4F/Ag8AWQDF1hrd9Z2HhEJXlUTkrPmnsWti26l2BZ79xnUaRBzJszh1ORTHYxUQoUxpnJXrSXB3VWrTGxkBNMnDaFLq/qPsVq+JYsp81P9GJWISMP5dAV2a+0SoA9wD/AJ7ml81wNfAA8Cx1prl/myThFpXhLjE3n+tOcxGEptaaVtN/a/kdfHv06XuC4ORSehqGJXrdyly7DFxUfYO3hkFxSxN6ewQcfMWZWmMSQi0qz4dAV2AGttNvCi5yEiIWZx+mIeWvIQlspdRqadMo3zep3nUFQSymKHDgOXC4qKKD10iMOrVxM7aJDTYTXZ+yvTKG7gTFpFJZb3v03j9jG9/RSViEjD+LRlRERCV35xPlOXT+W2RbdxoOCAt9xgAPjLD3+p1zokIr4WHh9H7ODB3tfBPqtWmeVbshp13Deb9/k4EhGRxlMyIiJN9nPWz1z+n8uZ/fPsSuW3DriVTy7+pM5pf0X8LX5UhdXYv2oZyUhOQeO6mzX2OBERf1AyIiKNVmpLeeOnN5j40UQ2H9zsLXeFufi/cf/HbQNvo2tC1yNO+ysSCPEVBrHnr11L8d69DkbjG/FRjetp3djjRET8QcmIiDTKztyd3PjZjTy/6nmKS8u/aY13xfOvC/7FiK4jvGV1rUMi4m+RvXsT0aV84oScJcE/xe+wnu3q3qkGsa5wrdouIs2GkhERabBPt3zKxR9ezIqdKyqVd4zpyNzz59KjVY9qx1RNSBZtXxSocEUwxlTuqtUC1hu5bEgyrnDT4OMWpO7igr8s4dutjRtzIiLiS2qrFZF6yynMYeryqczfPL9S+Zk9zqRvu76cd/R5JMYn1np8WUKyaPsirkq5yt/hilQSN3oUBz74AIDcJUuwJSWY8HCHo2q8TgnRXDq4G7NXpDX42J92HOLSv3/NhQOTeODs4+jSOtoPEYqI1E3JiIjUy/e7v+fBxQ+yI2eHtyw2IpaHhj3E+b3Ox5j6fUObGJ+oREQcETdiBEREQHExJQcPkr9mDTEDBzodVpM8NqEfm/fk1mtmrRO7t6Frmxj+s7q8i+S/f8jgs9Rd3D6mNzeM6klURPAmZyISnNRNS0SOqKi0iD9//2cmfTqpUiIyoOMA5pw/hwt6X1DvRETESeHx8cSeeKL3dU4LmFUr2hXOzMlDmTi0e61dtlzhholDuzP7xuG8csUgPrhlBP2SWnm35xWW8OyCXzjzha/4PHUX1mo8iYgEjtFNJzCMMWtTUlJS1q5d63QoIvW27dA2Hlz8IGv2rvGWhZtwbh5wMzf2v5GIMDWuSnDZ++qr7PnT8wBEn3ACPd9/z+GIfGdPdgHvf5vGN5v3kVNQTHxUBMOPbs9lJyXTMSGq0r4lpe7FD59d8AtZuZVXcR/dpyOPnpdC707xgQxfRIJYv379SE1NTbXW9mvosUpGAkTJiAQTay1zN8zlmZXPcLj4sLc8OSGZaaOmMaDjAAejE2m8/F9+YcsFF7pfGMMxS5cQ0a5xs1K1BAfzinhx0Xre/HpbpRm2IsIMk0YexR3jjqFVtMvBCEUkGDQlGVE3LRGpZH/+fu768i4e//rxSonIRb0v4oMJHygRkaAW1acPEZ06uV9YS+7Spc4G5LDWsS4em9CPT+4cxcm923vLi0stry3ZwunP/Zf3V6ZRqqmARcRPlIyIiNfSHUu5+MOL+SLtC29Z66jWvHDaCzxx8hPEueIcjE6k6YwxxI06xfs6Z3HwjxvxhT6dE5h1/TD+ftVgurWN8ZbvzSnk9/9czYV/XcqqbfsdjFBEWiolIyJCfnE+f1zxR25ZeAt7D5evTD0icQRzz5/LuB7jHIxOxLfiR432/py7ZCm2tNTBaJoPYwxnHd+Fhfecyr1n9CHaVf4nwur0g1zyt2Xc894P7DqU72CUItLSaPSpSIj7JesXHlj8ABsPbPSWucJc3D34bq487krCjL6zkJYlbuQICA+HkhJKsrLIX7uWmP79nQ6r2Yh2hfPbscdwyeBuTPvkZ+b/mOHdNvf7HSxYu5PfnH4Mk085qtJUwHuyC3hv5XaWb8mqcwC9iEgZDWAPEA1gl0CYlTqLsd3HHnHhQYDMnEwWbl9IqS3lpe9eoqi0yLutd5vePDP6Gfq07ePvcEUcs/XKqzi8ahUAHX77GzrefrvDETVfyzfv4/H5qazLPFSp/Kj2sTxyXgoje7Xnif+kMmdVOkUl1f+mcIUbLh2czGMTUoh2aR0TkZZIs2kFASUj4m+zUmfxzMpn6Bbfjenjp9eakGTmZHLtp9eSmZtZbdtVx13FXYPvIipc32JKy7b37/9gz4svAhAzcCBHvTvb2YCauZJSy+wV2/nTZ7+wP6+o0rbWMS4OHi6q5chyw3q2Y+bkoUpIRFqgkJ1NyxgTbYyZYoxZb4zJN8ZkGGOmG2O6NeAcbYwxVxhj3jHGpBpjco0x2caY5caYO40xmtNQgsLY7mPpFt+N9Jx0Ji+YTGZO9WQjMyeTiR9NrJaIdIzpyD/G/YP7h96vRERCQvzoUd6fD69eTcmBA84FEwTCwwxXDe/Bl787jWtH9CCswvqK9UlEAJZvyWLK/FQ/RSgiwSpokxFjTDSwCHgUiAfmAWnAdcB3xphe9TzV74C3gcuBPGA+sAIYALwIfGGMifVp8CJ+kBifyPTx02tNSDYd2MSF8y5kX/6+SseN7T6Wf57/T0Z2HRnokEUcE9W3L+EdOrhflJaSE+JT/NZXm9hIplxwPB/fOYrB3ds2+Pg5q9LYk13gh8hEJFgFbTICPASMBL4G+lhrL7fWDgPuBToC0+t5nhxgKtDdWnuStfbX1tqxQH9gO3AK8AefRy/iB7UlJAu3LeSSDy8hrzjPu29MRAxPjHyCF057gbbRDf+jQiSYmbAw4k8pn+I3d/ESB6MJPn27tOK0vh0bfFxRiXvldxGRMkGZjHi6Tv3W8/J2a21O2TZr7fPAamC0MWZwXeey1v7RWvuwtXZHlfINwAOelxN9E7mI/1VNSM6aexZ3//duSmyJd58TOpzAnAlzuOiYizDGHOFsIi1Xxa5aOUuWaIrfBlqxJatRx32zeV/dO4lIyAjKZAR3a0UbYJO19vsats/xPE9oYj0/ep6TmngekYBKjE9kysgpGAyltvwPLIPhlgG3MOPsGXRv1d3BCEWcFzdyJIS5fw2W7N1L/rp1DkcUXHIKiht1XHZ+444TkZYpWJORAZ7n72rZ/l2V/RrraM/zziaeRyRgSm0ps1JncevCW7FUni3v+dOe5/aBt+MK07wMIuFt2hBzwgne17lajb1B4qMat1TZ2oyDPLfgF9L359W9s4i0eMGajJR9pZtey/b0Kvs11p2e53n1PcAYs7amB1DfAfUijbYjZwc3fHYDz6x8hsLSQm+5wd0V60/f/qnGWbZEQlVcxa5aGjfSIMN6tmvUcUUllle+3Mjo//cl189YyRc/76KkVMsMiISqYE1G4j3PtX2tkltlvwYzxtwCjAMOAH9s7HlEAsFay9wNc7nkw0tYuXOltzzMhPHYiMdYcMmCOqf9FQlF8aNGe38+/MMPlBw6dIS9paLLhiTjCm/8mLNSC4t+3s3kGd8y+v99yStfbGB3dr4PIxSRYBCsyUjZ3a+2r1KaNCLXGHMq8JLn/JOttRn1PdZa26+mB7CpKTGJ1GZP3h5+88VveGzZY+QW5XrLYyJiePfcd7m0z6V1TvsrEqqi+6UQ3s7zDX9JCbnLljkbUBDplBDNpYPrvawXABcP6spD5/SlR/vKM+bvOHCY5z5bz8hpX3D729+xbNNetCizSGgI1mQk2/McV8v2srtcTi3ba2WMOQH4NxAJ3Gmt/VeDoxMJkE+3fMpFH17EV+lfVSpvF92Of5//b45rf5y3TAmJSHUmLIy4U072vs75SuNGGuKxCf3q3V1rWM92TL2oPzeN7sWX957GW9cP5ax+XQivsIJicanlozWZXPHqcsb+6X+8tngzB/IKj3BWEQl2jRt95rztnufavpLpVmW/evEslLgA90xdj1tr/9yo6ET87ED+AZ5e/jSfbv202rbEuERmnjWTxPjE6ts8CcnkBZNJz0ln0fZFXJVyVSBCFmm24keN5tCH8wH3IHZrraa8rqdoVzgzJw9lyvxU5qxKo6ikemuGK9xw6eBkHpuQQrQrHICwMMOoYzoy6piO7DqUz7sr0nh35XYyD5Z309q8N5enPlrHswt+4bwTkrhyeHdOTG5T67/NnuwC3lu5neVbssgpKCY+KoLhR7fnspOS6ZgQ5Z8LICJNZoKxGdQYMwb4AvfUvr1r2P4I8ATwhLX2sXqeMwlYAvQEXrLW3uW7iN0D21NSUlLWrl3ry9NKCPoq/SseW/YYew/v9ZbFRMRw35D7yC/OZ1z3cTUmIhVl5mQqERHxKN6/nw0jTwbP78Oe//4X0X37OhxV8NmTXcD736bxzeZ9jUoGiktK+fKXPby9fBv/W7+Hmv48SUlsxZXDu3PBwK7e2bzyi0qYMn8tc1al1zsZEhHf6tevH6mpqameoQkNEqzJSCSwG2gNDKq61ogx5kfgBGCotXZlDaeoer62wFfA8cAbwPXWxxdGyYg0VU5hDv9v5f/jXxsr9xwc1GkQT53yFMkJyQ5FJhL8tlx2OfmrVwPQ8Z576HDTjQ5HFNq278tj9srtvL8yjX251btpxUdFcOGJSVw6uBvTPv6Z5fVYgHFYz3bMnDxUCYmIHzQlGQnKMSPW2kLgFc/LV4wx3rEjxph7cCciSyomIsaY3xhjfjbGTKt4LmNMLPAx7kTkfeBGXyciIk21InMFl3x4SaVEJDIskt+d9Dumj5+uRESkieJHlU/xq/VGnNe9fSz3n9WXZQ+ezssTT6w2LiWnoJhZ32znwr8sq1ciArB8SxZT5qf6I1wRaYJgHTMC8BTuqXdHAhuMMYuBHsAwYB9wXZX9OwDHAlX7rzwNDAdKgGLg9Zr6o1prJ/kwdpF6OVx8mJe+e4m3171dqTylfQpTT5lKrzZavkbEF+JHj2LvX/4CQN7331OSk0N4fKNnhxcfiYoI5/wBSZw/IImNu7N5e/l25qxKb/Qq7nNWpXHPGX00hkSkGQnaZMRam+8ZO/IgcAVwIbAfmAk8Yq1Nq+ep2nqewz3nqc2kxkUq0jir96zm4SUPs/XQVm9ZhIngpgE3cUP/G7SKuogPRR9/POFt2lBy4AAUF5O7bBmtzjzT6bCkgt6dEnhsQj9+P74v81dn8Pzn69l5sGHrkhSVWN7/No3bx1QbbioiDgnKblplrLWHrbWPWmt7W2ujrLVdrLWTakpErLWPW2tN1RYOz/6mrkfA3pSEvKKSIl7+7mWu/uTqSolIr9a9mHXuLG4dcKsSEREfM+HhxJ1cPsWvumo1XzGR4Vx2UjLHdGpcy9U3m/f5OCIRaYqgTkZEWppfsn5h4kcTeXXNq5TaUgAMhuv6Xcd7E96jX/sGjwsTkXqKH10+biRn8RItutfM5RQ0rqvW2h0H+XhNZqOPFxHfCtpuWiItSXFpMTPWzuAvP/yF4tLyX5Dd4rvx9ClPM6jzIAejEwkNcaec4v25eOdOCjZsILpPHwcjkiMpm9q3obLyirjt7e+IDA9jRK/2jEvpzLjjOpHYOsbHEYpIfSgZEXHY1oNbeXjpw6zes7pS+eXHXs49g+8h1hXrUGQioSWifXui+/Uj3zMFe+7ixUpGmrFhPduxeMPeunesRWFJKf9bv4f/rd/DI/+G47u24ozjujAupRMpia0atPClFlwUaTwlIyI+NCt1FmO7j63XooMLty8kzITx4qoXyS8pH4TZKbYTT458kpFdR/o7XBGpIm70KG8ykvPVYtpff73DEUltLhuSzEuLNtS40GFtwsMMZ6Z0ZtmmfRw8XFRp2087DvHTjkO8sHA9Sa2jGZfSmTNSOjOsZ3siI2ru1X6kBRcXb9jLiwvXa8FFkTooGRHxkVmps3hm5TO8ve5tpo+fXmtCkpmTybWfXktmbma1bROOnsD9Q++ndVRrf4crIjWIHzWafX/7OwB5331HSU4u4fFxdRwlTuiUEM2lg7sxe0V9J8+Ey05KZtrF/SkuKeXbbfv5PHUXn6fuYntWXqX9Mg7m8+bX23jz623ER0Vw6rEdOeO4zow5thOtY90TiOQXlXDt9BVHXOekqMQye8V2Nu/J0YKLIrUIyhXYg5FWYG/5MnMymbxgMuk56XSL71ZjQpKRncGvP/o1+wv2VypvF92OR4Y/wrge4wIZsohUYYuLWX/yKZQePAhAt7+8QsLYsQ5HJbWpT0JQprYV2K21bNydw+frdrEwdRffpx2gtj+NwsMMQ49qx7iUzvywfT/zV1f/Uqk2E4d2Z9rF/eu9v0gwCbkV2EWao8T4RKaPn063+G6k56QzecFkMnPKf1Gt3buWC+ZdUC0RGdt9LHPPn6tERKQZMBERxI0c4X2d85Wm+G3Ool3hzJw8lIlDu+MKr3mMhyvcMHFo91pbJowxHNM5gdtO683c205m+UNjeeaS/ow7rjPRrsp/JpWUWr7evI8n/5PaoEQE3Asu7skuaNAxIqFALSMBopaR0FFTC8n/0v/HtOXTKKXUu1+CK4EHhz3IeUef16CBkiLiXwfm/ovMhx4CwJWURK9FC/V/NAjsyS7g/W/T+GbzPp8NIj9cWMKSjXtZmLqLRT/vYm9OYZNivG/8sVpwUVqkprSMKBkJECUjoaViQmIwWCr/PxuZNJIpI6fQJa6LQxGKSG2Kdu9m4+hTva+P/ug/RPXq5WBE0hyUllp+SD/A56m7mLl0K3lFJQ0+x6hjOvDW9cP8EJ2Is9RNS6SZSYxP5Pr+7ll4KiYiUeFRPDL8Ef4+7u9KRESaKVenTkQdd5z3tbpqCUBYmGFQ97bcf1Zfjk1MaNQ5UjMO8c9V6aRl5WlRTREPzaYl4mMH8g8wbcU0Pt7ycaVyg+Ef4/7B4C6DHYpMROorftQoCtatAyB38Ve0v26SswFJs9LYBRf35RZy7wc/ApDYOpohR7VjSM92DOvZjt4d4wkLa1x3QK1zIsFM3bQCRN20QsOi7Yt48usn2Ze/r1J5mAmj1JbWOsuWiDQved9+y7arrgbAuFz0Wf4NYbFagFTcXvliA899tt6n52wT6+KkHu0Y2rMtQ3u2p19SK1zhR+7AcqR1TsA9eF/rnEggNKWbllpGRHzgQP4Bpq6YyidbPqlUHhkWyd/G/Y3khGTvGJLJCyYrIRFp5mIGDiQsIYHS7GxsURG5y5eTMGaM02FJM9GoBReN4bwTElmz4yCb9+ZW234gr4iF63axcN0uAGJc4Qzq0YahR7VnSM+2nJjclpjI8oRC65xIS6FkRKSJFm1bxBPfPEFWfuVfCK2jWvPuue/SLaEbANPHT1dCIhIkTEQEcSNGkP3ZZwDkLl6sZES8GrXg4pBk7zoju7Pz+XbrflZsyWLFlizW7TxUbW2Tw0UlLN24j6Ub3S3trnBD/66tGdKzHUOPasfHazLrtb4KwPItWUyZn6p1TqRZUjetAFE3rZZnf/5+pi2fxidbP6m2rUtsF948+81qyUZ9FkYUkebhwJw5ZP7hEQBc3brR6/PPNMWvePliwcUyh/KLWLXNnZys3JLF6vSDFJaU1rhvY7nCDcseGKsxJOIXmk1LJMAWblvIhfMurJSIhBv3L5lu8d1qTESg+sKIi7YvCljMItIwcaNGeX8uSk+ncOtW54KRZscXCy6WaRXtYsyxnbj/rL7MuXUkqx8/k3dvGs69Z/Rh1DEdiI1seveqohLL+9/WvyVHJFDUTUukAfbn72fq8ql8uvXTSuUDOw7kyZOfZMmOJYztPvaIrR1lCcmi7Yu4KuUqf4csIo3k6tyZqD59KFjvHqicu3gxUT17OhyVNCfRrnCmXdyfe87o49MFF6Nd4Qw/uj3Dj24PQHFJKamZh7zdur74eTfFpQ3v2fLJmkx+PSSZ9vG+bR3RbF7SFOqmFSDqphX8Pt/2OU9981SlsSFR4VHcceIdXHnclYSHaWCgSEuz+7nn2Pfa6wDEnXIK3V971eGIROCivyzl+7QDjT4+sXU0/ZJac3zXVt7nLq2iG9wNUbN5SRnNpiXiR1n5WUxdPpUFWxdUKj+x04k8efKT9GjVw6HIRMTf4kaN9iYjeStXUpqfT1h0tMNRSaiLj27an2+ZB/PJPJjvnbkLoF1cJP2SWlVKUnq0i6117RPN5iW+omRE5Ag+2/oZTy9/ulJrSHR4NHcMuoMr+l6h1hCRFi72xIGExcVRmpuLLSggb8UK4kePdjosCXHDerZj8Ya9DT4uxhXO4aKSGrdl5RayeMPeSueNj4ogJbEV/Sq0oPTuGE9EeBhT5q/VbF7iE0pGRGqQlZ/F0988zWfbPqtUPqjTIJ44+Qm1hoiECBMZSeyI4eQsdE82kfPVYiUj4rjGrHPiCjf8777TKLGWtTsO8VPGQdZmHGLtjoNkHMyv8ZicgmJWbM1ixdbypCMyIozeHeP4eWd2g2KesyqNe87oozEkUo2SEZEqFmxdwNPfPM3+gv3esujwaO4cdCdXHHcFYUaT0ImEkvhRo8uTkcVfAQ87G5CEvMasc3Lp4GQ6tXJ3MUxsHcO4lM7ebVm5haz1JCc/7ThIasahGhdmBCgsLiU1s2GJCJTP5nX7mN4NPvZINHg++CkZEfHYd3gfU5dPrbE15MmTn6R7q+4ORSYiToofXWGK323bKdy2jcgeah0VZz02oR+b9+TWe52Txyak1Lq9XVwko47pyKhjOnrLcgqKWZfpTk7KkpSNu3MaNYtXmfk/ZnB630707BDX5PEjRxo8v3jDXl5cuF6D54OEkhER4NOtnzL1m6nVWkPuGnwXE/tOVGuISAhzJSYSdUxvCjZsBCBn8RLaKRkRh5WtczJlfipzVqX5fDar+KgIhhzVjiFHtfOW5ReVsH5XNre/8x1pWYcbHPPPO7M5+6XFhBlIbhdL747x9OoUX/7cKZ7WMa46z6PB8y2LkhFpUWalzqpznQ9wr4S+aPsizu55Nk8vf5rPt31eafvgzoN5YuQTag0REQDiThlVIRn5inZXXelwRCL+W+fkSPWd0K0NR7WPa1QyUqbUwrZ9eWzbl8ein3dX2tYxIYpeHePo7UlSendKoHeneDq3ivJOPazB8y2L1hkJEK0z4n+zUmfxzMpn6Bbfjenjp9eakGTmZDJ5wWTSc9KJCY/hcEn5DTUmIoY7B92p1hARqST366/Zft1kAEx0NH2Wf0NYlPqjS2h65YsNPPfZ+gYf1z4+kkOHixo08L5MfFQEvTrG0bVtDJ/+tJOG9BZzhRuWPTDWL2NINGbFrSnrjCgZCRAlI/5XMcmoLSHJzMlk0qeTyMjNqHb84M6DeXLkkyS3Sg5UyCISJEoLC1k/fAQ2Lw+A5NdeI/6Ukx2OSsQZu7PzOfmPXzR4Nq9lD4ylbayL7Vl5bNydw8Y9OWzcncOmPbls2p1DTkGx32K+b/yxPh08rwUfK9OihyJAYnwi08dP9yYkkxdMrpSQZGRnMPHjiZXWDAF3a8hdg+7i131/rdYQEalRWGQkccOGkfPllwDkLv5KyYiErMbO5lXWUnB0x3iO7hjPmRW2W2vZdajAnaTszmbTnlxvwrInu6DJMb/9zTaiXeF0bxfrfcRENi5J0JgV31LLSICoZSRwamohycrP4roF13G4uHIf15M6n8QTI59Qa4iI1Gn/7NnsnPIEAJE9e9Lrk48djkjEOfX5g7zMsJ7tmvQH+cG8IjbuyWHT7hye++wXdvsgOQH3+JSKyUn3drF0b+9+7hgfVevq8w/OXd2gRGzi0O4tfsyKumkFASUjgVUxIQk34ZTYyivOxkTEcPfgu7n82MvVGiIi9VKYvoNN48Z5X/da+DmR3bo5GJGIs9xdlfwzm1dtrn59eaNWn2+oqIgwb4KS7Hnu0T6W+KgIrnxteYOmOA6FMStKRoKAkpHAW7N3DVd+dCWWyp/xIV2GMGXkFJIT1BoiIg2z6ZxzKdy8GYAujz1K24kTHY5IxHl7sgsCMpsXNH7w/MDkNrSKcZGWlUf6/rxGDaJvipY+ZkVjRkQqsNYyb9M8/rjij9USkdsH3s5NJ9yk1hARaZT4UaeQ5UlGcr5arGREBHd3p9vH9Pb56uo1uWxIMi8t2tDgwfOvXnOSNzEqKbXsPJTP9n15bM/KZXtWHtuzDruf9+WyP6/I53H/c1U6R3eII7FNDEmto+lwhG5gdWlpY1aUjEiLkpGTwRNfP8HSjKWVyg0Gi2Xexnlc0OuCOtchERGpSdyo0WTNfBOA3OXLKS0sJCwy0uGoREJHUwfPA4SHGbq2iaFrmxhG9Gpfbf9D+UWkZeWRlpXH9iz3eijbPa+37strVNyb9+Zy69vfeV+7wg2dW0WT1DqGxDbRJHmSlMSy161jaBPr8q6tUlFLW2dFyYi0CKW2lA9++YDnVz1PXnH5jcJguGPQHZxz1Dlc/9n1Nc6yJSJSX7FDTsJER2Pz87F5eRxetYq4ESOcDkskpDw2oR+b9+TWe/D8YxNSGnT+VtEu+iW1pl9S62rbrnztG5Zu3Neg89WkqMSSvv8w6ftrXzwy2hXmTVYSW7uTlfjoCN7/Nr1Bdc1ZlcY9Z/RptuueqK+KBL3th7Zz/YLreWr5U5USkejwaGadM4sb+t9AUkIS08dPp1t8N29CkpmT6WDUIhKMwqKiiB021Ps656vFDkYjEpqiXeHMnDyUiUO74wqvuauTK9wwcWh3n3dRGnF09ZaU+ujZIZZjOyeQEF3/doD8olI2781l6cZ9zFmVzstfbGTqxz9T0pAVH3EnPu9/W/+WpEDTAPYA0QB23yspLWHWulm88v0r5JfkV9rWLrods8+ZTVJCUqXy+iyMKCJyJFmz3mbXU08BENm7F73+8x+HIxIJXYEcPA9NW/CxLJ6cgmIyDxwm42C++7ns54OHyTyQT8bBw+QXlfo07lHHdOCt64f59JwVaQC7hJxNBzbx6NJHWb13dbVtSXFJzDhrRo1JRtWFERdtX8RVKVcFImQRaSHiR49il+fnwo2bKMrIwJWUdMRjRMQ/Ajl4HnwzZiU+KoJjOidwTOeEGve31nIgr4gMT3KSebA8cfl83S5yC0pqPO5I/Lm6fVMpGZGgUlRaxBs/vcHff/w7RaXls120i27Hw8MeZnfebsZ2H3vE1o6yhESJiIg0RmT37kT26EHhtm0A5CxeQtvLL3M4KhEJFH+PWTHG0DYukrZxkdXGrTR2nZX4qOb7J7/GjEjQ+DnrZ6746Ar+/P2fKyUi5x19HvMumMeZR53JVSlX1avbVWJ8ohIREWm0uFGjvD/nLP7KwUhEJNCcHLMyrGe7Rh03vJFjXQKh+aZJIh6FJYX8/ce/88ZPb1Bsy5sZO8V04tERj3Jq8qkORicioSh+9Cj2z5oFQN6yr7GFhRhN8SsSMqJd4Uy7uD/3nNEnoGNWGrvOymUnNd+FnpWMSLO2es9qHl36KJsObqpUfskxl3DPSffQKrKVQ5GJSCiLHToUExWFLSigNC+PvO9/IK7CLFsiEhqCccxKc6NuWtIsHS4+zHMrn+PqT66ulIh0je/K/53xfzw+8nElIiLimLDoaGKHDPG+zlVXLREJkMcm9Kt3d63GjFkJNCUj0uys3LmSSz+8lJmpMym17qntDIYr+l7B3PPnMiJJC4yJiPPiR1cYN6L1RkQkQJwcs+IP6qYlzUZuUS4vrHqB9355r1L5Ua2OYsrIKQzqPMihyEREqnMPYp8GQMH69RTt2oWrc2dngxKRkODUmBV/UDIizcKyHct4/OvHycwtXxU9zIRxbb9ruW3AbURHRDsYnYhIdZFHHYUrOZmiNHff7dzFi2lz6aUORyUioSTQY1b8Qd20xKdmpc4iMyezzv0yczKZlTqLgwUHeWTpI9y88OZKiUjvNr15+5y3uWfwPUpERKRZMsYQP+oU72t11RIRaTi1jIjPzEqdxTMrn+HtdW8zffz0Wtf7yMzJ9K6A/soPr5BblOvdFmEiuPGEG7mx/424wl2BCl1EpFHiRo1i/zuzAchdtgxbVIRx6d4lIlJfahkRnxnbfSzd4ruRnpPO5AWTa2whyczJZNKnk0jPSQeolIgc1+443j3vXW4beJsSEREJCnHDhnmTj9KcHA7/+KPDEYmIBBclI+IzifGJTB8/vdaEJCM7g19/9GsycjMqHRcZFsmdg+7knXPf4dh2xwY6bBGRRguLja00xa+6aomINIySEfGp2hKSn/b+xAXzLiArP6vS/gM7DuSD8z/ghv43EBGmXoMiEnzcs2q55SxWMiIi0hBKRsTnqiYkZ889m4kfTSS/JN+7T0xEDA8MfYAZZ83g6NZHOxitiEjTVFxvpGDdOop273YwGhGR4KJkRPwiMT6RqaOmYjCU2JJK24Z2Gco/z/8nVx53JeFhzXshHhGRukQefTSupCTv69wlSx2MRkQkuCgZEZ8rtaW8s+4dbvzsRiy20rY7B93Ja2e+RnJCskPRiYj4ljGmSletrxyMRkQkuCgZEZ/acnALkz6dxLQV0ygoKfCWGwwA/1z/T3bm7nQqPBERv6jYVSt36TJscbGD0YiIBA8lI+ITRaVFvLbmNS798FK+3/29tzyMMO4fcj+fXvxpndP+iogEm6w336QoI4PYYcOhbIrfQ4c4vHpNpf2KMjLIevNNJ0IUEWnWlIxIk63bt44rP7qSl757icLSQm95TEQMs8+dzVUpV5GUkHTEaX9FRIJN1ptvsmvqNLZdO4nSQweJHTzYu61iV62ijAy2XTuJXVOnKSEREalCyYg0WkFJAS9/9zITP5rIuqx1lba1j27PvAvmkdIhxVtW1zokIiLBJGHcOFzJyRSlpbHt2knEDDjBuy3Xs95IWSJSlJaGKzmZhHHjnApXRKRZUjIijfL97u+59MNLeXXNq9Vmy0qKT2L2ubNJjE+sdlzVhGTR9kWBCllExKdcSUn0mDnDm5AcnPehd1v+2rUc/mltpUSkx8wZlWbdEhER0Cpz0iB5RXm89N1LzP55dqWZsrrGd+XxkY+zcf9GxnYfW2MiUqYsIVm0fRFXpVwViLBFRPyiLCEpSzoID4cS9xc0aTfdRElWlhIREZEjUDIi9bZsxzKmfD2FjNwMb5nBcOVxV/LbE39LrCuW4YnD63WuxPhEJSIi0iJUS0g8lIiIiNRNyYjU6WDBQZ5d+SzzNs2rVN6rdS+mnDyFAR0HOBSZiEjzUJaQbDzjTG/LCIAtLGTf9DdIOH0MsSedhImMdDBKEZHmR8mIHNHCbQt5evnT7D2811sWYSK4vv/13HTCTUSG6xeriIiXMZVeFu/axf5Zs9g/axZh8fHEnXIKCWNOI270aCLatnUmRhGRZkTJiNRo7+G9TF0+lc+3fV6pPKV9Ck+MfIJj2x3rUGQiIs1P2axZFBe7x42UloK1lfYpzckh+9NPyf70UwgLI2bQiSSMGUP8mDFE9uyJqZLIiIiEAiUjUom1lvmb5/PMimc4VHjIWx4VHsXtA2/n6pSriQjTx0ZEpEzV6Xt7zJwBwNarrqY4IwMTHY21FgoKyg8qLeXwt6s4/O0qdj/7HK4e3UkYczrxY8YQO+hEjGcBRRGRlk5T+4pXZk4mty66lYeXPFwpERnUaRBzJszhuuOvUyIiIlJBTYmIKykJV1ISR816C1dyMjY/H1fHjnR56ina/OpXRHTsWP0827aTNWMG26+9lvUnn8KOe3/HwY8+ouTQoRpqLVe2Anx94tSCiyLSHAV1MmKMiTbGTDHGrDfG5BtjMowx040x3RpxrjbGmBeNMduMMQWe55eMMW38ELpfzUqdVa/FBDNzMpmVOotSW8rsn2dz4bwLWbpjqXd7bEQsfxj2B9446w2Oan2UHyMWEQk+tSUiZSqtQ5Kezr5//IMOt95C7//9l6M++IAOt91K1HHHVTtv6aFDHProIzLu/R3rR57MtmsnsW/GDAq3bau0X8UV4I+UkPhzBXinkyEn6w/l9+50/aH83ptD/T5nrQ3KBxANLAUskAG8Byz3vN4N9GrAudoD6z3HbvKc6yfP6w1Aex/EuzYlJcX621tr37LHzzjenjXnLJuRneEuLMix9rFW7kdBjrXW2ozsDHvWnLO8+x4/4/hKj1s+v6X8+Kaooe6AUv3O1R/K793p+kP5vQew/n0zZ9rUY/vaDePOsIU7dtRad+GOHXbDuDNs6rF97b6ZM6udp3DHDrvv7bftthtutOuO729Tj+1b62PjOefaXc8+a3O//dYWbE/zntcbQ5X6K9ZdKc5m9P6Dsf5Qfu9O1x/K77051F+blJQUC6y1jfgbOZhbRh4CRgJfA32stZdba4cB9wIdgekNONcLwDHAXOBYz7mOB/4M9Aae92nkfjS2+1jv6uaTF0yusYUkMyeTyQsmk56TDuB9Bmgd1Zqpp0zlr2P/esSFC0VEQl27a66h80MP1rmOSFkLSeeHHqTdNdfUuL3dFVfQ/dX/o883X9P1zy/T+qKLCG/Xrtq+hZs2se+119l25VVs/dWviO7bl/D27SlKS3O3kGTu9O5blLnTryvAJ4wb5119vrbWmaqtRwnjxrWI+kP5vTtdfyi/9+ZQvz8EZTJijHEBv/W8vN1am1O2zVr7PLAaGG2MGVyPc3UBrgSKgNustcUVNt8H7AGuNMZ09lX8/lS2unmlhCS3/JdTZu5Orv7k6koJSJnxR41n3gXzmNBrgmZ1ERGph3bXXFOvP/BdSUk1JiJVhcXF0eqMM0iaNpVjFn9Fj9nv0P7GG4k6pne1fUsOHCD7888p2bcPgKK0NLZcfiW5e1wU5Yaz7fqb/ZaIlL0nbzc0B5IhJ+sP5ffudP2h/N6bQ/3+YGyVqQeDgTFmDPAFsMlaW+0ObYx5BHgCmGKtfbyOc12HuxVlkbW2WupojHkdmAxcZ62d0YSY16akpKSsXbu2sadokIqtH93iuzI9dQWFwMRefckuyq60b4eYDvxh+B8Y232s7wMpzIWpnv8ED2VAZJzv61D9zbP+UH7vTtcfyu/d6fr9WHdhWho5X/6X7C+/IG/lt+5phGtlAYOJiiJm0CAi2rUjLDaWsNgYTGwsYTGxntfusrDYWMJiPNsqPmJiMBFHnrik0rew3brS48TvANj2/SCK0nf4/Q8iJ+sP5ffudP2h/N6bQ/1V9evXj9TU1FRrbb+GHhusUyOVLfn9XS3bv6uyX1PPNbme52oeCnNJjGzF9NNfYfIXt5Oes4Ozk5MoAaiSiFzc4yzuOfEOWkcmwOEDPqg7r/LrogqvD2WCK7by9sgqr1V/8NYfyu/d6fpD+b07XX8A647skEC7X02g3a8mUJKdQ+7X35D95f/IWbac0kPZVfZ2t27bggLyvv660XUCmMhIdxITHeNOXGKivYlKWJSLsNgYYgf0I/vAforSd7BxRyf3gXYH4W1a0/aic8j932fuFvfIaDBhmDADJgzCwjBhYRBm3D+bij+7nzGmfB8ThgkPKz+2pACMoct9vyHzqeeq1R/RuROJD92N3Z9B4f4MiIyp8MZMpWdTS3n5c+VrS3E+AElPPEjGg1Oq1e3q0pmuTz2EKdhP8Zb9leumyrnq9Q9RZd+iwxig69Q/sOP3j1avP7EzXaf9AVN8kOLtB6t/FpuqKM9d/7Q/sOO+muvvNu0RwooPUbL9UC3vv5EKDxMGdJv2COn3PRLYuoOh/m5d6fHa33B1aO3+giTQXwg1ULC2jDwP3A28YK29p4btA4AfgO+stUfsqmWMmQtcBNxprX25hu0XAP8G5lprL6lHbLU1ffRKSUmJCkjLyOOtvT9ujojgom6JlFa5iXUtKuaxvfsYkV9Q9WgREQkithQO743kwLYYDm6KpUF/4IpIy2Isvc/bjSuupLzs8YN+rzYUW0biPc95tWzPrbJfoM7V7OwND6e0YoG1XJSdwwNZB4gNwkRUREQqM2Hgiishb2cUYMB47u3WEB5VQtveuZgwKC02nkeY92frLatcjlVCIyKBEazJSNldsra/phtyF/XluagtI/S0mKQ05FyN9pB7ZoXM3J089sXtkLOjLDgwhpUx0RwMCyO2pOQIJxERkWBQlBvOti/aU5QbgSuumB6nuwe0l5Ud3BpLj9P3Vf6mtA62pEryUmIoLTKUlngSmKLKZcWHw8hOi6G0OIzyX6cGE2aJaltIWJj7VxDWuLda98P9s/Ee4t3HUq2s/OeyY43351ILlBjcv7LL6/f+XLHYWyDSAlnDtu8G0mP6q7gSuzgdTb0EazJS1jm2tk5wZR0jc2rZ7q9zNQ+Rce4B7F/8hvScHXSLS2JCeirnZudyS9+TSM/LZPJxQ5h+2sskxvl5krDCPHjBs6jX3et831dc9Tff+kP5vTtdfyi/d6frD3DdRZk72XbTbyjK9fQT/8uzuGafBkCPd95n2+33UZS+g23fD6LH/71S7z9ODBDuedQ3htLiHbi6JtJj0A8AbPtuIEU7MimJOZpuDai7obzXIL3m+l3dutb53r1d1mt99u5Zaf+inbtIu+0uinZk4ErqQvLAHwFI+2EARRk7cSUlkfzK87i61PS7tgG9E2rpyVC0cxdpv72XooxMd/0DfsQC6T8MoChzJ66kRLq9/Fwt9Tdd0c5dpN/xO3f9iV3o5nn/1erv7Pv6i3Y5V3ezrv9H9+d+2w23BsVMWhC8ych2z3NtK613q7JfoM7VLFSeSasb009/hcSXBgEwfdzfPElKOpP/dxfTx0/373oi4a7yn2NaB34Qlep3rv5Qfu9O1x/K793p+gNYd1FGBttuvqPyzDkdyscMunr2ocebb3pn3Nl28x0+/+OkWgyv/Q3XG+7fNz2mv8q2G271W92+rL+h7STGU3fab+51JyJV637jNW/daXfe57f3nnbX/e4/RqvWP6O8/vS7H/Bb/el3P+BI/U7W3ezrf6PC5/7aSUGRkATlOiPAj57nQbVsLytfHeBzOa5aIjJ+Oolx5d8GJcZ1qb4OSQ0LI4qISPOWvXBhnWsJVF2TIHvhQp/VX3VhtR4zZ1RqfXAldqm+HkINC7QFY/2h/N6drj+U33tzqN8fgjUZWQocBHoZY06sYfulnuf/1ONcnwKlwChjTKeKG4wxUcAEz/ZPGh9u4CzavqhyIlJDq0fVhREXbV/kQKQiItIUvloBvrGcToacrD+U37vT9Yfye28O9fuFtTYoH8BTuDtcLgXiKpTf4ylfXGX/3wA/A9NqONcszzFzgIgK5S95yt/yQbxrU1JSbCC8tfYtm5GdUed+GdkZ9q21bwUgIhERaYn2zZxpC3fsqHO/wh077L6ZM1tU/aH83p2uP5Tfe3OovyYpKSkWWGsb8TdyUK4zAmCMiQb+CwwDMoHFQA/P633AcGvtxgr7Pw48Bsy01k6qcq4OwDdAL2AT8C3QDzje83q4tXZvE+MN6ArsIiIiIiKB0JR1RoK1mxbW2nxgDPAk7jVCLgSOAmYCJ1ZMROpxrr3AEODPQCTuRRBbA68AQ5uaiIiIiIiISHVB2zISbNQyIiIiIiItUUi2jIiIiIiISHBTMiIiIiIiIo5QMiIiIiIiIo5QMiIiIiIiIo5QMiIiIiIiIo5QMiIiIiIiIo5QMiIiIiIiIo5QMiIiIiIiIo5QMiIiIiIiIo5QMiIiIiIiIo4w1lqnYwgJxphDUVFRCb169XI6FBERERERn9m0aRMFBQXZ1tpWDT1WyUiAGGN2ArFAmkMhlGVBmxyqP5jp2jWerl3T6Po1nq5d4+naNZ6uXdPo+jWe09cuGciz1nZp6IFKRkKEMWYtgLW2n9OxBBtdu8bTtWsaXb/G07VrPF27xtO1axpdv8YL5munMSMiIiIiIuIIJSMiIiIiIuIIJSMiIiIiIuIIJSMiIiIiIuIIJSMiIiIiIuIIzaYlIiIiIiKOUMuIiIiIiIg4QsmIiIiIiIg4QsmIiIiIiIg4QsmIiIiIiIg4QsmIiIiIiIg4QsmIiIiIiIg4QsmIiIiIiIg4QsmIiIiIiIg4QslIM2aMiTbGTDHGrDfG5BtjMowx040x3RpxrjbGmBeNMduMMQWe55eMMW2OcEyYMeYuY8waY8xhY8weY8wHxpiUJr2xAPDFtfNcsyuMMe8YY1KNMbnGmGxjzHJjzJ3GGFctx80wxtgjPG7x3Tv1PV997owxW+u4Dn1rOS5oP3fgs8/epDquXdnjmirHBe1nzxgz2BjzgDFmrjFmhyfe/CacL2Tueb66dqF4z/Pl5y4U73k+/OyF1D3PGBNrjLnQGPO6MWa1MeaQ5//bj8aYR40x8Y04Z9De87QCezNljIkGFgEjgUxgMXAUMBTYA4yw1m6q57naA18DxwCbgW+Bfp7HRmC4tXZflWMM8D5wKXDAE0sHYDSQD4yx1i5vynv0F19dO2PMU8DDQCnwPe5r1RE4GYgClgDjrbV5VY6bAVwLLAB21nDqmdbaLxvx1vzOx5+7rUAPYGYtuzxorc2sckzQfu7Ap5+9U4AbatncGrjQ83Mva+3mCsfNIHg/e/8GLqhSXGCtjW7EuULtnvdvfHDtQvSe929897nbSujd8/6Nbz57IXXPM8bcALzqebkWSAVa4f7dkQD8DJxqrd1dz/MF9z3PWqtHM3wATwAWWAbEVyi/x1P+vwac603PMf8EIiqUv+wpn1nDMZM929YDnSuUX+Ip31jxXM3p4atrBzwAPA10rVJ+DLDNc66pNRw3w7PtNKevhcOfu63uW0yD6g/az52vr98R6rjVc64lNWwL5s/e/cAU4Dygs+d95DfyXKF2z/PJtQvRe54vP3eheM/z2fU7Qh0t7p4HXAP8FTimSnki8J3nPb3TgPMF9T3P8X8QPWr4RwEXsN/zYTixhu0/erYNrse5ugAlQGHFD5tnWxSwGyiuYdtaTx0X1nDOeZ5tlzh9rfx57eqoZ6LnPFtq2BasN0efXrtG/mIOys+dP67fEepZ6jnPzTVsC8rPXi3vs7F/UIfUPc+X166Oc7a4e56vr12o3fN8ff2OcM6QuOdVeE8jyq4jEFmP/YP+nqcxI83TKUAbYJO19vsats/xPE+ox7nOxj026Ctr7a6KG6y1BcB8INyzHwDGmJ5ACnAY+KiJ9QeaL6/dkfzoeU5q4nmak0BduxoF+ecOAnD9PNdoJO5fOu839jwtXKjd8wKlJd7zHKXPXd1C9J5X9n8tCmhfj/2D/p4XEYhKpMEGeJ6/q2X7d1X2a+q5Jlc5V9nPP1lri5pYf6D58todydGe55r6qJa52BhzCe6bwBZgvrX25ybW609+uXbGmPuAXkAB7m9i/mWt3XOE+oPxcweB+exd5Xn+yFq7/wj7Bdtnz5dC7Z4XKC3xnucXIXTPC4RQvOeV/V8rArLqsX/Q3/OUjDRP3T3P6bVsT6+yn6/P5cv6Ay1Qsd/peZ53hH1+W+X1M8aYvwF3WmuLm1i/P/jr2v2/Kq9fMMbcYa19PUD1B0og4r/S8/xWHfsF22fPl0LtnhcoLfGe5y+hcs8LhFC855X9X/vU07JRl6C/56mbVvNUNqVbXi3bc6vs5+tz+bL+QPN77J6pAsfhnn3ijzXs8j1wC9AHiMX9Lcftnv1vA55tbN1+5utr9yFwMe7ZZWKB44HncTc9v2aMudDP9QeaX+M3xgwFjsU9LqWmZnUI3s+eL4XaPc/vWvA9z9dC7Z7nV6F4zzPGnANcj7tV5JF6Hhb09zy1jDRPxvNs69jur3PVdUxz5strV/1gY04FXvKcf7K1NqPqPtbal6oUbQH+aoz5ClgF/NYY87y1Nq0psfiBT6+dtfaOKkVrgXuNMb8A/wCeAf7dgPqbO79+9ijvrvCetbawph2C+LPnS6F2z/OrFn7P86kQvOf5W0jd84wxxwGzcH8u7rPW/ljHId5DPc9Be89Ty0jzlO15jqtle6znOcdP56rrmLLy+tQfaL68dpUYY07A/YskEnfT778acry19ifc35yF4/6Wsbnx27Wr4jXcs3v08Qyiq2/9zflzB/797EUAl3te1tVdoZog+Oz5Uqjd8/wmBO55gdJS73l+E2r3PONeFPdToC3wfA1J1pEE/T1PyUjztN3zXNuKzd2q7Ofrc/my/kDzS+zGmF64F1VqAzxurf1zo6KDDZ7nxEYe708B+Xe31pYCZQv/VbwOwfy5A//GfybQCdhsrV3WiOOheX/2fCnU7nl+ESL3vIBowfc8fwqZe54xpgPwOe7xGW8Av2vgKYL+nqdkpHkqa5obVMv2svLVfjpX2THHG2NcTaw/0Hx57QAwxiThvlF0AV6y1k5pfHi09Tw3x2+6fH7tjqCm6xDMnzvw7/Ur664wqxHHlmnOnz1fCrV7ns+F0D0vkFriPc+fQuKeZ4xJAD4B+gJzgRutZ6GPBgj+e14gFjPRo2EP3E3iB6h78bQh9ThXIu7FcAqATlW2lS2GUwJ0qbItlboXw/mV09fKn9fOs39bYI3nmOmAaUJsUbi/ZbDAyU5fK39fuyPU0w8oxT1ALrLKtqD83Pnz+uEeQJjrOfaYRsbWrD97tcTcqMXTQu2e58tr5zk2ZO55vr52Rzhni7zn+ev6hco9zxPjF54YP6362WjAeYL+nuf4P4YetfzDwFOeD8JSIK5C+T2e8sVV9v8N8DMwrYZzzfIcMweIqFBeNijxrRqOucGzbX3FDzfumUIssBlwOX2d/HntcPez/NpzzHtAeD3qPha4oOq+QEfgX55z/dCUX/BBcu3GU8NK48AJFW6AL7Wkz50vr1+Vfa7xHPt1S/7s1fB+jvhHje55vr92oXjP8+G1C8l7nq+uX5V9Wvw9D/dYlrmeGL8CYutxTIu952k2rebrKdyDrkYCG4wxi3FPFzgM2AdcV2X/Drj/Y9bUN/IuYDhwCfCzMeZb3N/UHI+7H+vdNRwzHTgHuMhzzCJPHacC+cBVtuaFcpoDX127p3FftxKgGHjdmOqTUlhrJ1V4mYh7wOc+Y8zPwA7c/V4HAwm45+6+zHr+xzdDvrp2I4DHjDHbcH/G9gA9cTf9RgD/Ax6sof5g/tyBb//flinrrlDXIM6g/uwZY86l+lSWkcaYbyq8ftJaWzbFp+55Hj68diF3z/PhtQvJe56P/9+WCYV73m9w/5sD7MU9A1hN+/3OWrvX83PLvec5nR3qUfsDiAGeADbibn7bCcwAkmvY93HcmeyMWs7VFngZd7Nlgef5z0C7I9Qfjvsb3Z+Aw7j/w/wT6Of0tQnEtfPsb+t6VDkmCXgB97eLmUAh7lkrVnnqaev0tQnQtRsBvI67v+le3HOm7wO+xP1tTK3fuAbz585X16/C9kTcfxQWAu3rqDeoP3vApHr8f5vUgGsXMvc8X107QvCe58NrF5L3PD/8vw2Je16F61DX46gGXLugvecZTzAiIiIiIiIBpdm0RERERETEEUpGRERERETEEUpGRERERETEEUpGRERERETEEUpGRERERETEEUpGRERERETEEUpGRERERETEEUpGRERERETEEUpGRERERETEEUpGRERERETEEUpGRERERETEEUpGRERERETEEUpGRERERETEEUpGREQkZBhjNhhjNlUpe80YY40x9zoVl4hIqFIyIiIiIcEY0xfoDXxYZdNAz/P3AQ1IRESUjIiISMg43/M8v6zAGBMB9PO8/CHQAYmIhDolIyIiEiomAAeBxRXK+gLRwHZrbZYjUYmIhDAlIyIiEnSMMRcZY/5rjMk2xuw1xrxljOlojPmjZ/zHJVX27wCMAD6x1hZV2DTQ8/xDYCIXEZGKIpwOQEREpCGMMc8DdwO5wH8BC1wADAC2e3arOv7jHCCcCl20PAbUsr+IiASAkhEREQkaxpjJuBORlcCvrLXbPOVHA98Cx+PuirWlyqHnA8XAJ1XKB3qef/BPxCIiciTqpiUiIkHBGNMa+BOwHzi/LBEBsNZuBmYABvjBWmsrHBcJnAkstdbur3JatYyIiDhIyYiIiASLm4A2wJ+ttTtr2F62fkjVxGIMkECVKX2NMV2BjsD+iomNiIgEjpIREREJFhM8z+/Xsj3G81w1GSk7rup4kYGe5x+aFJWIiDSakhEREQkWA4Eia+3aWraneJ5/qFI+AfjFWruhhvPVtL+IiASIkhEREWn2jDEu3F2tcmvZHgmcBxQA6yqUDwC6U71VBDReRETEcUpGRESk2fOsDZIFtDHGtKthl3txj//4qco6ImWrrn9Y/RC1jIiIOE3JiIiIBItVnufHKhYaY64GnvS8rGm8SBawrMoxcUAvqrSkiIhIYCkZERGRYPEE7gUO7zDGrDbGvGuMWQ28DvzHs8+Ksp2NMYnAScDH1tqSKucagPt34BprbbH/QxcRkZooGRERkaBgrV0CXAKsBo7FPWXvWmAQUOjZbVGFQ87Dve5ITeNFhnuel9WwTUREAsRUWBdKREQk6Bhj+uIe9/G1tXZMhfL5wHigg7X2UJVjlgAnA+OstRUTGBERCSC1jIiISLNnjAnzzIxVtfwU4BMgHPh9lc2Lgd/WkIhcjzsRWQ186Z+IRUSkPtQyIiIizZ6n9WOd57EBKAb6AMfj7qJ1g7X2rSMc3xl4DuiLexxJNnCqtVbT+oqIOEgtIyIiEgzygTdxr7I+DjgXiAVeA048UiLiMRKYCCR5zjNYiYiIiPPUMiIiIiIiIo5Qy4iIiIiIiDhCyYiIiIiIiDhCyYiIiIiIiDhCyYiIiIiIiDhCyYiIiIiIiDhCyYiIiIiIiDhCyYiIiIiIiDhCyYiIiIiIiDhCyYiIiIiIiDhCyYiIiIiIiDhCyYiIiIiIiDhCyYiIiIiIiDhCyYiIiIiIiDhCyYiIiIiIiDhCyYiIiIiIiDji/wNhR4mxHx76lQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 900x600 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.plot(gs, entropies, marker='o', label='max entropy $S$')\n",
    "# plot expecation values of sigmax and sigmaz as well\n",
    "plt.plot(gs, vals_X, marker='+', label=r'$<\\sigma^x>$')\n",
    "plt.plot(gs, vals_Z, marker='x', label=r'$<\\sigma^z>$')\n",
    "plt.plot(gs, corrs_XX, marker='x', label=r'$<\\sigma^X_0 \\sigma^X_{20}>$')\n",
    "\n",
    "\n",
    "plt.xlabel('$g/J$')\n",
    "plt.ylabel('expectation values')\n",
    "plt.legend(loc='best')\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "1dfbc061",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "id": "ad7f61d0",
   "metadata": {},
   "source": [
    "### Advanced exercises - if you're an expert and have time left ;-)\n",
    "\n",
    "These examples only scratch on the surface of what you can do with TeNPy.\n",
    "- There are plenty of [more examples](https://tenpy.readthedocs.io/en/latest/examples.html) in the documentation.  Take a look at them!\n",
    "- Try to learn how to define your own model from the TeNPy documentation. Define a model for the XX Chain.\n",
    "- Look at the documentation how to run TEBD and TDVP and reproduce the time-evolution plot for S(t) from the toy    code notebook.\n",
    "- Learn how to save and load data in TeNPy."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "ea502580",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.1"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}