diff --git a/test.ipynb b/test.ipynb
index 1df7eec..8a76042 100644
--- a/test.ipynb
+++ b/test.ipynb
@@ -2,17 +2,9 @@
  "cells": [
   {
    "cell_type": "code",
-   "execution_count": 1,
+   "execution_count": null,
    "metadata": {},
-   "outputs": [
-    {
-     "name": "stderr",
-     "output_type": "stream",
-     "text": [
-      "[Mac:57824] shmem: mmap: an error occurred while determining whether or not /var/folders/yh/dx7xl94n3g52ts3td8qcxjcc0000gn/T//ompi.Mac.501/jf.0/1921449984/sm_segment.Mac.501.72870000.0 could be created.\n"
-     ]
-    }
-   ],
+   "outputs": [],
    "source": [
     "import os\n",
     "from sys import stdout\n",
@@ -39,20 +31,12 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 2,
+   "execution_count": null,
    "metadata": {},
-   "outputs": [
-    {
-     "name": "stdout",
-     "output_type": "stream",
-     "text": [
-      "Number of nodes in the parallel cluster:  1\n"
-     ]
-    }
-   ],
+   "outputs": [],
    "source": [
     "# this cell mimicks an input file\n",
-    "fdf = sisl.get_sile(\"./lat3_791/Fe3GeTe2.fdf\")\n",
+    "fdf = sisl.get_sile(\"./Jij_for_Marci_6p45ang/CrBr.fdf\")  # ./lat3_791/Fe3GeTe2.fdf\n",
     "# this information needs to be given at the input!!\n",
     "scf_xcf_orientation = np.array([0, 0, 1])  # z\n",
     "# list of reference directions for around which we calculate the derivatives\n",
@@ -65,7 +49,8 @@
     "    dict(o=np.array([0, 0, 1]), vw=[np.array([1, 0, 0]), np.array([0, 1, 0])]),\n",
     "]\n",
     "\n",
-    "# human readable definition of magnetic entities\n",
+    "\"\"\"\n",
+    "# human readable definition of magnetic entities ./lat3_791/Fe3GeTe2.fdf\n",
     "magnetic_entities = [\n",
     "    dict(atom=3, l=2),\n",
     "    dict(atom=4, l=2),\n",
@@ -74,8 +59,7 @@
     "        atom=[3, 4],\n",
     "    ),\n",
     "]\n",
-    "\n",
-    "# pair information\n",
+    "# pair information ./lat3_791/Fe3GeTe2.fdf\n",
     "pairs = [\n",
     "    dict(ai=0, aj=1, Ruc=np.array([0, 0, 0])),  # isotropic should be -82 meV\n",
     "    dict(\n",
@@ -84,6 +68,27 @@
     "    dict(ai=1, aj=2, Ruc=np.array([0, 0, 0])),\n",
     "    dict(ai=0, aj=2, Ruc=np.array([-1, 0, 0])),\n",
     "    dict(ai=1, aj=2, Ruc=np.array([-1, 0, 0])),\n",
+    "] \"\"\"\n",
+    "\n",
+    "# human readable definition of magnetic entities ./Jij_for_Marci_6p45ang/CrBr.fdf\n",
+    "magnetic_entities = [\n",
+    "    dict(atom=0, l=2),\n",
+    "    dict(atom=1, l=2),\n",
+    "    dict(atom=2, l=2),\n",
+    "]\n",
+    "# pair information ./Jij_for_Marci_6p45ang/CrBr.fdf\n",
+    "pairs = [\n",
+    "    dict(ai=0, aj=1, Ruc=np.array([0, 0, 0])),\n",
+    "    dict(ai=0, aj=2, Ruc=np.array([0, 0, 0])),\n",
+    "    dict(ai=1, aj=2, Ruc=np.array([0, 0, 0])),\n",
+    "    dict(ai=0, aj=1, Ruc=np.array([1, 0, 0])),\n",
+    "    dict(ai=0, aj=2, Ruc=np.array([1, 0, 0])),\n",
+    "    dict(ai=0, aj=1, Ruc=np.array([-1, 0, 0])),\n",
+    "    dict(ai=0, aj=2, Ruc=np.array([-1, 0, 0])),\n",
+    "    dict(ai=0, aj=1, Ruc=np.array([0, 1, 0])),\n",
+    "    dict(ai=0, aj=2, Ruc=np.array([0, 1, 0])),\n",
+    "    dict(ai=0, aj=1, Ruc=np.array([0, 1, 0])),\n",
+    "    dict(ai=0, aj=2, Ruc=np.array([0, 1, 0])),\n",
     "]\n",
     "\n",
     "\n",
@@ -91,7 +96,7 @@
     "kset = 20\n",
     "kdirs = \"xy\"\n",
     "ebot = -30\n",
-    "eset = 50\n",
+    "eset = 100\n",
     "esetp = 10000\n",
     "\n",
     "\n",
@@ -131,7 +136,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 3,
+   "execution_count": null,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -224,7 +229,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 4,
+   "execution_count": null,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -286,17 +291,9 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 5,
+   "execution_count": null,
    "metadata": {},
-   "outputs": [
-    {
-     "name": "stdout",
-     "output_type": "stream",
-     "text": [
-      "k loop:   0%|          | 0/400 [00:00<?, ?it/s]"
-     ]
-    }
-   ],
+   "outputs": [],
    "source": [
     "kset = make_kset(dirs=kdirs, NUMK=kset)  # generate k space sampling\n",
     "wkset = np.ones(len(kset)) / len(kset)  # generate weights for k points\n",
@@ -308,7 +305,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 6,
+   "execution_count": null,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -355,20 +352,9 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 7,
+   "execution_count": null,
    "metadata": {},
-   "outputs": [
-    {
-     "name": "stdout",
-     "output_type": "stream",
-     "text": [
-      "Number of magnetic entities being calculated:  4\n",
-      "We have to calculate the Greens function for three reference direction and we are going to calculate 15 energy integrals per site.\n",
-      "The shape of the Hamiltonian and the Greens function is 84x84.\n",
-      "k loop: 100%|██████████| 400/400 [01:51<00:00,  3.59it/s]\n"
-     ]
-    }
-   ],
+   "outputs": [],
    "source": [
     "if rank == root_node:\n",
     "    print(\"Number of magnetic entities being calculated: \", len(magnetic_entities))\n",
@@ -396,6 +382,11 @@
     "        HK, SK = hsk(H, ss, dh.sc_off, k)\n",
     "        Gk = inv(SK * eran.reshape(eset, 1, 1) - HK)\n",
     "\n",
+    "        # solve Greens function sequentially for the energies, because of memory bound\n",
+    "        # Gk = np.zeros(shape=(eset, HK.shape[0], HK.shape[1]), dtype=\"complex128\")\n",
+    "        # for j in range(eset):\n",
+    "        #    Gk[j] = inv(SK * eran[j] - HK)\n",
+    "\n",
     "        # store the Greens function slice of the magnetic entities (for now) based on the on-site projections\n",
     "        for mag_ent in magnetic_entities:\n",
     "            mag_ent[\"Gii_tmp\"][i] += (\n",
@@ -429,180 +420,9 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 8,
+   "execution_count": null,
    "metadata": {},
-   "outputs": [
-    {
-     "name": "stdout",
-     "output_type": "stream",
-     "text": [
-      "##################################################################### GROGU OUTPUT #############################################################################\n",
-      "================================================================================================================================================================\n",
-      "Input file: \n",
-      "Not yet specified.\n",
-      "Number of nodes in the parallel cluster:  1\n",
-      "================================================================================================================================================================\n",
-      "Cell [Ang]: \n",
-      "[[ 3.79100000e+00  0.00000000e+00  0.00000000e+00]\n",
-      " [-1.89550000e+00  3.28310231e+00  0.00000000e+00]\n",
-      " [ 1.25954923e-15  2.18160327e-15  2.05700000e+01]]\n",
-      "================================================================================================================================================================\n",
-      "DFT axis: \n",
-      "[0 0 1]\n",
-      "Quantization axis and perpendicular rotation directions:\n",
-      "[1 0 0]  --»  [array([0, 1, 0]), array([0, 0, 1])]\n",
-      "[0 1 0]  --»  [array([1, 0, 0]), array([0, 0, 1])]\n",
-      "[0 0 1]  --»  [array([1, 0, 0]), array([0, 1, 0])]\n",
-      "================================================================================================================================================================\n",
-      "number of k points:  20\n",
-      "k point directions:  xy\n",
-      "================================================================================================================================================================\n",
-      "Parameters for the contour integral:\n",
-      "Ebot:  -30\n",
-      "Eset:  50\n",
-      "Esetp:  10000\n",
-      "================================================================================================================================================================\n",
-      "Atomic informations: \n",
-      "----------------------------------------------------------------------------------------------------------------------------------------------------------------\n",
-      "[atom index]Element(orbitals)        x [Ang]        y [Ang]        z [Ang]        Sx        Sy        Sz        Q        Lx        Ly        Lz        Jx        Jy        Jz\n",
-      "----------------------------------------------------------------------------------------------------------------------------------------------------------------\n",
-      "[3]Fe(2)        -7.339158738013707e-06        4.149278510690423e-06        11.657585837928032\n",
-      "\n",
-      "[4]Fe(2)        -7.326987662162937e-06        4.158274523275774e-06        8.912422537596708\n",
-      "\n",
-      "[5]Fe(2)        1.8954667088117545        1.0943913231921656        10.285002698393109\n",
-      "\n",
-      "[3]Fe(all)        -7.339158738013707e-06        4.149278510690423e-06        11.657585837928032\n",
-      "[4]Fe(all)        -7.326987662162937e-06        4.158274523275774e-06        8.912422537596708\n",
-      "\n",
-      "================================================================================================================================================================\n",
-      "Exchange [meV]\n",
-      "----------------------------------------------------------------------------------------------------------------------------------------------------------------\n",
-      "Magnetic entity1        Magnetic entity2       [i  j  k]       d [Ang]\n",
-      "----------------------------------------------------------------------------------------------------------------------------------------------------------------\n",
-      "[3]Fe(2)        {[3]Fe(all)--[4]Fe(all)}             [0 0 0]    d [Ang] Not yet.\n",
-      "Isotropic:  -7405.318419294225\n",
-      "DMI:  [-9.31336133e-01 -7.86063119e-04  1.83575526e-06]\n",
-      "Symmetric-anisotropy:  [ 2.65595540e+00 -2.44949089e+00 -1.40962466e-03 -6.58426045e-04\n",
-      " -6.97969537e-02]\n",
-      "Energies for debugging: \n",
-      "array([[-7.40776791e+00, -8.61539179e-04,  1.00113309e-03,\n",
-      "        -7.40552488e+00],\n",
-      "       [-7.40477088e+00,  6.58426045e-07, -9.13700193e-07,\n",
-      "        -7.40266246e+00],\n",
-      "       [-7.40364830e+00,  1.41146041e-06,  1.40778890e-06,\n",
-      "        -7.40364662e+00]])\n",
-      "J_ii for debugging: (check if this is the same as in calculate_exchange_tensor)\n",
-      "array([-7.40266246, -7.4036483 , -7.40776791])\n",
-      "Test J_xx = E(y,z) = E(z,y)\n",
-      "-7.402662463889994 -7.403646622630278\n",
-      "\n",
-      "[3]Fe(2)        [4]Fe(2)             [0 0 0]    d [Ang] Not yet.\n",
-      "Isotropic:  -63.511236909297565\n",
-      "DMI:  [-9.32966923e-01 -8.92579299e-04 -2.04258659e-06]\n",
-      "Symmetric-anisotropy:  [-3.33948782e+00  1.29796666e+00  6.90431275e-04 -8.11057566e-04\n",
-      " -5.49031203e-06]\n",
-      "Energies for debugging: \n",
-      "array([[-6.22132703e-02, -9.32961433e-04,  9.32972414e-04,\n",
-      "        -6.14697157e-02],\n",
-      "       [-6.75757803e-02,  8.11057566e-07, -9.74101032e-07,\n",
-      "        -6.68507247e-02],\n",
-      "       [-5.36159327e-02, -6.92473862e-07, -6.88388688e-07,\n",
-      "        -5.36168088e-02]])\n",
-      "J_ii for debugging: (check if this is the same as in calculate_exchange_tensor)\n",
-      "array([-0.06685072, -0.05361593, -0.06221327])\n",
-      "Test J_xx = E(y,z) = E(z,y)\n",
-      "-0.0668507247298092 -0.05361680878933305\n",
-      "\n",
-      "[4]Fe(2)        [3]Fe(2)             [0 0 0]    d [Ang] Not yet.\n",
-      "Isotropic:  -63.511236909297594\n",
-      "DMI:  [9.32966923e-01 8.92579299e-04 2.04258658e-06]\n",
-      "Symmetric-anisotropy:  [-3.33948782e+00  1.29796666e+00  6.90431275e-04  9.74101032e-04\n",
-      " -5.49031203e-06]\n",
-      "Energies for debugging: \n",
-      "array([[-6.22132703e-02,  9.32972414e-04, -9.32961433e-04,\n",
-      "        -6.14697157e-02],\n",
-      "       [-6.75757803e-02, -9.74101032e-07,  8.11057566e-07,\n",
-      "        -6.68507247e-02],\n",
-      "       [-5.36159327e-02, -6.88388689e-07, -6.92473862e-07,\n",
-      "        -5.36168088e-02]])\n",
-      "J_ii for debugging: (check if this is the same as in calculate_exchange_tensor)\n",
-      "array([-0.06685072, -0.05361593, -0.06221327])\n",
-      "Test J_xx = E(y,z) = E(z,y)\n",
-      "-0.06685072472980924 -0.05361680878933306\n",
-      "\n",
-      "[3]Fe(2)        [5]Fe(2)             [0 0 0]    d [Ang] Not yet.\n",
-      "Isotropic:  -60.982930301608\n",
-      "DMI:  [ 3.78506176e+00 -6.13838308e+00  3.59037036e-03]\n",
-      "Symmetric-anisotropy:  [ 0.08076008  0.11287844  0.07106945 -6.23149871 -0.0424978 ]\n",
-      "Energies for debugging: \n",
-      "array([[-6.08700519e-02,  3.82755956e-03, -3.74256396e-03,\n",
-      "        -6.11765688e-02],\n",
-      "       [-6.06451611e-02,  6.23149871e-03, -6.04526746e-03,\n",
-      "        -6.09021702e-02],\n",
-      "       [-5.98973147e-02, -6.74790753e-05, -7.46598160e-05,\n",
-      "        -5.98149915e-02]])\n",
-      "J_ii for debugging: (check if this is the same as in calculate_exchange_tensor)\n",
-      "array([-0.06090217, -0.05989731, -0.06087005])\n",
-      "Test J_xx = E(y,z) = E(z,y)\n",
-      "-0.06090217021953539 -0.059814991464384966\n",
-      "\n",
-      "[4]Fe(2)        [5]Fe(2)             [0 0 0]    d [Ang] Not yet.\n",
-      "Isotropic:  -60.97738429102214\n",
-      "DMI:  [-3.79945963e+00  6.15244494e+00  3.58990840e-03]\n",
-      "Symmetric-anisotropy:  [0.09125387 0.11162376 0.07106826 6.24185638 0.03636701]\n",
-      "Energies for debugging: \n",
-      "array([[-6.08657605e-02, -3.83582664e-03,  3.76309262e-03,\n",
-      "        -6.11802619e-02],\n",
-      "       [-6.06341379e-02, -6.24185638e-03,  6.06303351e-03,\n",
-      "        -6.08861304e-02],\n",
-      "       [-5.98973587e-02, -6.74783547e-05, -7.46581715e-05,\n",
-      "        -5.98150413e-02]])\n",
-      "J_ii for debugging: (check if this is the same as in calculate_exchange_tensor)\n",
-      "array([-0.06088613, -0.05989736, -0.06086576])\n",
-      "Test J_xx = E(y,z) = E(z,y)\n",
-      "-0.06088613041916854 -0.05981504131171201\n",
-      "\n",
-      "[3]Fe(2)        [5]Fe(2)             [-1  0  0]    d [Ang] Not yet.\n",
-      "Isotropic:  -6.828322238590878\n",
-      "DMI:  [5.95251705 7.64859703 6.50501652]\n",
-      "Symmetric-anisotropy:  [-0.45523613  0.13725614 -0.031302    7.69961304  0.03239586]\n",
-      "Energies for debugging: \n",
-      "array([[-0.00669107,  0.00592012, -0.00598491, -0.00651034],\n",
-      "       [-0.00748301, -0.00769961,  0.00759758, -0.00728356],\n",
-      "       [-0.00590136,  0.00653632, -0.00647371, -0.00567313]])\n",
-      "J_ii for debugging: (check if this is the same as in calculate_exchange_tensor)\n",
-      "array([-0.00728356, -0.00590136, -0.00669107])\n",
-      "Test J_xx = E(y,z) = E(z,y)\n",
-      "-0.007283558363709989 -0.005673134248977463\n",
-      "\n",
-      "[4]Fe(2)        [5]Fe(2)             [-1  0  0]    d [Ang] Not yet.\n",
-      "Isotropic:  -6.332622110203909\n",
-      "DMI:  [6.19414647 4.23019689 6.50504332]\n",
-      "Symmetric-anisotropy:  [ 0.53646663 -0.33927457 -0.03129943  4.24610256 -0.09833472]\n",
-      "Energies for debugging: \n",
-      "array([[-0.0066719 ,  0.00629248, -0.00609581, -0.00652981],\n",
-      "       [-0.00599877, -0.0042461 ,  0.00421429, -0.00579616],\n",
-      "       [-0.00590121,  0.00653634, -0.00647374, -0.00567299]])\n",
-      "J_ii for debugging: (check if this is the same as in calculate_exchange_tensor)\n",
-      "array([-0.00579616, -0.00590121, -0.0066719 ])\n",
-      "Test J_xx = E(y,z) = E(z,y)\n",
-      "-0.005796155475608501 -0.00567298926736884\n",
-      "\n",
-      "================================================================================================================================================================\n",
-      "Runtime information: \n",
-      "Total runtime: 112.18818 s\n",
-      "----------------------------------------------------------------------------------------------------------------------------------------------------------------\n",
-      "Initial setup: 0.13208074999999986 s\n",
-      "Hamiltonian conversion and XC field extraction: 0.530 s\n",
-      "Pair and site datastructure creatrions: 0.010 s\n",
-      "k set cration and distribution: 0.014 s\n",
-      "Rotating XC potential: 0.220 s\n",
-      "Greens function inversion: 111.167 s\n",
-      "Calculate energies and magnetic components: 0.115 s\n"
-     ]
-    }
-   ],
+   "outputs": [],
    "source": [
     "if rank == root_node:\n",
     "    # iterate over the magnetic entities\n",
@@ -644,18 +464,251 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 9,
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "##################################################################### GROGU OUTPUT #############################################################################\n",
+    "================================================================================================================================================================\n",
+    "Input file: \n",
+    "Not yet specified.\n",
+    "Number of nodes in the parallel cluster:  1\n",
+    "================================================================================================================================================================\n",
+    "Cell [Ang]: \n",
+    "[[ 6.47        0.          0.        ]\n",
+    " [-3.235       5.60318436  0.        ]\n",
+    " [ 0.          0.         29.999449  ]]\n",
+    "================================================================================================================================================================\n",
+    "DFT axis: \n",
+    "[0 0 1]\n",
+    "Quantization axis and perpendicular rotation directions:\n",
+    "[1 0 0]  --»  [array([0, 1, 0]), array([0, 0, 1])]\n",
+    "[0 1 0]  --»  [array([1, 0, 0]), array([0, 0, 1])]\n",
+    "[0 0 1]  --»  [array([1, 0, 0]), array([0, 1, 0])]\n",
+    "================================================================================================================================================================\n",
+    "number of k points:  4\n",
+    "k point directions:  xy\n",
+    "================================================================================================================================================================\n",
+    "Parameters for the contour integral:\n",
+    "Ebot:  -30\n",
+    "Eset:  100\n",
+    "Esetp:  10000\n",
+    "================================================================================================================================================================\n",
+    "Atomic informations: \n",
+    "----------------------------------------------------------------------------------------------------------------------------------------------------------------\n",
+    "[atom index]Element(orbitals)        x [Ang]        y [Ang]        z [Ang]        Sx        Sy        Sz        Q        Lx        Ly        Lz        Jx        Jy        Jz\n",
+    "----------------------------------------------------------------------------------------------------------------------------------------------------------------\n",
+    "[0]Cr(2)        0.025224941687312104        -0.014371429816522993        -0.020792425897793128\n",
+    "\n",
+    "[1]Cr(2)        3.2508367285399955        1.8478445793345957        -0.02082986959108256\n",
+    "\n",
+    "[2]Br(2)        2.1053343349045215        0.10721337338160676        -1.471910834765145\n",
+    "\n",
+    "================================================================================================================================================================\n",
+    "Exchange [meV]\n",
+    "----------------------------------------------------------------------------------------------------------------------------------------------------------------\n",
+    "Magnetic entity1        Magnetic entity2       [i  j  k]       d [Ang]\n",
+    "----------------------------------------------------------------------------------------------------------------------------------------------------------------\n",
+    "[0]Cr(2)        [1]Cr(2)             [0 0 0]    d [Ang] Not yet.\n",
+    "Isotropic:  15.91078131425419\n",
+    "DMI:  [-3.16312278e-05  1.94021443e-05 -1.65972633e-04]\n",
+    "Symmetric-anisotropy:  [-0.03896181 -0.01392426 -0.07639227  0.05293223 -0.08712209]\n",
+    "Energies for debugging: \n",
+    "array([[ 1.58968571e-02,  8.70904601e-05,  8.71537225e-05,\n",
+    "         1.59636674e-02],\n",
+    "       [ 1.58999738e-02, -5.29322330e-05, -5.28934287e-05,\n",
+    "         1.58718195e-02],\n",
+    "       [ 1.59589510e-02,  7.62263019e-05,  7.65582472e-05,\n",
+    "         1.58707914e-02]])\n",
+    "J_ii for debugging: (check if this is the same as in calculate_exchange_tensor)\n",
+    "array([0.01587182, 0.01595895, 0.01589686])\n",
+    "Test J_xx = E(y,z) = E(z,y)\n",
+    "0.015871819503030533 0.015870791381430095\n",
+    "\n",
+    "[0]Cr(2)        [2]Br(2)             [0 0 0]    d [Ang] Not yet.\n",
+    "Isotropic:  0.14694676346924165\n",
+    "DMI:  [ 8.01967166e-05  8.29086493e-04 -9.27235617e-05]\n",
+    "Symmetric-anisotropy:  [ 2.66339812e-05 -4.77282294e-06 -6.60708923e-06  8.03002228e-04\n",
+    " -2.98623982e-06]\n",
+    "Energies for debugging: \n",
+    "array([[ 1.46941991e-04,  8.31829565e-08, -7.72104768e-08,\n",
+    "         1.46924902e-04],\n",
+    "       [ 1.46982775e-04, -8.03002228e-07,  8.55170759e-07,\n",
+    "         1.46973397e-04],\n",
+    "       [ 1.46942514e-04, -8.61164724e-08,  9.93306509e-08,\n",
+    "         1.46901824e-04]])\n",
+    "J_ii for debugging: (check if this is the same as in calculate_exchange_tensor)\n",
+    "array([0.00014697, 0.00014694, 0.00014694])\n",
+    "Test J_xx = E(y,z) = E(z,y)\n",
+    "0.00014697339745047312 0.00014690182371456808\n",
+    "\n",
+    "[1]Cr(2)        [2]Br(2)             [0 0 0]    d [Ang] Not yet.\n",
+    "Isotropic:  0.146966596490555\n",
+    "DMI:  [ 7.83469069e-04 -3.25183011e-04  9.29353203e-05]\n",
+    "Symmetric-anisotropy:  [-1.73027923e-05  1.26068464e-05 -2.08606748e-05 -2.97793770e-04\n",
+    "  3.14336799e-07]\n",
+    "Energies for debugging: \n",
+    "array([[ 1.46979203e-04,  7.83154732e-07, -7.83783406e-07,\n",
+    "         1.46971292e-04],\n",
+    "       [ 1.46912608e-04,  2.97793770e-07, -3.52572253e-07,\n",
+    "         1.46949294e-04],\n",
+    "       [ 1.46905570e-04,  1.13795995e-07, -7.20746454e-08,\n",
+    "         1.46914411e-04]])\n",
+    "J_ii for debugging: (check if this is the same as in calculate_exchange_tensor)\n",
+    "array([0.00014695, 0.00014691, 0.00014698])\n",
+    "Test J_xx = E(y,z) = E(z,y)\n",
+    "0.00014694929369829223 0.00014691441062702917\n",
+    "\n",
+    "[0]Cr(2)        [1]Cr(2)             [1 0 0]    d [Ang] Not yet.\n",
+    "Isotropic:  0.09111814992141927\n",
+    "DMI:  [0.03377305 0.0340363  0.00964244]\n",
+    "Symmetric-anisotropy:  [-0.00291957  0.0031014  -0.0021492   0.03416682  0.0051237 ]\n",
+    "Energies for debugging: \n",
+    "array([[ 9.42195473e-05,  2.86493440e-05, -3.88967487e-05,\n",
+    "         9.09363256e-05],\n",
+    "       [ 9.41093226e-05, -3.41668178e-05,  3.39057913e-05,\n",
+    "         8.81985768e-05],\n",
+    "       [ 9.14933866e-05,  1.17916416e-05, -7.49323868e-06,\n",
+    "         8.86473282e-05]])\n",
+    "J_ii for debugging: (check if this is the same as in calculate_exchange_tensor)\n",
+    "array([8.81985768e-05, 9.14933866e-05, 9.42195473e-05])\n",
+    "Test J_xx = E(y,z) = E(z,y)\n",
+    "8.819857682132204e-05 8.8647328224418e-05\n",
+    "\n",
+    "[0]Cr(2)        [2]Br(2)             [1 0 0]    d [Ang] Not yet.\n",
+    "Isotropic:  0.00023346314999400802\n",
+    "DMI:  [-0.0002294  -0.0002349  -0.00020879]\n",
+    "Symmetric-anisotropy:  [-1.12556482e-05  1.43988058e-05 -1.55475448e-06 -2.34959571e-04\n",
+    "  6.19663196e-06]\n",
+    "Energies for debugging: \n",
+    "array([[ 2.47861956e-07, -2.35597865e-07,  2.23204601e-07,\n",
+    "         2.30319992e-07],\n",
+    "       [ 2.10261125e-07,  2.34959571e-07, -2.34837442e-07,\n",
+    "         2.22207502e-07],\n",
+    "       [ 2.12323581e-07, -2.07238961e-07,  2.10348470e-07,\n",
+    "         2.07858965e-07]])\n",
+    "J_ii for debugging: (check if this is the same as in calculate_exchange_tensor)\n",
+    "array([2.22207502e-07, 2.12323581e-07, 2.47861956e-07])\n",
+    "Test J_xx = E(y,z) = E(z,y)\n",
+    "2.2220750178916023e-07 2.078589648170586e-07\n",
+    "\n",
+    "[0]Cr(2)        [1]Cr(2)             [-1  0  0]    d [Ang] Not yet.\n",
+    "Isotropic:  0.09111258017019819\n",
+    "DMI:  [-0.03378783 -0.03403353 -0.0096277 ]\n",
+    "Symmetric-anisotropy:  [-0.00291741  0.00310155 -0.00214665 -0.03390027  0.00512057]\n",
+    "Energies for debugging: \n",
+    "array([[ 9.42141309e-05, -3.89083994e-05,  2.86672518e-05,\n",
+    "         9.09284386e-05],\n",
+    "       [ 9.41042657e-05,  3.39002650e-05, -3.41667971e-05,\n",
+    "         8.81951711e-05],\n",
+    "       [ 9.14629700e-05, -7.48105019e-06,  1.17743519e-05,\n",
+    "         8.86174723e-05]])\n",
+    "J_ii for debugging: (check if this is the same as in calculate_exchange_tensor)\n",
+    "array([8.81951711e-05, 9.14629700e-05, 9.42141309e-05])\n",
+    "Test J_xx = E(y,z) = E(z,y)\n",
+    "8.819517110814147e-05 8.861747229207156e-05\n",
+    "\n",
+    "[0]Cr(2)        [2]Br(2)             [-1  0  0]    d [Ang] Not yet.\n",
+    "Isotropic:  0.0002516119318884681\n",
+    "DMI:  [0.0002271  0.00022722 0.00023639]\n",
+    "Symmetric-anisotropy:  [ 8.12786238e-06 -1.33724299e-05  2.28353801e-06  2.31211900e-04\n",
+    " -2.16488405e-06]\n",
+    "Energies for debugging: \n",
+    "array([[ 2.38239502e-07,  2.29269073e-07, -2.24939305e-07,\n",
+    "         2.56856499e-07],\n",
+    "       [ 2.77408158e-07, -2.31211900e-07,  2.23219979e-07,\n",
+    "         2.59739794e-07],\n",
+    "       [ 2.73816864e-07,  2.34111415e-07, -2.38678491e-07,\n",
+    "         2.71887108e-07]])\n",
+    "J_ii for debugging: (check if this is the same as in calculate_exchange_tensor)\n",
+    "array([2.59739794e-07, 2.73816864e-07, 2.38239502e-07])\n",
+    "Test J_xx = E(y,z) = E(z,y)\n",
+    "2.5973979426862243e-07 2.718871078150382e-07\n",
+    "\n",
+    "[0]Cr(2)        [1]Cr(2)             [0 1 0]    d [Ang] Not yet.\n",
+    "Isotropic:  -0.13657014924460614\n",
+    "DMI:  [ 0.00026859 -0.00029338 -0.0001818 ]\n",
+    "Symmetric-anisotropy:  [-0.00104734  0.00154555 -0.00046486 -0.00072287  0.00087507]\n",
+    "Energies for debugging: \n",
+    "array([[-1.35024601e-04, -6.06475771e-07, -1.14365755e-06,\n",
+    "        -1.37068359e-04],\n",
+    "       [-1.34947686e-04,  7.22871593e-07,  1.36112752e-07,\n",
+    "        -1.37617487e-04],\n",
+    "       [-1.36691877e-04,  2.83065496e-07,  6.46663000e-07,\n",
+    "        -1.37228934e-04]])\n",
+    "J_ii for debugging: (check if this is the same as in calculate_exchange_tensor)\n",
+    "array([-0.00013762, -0.00013669, -0.00013502])\n",
+    "Test J_xx = E(y,z) = E(z,y)\n",
+    "-0.00013761748740593208 -0.00013722893403347796\n",
+    "\n",
+    "[0]Cr(2)        [2]Br(2)             [0 1 0]    d [Ang] Not yet.\n",
+    "Isotropic:  -0.00010028627777987032\n",
+    "DMI:  [-8.54009417e-05 -1.21477411e-04 -8.32475011e-05]\n",
+    "Symmetric-anisotropy:  [ 1.89621993e-05 -7.40845210e-06  6.91965022e-07 -1.13803644e-04\n",
+    " -7.00808812e-06]\n",
+    "Energies for debugging: \n",
+    "array([[-1.07694730e-07, -7.83928536e-08,  9.24090298e-08,\n",
+    "        -1.11840025e-07],\n",
+    "       [-8.13196124e-08,  1.13803644e-07, -1.29151178e-07,\n",
+    "        -8.13240785e-08],\n",
+    "       [-1.76542958e-07, -8.39394661e-08,  8.25555360e-08,\n",
+    "        -1.76140161e-07]])\n",
+    "J_ii for debugging: (check if this is the same as in calculate_exchange_tensor)\n",
+    "array([-8.13240785e-08, -1.76542958e-07, -1.07694730e-07])\n",
+    "Test J_xx = E(y,z) = E(z,y)\n",
+    "-8.132407846090182e-08 -1.7614016112304844e-07\n",
+    "\n",
+    "[0]Cr(2)        [1]Cr(2)             [0 1 0]    d [Ang] Not yet.\n",
+    "Isotropic:  -0.13657014924460614\n",
+    "DMI:  [ 0.00026859 -0.00029338 -0.0001818 ]\n",
+    "Symmetric-anisotropy:  [-0.00104734  0.00154555 -0.00046486 -0.00072287  0.00087507]\n",
+    "Energies for debugging: \n",
+    "array([[-1.35024601e-04, -6.06475771e-07, -1.14365755e-06,\n",
+    "        -1.37068359e-04],\n",
+    "       [-1.34947686e-04,  7.22871593e-07,  1.36112752e-07,\n",
+    "        -1.37617487e-04],\n",
+    "       [-1.36691877e-04,  2.83065496e-07,  6.46663000e-07,\n",
+    "        -1.37228934e-04]])\n",
+    "J_ii for debugging: (check if this is the same as in calculate_exchange_tensor)\n",
+    "array([-0.00013762, -0.00013669, -0.00013502])\n",
+    "Test J_xx = E(y,z) = E(z,y)\n",
+    "-0.00013761748740593208 -0.00013722893403347796\n",
+    "\n",
+    "[0]Cr(2)        [2]Br(2)             [0 1 0]    d [Ang] Not yet.\n",
+    "Isotropic:  -0.00010028627777987032\n",
+    "DMI:  [-8.54009417e-05 -1.21477411e-04 -8.32475011e-05]\n",
+    "Symmetric-anisotropy:  [ 1.89621993e-05 -7.40845210e-06  6.91965022e-07 -1.13803644e-04\n",
+    " -7.00808812e-06]\n",
+    "Energies for debugging: \n",
+    "array([[-1.07694730e-07, -7.83928536e-08,  9.24090298e-08,\n",
+    "        -1.11840025e-07],\n",
+    "       [-8.13196124e-08,  1.13803644e-07, -1.29151178e-07,\n",
+    "        -8.13240785e-08],\n",
+    "       [-1.76542958e-07, -8.39394661e-08,  8.25555360e-08,\n",
+    "        -1.76140161e-07]])\n",
+    "J_ii for debugging: (check if this is the same as in calculate_exchange_tensor)\n",
+    "array([-8.13240785e-08, -1.76542958e-07, -1.07694730e-07])\n",
+    "Test J_xx = E(y,z) = E(z,y)\n",
+    "-8.132407846090182e-08 -1.7614016112304844e-07\n",
+    "\n",
+    "================================================================================================================================================================\n",
+    "Runtime information: \n",
+    "Total runtime: 16.641370917 s\n",
+    "----------------------------------------------------------------------------------------------------------------------------------------------------------------\n",
+    "Initial setup: 0.07859033299999996 s\n",
+    "Hamiltonian conversion and XC field extraction: 0.683 s\n",
+    "Pair and site datastructure creatrions: 0.012 s\n",
+    "k set cration and distribution: 0.023 s\n",
+    "Rotating XC potential: 0.243 s\n",
+    "Greens function inversion: 15.524 s\n",
+    "Calculate energies and magnetic components: 0.078 s"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
    "metadata": {},
-   "outputs": [
-    {
-     "ename": "SyntaxError",
-     "evalue": "invalid syntax (3105939143.py, line 1)",
-     "output_type": "error",
-     "traceback": [
-      "\u001b[0;36m  Cell \u001b[0;32mIn[9], line 1\u001b[0;36m\u001b[0m\n\u001b[0;31m    ========================================\u001b[0m\n\u001b[0m    ^\u001b[0m\n\u001b[0;31mSyntaxError\u001b[0m\u001b[0;31m:\u001b[0m invalid syntax\n"
-     ]
-    }
-   ],
+   "outputs": [],
    "source": [
     "========================================\n",
     " \n",
@@ -699,27 +752,9 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 27,
+   "execution_count": null,
    "metadata": {},
-   "outputs": [
-    {
-     "name": "stdout",
-     "output_type": "stream",
-     "text": [
-      "xyz[-3:]: red, green, blue\n"
-     ]
-    },
-    {
-     "data": {
-      "image/png": "iVBORw0KGgoAAAANSUhEUgAABNoAAAHACAYAAAB0/gUQAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAABCsElEQVR4nO3df5hWdZ0//ucMyAwWM2rK8GMnf5b4E0iTMFu1Rln1YmU/5aLtAktqW9l+VLYfUiqWJqlpbEayaUbq9vG35ie9IGUjM1lNlNZMTRSElEH9qjOICjJzf//g47QToAOcmXtmeDyu61x5n/u87/v1PkznNfdzzn1ORalUKgUAAAAA2CqV5S4AAAAAAHoDQRsAAAAAFEDQBgAAAAAFELQBAAAAQAEEbQAAAABQAEEbAAAAABRA0AYAAAAABRC0AQAAAEAB+pa7gO6otbU1zz//fAYMGJCKiopylwPQ45VKpaxatSpDhgxJZaW/8egzAMXSZzak1wAUq6O9RtC2Ec8//3zq6+vLXQZAr7N8+fL81V/9VbnLKDt9BqBz6DN/ptcAdI536zWCto0YMGBAkvU7r6ampszVAPR8zc3Nqa+vbzu+buv0GYBi6TMb0msAitXRXiNo24i3T62uqanRlAAK5Ksr6+kzAJ1Dn/kzvQagc7xbr3EBAwAAAAAogKANAAAAAAogaAMAAACAAgjaAAAAAKAAgjYAAAAAKICgDQAAAAAKIGgDAAAAgAII2gAAAACgAII2AAAAAChA33IX0KusWJHMnp388Y/JgAHJCSckhx2WVFSUuzIAeoM330xuvDG59971veXww5NPfSqpri53ZQAA0C2tW5fccUcyd+76//7wh5N/+If1sU1nELQV5fLLkzPPTEqlpLLyz+uOOCK5/faktrac1QHQ0/32t8lxxyUvvpj0/X/t+6qrki99KbnrruRDHypvfQAA0M0sXpyMGZM888yff4X+8Y+TL385ueWW5Oiji39PXx0twq23Jv/7fyctLUlr6/qIdN269c/9+tfJ+PHlrQ+Anm3FiuSoo5KXX17/+H/2mZdeShoakhdeKF99AADQzbz+evLxjyfPPrv+8du/QpdK658bOzb5wx+Kf19B29YqlZJvfGPTXw9taVl/fuKiRV1aFgC9yKxZyWuvre8pf6mlJWlqWn92GwAAkCS54YZk+fKN/wrd2rp++bd/K/59BW1b67nnkv/+7/WB26b07bv+66MAsCVuumnjvyG8rbV1/TYAAECS9V8+rHyH1GvduvWXPy6aoG1rvf76u29TUZG88Ubn1wJA77R6dTHbAADANmL16vV/j34nb75Z/PsK2rZWfX3ynve88zZvvZXsv3/X1ANA7zNy5J+v3roxffuu3wYAAEiSDB/+zr9CV1Z2TlQjaNta/fsnJ5+c9Omz8ecrKpIddkhOOKFLywKgF/nCF/5884ONWbcu+fznu64eAADo5v75n9/5V+jW1uSLXyz+fQVtRfjmN5NhwzYM2/r0Wb/8x38k1dXlqQ2Anu+oo5LTTlv/3//z5jtvX3RiypTkiCO6vCwAAOiuhg1LLr54/X//z2u1VVSsX/7X/0r+8R+Lf19BWxFqa5Pf/CY566xkp53Wr6usXH+v2PvvT449trz1AdCzVVQkl1+e/PjHyb77/nn9/vsn116bfOc75asNAAC6qS9/OfnZz5KPfOTP63bfff3dRm+8cdNfTtwaFaXSO90uc9vU3Nyc2traNDU1paamZvMGt7YmTU3J9tsnVVWdUyBAD7NVx9VeaKv3R3Pz+vBtwIDiiwPogfSZDdknAO2tXr3+q6Q1Ne2/JNJRHT2uvsNl4dgilZXJjjuWuwoAejMfmAAAYLO8230si+KrowAAAABQAEEbAAAAABRA0AYAAAAABRC0AQAAAEABBG0AAAAAUABBGwAAAAAUQNAGAAAAAAUQtAEAAABAAQRtAAAAAFCAsgZt9957b8aOHZshQ4akoqIit99+e7vnb7311hx99NF53/vel4qKiixatOhdX3P27NmpqKhot1RXV3fOBADo1vQZAACgK5U1aFu9enWGDx+emTNnbvL5ww47LBdddNFmvW5NTU1WrFjRtjz77LNFlAtAD6PPAAAAXalvOd/8mGOOyTHHHLPJ5ydMmJAkWbp06Wa9bkVFRQYNGrQ1pQHQC+gzAABAV+qV12h77bXXsuuuu6a+vj7HH398HnvssXfcfs2aNWlubm63AMCm6DMAAMDG9Lqgbe+9987VV1+dn/3sZ7nuuuvS2tqaQw89NH/60582OWb69Ompra1tW+rr67uwYgB6En0GAADYlF4XtI0ePToTJ07MiBEjcvjhh+fWW2/NLrvskn//93/f5JipU6emqampbVm+fHkXVgxAT6LPAAAAm1LWa7R1he222y4jR47M4sWLN7lNVVVVqqqqurAqAHoLfQYAAHhbrzuj7S+1tLTk0UcfzeDBg8tdCgC9kD4DAAC8raxntL322mvtzgBYsmRJFi1alJ122invf//78/LLL2fZsmV5/vnnkyRPPvlkkmTQoEFtd3ubOHFihg4dmunTpydJvvnNb+YjH/lI9tprr7z66qu55JJL8uyzz+aUU07p4tkBUG76DAAA0JXKGrQ99NBDOfLII9seT5kyJUkyadKkzJ49O3fccUcmT57c9vyJJ56YJJk2bVrOO++8JMmyZctSWfnnE/NeeeWVnHrqqWlsbMyOO+6Ygw46KPfff3/23XffLpgRAN2JPgMAAHSlilKpVCp3Ed1Nc3Nzamtr09TUlJqamnKXA9DjOa62Z38AFKu7H1fvvffeXHLJJVm4cGFWrFiR2267LePGjXvHMfPnz8+UKVPy2GOPpb6+PmeffXb+6Z/+qcPv2d33CUBP09Hjaq+/RhsAAEA5rV69OsOHD8/MmTM7tP2SJUty3HHH5cgjj8yiRYtyxhln5JRTTsncuXM7uVIAtlavv+soAABAOR1zzDE55phjOrz9rFmzsvvuu+fSSy9Nkuyzzz6577778t3vfjdjxozprDIBKIAz2gAAALqRBQsWpKGhod26MWPGZMGCBZscs2bNmjQ3N7dbAOh6gjYAAIBupLGxMXV1de3W1dXVpbm5OW+88cZGx0yfPj21tbVtS319fVeUCsBfELQBAAD0cFOnTk1TU1Pbsnz58nKXBLBNco02AACAbmTQoEFZuXJlu3UrV65MTU1N+vfvv9ExVVVVqaqq6oryAHgHzmgDAADoRkaPHp158+a1W3f33Xdn9OjRZaoIgI4StAEAAHSi1157LYsWLcqiRYuSJEuWLMmiRYuybNmyJOu/9jlx4sS27T/3uc/lmWeeyVe+8pU88cQT+cEPfpAbb7wxZ555ZjnKB2AzCNoAAAA60UMPPZSRI0dm5MiRSZIpU6Zk5MiROffcc5MkK1asaAvdkmT33XfPnXfembvvvjvDhw/PpZdemquuuipjxowpS/0AdJxrtAEAAHSiI444IqVSaZPPz549e6NjHnnkkU6sCoDO4Iw2AAAAACiAoA0AAAAACiBoAwAAAIACCNoAAAAAoACCNgAAAAAogKANAAAAAAogaAMAAACAAgjaAAAAAKAAgjYAAAAAKICgDQAAAAAKIGgDAAAAgAII2gAAAACgAII2AAAAACiAoA0AAAAACiBoAwAAAIACCNoAAAAAoACCNgAAAAAogKANAAAAAAogaAMAAACAAgjaAAAAAKAAgjYAAAAAKICgDQAAAAAKIGgDAAAAgAII2gAAAACgAII2AAAAACiAoA0AAAAACiBoAwAAAIACCNoAAAAAoACCNgAAAAAogKANAAAAAAogaAMAAACAAgjaAAAAAKAAgjYAAAAAKICgDQAAAAAKIGgDAAAAgAII2gAAAACgAII2AAAAACiAoA0AAAAACiBoAwAAAIACCNoAAAAAoACCNgAAAAAogKANAAAAAAogaAMAAACAApQ1aLv33nszduzYDBkyJBUVFbn99tvbPX/rrbfm6KOPzvve975UVFRk0aJFHXrdm266KcOGDUt1dXUOOOCA3HXXXcUXD0C3p88AAABdqaxB2+rVqzN8+PDMnDlzk88fdthhueiiizr8mvfff39OOumknHzyyXnkkUcybty4jBs3Lr///e+LKhuAHkKfAQAAulJFqVQqlbuIJKmoqMhtt92WcePGbfDc0qVLs/vuu+eRRx7JiBEj3vF1xo8fn9WrV+fnP/9527qPfOQjGTFiRGbNmtWhWpqbm1NbW5umpqbU1NRszjQA2IjucFzVZwB6L8fVDdknAMXq6HG1112jbcGCBWloaGi3bsyYMVmwYMEmx6xZsybNzc3tFgDYGH0GAADYlF4XtDU2Nqaurq7durq6ujQ2Nm5yzPTp01NbW9u21NfXd3aZAPRQ+gwAALApvS5o2xJTp05NU1NT27J8+fJylwRAL6LPAADAtqFvuQso2qBBg7Jy5cp261auXJlBgwZtckxVVVWqqqo6uzQAegF9BgAA2JRed0bb6NGjM2/evHbr7r777owePbpMFQHQm+gzAADAppT1jLbXXnstixcvbnu8ZMmSLFq0KDvttFPe//735+WXX86yZcvy/PPPJ0mefPLJJOvPJnj7zIGJEydm6NChmT59epLk9NNPz+GHH55LL700xx13XK6//vo89NBD+eEPf9jFswOg3PQZAACgK5X1jLaHHnooI0eOzMiRI5MkU6ZMyciRI3PuuecmSe64446MHDkyxx13XJLkxBNPzMiRIzNr1qy211i2bFlWrFjR9vjQQw/NT3/60/zwhz/M8OHDc/PNN+f222/P/vvv34UzA6A70GcAAICuVFEqlUrlLqK7aW5uTm1tbZqamlJTU1PucgB6PMfV9uwPgGL1hOPqzJkzc8kll6SxsTHDhw/P5ZdfnkMOOWST28+YMSNXXHFFli1blp133jmf+tSnMn369FRXV3fo/XrCPgHoSTp6XO1112gDAADoTm644YZMmTIl06ZNy8MPP5zhw4dnzJgxeeGFFza6/U9/+tOcddZZmTZtWh5//PH86Ec/yg033JCvfe1rXVw5AJtL0AYAANCJLrvsspx66qmZPHly9t1338yaNSvbb799rr766o1uf//99+ejH/1oPv3pT2e33XbL0UcfnZNOOikPPvhgF1cOwOYStAEAAHSStWvXZuHChWloaGhbV1lZmYaGhixYsGCjYw499NAsXLiwLVh75plnctddd+XYY4/d5PusWbMmzc3N7RYAul5Z7zoKAADQm7300ktpaWlJXV1du/V1dXV54oknNjrm05/+dF566aUcdthhKZVKWbduXT73uc+941dHp0+fnm984xuF1g7A5nNGGwAAQDcyf/78XHjhhfnBD36Qhx9+OLfeemvuvPPOnH/++ZscM3Xq1DQ1NbUty5cv78KKAXibM9oAAAA6yc4775w+ffpk5cqV7davXLkygwYN2uiYc845JxMmTMgpp5ySJDnggAOyevXqfPazn83Xv/71VFZueL5EVVVVqqqqip8AAJvFGW0AAACdpF+/fjnooIMyb968tnWtra2ZN29eRo8evdExr7/++gZhWp8+fZIkpVKp84oFYKs5ow0AAKATTZkyJZMmTcrBBx+cQw45JDNmzMjq1aszefLkJMnEiRMzdOjQTJ8+PUkyduzYXHbZZRk5cmRGjRqVxYsX55xzzsnYsWPbAjcAuidBGwAAQCcaP358XnzxxZx77rlpbGzMiBEjMmfOnLYbJCxbtqzdGWxnn312KioqcvbZZ+e5557LLrvskrFjx+Zb3/pWuaYAQAdVlJx7vIHm5ubU1tamqakpNTU15S4HoMdzXG3P/gAoluPqhuwTgGJ19LjqGm0AAAAAUABBGwAAAAAUQNAGAAAAAAUQtAEAAABAAQRtAAAAAFAAQRsAAAAAFEDQBgAAAAAFELQBAAAAQAEEbQAAAABQAEEbAAAAABRA0AYAAAAABRC0AQAAAEABBG0AAAAAUABBGwAAAAAUQNAGAAAAAAUQtAEAAABAAQRtAAAAAFAAQRsAAAAAFEDQBgAAAAAFELQBAAAAQAEEbQAAAABQAEEbAAAAABRA0AYAAAAABRC0AQAAAEABBG0AAAAAUABBGwAAAAAUQNAGAAAAAAUQtAEAAABAAQRtAAAAAFAAQRsAAAAAFEDQBgAAAAAFELQBAAAAQAEEbQAAAABQAEEbAAAAABRA0AYAAAAABRC0AQAAAEABBG0AAAAAUABBGwAAAAAUQNAGAAAAAAUQtAEAAABAAQRtAAAAAFAAQRsAAAAAFEDQBgAAAAAFKGvQdu+992bs2LEZMmRIKioqcvvtt7d7vlQq5dxzz83gwYPTv3//NDQ05KmnnnrH1zzvvPNSUVHRbhk2bFgnzgKA7kyvAQAAukpZg7bVq1dn+PDhmTlz5kafv/jii/O9730vs2bNygMPPJD3vOc9GTNmTN588813fN399tsvK1asaFvuu+++zigfgB5ArwEAALpK33K++THHHJNjjjlmo8+VSqXMmDEjZ599do4//vgkyTXXXJO6urrcfvvtOfHEEzf5un379s2gQYM6pWYAeha9BgAA6Crd9hptS5YsSWNjYxoaGtrW1dbWZtSoUVmwYME7jn3qqacyZMiQ7LHHHvmHf/iHLFu27B23X7NmTZqbm9stAPR+XdVr9BkAANg2dNugrbGxMUlSV1fXbn1dXV3bcxszatSozJ49O3PmzMkVV1yRJUuW5GMf+1hWrVq1yTHTp09PbW1t21JfX1/MJADo1rqq1+gzAACwbei2QduWOuaYY3LCCSfkwAMPzJgxY3LXXXfl1VdfzY033rjJMVOnTk1TU1Pbsnz58i6sGICeZnN7jT4DAADbhm4btL193ZuVK1e2W79y5crNuibODjvskA9+8INZvHjxJrepqqpKTU1NuwWA3q+reo0+AwAA24ZuG7TtvvvuGTRoUObNm9e2rrm5OQ888EBGjx7d4dd57bXX8vTTT2fw4MGdUSYAPZheAwAAFKmsQdtrr72WRYsWZdGiRUnWX5R60aJFWbZsWSoqKnLGGWfkggsuyB133JFHH300EydOzJAhQzJu3Li21/jEJz6R73//+22Pv/SlL+VXv/pVli5dmvvvvz9/93d/lz59+uSkk07q4tkB0B3oNQAAQFfpW843f+ihh3LkkUe2PZ4yZUqSZNKkSZk9e3a+8pWvZPXq1fnsZz+bV199NYcddljmzJmT6urqtjFPP/10XnrppbbHf/rTn3LSSSfl//v//r/ssssuOeyww/Jf//Vf2WWXXbpuYgB0G3oNAADQVSpKpVKp3EV0N83NzamtrU1TU5Pr6AAUwHG1PfsDoFg94bg6c+bMXHLJJWlsbMzw4cNz+eWX55BDDtnk9q+++mq+/vWv59Zbb83LL7+cXXfdNTNmzMixxx7boffrCfsEoCfp6HG1rGe0AQAA9HY33HBDpkyZklmzZmXUqFGZMWNGxowZkyeffDIDBw7cYPu1a9fmqKOOysCBA3PzzTdn6NChefbZZ7PDDjt0ffEAbBZBGwAAQCe67LLLcuqpp2by5MlJklmzZuXOO+/M1VdfnbPOOmuD7a+++uq8/PLLuf/++7PddtslSXbbbbeuLBmALdRt7zoKAADQ061duzYLFy5MQ0ND27rKyso0NDRkwYIFGx1zxx13ZPTo0TnttNNSV1eX/fffPxdeeGFaWlo2+T5r1qxJc3NzuwWAridoAwAA6CQvvfRSWlpaUldX1259XV1dGhsbNzrmmWeeyc0335yWlpbcddddOeecc3LppZfmggsu2OT7TJ8+PbW1tW1LfX19ofMAoGMEbQAAAN1Ia2trBg4cmB/+8Ic56KCDMn78+Hz961/PrFmzNjlm6tSpaWpqaluWL1/ehRUD8DbXaAMAAOgkO++8c/r06ZOVK1e2W79y5coMGjRoo2MGDx6c7bbbLn369Glbt88++6SxsTFr165Nv379NhhTVVWVqqqqYosHYLM5ow0AAKCT9OvXLwcddFDmzZvXtq61tTXz5s3L6NGjNzrmox/9aBYvXpzW1ta2dX/84x8zePDgjYZsAHQfgjYAAIBONGXKlFx55ZX5yU9+kscffzyf//zns3r16ra7kE6cODFTp05t2/7zn/98Xn755Zx++un54x//mDvvvDMXXnhhTjvttHJNAYAO8tVRAACATjR+/Pi8+OKLOffcc9PY2JgRI0Zkzpw5bTdIWLZsWSor/3wORH19febOnZszzzwzBx54YIYOHZrTTz89X/3qV8s1BQA6qKJUKpXKXUR309zcnNra2jQ1NaWmpqbc5QD0eI6r7dkfAMVyXN2QfQJQrI4eV311FAAAAAAKIGgDAAAAgAII2gAAAACgAII2AAAAACiAoA0AAAAACtC33AUA0L09+WTyk58kzz+fDBqUTJiQ7LdfuasCoLd45ZXkmmuS3/0uqapKxo5NxoxJ+vQpd2UAsPkEbQBsVGtrcvrpyfe/n/Ttm5RKSUVFctFFycknJ7NmrV8PAFvqlluSf/zHZM2aPwdrs2Yl+++fzJmTDB1a3voAYHP56igAG/Wtb60P2ZJk3bqkpWX9/ybJ1VcnX/96+WoDoOd74IFk/Pj1IVuptL7HvN1nnngiOfroPz8GgJ5ii4K2j3/84/nGN76xwfpXXnklH//4x7e6KADK6/XXk0su2fTzpVLyve8lTU2d8/76DEDv9+1vrz9TulTa8Ll165I//CG5666ur+ttkyZNyr333lu+AgDokbYoaJs/f36+//3vZ9y4cVm9enXb+rVr1+ZXv/pVYcUBUB733pusWvXO27z5ZnL33Z3z/voMQO/W0pL83//7zmes9e2b3HZb19X0l5qamtLQ0JAPfOADufDCC/Pcc8+VrxgAeowt/uroPffck8bGxnzkIx/J0qVLCywJgHJ7/fWObffGG51Xgz4D0Hu9fUmCd9La2rl95t3cfvvtee655/L5z38+N9xwQ3bbbbccc8wxufnmm/PWW2+VrzAAurUtDtoGDx6cX/3qVznggAPy4Q9/OPPnzy+wLADKaf/9i91uS+gzAL1XVVWy++7rvzr6Tjqzz3TELrvskilTpuR3v/tdHnjggey1116ZMGFChgwZkjPPPDNPPfVUeQsEoNvZoqCt4v91xKqqqvz0pz/N6aefnr/5m7/JD37wg0KLA6A8PvjB5PDD/3wHuL/Up09y0EHJyJGd8/76DEDv9y//8s7PV1Ymn/lM19TyblasWJG77747d999d/r06ZNjjz02jz76aPbdd99897vfLXd5AHQjfbdkUOkvrlh69tlnZ5999smkSZMKKQqA8rvqqmT06OTVV9tfQ6dv3+S9701+8pPOe299BqD3O+205Oc/T+bPX/810bf16bP+a6VXXJEMGVK28vLWW2/ljjvuyI9//OP84he/yIEHHpgzzjgjn/70p1NTU5Mkue222/KZz3wmZ555ZvkKBaBb2aKgbcmSJdlll13arfvkJz+ZYcOG5aGHHiqkMADKa6+9kocfTqZPT2bPXn+dnKqq5B//Mfna15I99ui899ZnAHq/fv3W31X03/4tufzy5E9/Wr/+r/86mTo1Oeqo8tY3ePDgtLa25qSTTsqDDz6YESNGbLDNkUcemR122KHLawOg+6oo/eVpA6S5uTm1tbVpampq+2sVwLZs3bqkuTkZMCDZbrvNH++42p79AdBeqZQ0Na0P37bffvPHd8Zx9dprr80JJ5yQ6urqQl6vq+k1AMXq6HF1i85oA2Db0rdvstNO5a4CgN6qoiLpbieGTZgwodwlANADbfFdRwEAAACAPxO0AQAAAEABBG0AAAAAUABBGwAAAAAUQNAGAAAAAAUQtAEAAABAAQRtAAAAAFAAQRsAAAAAFEDQBgAAAAAFELQBAAAAQAEEbQAAAABQAEEbAAAAABRA0AYAAAAABRC0AQAAAEABBG0AAAAAUABBGwAAAAAUQNAGAAAAAAUQtAEAAABAAQRtAAAAAFAAQRsAAAAAFEDQBgAAAAAFELQBAAAAQAEEbQAAAABQAEEbAAAAABRA0AYAAAAABShr0Hbvvfdm7NixGTJkSCoqKnL77be3e75UKuXcc8/N4MGD079//zQ0NOSpp55619edOXNmdtttt1RXV2fUqFF58MEHO2kGAHR3eg0AANBVyhq0rV69OsOHD8/MmTM3+vzFF1+c733ve5k1a1YeeOCBvOc978mYMWPy5ptvbvI1b7jhhkyZMiXTpk3Lww8/nOHDh2fMmDF54YUXOmsaAHRjeg0AANBVKkqlUqncRSRJRUVFbrvttowbNy7J+jMMhgwZkn/913/Nl770pSRJU1NT6urqMnv27Jx44okbfZ1Ro0blwx/+cL7//e8nSVpbW1NfX59/+Zd/yVlnndWhWpqbm1NbW5umpqbU1NRs/eQAtnHd5bjaXXpNd9kfAL2F4+qG7BOAYnX0uNptr9G2ZMmSNDY2pqGhoW1dbW1tRo0alQULFmx0zNq1a7Nw4cJ2YyorK9PQ0LDJMQBsu/QaAACgSH3LXcCmNDY2Jknq6urara+rq2t77i+99NJLaWlp2eiYJ554YpPvtWbNmqxZs6btcXNz85aWDUAP0lW9Rp8BAIBtQ7c9o60rTZ8+PbW1tW1LfX19uUsCoBfRZwAAYNvQbYO2QYMGJUlWrlzZbv3KlSvbnvtLO++8c/r06bNZY5Jk6tSpaWpqaluWL1++ldUD0BN0Va/RZwAAYNvQbYO23XffPYMGDcq8efPa1jU3N+eBBx7I6NGjNzqmX79+Oeigg9qNaW1tzbx58zY5JkmqqqpSU1PTbgGg9+uqXqPPAADAtqGs12h77bXXsnjx4rbHS5YsyaJFi7LTTjvl/e9/f84444xccMEF+cAHPpDdd98955xzToYMGdJ2t7gk+cQnPpG/+7u/yxe/+MUkyZQpUzJp0qQcfPDBOeSQQzJjxoysXr06kydP7urpAdAN6DUAAEBXKWvQ9tBDD+XII49sezxlypQkyaRJkzJ79ux85StfyerVq/PZz342r776ag477LDMmTMn1dXVbWOefvrpvPTSS22Px48fnxdffDHnnntuGhsbM2LEiMyZM2eDi1YDsG3QawAAgK5SUSqVSuUuortpbm5ObW1tmpqafL0HoACOq+3ZHwDFclzdkH0CUKyOHle77TXaAAAAeouZM2dmt912S3V1dUaNGpUHH3ywQ+Ouv/76VFRUtLukAQDdl6ANAACgE91www2ZMmVKpk2blocffjjDhw/PmDFj8sILL7zjuKVLl+ZLX/pSPvaxj3VRpQBsLUEbAABAJ7rsssty6qmnZvLkydl3330za9asbL/99rn66qs3OaalpSX/8A//kG984xvZY489urBaALaGoA0AAKCTrF27NgsXLkxDQ0PbusrKyjQ0NGTBggWbHPfNb34zAwcOzMknn9yh91mzZk2am5vbLQB0PUEbAABAJ3nppZfS0tKywZ2p6+rq0tjYuNEx9913X370ox/lyiuv7PD7TJ8+PbW1tW1LfX39VtUNwJYRtAEAAHQTq1atyoQJE3LllVdm55137vC4qVOnpqmpqW1Zvnx5J1YJwKb0LXcBAAAAvdXOO++cPn36ZOXKle3Wr1y5MoMGDdpg+6effjpLly7N2LFj29a1trYmSfr27Zsnn3wye+655wbjqqqqUlVVVXD1AGwuZ7QBAAB0kn79+uWggw7KvHnz2ta1trZm3rx5GT169AbbDxs2LI8++mgWLVrUtvzt3/5tjjzyyCxatMhXQgG6OWe0AQAAdKIpU6Zk0qRJOfjgg3PIIYdkxowZWb16dSZPnpwkmThxYoYOHZrp06enuro6+++/f7vxO+ywQ5JssB6A7kfQBgAA0InGjx+fF198Meeee24aGxszYsSIzJkzp+0GCcuWLUtlpS8bAfQGFaVSqVTuIrqb5ubm1NbWpqmpKTU1NeUuB6DHc1xtz/4AKJbj6obsE4BidfS46s8mAAAAAFAAQRsAAAAAFEDQBgAAAAAFELQBAAAAQAEEbQAAAABQAEEbAAAAABRA0AYAAAAABRC0AQAAAEABBG0AAAAAUABBGwAAAAAUQNAGAAAAAAUQtAEAAABAAQRtAAAAAFAAQRsAAAAAFEDQBgAAAAAFELQBAAAAQAEEbQAAAABQAEEbAAAAABRA0AYAAAAABRC0AQAAAEABBG0AAAAAUABBGwAAAAAUQNAGAAAAAAUQtAEAAABAAQRtAAAAAFAAQRsAAAAAFEDQBgAAAAAFELQBAAAAQAEEbQAAAABQAEEbAAAAABRA0AYAAAAABRC0AQAAAEABBG0AAAAAUABBGwAAAAAUQNAGAAAAAAUQtAEAAABAAQRtAAAAAFAAQRsAAAAAFEDQBgAAAAAFELQBAAAAQAEEbQAAAABQAEEbAAAAABSg2wdtq1atyhlnnJFdd901/fv3z6GHHprf/va3m9x+/vz5qaio2GBpbGzswqoB6En0GgAAoAh9y13AuznllFPy+9//Ptdee22GDBmS6667Lg0NDfnDH/6QoUOHbnLck08+mZqamrbHAwcO7IpyAeiB9BoAAKAI3fqMtjfeeCO33HJLLr744vz1X/919tprr5x33nnZa6+9csUVV7zj2IEDB2bQoEFtS2Vlt54qAGWi1wAAAEXp1p8I1q1bl5aWllRXV7db379//9x3333vOHbEiBEZPHhwjjrqqPzmN795x23XrFmT5ubmdgsA24au6DX6DAAAbBu6ddA2YMCAjB49Oueff36ef/75tLS05LrrrsuCBQuyYsWKjY4ZPHhwZs2alVtuuSW33HJL6uvrc8QRR+Thhx/e5PtMnz49tbW1bUt9fX1nTQmAbqYreo0+AwAA24aKUqlUKncR7+Tpp5/OZz7zmdx7773p06dPPvShD+WDH/xgFi5cmMcff7xDr3H44Yfn/e9/f6699tqNPr9mzZqsWbOm7XFzc3Pq6+vT1NTU7to7AGyZ5ubm1NbWdtvjamf3Gn0GoHN19z5TDvYJQLE6elzt1me0Jcmee+6ZX/3qV3nttdeyfPnyPPjgg3nrrbeyxx57dPg1DjnkkCxevHiTz1dVVaWmpqbdAsC2o7N7jT4DAADbhm4ftL3tPe95TwYPHpxXXnklc+fOzfHHH9/hsYsWLcrgwYM7sToAegO9BgAA2Bp9y13Au5k7d25KpVL23nvvLF68OF/+8pczbNiwTJ48OUkyderUPPfcc7nmmmuSJDNmzMjuu++e/fbbL2+++Wauuuqq/Od//md+8YtflHMaAHRjeg0AAFCEbh+0NTU1ZerUqfnTn/6UnXbaKZ/85CfzrW99K9ttt12SZMWKFVm2bFnb9mvXrs2//uu/5rnnnsv222+fAw88MPfcc0+OPPLIck0BgG5OrwEAAIrQ7W+GUA4uHApQLMfV9uwPgGI5rm7IPgEoVq+5GQIAAEBPN3PmzOy2226prq7OqFGj8uCDD25y2yuvvDIf+9jHsuOOO2bHHXdMQ0PDO24PQPchaAMAAOhEN9xwQ6ZMmZJp06bl4YcfzvDhwzNmzJi88MILG91+/vz5Oemkk/LLX/4yCxYsSH19fY4++ug899xzXVw5AJtL0AYAANCJLrvsspx66qmZPHly9t1338yaNSvbb799rr766o1u/x//8R/5whe+kBEjRmTYsGG56qqr0tramnnz5nVx5QBsLkEbAABAJ1m7dm0WLlyYhoaGtnWVlZVpaGjIggULOvQar7/+et56663stNNOm9xmzZo1aW5ubrcA0PUEbQAAAJ3kpZdeSktLS+rq6tqtr6urS2NjY4de46tf/WqGDBnSLqz7S9OnT09tbW3bUl9fv1V1A7BlBG0AAADd1Le//e1cf/31ue2221JdXb3J7aZOnZqmpqa2Zfny5V1YJQBv61vuAgAAAHqrnXfeOX369MnKlSvbrV+5cmUGDRr0jmO/853v5Nvf/nbuueeeHHjgge+4bVVVVaqqqra6XgC2jjPaAAAAOkm/fv1y0EEHtbuRwds3Nhg9evQmx1188cU5//zzM2fOnBx88MFdUSoABXBGGwAAQCeaMmVKJk2alIMPPjiHHHJIZsyYkdWrV2fy5MlJkokTJ2bo0KGZPn16kuSiiy7Kueeem5/+9KfZbbfd2q7l9t73vjfvfe97yzYPAN6doA0AAKATjR8/Pi+++GLOPffcNDY2ZsSIEZkzZ07bDRKWLVuWyso/f9noiiuuyNq1a/OpT32q3etMmzYt5513XleWDsBmErQBAAB0si9+8Yv54he/uNHn5s+f3+7x0qVLO78gADqFa7QBAAAAQAEEbQAAAABQAEEbAAAAABRA0AYAAAAABRC0AQAAAEABBG0AAAAAUABBGwAAAAAUQNAGAAAAAAUQtAEAAABAAQRtAAAAAFAAQRsAAAAAFEDQBgAAAAAFELQBAAAAQAEEbQAAAABQAEEbAAAAABRA0AYAAAAABRC0AQAAAEABBG0AAAAAUABBGwAAAAAUQNAGAAAAAAUQtAEAAABAAQRtAAAAAFAAQRsAAAAAFEDQBgAAAAAFELQBAAAAQAEEbQAAAABQAEEbAAAAABRA0AYAAAAABRC0AQAAAEABBG0AAAAAUABBGwAAAAAUQNAGAAAAAAUQtAEAAABAAQRtAAAAAFAAQRsAAAAAFEDQBgAAAAAFELQBAAAAQAEEbQAAAABQAEEbAAAAABRA0AYAAAAABRC0AQAAAEABBG0AAAAAUIBuH7StWrUqZ5xxRnbdddf0798/hx56aH7729++45j58+fnQx/6UKqqqrLXXntl9uzZXVMsAD2SXgMAABSh2wdtp5xySu6+++5ce+21efTRR3P00UenoaEhzz333Ea3X7JkSY477rgceeSRWbRoUc4444yccsopmTt3bhdXDkBPodcAAABFqCiVSqVyF7Epb7zxRgYMGJCf/exnOe6449rWH3TQQTnmmGNywQUXbDDmq1/9au688878/ve/b1t34okn5tVXX82cOXM69L7Nzc2pra1NU1NTampqtn4iANu47nxcLUev6c77A6AnclzdkH0CUKyOHle79Rlt69atS0tLS6qrq9ut79+/f+67776NjlmwYEEaGhrarRszZkwWLFiwyfdZs2ZNmpub2y0AbBu6otfoMwAAsG3o1kHbgAEDMnr06Jx//vl5/vnn09LSkuuuuy4LFizIihUrNjqmsbExdXV17dbV1dWlubk5b7zxxkbHTJ8+PbW1tW1LfX194XMBoHvqil6jzwAAwLahWwdtSXLttdemVCpl6NChqaqqyve+972cdNJJqawsrvSpU6emqampbVm+fHlhrw1A99fZvUafAQCAbUPfchfwbvbcc8/86le/yurVq9Pc3JzBgwdn/Pjx2WOPPTa6/aBBg7Jy5cp261auXJmampr0799/o2OqqqpSVVVVeO0A9Ayd3Wv0GQAA2DZ0+zPa3vae97wngwcPziuvvJK5c+fm+OOP3+h2o0ePzrx589qtu/vuuzN69OiuKBOAHkyvAQAAtka3D9rmzp2bOXPmZMmSJbn77rtz5JFHZtiwYZk8eXKS9V/HmThxYtv2n/vc5/LMM8/kK1/5Sp544on84Ac/yI033pgzzzyzXFMAoJvTawAAgCJ0+6Ctqakpp512WoYNG5aJEyfmsMMOy9y5c7PddtslSVasWJFly5a1bb/77rvnzjvvzN13353hw4fn0ksvzVVXXZUxY8aUawoAdHN6DQAAUISKUqlUKncR3U1zc3Nqa2vT1NSUmpqaDo9bsWpFZi+anT++/McM6DcgJ+x7Qg57/2GpqKjoxGoBur8tPa72Vlu6P95c92ZufOzG3PvsvalIRQ7f7fB8at9PpbpvdSdWC9D96TMb2pJ9Umptzb3/9/LcvOBHeW3dGxm24175p09fkrrd9+/kagG6v44eV7v9GW09xeUPXJ7679bn7F+enev++7pc8dAV+evZf52PX/PxNL3ZVO7yAOjhfvvcb/P+774/k26flJ/87ieZ/bvZmXDbhOw2Y7c8vOLhcpcHwLuYOXNmdtttt1RXV2fUqFF58MEH33H7m266KcOGDUt1dXUOOOCA3HXXXZ1a3ysrluSvp+yQIxadkVn9Hs112y/O19bOyV/9+IDMuuzTnfreAL2JoK0Atz5+a/73nP+dllJLWkutWde6Luta1yVJfv3srzP+5vFlrhCAnmzFqhU56tqj8vIbLydJuz7z0usvpeGahryw+oVylgjAO7jhhhsyZcqUTJs2LQ8//HCGDx+eMWPG5IUXNn7svv/++3PSSSfl5JNPziOPPJJx48Zl3Lhx+f3vf98p9ZVaW/PJiz6UBbWrkiTr+qxfWivX/+/nV/2f/Py6czvlvQF6G0HbViqVSvnGr76Rimz866EtpZbMfXpuFjUu6trCAOg1Zj00K6+tfS0tpZYNnmsptaRpTVOueviqMlQGQEdcdtllOfXUUzN58uTsu+++mTVrVrbffvtcffXVG93+3/7t3/I3f/M3+fKXv5x99tkn559/fj70oQ/l+9//fqfU99t7fpJf7vhqWjbx6bCyNblg4WWd8t4AvY2gbSs9t+q5/PfK/04pm77UXd+Kvrn9idu7rigAepWb/nDTRkO2t7WWWnPTYzd1YUUAdNTatWuzcOHCNDQ0tK2rrKxMQ0NDFixYsNExCxYsaLd9kowZM2aT2yfJmjVr0tzc3G7pqNt/fWX6brrNpLUyeWCH1Xlh6WMdfk2AbZWgbSu9/tbr77pNRUVF3njrjS6oBoDeaPVbqwvZBoCu99JLL6WlpSV1dXXt1tfV1aWxsXGjYxobGzdr+ySZPn16amtr25b6+voO1/j6ujc28f2c9t5Y/WqHXxNgWyVo20r1NfV5z3bvecdt3mp9K/sPdKceALbMyEEj07ey7yaf71vZNyMHjezCigDobqZOnZqmpqa2Zfny5R0ee8Dg4XnrXT4Z1q5JBu8xfCurBOj9BG1bqf92/XPyyJPTp6LPRp+vSEV2qN4hJ+x3QhdXBkBv8YUPf6Ht5gcbs651XT7/4c93YUUAdNTOO++cPn36ZOXKle3Wr1y5MoMGDdromEGDBm3W9klSVVWVmpqadktHnTjh4gxYm1Rs4mo4fVqTUysOTr/+7+3wawJsqwRtBfjmkd/MsJ2HbRC29anokz6VffIf/+s/Ut23ukzVAdDTHbXHUTntw6clSbub71T+vzY+5SNTcsRuR5SjNADeRb9+/XLQQQdl3rx5betaW1szb968jB49eqNjRo8e3W77JLn77rs3uf3Wes+OA3PdPl9LZSnp8xfXauvTmuzfXJ1zp/ysU94boLcRtBWgtro2v/nMb3LWYWdlp/47JUkqKyozdu+xuf8z9+fYDxxb5goB6MkqKipy+TGX58fH/zj77rJv2/r96/bPtX93bb5z9HfKWB0A72bKlCm58sor85Of/CSPP/54Pv/5z2f16tWZPHlykmTixImZOnVq2/ann3565syZk0svvTRPPPFEzjvvvDz00EP54he/2Gk1/u2Eb+W+0Vfm2FV1qWxdv27nNyrytXwsvz776Qx435BOe2+A3qSiVCpt+naZ26jm5ubU1tamqalps065Ttbf+a3pzaZsv932qepb1UkVAvQsW3Nc7Y22dn80r2lORSoyoGpAJ1QH0PP0hD7z/e9/P5dcckkaGxszYsSIfO9738uoUaOSJEcccUR22223zJ49u237m266KWeffXaWLl2aD3zgA7n44otz7LEd/wP+1uyTN197NW+seiW1A+tT2WfT1wgF2JZ09LgqaNuIntCoAXoSx9X27A+AYjmubsg+AShWR4+rvjoKAAAAAAUQtAEAAABAAQRtAAAAAFAAQRsAAAAAFEDQBgAAAAAFELQBAAAAQAEEbQAAAABQAEEbAAAAABRA0AYAAAAABRC0AQAAAEAB+pa7gO6oVColSZqbm8tcCUDv8Pbx9O3j67ZOnwEolj6zIb0GoFgd7TWCto1YtWpVkqS+vr7MlQD0LqtWrUptbW25yyg7fQagc+gzf6bXAHSOd+s1FSV/9tlAa2trnn/++QwYMCAVFRWbPb65uTn19fVZvnx5ampqOqHCbYP9WAz7sRj249YplUpZtWpVhgwZkspKVy3QZ96dOfZ8vX1+iTl2J/rMhram1/SUf/fuzn4shv1YDPtx63W01zijbSMqKyvzV3/1V1v9OjU1NX6AC2A/FsN+LIb9uOWcYfBn+kzHmWPP19vnl5hjd6HPtFdEr+kJ/+49gf1YDPuxGPbj1ulIr/HnHgAAAAAogKANAAAAAAogaOsEVVVVmTZtWqqqqspdSo9mPxbDfiyG/Uh3si38PJpjz9fb55eYI72Xf/di2I/FsB+LYT92HTdDAAAAAIACOKMNAAAAAAogaAMAAACAAgjaAAAAAKAAgjYAAAAAKICgbQvNnDkzu+22W6qrqzNq1Kg8+OCD77j9TTfdlGHDhqW6ujoHHHBA7rrrri6qtHvbnP04e/bsVFRUtFuqq6u7sNru5957783YsWMzZMiQVFRU5Pbbb3/XMfPnz8+HPvShVFVVZa+99srs2bM7vc7ubnP34/z58zf4WayoqEhjY2PXFMw2YVvoM5szxyuvvDIf+9jHsuOOO2bHHXdMQ0PDu+6T7mBz/x3fdv3116eioiLjxo3r3AK30ubO79VXX81pp52WwYMHp6qqKh/84Ae7/c/q5s5xxowZ2XvvvdO/f//U19fnzDPPzJtvvtlF1W4+v0tsu7aFPtMVfJ7Zeo5DxfCZpvsQtG2BG264IVOmTMm0adPy8MMPZ/jw4RkzZkxeeOGFjW5///3356STTsrJJ5+cRx55JOPGjcu4cePy+9//vosr7142dz8mSU1NTVasWNG2PPvss11YcfezevXqDB8+PDNnzuzQ9kuWLMlxxx2XI488MosWLcoZZ5yRU045JXPnzu3kSru3zd2Pb3vyySfb/TwOHDiwkypkW7Mt9JnNneP8+fNz0kkn5Ze//GUWLFiQ+vr6HH300Xnuuee6uPKO25I+lyRLly7Nl770pXzsYx/rokq3zObOb+3atTnqqKOydOnS3HzzzXnyySdz5ZVXZujQoV1cecdt7hx/+tOf5qyzzsq0adPy+OOP50c/+lFuuOGGfO1rX+viyjvO7xLbpm2hz3QFn2eK4ThUDJ9pupESm+2QQw4pnXbaaW2PW1paSkOGDClNnz59o9v//d//fem4445rt27UqFGlf/7nf+7UOru7zd2PP/7xj0u1tbVdVF3Pk6R02223veM2X/nKV0r77bdfu3Xjx48vjRkzphMr61k6sh9/+ctflpKUXnnllS6piW3PttBnNneOf2ndunWlAQMGlH7yk590VolbbUvmuG7dutKhhx5auuqqq0qTJk0qHX/88V1Q6ZbZ3PldccUVpT322KO0du3aripxq23uHE877bTSxz/+8XbrpkyZUvroRz/aqXUWxe8S245toc90BZ9niuc4VAyfacrLGW2bae3atVm4cGEaGhra1lVWVqahoSELFizY6JgFCxa02z5JxowZs8nttwVbsh+T5LXXXsuuu+6a+vr6HH/88Xnssce6otxew89isUaMGJHBgwfnqKOOym9+85tyl0MvsS30mS3tAf/T66+/nrfeeis77bRTZ5W5VbZ0jt/85jczcODAnHzyyV1R5hbbkvndcccdGT16dE477bTU1dVl//33z4UXXpiWlpauKnuzbMkcDz300CxcuLDtq2PPPPNM7rrrrhx77LFdUnNX6GnHGza0LfSZruDzTPn4eSyWzzTFE7RtppdeeiktLS2pq6trt76urm6T32VubGzcrO23BVuyH/fee+9cffXV+dnPfpbrrrsura2tOfTQQ/OnP/2pK0ruFTb1s9jc3Jw33nijTFX1PIMHD86sWbNyyy235JZbbkl9fX2OOOKIPPzww+UujV5gW+gzWzLHv/TVr341Q4YM2eAX7e5iS+Z433335Uc/+lGuvPLKrihxq2zJ/J555pncfPPNaWlpyV133ZVzzjknl156aS644IKuKHmzbckcP/3pT+eb3/xmDjvssGy33XbZc889c8QRR3Trr45uLr9L9HzbQp/pCj7PlI/jUDF8puk8fctdAHTU6NGjM3r06LbHhx56aPbZZ5/8+7//e84///wyVsa2Zu+9987ee+/d9vjQQw/N008/ne9+97u59tpry1gZbBu+/e1v5/rrr8/8+fN7zUWkV61alQkTJuTKK6/MzjvvXO5yOkVra2sGDhyYH/7wh+nTp08OOuigPPfcc7nkkksybdq0cpdXiPnz5+fCCy/MD37wg4waNSqLFy/O6aefnvPPPz/nnHNOucsDysznGboTn2k6j6BtM+28887p06dPVq5c2W79ypUrM2jQoI2OGTRo0GZtvy3Ykv34l7bbbruMHDkyixcv7owSe6VN/SzW1NSkf//+ZaqqdzjkkENy3333lbsMeoFtoc9sTQ/4zne+k29/+9u55557cuCBB3ZmmVtlc+f49NNPZ+nSpRk7dmzbutbW1iRJ37598+STT2bPPffs3KI3w5b8Gw4ePDjbbbdd+vTp07Zun332SWNjY9auXZt+/fp1as2ba0vmeM4552TChAk55ZRTkiQHHHBAVq9enc9+9rP5+te/nsrKnv9lEr9L9HzbQp/pCj7PlI/jUOfxmaYYPb/bd7F+/frloIMOyrx589rWtba2Zt68ee3+OvE/jR49ut32SXL33XdvcvttwZbsx7/U0tKSRx99NIMHD+6sMnsdP4udZ9GiRX4WKcS20Ge2tAdcfPHFOf/88zNnzpwcfPDBXVHqFtvcOQ4bNiyPPvpoFi1a1Lb87d/+bdsd1err67uy/He1Jf+GH/3oR7N48eK2ADFJ/vjHP2bw4MHdLmRLtmyOr7/++gZh2tvBYqlU6rxiu1BPO96woW2hz3QFn2fKx89j5/GZpiDlvhtDT3T99deXqqqqSrNnzy794Q9/KH32s58t7bDDDqXGxsZSqVQqTZgwoXTWWWe1bf+b3/ym1Ldv39J3vvOd0uOPP16aNm1aabvttis9+uij5ZpCt7C5+/Eb3/hGae7cuaWnn366tHDhwtKJJ55Yqq6uLj322GPlmkLZrVq1qvTII4+UHnnkkVKS0mWXXVZ65JFHSs8++2ypVCqVzjrrrNKECRPatn/mmWdK22+/fenLX/5y6fHHHy/NnDmz1KdPn9KcOXPKNYVuYXP343e/+93S7bffXnrqqadKjz76aOn0008vVVZWlu65555yTYFeZlvoM5s7x29/+9ulfv36lW6++ebSihUr2pZVq1aVawrvanPn+Je6+11HN3d+y5YtKw0YMKD0xS9+sfTkk0+Wfv7zn5cGDhxYuuCCC8o1hXe1uXOcNm1aacCAAaX/83/+T+mZZ54p/eIXvyjtueeepb//+78v1xTeld8ltk3bQp/pCj7PFMNxqBg+03QfgrYtdPnll5fe//73l/r161c65JBDSv/1X//V9tzhhx9emjRpUrvtb7zxxtIHP/jBUr9+/Ur77bdf6c477+ziirunzdmPZ5xxRtu2dXV1pWOPPbb08MMPl6Hq7uPtWzL/5fL2fps0aVLp8MMP32DMiBEjSv369SvtsccepR//+MddXnd3s7n78aKLLirtueeeperq6tJOO+1UOuKII0r/+Z//WZ7i6bW2hT6zOXPcddddN/r/02nTpnV94Zthc/8d/6fuHrSVSps/v/vvv780atSoUlVVVWmPPfYofetb3yqtW7eui6vePJszx7feeqt03nnntfWI+vr60he+8IXSK6+80vWFd5DfJbZd20Kf6Qo+z2w9x6Fi+EzTfVSUSr3kPHYAAAAAKCPXaAMAAACAAgjaAAAAAKAAgjYAAAAAKICgDQAAAAAKIGgDAAAAgAII2gAAAACgAII2AAAAACiAoA0AAAAACiBoAwAAAIACCNoAAAAAoACCNuhhXnzxxQwaNCgXXnhh27r7778//fr1y7x588pYGQC9xTXXXJP3ve99WbNmTbv148aNy4QJE8pUFQC9hc809GYVpVKpVO4igM1z1113Zdy4cbn//vuz9957Z8SIETn++ONz2WWXlbs0AHqBN954I4MHD86VV16ZE044IUnywgsvZOjQofnFL36RI488sswVAtDT+UxDbyVogx7qtNNOyz333JODDz44jz76aH7729+mqqqq3GUB0Et84QtfyNKlS3PXXXclSS677LLMnDkzixcvTkVFRZmrA6A38JmG3kjQBj3UG2+8kf333z/Lly/PwoULc8ABB5S7JAB6kUceeSQf/vCH8+yzz2bo0KE58MADc8IJJ+Scc84pd2kA9BI+09AbuUYb9FBPP/10nn/++bS2tmbp0qXlLgeAXmbkyJEZPnx4rrnmmixcuDCPPfZY/umf/qncZQHQi/hMQ2/kjDbogdauXZtDDjkkI0aMyN57750ZM2bk0UcfzcCBA8tdGgC9yBVXXJEZM2bkqKOOylNPPZW5c+eWuyQAegmfaeitBG3QA335y1/OzTffnN/97nd573vfm8MPPzy1tbX5+c9/Xu7SAOhFmpqaMmTIkKxbty7XXHNNxo8fX+6SAOglfKaht/LVUehh5s+fnxkzZuTaa69NTU1NKisrc+211+bXv/51rrjiinKXB0AvUltbm09+8pN573vfm3HjxpW7HAB6CZ9p6M2c0QYAwCZ94hOfyH777Zfvfe975S4FAKDbE7QBALCBV155JfPnz8+nPvWp/OEPf8jee+9d7pIAALq9vuUuAACA7mfkyJF55ZVXctFFFwnZAAA6yBltAAAAAFAAN0MAAAAAgAII2gAAAACgAII2AAAAACiAoA0AAAAACiBoAwAAAIACCNoAAAAAoACCNgAAAAAogKANAAAAAAogaAMAAACAAvz/e/gHozxT/KgAAAAASUVORK5CYII=",
-      "text/plain": [
-       "<Figure size 1500x500 with 3 Axes>"
-      ]
-     },
-     "metadata": {},
-     "output_type": "display_data"
-    }
-   ],
+   "outputs": [],
    "source": [
     "import matplotlib.pyplot as plt\n",
     "\n",