diff --git a/test.ipynb b/test.ipynb
index c267153..9dc5c96 100644
--- a/test.ipynb
+++ b/test.ipynb
@@ -2,23 +2,16 @@
  "cells": [
   {
    "cell_type": "code",
-   "execution_count": 1,
+   "execution_count": 41,
    "metadata": {},
    "outputs": [
-    {
-     "name": "stderr",
-     "output_type": "stream",
-     "text": [
-      "[Mac:27855] shmem: mmap: an error occurred while determining whether or not /var/folders/yh/dx7xl94n3g52ts3td8qcxjcc0000gn/T//ompi.Mac.501/jf.0/204800000/sm_segment.Mac.501.c350000.0 could be created.\n"
-     ]
-    },
     {
      "data": {
       "text/plain": [
        "'0.14.3'"
       ]
      },
-     "execution_count": 1,
+     "execution_count": 41,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -54,7 +47,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 2,
+   "execution_count": 42,
    "metadata": {},
    "outputs": [
     {
@@ -81,13 +74,13 @@
       "[0 0 1]  --»  [array([1, 0, 0]), array([0, 1, 0])]\n",
       "================================================================================================================================================================\n",
       "Parameters for the contour integral:\n",
-      "Number of k points:  20\n",
-      "k point directions:  xy\n",
+      "Number of k points:  1000\n",
+      "k point directions:  x\n",
       "Ebot:  -13\n",
-      "Eset:  100\n",
-      "Esetp:  100\n",
+      "Eset:  500\n",
+      "Esetp:  10000\n",
       "================================================================================================================================================================\n",
-      "Setup done. Elapsed time: 1.811041458 s\n",
+      "Setup done. Elapsed time: 5447.39860525 s\n",
       "================================================================================================================================================================\n"
      ]
     }
@@ -125,7 +118,7 @@
     "    # isotropic should be -82 meV\n",
     "    dict(ai=0, aj=1, Ruc=np.array([0, 0, 0])),\n",
     "    dict(ai=1, aj=2, Ruc=np.array([0, 0, 0])),\n",
-    "    dict(ai=0, aj=2, Ruc=np.array([0, 0, 0])),\n",
+    "    #    dict(ai=0, aj=2, Ruc=np.array([0, 0, 0])),\n",
     "]\n",
     "\n",
     "\"\"\"\n",
@@ -152,11 +145,11 @@
     "\"\"\"\n",
     "\n",
     "# Brilloun zone sampling and Green function contour integral\n",
-    "kset = 20\n",
-    "kdirs = \"xy\"\n",
+    "kset = 1000\n",
+    "kdirs = \"x\"\n",
     "ebot = -13\n",
-    "eset = 100\n",
-    "esetp = 100\n",
+    "eset = 500\n",
+    "esetp = 10000\n",
     "\n",
     "\n",
     "# MPI parameters\n",
@@ -197,33 +190,33 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 3,
+   "execution_count": 43,
    "metadata": {},
    "outputs": [
     {
-     "name": "stderr",
+     "name": "stdout",
      "output_type": "stream",
      "text": [
-      "/Users/danielpozsar/Documents/oktatás/elte/phd/grogu_project/.venv/lib/python3.9/site-packages/matplotlib/cbook.py:1762: ComplexWarning: Casting complex values to real discards the imaginary part\n",
-      "  return math.isfinite(val)\n",
-      "/Users/danielpozsar/Documents/oktatás/elte/phd/grogu_project/.venv/lib/python3.9/site-packages/matplotlib/cbook.py:1398: ComplexWarning: Casting complex values to real discards the imaginary part\n",
-      "  return np.asarray(x, float)\n"
+      "-12.806739\n",
+      "-0.01254111\n",
+      "xyz[-3:]: red, green, blue\n"
      ]
     },
     {
-     "name": "stdout",
+     "name": "stderr",
      "output_type": "stream",
      "text": [
-      "-12.806739\n",
-      "-0.01254111\n",
-      "xyz[-3:]: red, green, blue\n"
+      "/Users/danielpozsar/Documents/oktatás/elte/phd/grogu_project/.venv/lib/python3.9/site-packages/matplotlib/cbook.py:1762: ComplexWarning: Casting complex values to real discards the imaginary part\n",
+      "  return math.isfinite(val)\n",
+      "/Users/danielpozsar/Documents/oktatás/elte/phd/grogu_project/.venv/lib/python3.9/site-packages/matplotlib/cbook.py:1398: ComplexWarning: Casting complex values to real discards the imaginary part\n",
+      "  return np.asarray(x, float)\n"
      ]
     },
     {
      "data": {
-      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjcAAAGiCAYAAAAba+fDAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA3T0lEQVR4nO3deXhU9aH/8c9km5BtQvYEEkgAQfZFjHErKrJorSi1bq1g3YWqoK3SX0XBalx6tddeqn1uRbRWqfa6VNuqiCxawhaIC0sEDCRAFkhIJguZJDPn9wdmdCRAApk5M5P363nmeZgzZ5JPjieZj9/zPedYDMMwBAAAECRCzA4AAADQnSg3AAAgqFBuAABAUKHcAACAoEK5AQAAQYVyAwAAggrlBgAABBXKDQAACCqUGwAAEFQoNwAAIKh4tdzk5+dr/Pjxio2NVUpKiqZNm6bi4mKPdZqbmzVr1iwlJiYqJiZG06dPV2VlpTdjAQCAIObVcrNq1SrNmjVLa9eu1bJly9Ta2qpJkyapsbHRvc6cOXP07rvv6o033tCqVau0f/9+XXnlld6MBQAAgpjFlzfOPHDggFJSUrRq1Sqdf/75qqurU3Jysl599VX9+Mc/liRt375dp59+ugoKCnTWWWf5KhoAAAgSYb78ZnV1dZKkhIQESVJhYaFaW1s1ceJE9zpDhgxRVlbWMcuNw+GQw+FwP3e5XKqpqVFiYqIsFouXfwIAANAdDMNQfX29MjIyFBLSvQeSfFZuXC6X7rnnHp1zzjkaPny4JKmiokIRERGKj4/3WDc1NVUVFRUdfp38/HwtWLDA23EBAIAPlJWVqW/fvt36NX1WbmbNmqUvv/xSn3766Sl9nXnz5mnu3Lnu53V1dcrKylJZWZni4uJONSYAdAuXy1De48vV6HDqrTvP1qDUWLMjAX7FbrcrMzNTsbHd/7vhk3Ize/Zsvffee1q9erVHO0tLS1NLS4tqa2s9Rm8qKyuVlpbW4deyWq2yWq1HLY+Li6PcAPAbpdVNOiyrIqNCNDInXeGhXHkD6Ig3ppR49bfNMAzNnj1bb731lj7++GNlZ2d7vD5u3DiFh4dr+fLl7mXFxcUqLS1VXl6eN6MBgFdtLbdLkgalxlBsAB/z6sjNrFmz9Oqrr+qdd95RbGysex6NzWZTr169ZLPZdNNNN2nu3LlKSEhQXFycfvGLXygvL48zpQAEtG3flJvT0xlRBnzNq+XmueeekyRNmDDBY/mLL76omTNnSpKeeeYZhYSEaPr06XI4HJo8ebL++Mc/ejMWAHiVo82p9SU1kig3gBl8ep0bb7Db7bLZbKqrq2PODQBTtTpd+tcX5frdh8UqqzksSfrH7HM0sm+8ucEAP+TNz2+fXucGAIJNq9Ol9z7frw+3VOrTHQdV72iTJKXEWvXA1CEUG8AElBsAOAkul6F3P9+vZ5Z9pd3VTe7lidERmnF2f918XraiIvgTC5iB3zwAOAm3v1KoD7ceuclvYnSErj+rny4ckqKRfWwKCeFq6YCZKDcA0EWtTpc+3l4lSbr7okG69fwcRVv5cwr4C34bAaCL9tceVpvLUERYiO6+aBAjNYCf4cpSANBFe76ZY5OVEEWxAfwQ5QYAumhPzZFy0z8xyuQkADpCuQGALiqtbpQkZSVEm5wEQEcoNwDQRe2HpfoxcgP4JcoNAHRR6TeHpbIoN4BfotwAQBcYhvHtyE0C5QbwR5QbAOiCA/UOHW51KsQi9e1NuQH8EeUGALqg/UypjPheigjjTyjgj/jNBIAuYDIx4P8oNwDQBZwGDvg/yg0AdEH7YSlGbgD/RbkBgC7gTCnA/1FuAKAL9nxzWKpfIoelAH9FuQGATrI3t+pQU6skLuAH+DPKDQB0Uuk3h6SSYiIUYw0zOQ2AY6HcAEAntc+3yWK+DeDXKDcA0El7aphvAwQCyg0AdEJVfbPe2rRPEqeBA/6Og8YAcAJlNU366QvrtKe6ScmxVv14XF+zIwE4DsoNABzH3kNNuur5AlXYm5WZ0Et/veksbpgJ+DnKDQAcQ5vTpXuWFqnC3qxBKTF65eZcpcZFmh0LwAkw5wYAjmHRil3auOeQYqxhWjxzPMUGCBCUGwDoQOGeQ3r24x2SpEemDVMmp38DAYNyAwDfU1HXrLte2yyny9DlozN0xRgmEAOBhHIDAN9R3eDQ9X9eq321h9U/MUqPTBtudiQAXUS5AYBv1B1u1c9eWK9dBxqVbovUKzfnKi4y3OxYALqIcgMAkgzD0F2vbdbWcruSYiL015tzOeUbCFCUGwCQtHRDmVZ9dUDWsBC9/PNc5STHmB0JwEmi3ADo8fYeatJv39sqSfrl5MEamhFnciIAp4JyA6BHc7kM/ervn6uxxakz+vXWjedkmx0JwCmi3ADo0f5euFdrdlUrMjxET101SqEhFrMjAThFlBsAPVaDo01PflAsSZp78WnKToo2ORGA7kC5AdBjPbdypw42ONQ/MUozz+ZwFBAsvFpuVq9ercsuu0wZGRmyWCx6++23PV6fOXOmLBaLx2PKlCnejAQAko5MIv7fT0okSfMuOV0RYfy/HhAsvPrb3NjYqFGjRmnRokXHXGfKlCkqLy93P1577TVvRgIASdKT7xerpc2l3OwETRqaanYcAN0ozJtffOrUqZo6depx17FarUpLS/NmDADwsOtAg/7x2X5ZLNKDPxwqi4VJxEAwMX0cduXKlUpJSdHgwYN1xx13qLq6+rjrOxwO2e12jwcAdMUbG/dKki4cnKLhfWwmpwHQ3UwtN1OmTNHLL7+s5cuX64knntCqVas0depUOZ3OY74nPz9fNpvN/cjMzPRhYgCBrs3p0pubjpSbq87gbt9AMPLqYakTueaaa9z/HjFihEaOHKkBAwZo5cqVuuiiizp8z7x58zR37lz3c7vdTsEB0Gmf7DioqnqHEqIjdOEQ5toAwcj0w1LflZOTo6SkJO3cufOY61itVsXFxXk8AKCz3igskyRdPjqDM6SAIOVXv9l79+5VdXW10tPTzY4CIAjVNLZo2dZKSdJV4xjxBYKVVw9LNTQ0eIzClJSUqKioSAkJCUpISNCCBQs0ffp0paWladeuXfrVr36lgQMHavLkyd6MBaCH+r/CvWp1GhreJ46bYwJBzKvlZuPGjbrgggvcz9vnysyYMUPPPfecPv/8c7300kuqra1VRkaGJk2apEceeURWq9WbsQAEuQZHmz7dcVBV9c1qbnWqyu7Qyq8OaGdVgyRGbYBg59VyM2HCBBmGcczXP/jgA29+ewBBrL65VfXNbTIkHW5xqrSmUSUHm1Swq1qrdxxQS5vrqPeEhVg0eViafnIG5QYIZqaeLQUAXVVlb9Z/L9+hpRvK5HQd+3+e+idGaUhanKIiQhVtDdOZ2Qn6weBkxUWG+zAtADNQbgD4tS/31emVtXvU0uZSi9Ol5duqdLj1yLWwIkJDZLFI4aEhykyIUnZSlE5Pi9OkYWk6LTWGKw8DPRTlBoBfW/DuFm3Yfchj2ZiseD0wZYhycxJNSgXAn1FuAPitwy1OFZXVSpLmTDxN0dZQDUiJ0YTTkhmVAXBMlBsAfmtT6SG1Og1l2CJ110UDKTQAOsWvLuIHAN+17usjN9LNzUmk2ADoNMoNAL+1tqRGkpSbnWByEgCBhHIDwC81t3473+ZMyg2ALqDcAPBLRWW1amlzKTnWquykaLPjAAgglBsAfmnd198ekmK+DYCuoNwA8EvrSr6dTAwAXUG5AeB3Wtpc2lR65MJ9ZzHfBkAXUW4A+J2Ne2rU3OpSYnSEBqbEmB0HQICh3ADwK+9/Wa7b/lIoSTpnYBLzbQB0GVcoBmCaspomrSyuUnFlvVrbDFU3OvTRtipJ0ujMeP36ktNNTgggEFFuAPhUS5tLf9tQqiVrdmvXgcYO17ntBzm6b9JghYcyuAyg6yg3AHziQL1DH26t0B9X7NK+2sOSpNAQi8b1663x/XurV3iowkJDNL5/b43rxyRiACePcgPAqz7YUqE/fLxDX+6zu5elxFo1+8KBmjamj+Iiw01MByAYUW4AeE1tU4t+8epmtThdkqThfeI0bXQfXZ/bT70iQk1OByBYUW4AeM07RfvV4nRpSFqs/nJTrpJjrWZHAtADMFsPgNe8vrFMknTN+EyKDQCfodwA8Iot++u0Zb9dEaEhunx0H7PjAOhBKDcAvOKNjXslSRcPTVXv6AiT0wDoSSg3ALqdo82pt4v2SZKuOqOvyWkA9DSUGwDd7t9fVKi2qVVpcZE6b1Cy2XEA9DCcLQWg27hchhb/p0RPvL9d0pFRm9AQ7g0FwLcoNwC6RXOrU3e8UqgVxQckSZOHpeqOCQNMTgWgJ6LcAOgWS9eXakXxAUWEhWj+D4fq+tws7ugNwBSUGwCnzDAM/e2bs6PmTR2in57Vz+REAHoyJhQDOGVb9tu1rfzINW2uGMM1bQCYi3ID4JS1X4l40rBUxUdxTRsA5qLcADglza1Ovb35yDVtfnJGpslpAIByA+AUfbClQvbmNvWJ76VzBiaZHQcAmFAMoOsOtzhVuOeQSmua9Je1eyRJ08dxTRsA/oFyA6DTGh1t+svaPfrzJ1/rYEOLe3mIRbpqHLdZAOAfKDcAjuuDLRVa/dUB7TrQoC377apvbpMkpcVFamhGnDJ799J5g5KVmRBlclIAOIJyA+CYKu3Nuv2VQhnGt8tykqI164KBunx0hsJCmbYHwP9QbgAcU8nBRhmGlBRj1a8vGaKBKTEalmFjbg0Av+bV/+1avXq1LrvsMmVkZMhisejtt9/2eN0wDM2fP1/p6enq1auXJk6cqB07dngzEoAu2HfosCRpcFqMrhzbVyP7xlNsAPg9r5abxsZGjRo1SosWLerw9SeffFLPPvusnn/+ea1bt07R0dGaPHmympubvRkLQCft/abc9I1nPg2AwOHVw1JTp07V1KlTO3zNMAz9/ve/129+8xtdfvnlkqSXX35Zqampevvtt3XNNdd4MxqATthX2yRJ6tO7l8lJAKDzTJsNWFJSooqKCk2cONG9zGazKTc3VwUFBcd8n8PhkN1u93gA8I59tUdGbvrEU24ABA7Tyk1FRYUkKTU11WN5amqq+7WO5Ofny2azuR+ZmVzuHfAW92EpRm4ABJCAO49z3rx5qqurcz/KysrMjgQEJZfLUHntkflvHJYCEEhMKzdpaWmSpMrKSo/llZWV7tc6YrVaFRcX5/EA0P0ONDjU4nQpNMSitLhIs+MAQKeZVm6ys7OVlpam5cuXu5fZ7XatW7dOeXl5ZsUC8I32Q1JpcZFcrA9AQPHq2VINDQ3auXOn+3lJSYmKioqUkJCgrKws3XPPPfrtb3+rQYMGKTs7Ww8++KAyMjI0bdo0b8YC0Al7D3GmFIDA5NVys3HjRl1wwQXu53PnzpUkzZgxQ0uWLNGvfvUrNTY26tZbb1Vtba3OPfdcvf/++4qMZAgcMFv7mVJ9OVMKQIDxarmZMGGCjO/elOZ7LBaLFi5cqIULF3ozBoCT0H51YkZuAAQaDqQD6BDXuAEQqCg3ADr07TVuuPUCgMBCuQFwFMMwOCwFIGBRbgAc5VBTqw63OiVJ6TYm+AMILJQbAEdpH7VJjrUqMjzU5DQA0DWUGwBHab/GDfeUAhCIKDcAjsKZUgACGeUGwFH2MpkYQACj3AA4ivuwFCM3AAKQV69QDCAw1Da1aE91kz7fV6e/F+7VZ2W1kqS+CVzjBkDgodwAPZRhGHqnaL9+92Gx+zBUu7AQi6aOSFdeTqJJ6QDg5FFugB6ktqlFew8dVmlNk174tESFew65X0uJtap/UrQmDU3VtDF9lBRjNTEpAJw8yg0QxFqdLv35kxJt2F2jL/fVqare4fF6r/BQzb5woGac3V8xVv4cAAgO/DUDgthbm/fpife3eyxLirGqT+9eGp4Rp9kXDlS6jUnDAIIL5QYIYu8U7ZMkXTmmj64/K0tD0uIUzQgNgCDHXzkgSFXam7VmV7Ukac7FpymTM58A9BBc5wYIUu9+tl+GIY3r15tiA6BHodwAQeqdov2SpMtHZ5icBAB8i3IDBKFdBxr0xb46hYZYdOmIdLPjAIBPUW6AIPSPb0ZtzhuUpESuVwOgh2FCMRAEDMPQ3kOHtan0kDbtOaR/fHak3Ewb3cfkZADge5QbIIBtLj2k//3ka60vOaSDDZ4X6EuIjtDFQ1NNSgYA5qHcAAGorKZJT35QrHe/GaGRpPBQi4Zm2DQ2K15js3rr7AGJXNMGQI/EXz4gQLhchv6z66BeWbtHH22rktNlyGKRpo/tq2vGZ2p4H5siw0PNjgkApqPcAAFgy/46/frNL/TZ3jr3snMHJmneJUM0LMNmYjIA8D+UG8CPNbc69cxHX+nPn5TI6TIUHRGq6eP66vrcfhqcFmt2PADwS5QbwE+1OV2a/eomfbStSpJ06Yh0PXTZUKXERZqcDAD8G+UG8EOGYWj+P7boo21VsoaF6Nlrx2jysDSzYwFAQKDcAH6mzenSohW79Oq6Ulks0n9fQ7EBgK6g3AB+wOUy9MnOg/pH0X4t316p2qZWSdLDlw3TlOEUGwDoCsoNYKL65la993m5Xvi0RDurGtzLe0eF6/YfDNCMs/ubFw4AAhTlBvAxp8vQ3wvL9M8vKlSw66BanYYkKcYapulj++iSEeka16+3wkK59RsAnAzKDeBDza1O3b10sz7YUulelpMUretys3T1+EzFRoabmA4AggPlBvCR2qYW3fzSRm3cc0gRoSGafeFAXTIiXQNTYsyOBgBBhXIDeFl53WG9vmGvXltfqgp7s2Ijw/S/N5yhs3ISzY4GAEGJcgN40csFu/XwP7bIdWRajfrE99LimeO5ujAAeBHlBvCSRkebnvqgWC5DOrN/gq7LzdKU4Wnc3BIAvIxyA3jJW5v3qb65TdlJ0Vp661kKCbGYHQkAegTONQW8wDAMvVywW5L0s7P6UWwAwIdMLzcPP/ywLBaLx2PIkCFmxwJOScHX1fqqskFREaH68Rl9zY4DAD2KXxyWGjZsmD766CP387Awv4gFnLSX1uyWJF05to/iuHYNAPiUX7SIsLAwpaV17v45DodDDofD/dxut3srFtBpza1ObdpzSKU1Taqqd2jZ1iMX6ZuR19/cYADQA/lFudmxY4cyMjIUGRmpvLw85efnKysrq8N18/PztWDBAh8nBDwdamxRUVmtNpfVakNJjQpLD6mlzeWxztkDEjUolVO+AcDXLIZhGGYG+Pe//62GhgYNHjxY5eXlWrBggfbt26cvv/xSsbFHfzB0NHKTmZmpuro6xcXF+TI6eqgPt1Tojr9uktPl+auTGmfVsAybkmIilBIbqavHZyozIcqklADg3+x2u2w2m1c+v00fuZk6dar73yNHjlRubq769eun119/XTfddNNR61utVlmtVl9GBDys3nFATpeh1DirzhmQpDFZ8cobkKQBydGyWDgrCgDMZnq5+b74+Hiddtpp2rlzp9lRgA5V1B0ZOfzFhYP007P6mZwGAPB9pp8K/n0NDQ3atWuX0tPTzY4CdKjS3ixJSouLNDkJAKAjppeb++67T6tWrdLu3bu1Zs0aXXHFFQoNDdW1115rdjSgQxXt5cZGuQEAf2T6Yam9e/fq2muvVXV1tZKTk3Xuuedq7dq1Sk5ONjsacJRWp0sHG44clkpl5AYA/JLp5Wbp0qVmRwA6rareIcOQwkMtSoyOMDsOAKADph+WAgJJRd2RQ1IpsZHcLwoA/BTlBuiCSubbAIDfo9wAXdA+csOZUgDgvyg3QBe0j9wwmRgA/BflBuiC8vaRGxtXyQYAf0W5AbqggpEbAPB7lBugC9oPS6XbepmcBABwLJQboJMMw2BCMQAEAMoN0El1h1vlaHNJklLimHMDAP6KcgN0Uvt8m95R4YoMDzU5DQDgWCg3QCe1nynFZGIA8G+UG6CTKuu4OjEABALKDdBJFe4zpSg3AODPKDdAJ3F1YgAIDJQboJM4DRwAAgPlBuikCrtDkpTKYSkA8GuUG6CTKuoOS2LkBgD8HeUG6ITmVqcONbVKotwAgL8LMzsA4M8cbU4VldZqc1mtJMkaFqL4qHBzQwEAjotyAxzDocYWXfWnAu2sanAv658YLYvFYmIqAMCJUG6ADrS0uXT7K4XaWdUgW69wjc2K18CUGF0xpq/Z0QAAJ0C5Ab7HMAz95u0vtK6kRjHWML1+W54Gp8WaHQsA0EmUG/R4DY42ldce1t5Dh7VxT43+s7NaRWW1CrFIf7huDMUGAAIM5QY9UpvTpX98tl/PrdylHd+ZU9POYpEW/GiYLhicYkI6AMCpoNygxyiuqNdne2u1o7JeH2ypVGlNk/u12MgwZdh6aXgfm3JzEpSXk6jMhCgT0wIAThblBj3CH1fu1JPvF3ssS4iO0M3nZeu6M7MUHxVhUjIAQHej3CDovVO0z11scrMTdHp6nIZmxOmHI9MVFcGvAAAEG/6yI6ht2F2jX77xuSTppnOz9eAPh5qcCADgbdx+AUHr/S8rdNOSDWpxujR5WKp+fcnpZkcCAPgAIzcIOo2ONv32n1v12voySdIZ/Xrr91ePUWgIVxYGgJ6AcoOA5XIZsje3qrapVYeaWlRa06QPt1ZqxfYqNbU4ZbFIt50/QHMvPk0RYQxSAkBPQblBQFj11QHl/2ubGlva5HQaamp1qu5wqwyj4/X7JUYp/8oROntAkm+DAgBMR7lBQHh+5S5tr6jv8LXoiFDFR0UoITpCZw9M1CXD0zWyr40bXAJAD0W5gd9ramnTxj01kqQ//Wyc0m2RigwPVXxUuGy9wmUNCzU5IQDAn1Bu4PfWldSo1WmoT3wvTRqayogMAOC4mGUJv/fpjoOSpHMHJlFsAAAnRLmB32svN+edxuRgAMCJUW7g16rszSqurJfFIp3DmU8AgE7wi3KzaNEi9e/fX5GRkcrNzdX69evNjgQ/8ck3ozbDM2zqHc3NLQEAJ2Z6ufnb3/6muXPn6qGHHtKmTZs0atQoTZ48WVVVVWZHgx/4dOc3820GMWoDAOgc08vN008/rVtuuUU33nijhg4dqueff15RUVFavHix2dFgMsMw3OXmvIGUGwBA55hablpaWlRYWKiJEye6l4WEhGjixIkqKCjo8D0Oh0N2u93jgeC0cc8hHah3KDI8ROP69zY7DgAgQJhabg4ePCin06nU1FSP5ampqaqoqOjwPfn5+bLZbO5HZmamL6LCx9aX1OjnSzZIki4cksKF+gAAnWb6Yamumjdvnurq6tyPsrIysyOhm73/ZYV+9sI61Te3aXz/3sq/YqTZkQAAAcTUKxQnJSUpNDRUlZWVHssrKyuVlpbW4XusVqusVqsv4sHHquqb9eg/t+mdov2SpImnp+p/rhujyHBGbQAAnWdquYmIiNC4ceO0fPlyTZs2TZLkcrm0fPlyzZ4928xo8JED9Q5t2F2jdV9X683N+1Tf3CaLRfr5OdmaN3WIwkIDbnARAGAy0+8tNXfuXM2YMUNnnHGGzjzzTP3+979XY2OjbrzxRrOjoZu1tLm0o6pe28rrVbjnkNaVVOvrA40e64zoY9OjVwzXyL7x5oQEAAQ808vN1VdfrQMHDmj+/PmqqKjQ6NGj9f777x81yRiBq6q+WY+8t03//qJcbS7jqNeHpMXqzOwEnT0gURcPTVNoCPePAgCcPIthGEd/2gQQu90um82muro6xcXFmR2nxzMMQ9WNLWp0tKnV6VLhnkN67F/bVXe4VZIUFxmmIelxGtnHptycRI3v31vxUVx5GAB6Gm9+fps+coPA19zq1O8+KNby7VXaV3tYLW2uo9YZ3idOj10xQiP62LizNwDAqyg3OCW7Dzbqzr9u0tZyz4spRkWEKiIsRDHWMP30rH66+dxsJgcDAHyCcoOTtuqrA5r1101qcLQpITpCCy8fplF945UaF6mIMIoMAMAclBuclKKyWt3+l0IdbnXqzP4JevbaMUqzRZodCwAAyg26bvfBRv18yQYdbnXq/NOS9cKMMxTOIScAgJ/gEwldUrinRjNeXK+axhaN6GPTc9ePpdgAAPwKIzc4oTanSxt2H9Lzq3Zp1VcHJEmZCb20eOZ4RVvZhQAA/oVPJnSopc2lj7dX6p9fVGj1Vwfc16kJDbHoqnF9Nefi05Qcyz2+AAD+h3IDD4caW/SHj3fqrc17daip1b08PipcU4al6Y4JA9QvMdrEhAAAHB/lBm4bdtfortc2q7yuWZKUEmvVFWP6aNKwVI3O7M1tEQAAAYFyAxmGoedXfa3ffVgsp8tQTlK0fvPD03X+oGQuvAcACDiUmx6u0dGmX/79M/3riwpJ0hVj+ui304YzURgAELD4BOvBiivqdffSzdpeUa/wUIsWXj5c14zP5N5PAICARrnpARxtTn3y1UFV2JvlaHPpQL1Dy7dVakdVgyQpKcaqP/1srMb1SzA5KQAAp45yE8T21x7Wi/8p0f9t2qeaxpajXg8PtegHp6XokWnDlG7rZUJCAAC6H+UmSDU62nTlH9eown7kzKfUOKtG9Y1XZHiooq2hys1O1AVDUmTrFW5yUgAAuhflJkj9+ZMSVdib1Se+lxb8aJgmDObMJwBAz0C5CUIH6h360+pdkqRfX3K6Jg5NNTkRAAC+w//KB6H/Xv6VmlqcGpUZr0tGpJkdBwAAn2LkJkg0tbRp76HD2lXVoNfWl0mS5k0dwmndAIAeh3ITBLaV23XlH9focKvTvWzi6Sk6KyfRxFQAAJiDchMEXi7YrcOtTkVHhCozIUpZCVF68IdDzY4FAIApKDcBrqmlTe9+Vi5JemHmeEZrAAA9HhOKA9y/v6hQg6NNWQlRys3mCsMAAFBuAtwbhUcmD181ri+ThwEAEOUmoO2pbtTar2tksUjTx/U1Ow4AAH6BchPA/l64V5J03qBkZcRzbygAACQmFAcUp8vQ3wvL9M8vKrS/9rD2VDdKOnJICgAAHEG5CQCGYeiTHQf12L+2aXtFvcdr2UnRupjbKwAA4Ea58VMtbS5tLbfrwy0V+ucX5dpT3SRJiosM020/GKDRmfFKt0UqMyFK4dwQEwAAN8qNn1mxvUq/X75D2/bb1eJ0uZdHhofomvFZuvuiQeodHWFiQgAA/Bvlxo/UNbVq9qub1Nhy5DYKtl7hystJ1KUj03XhkBRFW/nPBQDAifBp6Uf+sna3GlucOi01Rv97wxnKSoji2jUAAHQR5cZPNLc69eJ/dkuS7pwwUP0So80NBABAgGImqp94Y2OZqhtb1Ld3L/1wZLrZcQAACFiUGz/Q5nTpT6u/liTden6Owjj7CQCAk8anqB9Ysma39h46rMToCF01LtPsOAAABDTm3JjIMAwtWrFTv/vwK0nSbT/IUa+IUJNTAQAQ2Ewduenfv78sFovH4/HHHzczktfZm1v1/pcVenVdqea+/pm72My6YIBuOS/H5HQAAAQ+00duFi5cqFtuucX9PDY21sQ03nWosUWXPvuJ9tc1u5dZLNJDPxyqmedkm5gMAIDgYXq5iY2NVVpamtkxvM4wDP2/t7/Q/rpmJcdaNaqvTQnREfrhyAydf1qy2fEAAAgaFsMwDLO+ef/+/dXc3KzW1lZlZWXpuuuu05w5cxQWduzO5XA45HA43M/tdrsyMzNVV1enuLg4X8Q+KW9u2qu5r3+msBCL3rrzHI3oazM7EgAAprHb7bLZbF75/DZ15Oauu+7S2LFjlZCQoDVr1mjevHkqLy/X008/fcz35Ofna8GCBT5MeerKapr00DtbJEn3TBxEsQEAwIu6feTmgQce0BNPPHHcdbZt26YhQ4YctXzx4sW67bbb1NDQIKvV2uF7A23kZu3X1brrtc2qqndoXL/e+tutZ3EdGwBAj+fNkZtuLzcHDhxQdXX1cdfJyclRRMTRd7besmWLhg8fru3bt2vw4MGd+n7e3DinatGKnfqvD4vlMqSBKTFacuN49e0dZXYsAABMF1CHpZKTk5WcfHITZIuKihQSEqKUlJRuTuV7G3fX6KkPiiVJV47to99OG66oCNPnbwMAEPRM+7QtKCjQunXrdMEFFyg2NlYFBQWaM2eOfvrTn6p3795mxeo2u6ubJEl5OYl6+iejzQ0DAEAPYlq5sVqtWrp0qR5++GE5HA5lZ2drzpw5mjt3rlmRulVDc6skKSHm6MNvAADAe0wrN2PHjtXatWvN+vZe1+BokyTFWjkUBQCAL3HajpfUf1NuYig3AAD4FOXGSxqavyk3kZQbAAB8iXLjJY2M3AAAYArKjZe459wwcgMAgE9Rbrykvv2wlDXc5CQAAPQslBsvaR+5Yc4NAAC+Rbnxkgbm3AAAYArKjZe4z5ai3AAA4FOUGy+p57AUAACmoNx4gaPNqZY2lyRGbgAA8DXKjRc0Opzuf1NuAADwLcqNF7TPt4mKCFVoiMXkNAAA9CyUGy+odxy5IzijNgAA+B7lxgvaD0sxmRgAAN+j3HhBwzcjN7GM3AAA4HOUGy+o547gAACYhnLjBVydGAAA81BuvKD9bKloyg0AAD5HufGC9pEb5twAAOB7lBsvYM4NAADmodx4wbdzbsJNTgIAQM9DufGCBkZuAAAwDeXGC5hzAwCAeSg3XsCp4AAAmIdy4wXucsNhKQAAfI5y4wXuOTeM3AAA4HOUGy/gsBQAAOah3HQzl8vgsBQAACai3HSzxpY2978ZuQEAwPcoN92sfdQmPNQiaxibFwAAX+PTt5t9dzKxxWIxOQ0AAD0P5aab1TPfBgAAU1Fuutm3IzfcVwoAADNQbroZt14AAMBclJtuxmngAACYi3LTzbg6MQAA5qLcdLP2kZtoyg0AAKag3HQz95wbDksBAGAKr5WbRx99VGeffbaioqIUHx/f4TqlpaW69NJLFRUVpZSUFP3yl79UW1tbh+sGinoOSwEAYCqvfQK3tLToqquuUl5enl544YWjXnc6nbr00kuVlpamNWvWqLy8XDfccIPCw8P12GOPeSuW13HTTAAAzOW1kZsFCxZozpw5GjFiRIevf/jhh9q6dateeeUVjR49WlOnTtUjjzyiRYsWqaWlxVuxvK6huVUSZ0sBAGAW0+bcFBQUaMSIEUpNTXUvmzx5sux2u7Zs2XLM9zkcDtntdo+HP+E6NwAAmMu0clNRUeFRbCS5n1dUVBzzffn5+bLZbO5HZmamV3N2lXvODSM3AACYokvl5oEHHpDFYjnuY/v27d7KKkmaN2+e6urq3I+ysjKvfr+uYs4NAADm6tIn8L333quZM2ced52cnJxOfa20tDStX7/eY1llZaX7tWOxWq2yWq2d+h5m4FRwAADM1aVP4OTkZCUnJ3fLN87Ly9Ojjz6qqqoqpaSkSJKWLVumuLg4DR06tFu+hy852pxaX1LjvkIxF/EDAMAcXvsELi0tVU1NjUpLS+V0OlVUVCRJGjhwoGJiYjRp0iQNHTpUP/vZz/Tkk0+qoqJCv/nNbzRr1iy/HpnpyOsbyvTwu1vU1OKUJEWEhah3VITJqQAA6Jm8Vm7mz5+vl156yf18zJgxkqQVK1ZowoQJCg0N1Xvvvac77rhDeXl5io6O1owZM7Rw4UJvRfKal9fuVlOLU8mxVl0wOFk/HpepyPBQs2MBANAjWQzDMMwOcSrsdrtsNpvq6uoUFxfn8+/vchka+tD7am51acV9E5SdFO3zDAAABBpvfn5zb6lTVHaoSc2tLlnDQpSVEGV2HAAAejzKzSkqrqiXJA1IjlFoiMXkNAAAgHJzinZUNUiSTkuNMTkJAACQKDen7KvKIyM3g1JjTU4CAAAkys0p+6qyfeSGcgMAgD+g3JwCp8vQrgNHys1gyg0AAH6BcnMK9lQ3qqXNpV7hoerbu5fZcQAAgCg3p6R9vs3AlBiFcKYUAAB+gXJzCtrn2wziTCkAAPwG5eYUtI/cMJkYAAD/Qbk5BTsqmUwMAIC/odycpFanS18f5LAUAAD+hnJzkvZUN6rVaSg6IlR94jlTCgAAfxFmdoBAc6ixRSu/qtK7n5VLkgamxspi4UwpAAD8BeWmCyrtzbr02U91sMHhXnZWToKJiQAAwPdRbrrgfz7eqYMNDmXYIvWj0X107sAknT0g0exYAADgOyg3nVRW06SlG0olSf/1k9HKo9QAAOCXmFDcSX/4eIdanYbOGZhIsQEAwI9Rbjrh6wMN+r9N+yRJ904abHIaAABwPByWOobNpYe0Zle1tlfUq3B3jZwuQxcNSdHYrN5mRwMAAMdBuTmGdz8r1+L/lLifx0WG6b7JjNoAAODvKDfHkDcgUQcbHDo9PU5D0mM1um+8ekdHmB0LAACcAOXmGC4emqqLh6aaHQMAAHQRE4oBAEBQodwAAICgQrkBAABBhXIDAACCCuUGAAAEFcoNAAAIKpQbAAAQVCg3AAAgqFBuAABAUKHcAACAoEK5AQAAQYVyAwAAggrlBgAABBXKDQAACCqUGwAAEFS8Vm4effRRnX322YqKilJ8fHyH61gslqMeS5cu9VYkAADQA4R56wu3tLToqquuUl5enl544YVjrvfiiy9qypQp7ufHKkIAAACd4bVys2DBAknSkiVLjrtefHy80tLSvBUDAAD0MKbPuZk1a5aSkpJ05plnavHixTIM47jrOxwO2e12jwcAAEA7r43cdMbChQt14YUXKioqSh9++KHuvPNONTQ06K677jrme/Lz892jQgAAAN9nMU40VPIdDzzwgJ544onjrrNt2zYNGTLE/XzJkiW65557VFtbe8KvP3/+fL344osqKys75joOh0MOh8P93G63KzMzU3V1dYqLizvxDwEAAExnt9tls9m88vndpZGbe++9VzNnzjzuOjk5OScdJjc3V4888ogcDoesVmuH61it1mO+BgAA0KVyk5ycrOTkZG9lUVFRkXr37k15AQAAJ81rc25KS0tVU1Oj0tJSOZ1OFRUVSZIGDhyomJgYvfvuu6qsrNRZZ52lyMhILVu2TI899pjuu+8+b0UCAAA9gNfKzfz58/XSSy+5n48ZM0aStGLFCk2YMEHh4eFatGiR5syZI8MwNHDgQD399NO65ZZbvBUJAAD0AF2aUOyPvDkhCQAAeIc3P79Nv84NAABAd6LcAACAoEK5AQAAQYVyAwAAggrlBgAABBXKDQAACCqUGwAAEFQoNwAAIKhQbgAAQFCh3AAAgKBCuQEAAEGFcgMAAIIK5QYAAAQVyg0AAAgqlBsAABBUKDcAACCoUG4AAEBQCTM7wKkyDEOSZLfbTU4CAAA6q/1zu/1zvDsFfLmpr6+XJGVmZpqcBAAAdFV9fb1sNlu3fk2L4Y3K5EMul0v79+9XbGysLBZLt35tu92uzMxMlZWVKS4urlu/dqBhW3hie3yLbeGJ7fEttoUntse32rfF1q1bNXjwYIWEdO8smYAfuQkJCVHfvn29+j3i4uJ6/I7Yjm3hie3xLbaFJ7bHt9gWntge3+rTp0+3FxuJCcUAACDIUG4AAEBQodwch9Vq1UMPPSSr1Wp2FNOxLTyxPb7FtvDE9vgW28IT2+Nb3t4WAT+hGAAA4LsYuQEAAEGFcgMAAIIK5QYAAAQVyg0AAAgqlBsAABBUKDfHsGjRIvXv31+RkZHKzc3V+vXrzY7kdfn5+Ro/frxiY2OVkpKiadOmqbi42GOdCRMmyGKxeDxuv/12kxJ718MPP3zUzzpkyBD3683NzZo1a5YSExMVExOj6dOnq7Ky0sTE3tO/f/+jtoXFYtGsWbMkBf9+sXr1al122WXKyMiQxWLR22+/7fG6YRiaP3++0tPT1atXL02cOFE7duzwWKempkbXX3+94uLiFB8fr5tuukkNDQ0+/Cm6x/G2RWtrq+6//36NGDFC0dHRysjI0A033KD9+/d7fI2O9qfHH3/cxz9J9zjRvjFz5syjftYpU6Z4rNMT9g1JHf4NsVgseuqpp9zrdNe+QbnpwN/+9jfNnTtXDz30kDZt2qRRo0Zp8uTJqqqqMjuaV61atUqzZs3S2rVrtWzZMrW2tmrSpElqbGz0WO+WW25ReXm5+/Hkk0+alNj7hg0b5vGzfvrpp+7X5syZo3fffVdvvPGGVq1apf379+vKK680Ma33bNiwwWM7LFu2TJJ01VVXudcJ5v2isbFRo0aN0qJFizp8/cknn9Szzz6r559/XuvWrVN0dLQmT56s5uZm9zrXX3+9tmzZomXLlum9997T6tWrdeutt/rqR+g2x9sWTU1N2rRpkx588EFt2rRJb775poqLi/WjH/3oqHUXLlzosb/84he/8EX8bneifUOSpkyZ4vGzvvbaax6v94R9Q5LHNigvL9fixYtlsVg0ffp0j/W6Zd8wcJQzzzzTmDVrlvu50+k0MjIyjPz8fBNT+V5VVZUhyVi1apV72Q9+8APj7rvvNi+UDz300EPGqFGjOnyttrbWCA8PN9544w33sm3bthmSjIKCAh8lNM/dd99tDBgwwHC5XIZh9Kz9QpLx1ltvuZ+7XC4jLS3NeOqpp9zLamtrDavVarz22muGYRjG1q1bDUnGhg0b3Ov8+9//NiwWi7Fv3z6fZe9u398WHVm/fr0hydizZ497Wb9+/YxnnnnGu+FM0NH2mDFjhnH55Zcf8z09ed+4/PLLjQsvvNBjWXftG4zcfE9LS4sKCws1ceJE97KQkBBNnDhRBQUFJibzvbq6OklSQkKCx/K//vWvSkpK0vDhwzVv3jw1NTWZEc8nduzYoYyMDOXk5Oj6669XaWmpJKmwsFCtra0e+8mQIUOUlZUV9PtJS0uLXnnlFf385z+XxWJxL+9J+8V3lZSUqKKiwmNfsNlsys3Nde8LBQUFio+P1xlnnOFeZ+LEiQoJCdG6det8ntmX6urqZLFYFB8f77H88ccfV2JiosaMGaOnnnpKbW1t5gT0gZUrVyolJUWDBw/WHXfcoerqavdrPXXfqKys1D//+U/ddNNNR73WHftGwN8VvLsdPHhQTqdTqampHstTU1O1fft2k1L5nsvl0j333KNzzjlHw4cPdy+/7rrr1K9fP2VkZOjzzz/X/fffr+LiYr355psmpvWO3NxcLVmyRIMHD1Z5ebkWLFig8847T19++aUqKioUERFx1B/s1NRUVVRUmBPYR95++23V1tZq5syZ7mU9ab/4vvb/3h39zWh/raKiQikpKR6vh4WFKSEhIaj3l+bmZt1///269tprPe6Cfdddd2ns2LFKSEjQmjVrNG/ePJWXl+vpp582Ma13TJkyRVdeeaWys7O1a9cu/frXv9bUqVNVUFCg0NDQHrtvvPTSS4qNjT3qUH537RuUG3Ro1qxZ+vLLLz3mmEjyOA48YsQIpaen66KLLtKuXbs0YMAAX8f0qqlTp7r/PXLkSOXm5qpfv356/fXX1atXLxOTmeuFF17Q1KlTlZGR4V7Wk/YLdE5ra6t+8pOfyDAMPffccx6vzZ071/3vkSNHKiIiQrfddpvy8/OD7r5L11xzjfvfI0aM0MiRIzVgwACtXLlSF110kYnJzLV48WJdf/31ioyM9FjeXfsGh6W+JykpSaGhoUed9VJZWam0tDSTUvnW7Nmz9d5772nFihXq27fvcdfNzc2VJO3cudMX0UwVHx+v0047TTt37lRaWppaWlpUW1vrsU6w7yd79uzRRx99pJtvvvm46/Wk/aL9v/fx/makpaUddUJCW1ubampqgnJ/aS82e/bs0bJlyzxGbTqSm5urtrY27d692zcBTZSTk6OkpCT370ZP2zck6ZNPPlFxcfEJ/45IJ79vUG6+JyIiQuPGjdPy5cvdy1wul5YvX668vDwTk3mfYRiaPXu23nrrLX388cfKzs4+4XuKiookSenp6V5OZ76Ghgbt2rVL6enpGjdunMLDwz32k+LiYpWWlgb1fvLiiy8qJSVFl1566XHX60n7RXZ2ttLS0jz2BbvdrnXr1rn3hby8PNXW1qqwsNC9zscffyyXy+UugsGivdjs2LFDH330kRITE0/4nqKiIoWEhBx1eCYY7d27V9XV1e7fjZ60b7R74YUXNG7cOI0aNeqE6570vnHKU5KD0NKlSw2r1WosWbLE2Lp1q3Hrrbca8fHxRkVFhdnRvOqOO+4wbDabsXLlSqO8vNz9aGpqMgzDMHbu3GksXLjQ2Lhxo1FSUmK88847Rk5OjnH++eebnNw77r33XmPlypVGSUmJ8Z///MeYOHGikZSUZFRVVRmGYRi33367kZWVZXz88cfGxo0bjby8PCMvL8/k1N7jdDqNrKws4/777/dY3hP2i/r6emPz5s3G5s2bDUnG008/bWzevNl9BtDjjz9uxMfHG++8847x+eefG5dffrmRnZ1tHD582P01pkyZYowZM8ZYt26d8emnnxqDBg0yrr32WrN+pJN2vG3R0tJi/OhHPzL69u1rFBUVefwdcTgchmEYxpo1a4xnnnnGKCoqMnbt2mW88sorRnJysnHDDTeY/JOdnONtj/r6euO+++4zCgoKjJKSEuOjjz4yxo4dawwaNMhobm52f42esG+0q6urM6KiooznnnvuqPd3575BuTmGP/zhD0ZWVpYRERFhnHnmmcbatWvNjuR1kjp8vPjii4ZhGEZpaalx/vnnGwkJCYbVajUGDhxo/PKXvzTq6urMDe4lV199tZGenm5EREQYffr0Ma6++mpj586d7tcPHz5s3HnnnUbv3r2NqKgo44orrjDKy8tNTOxdH3zwgSHJKC4u9ljeE/aLFStWdPi7MWPGDMMwjpwO/uCDDxqpqamG1Wo1LrrooqO2U3V1tXHttdcaMTExRlxcnHHjjTca9fX1Jvw0p+Z426KkpOSYf0dWrFhhGIZhFBYWGrm5uYbNZjMiIyON008/3Xjsscc8PuwDyfG2R1NTkzFp0iQjOTnZCA8PN/r162fccsstR/2Pck/YN9r96U9/Mnr16mXU1tYe9f7u3DcshmEYXRvrAQAA8F/MuQEAAEGFcgMAAIIK5QYAAAQVyg0AAAgqlBsAABBUKDcAACCoUG4AAEBQodwAAICgQrkBAABBhXIDAACCCuUGAAAElf8P7Y3kKRa5/nAAAAAASUVORK5CYII=",
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAABLwAAAGsCAYAAADXMb4GAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAAB3r0lEQVR4nO3deXxU9b3/8ffsM9lJgCwQdhAXXFFErNWKoLVWr3axtS1Vq+0VW5V7W6W32p9Wpdpby9Vaqf1Ze71XbeuvVm1rrUjdRUAQFUH2HZIAIZlkktnP74+ZTGayDEmY5Ewmr+fjwSMz55w5+c7XEU7e+Xw/x2IYhiEAAAAAAAAgR1jNHgAAAAAAAACQSQReAAAAAAAAyCkEXgAAAAAAAMgpBF4AAAAAAADIKQReAAAAAAAAyCkEXgAAAAAAAMgpBF4AAAAAAADIKXazB9BRNBrVvn37VFhYKIvFYvZwAADAIGEYhpqamlRVVSWrld/pZSOu8wAAQF/05Tov6wKvffv2qbq62uxhAACAQWr37t0aPXq02cNAF7jOAwAAR6M313lZF3gVFhZKir2JoqKijJ8/FArp5Zdf1pw5c+RwODJ+/sGO+UmP+UmP+UmP+UmP+UmP+UkvFArpueee07e+9a3EtQSyT39f5wEAgNzk9XpVXV3dq+u8rAu82srbi4qK+i3wysvLU1FRET8wdIH5SY/5SY/5SY/5SY/5SY/5Sa9tfiSxVC6L9fd1HgAAyG29uc6jwQUAAAAAAAByCoEXAAAAAAAAcgqBFwAAAAAAAHIKgRcAAAAAAAByCoEXAAAAAAAAcgqBFwAAAAAAAHIKgRcAAAAAAAByCoEXAAAAAAAAcgqBFwAAAAAAAHIKgRcAAAAAAAByCoEXAAAAAAAAcgqBFwAAAAAAAHIKgRcAAAAAAAByit3sAQAAgNxnGIaWbzskp82q6eNKzR4OctSWumbVNPo1pjRPY8ryzB4OAAAwERVeAACgX0Sjhmq9fv3fN7fp/J+/rq/+ZoX+8+WNZg8LOex372zX1x5boT+t2WP2UAAAgMmo8AIAABnhC4T10roaPf/BPm2qadLB5oDCUSOxP99p08QRBQpHorLb+J0bAAAA+g+BFwAAOGoPLdusX722Va2hSKd900YV6ytnjNHnT65SgYtLDwAAAPQ/rjoBAMBRe/TNbWoNRTR+eL7+5ZRROmfKCJUXuVSW75LTTjUXAAAABhaBFwAAOGr+eGXXk9+aoaoSj8mjAQAAwFDXq1+5RiIR3X777Ro/frw8Ho8mTpyon/zkJzKM9v4chmHojjvuUGVlpTwej2bPnq3NmzdnfOAAACA7RKOGQpHYtYCLai4AAABkgV5dld5333165JFH9Mtf/lIbNmzQfffdp/vvv18PPfRQ4pj7779fDz74oJYsWaIVK1YoPz9fc+fOld/vz/jgAQCA+YKRaOKxy2EzcSQAAABATK+WNL7zzju69NJLdfHFF0uSxo0bp6efflorV66UFKvuWrx4sX70ox/p0ksvlSQ98cQTKi8v13PPPacrr7yy0zkDgYACgUDiudfrlSSFQiGFQqG+vas02s7ZH+fOBcxPesxPesxPesxPesxPetk8P82t7WOyGhGFQkaao/tHNs4LAAAAzNOrwOuss87So48+qk2bNmnKlCn64IMP9NZbb+mBBx6QJG3fvl01NTWaPXt24jXFxcWaMWOGli9f3mXgtWjRIt15552dtr/88svKy8vr7fvpsaVLl/bbuXMB85Me85Me85Me85Me85NeNs5PY1CS7LLI0MsvvSSLxewRAQAAYKjrVeB12223yev1aurUqbLZbIpEIrrnnnt01VVXSZJqamokSeXl5SmvKy8vT+zraOHChVqwYEHiudfrVXV1tebMmaOioqJevZmeCIVCWrp0qS644AI5HI6Mn3+wY37SY37SY37SY37SY37Sy+b52XO4VVr9plwOmy6+eK4pYwiFQnr++edN+d4AAADIPr0KvP74xz/qySef1FNPPaXjjz9ea9eu1c0336yqqirNmzevTwNwuVxyuVydtjscjn69oO/v8w92zE96zE96zE96zE96zE962Tg/EcVaE7jstqwbGwAAAIamXgVe3//+93XbbbclliZOmzZNO3fu1KJFizRv3jxVVFRIkmpra1VZWZl4XW1trU4++eTMjRoAAGSNQDgiiTs0AgAAIHv06sq0paVFVmvqS2w2m6LR2N2Zxo8fr4qKCi1btiyx3+v1asWKFZo5c2YGhgsAALJNIBy7DnA5CLwAAACQHXpV4XXJJZfonnvu0ZgxY3T88cfr/fff1wMPPKBrrrlGkmSxWHTzzTfr7rvv1uTJkzV+/Hjdfvvtqqqq0mWXXdYf4wcAACYLhGKBl9NG4AUAAIDs0KvA66GHHtLtt9+uG264QXV1daqqqtK3v/1t3XHHHYljfvCDH8jn8+n6669XQ0ODzj77bL300ktyu90ZHzwAADBfMBKv8LLbTB4JAAAAENOrX8UWFhZq8eLF2rlzp1pbW7V161bdfffdcjqdiWMsFovuuusu1dTUyO/365VXXtGUKVMyPnAAAJAdAqF4Dy+WNA5qkUhEt99+u8aPHy+Px6OJEyfqJz/5iQzDSBxjGIbuuOMOVVZWyuPxaPbs2dq8ebOJowYAAOgaV6YAAOCoJHp40bR+ULvvvvv0yCOP6Je//KU2bNig++67T/fff78eeuihxDH333+/HnzwQS1ZskQrVqxQfn6+5s6dK7/fb+LIAQAAOuvVkkYAAICO2gIvJ0saB7V33nlHl156qS6++GJJ0rhx4/T0009r5cqVkmLVXYsXL9aPfvQjXXrppZKkJ554QuXl5XruuecSd/FOFggEFAgEEs+9Xu8AvBMAAAAqvAAAwFEKUuGVE8466ywtW7ZMmzZtkiR98MEHeuutt3TRRRdJkrZv366amhrNnj078Zri4mLNmDFDy5cv7/KcixYtUnFxceJPdXV1/78RAAAAUeEFAACOUiAc7+FF4DWo3XbbbfJ6vZo6dapsNpsikYjuueceXXXVVZKkmpoaSVJ5eXnK68rLyxP7Olq4cKEWLFiQeO71egm9AADAgCDwAgAAR6W9hxdLGgezP/7xj3ryySf11FNP6fjjj9fatWt18803q6qqSvPmzevTOV0ul1wuV4ZHCgAAcGQEXgAA4KgEQvHAi7s0Dmrf//73ddtttyV6cU2bNk07d+7UokWLNG/ePFVUVEiSamtrVVlZmXhdbW2tTj75ZDOGDAAA0C2uTAEAwFFpW9LotHFZMZi1tLTIak39b2iz2RSNxgLN8ePHq6KiQsuWLUvs93q9WrFihWbOnDmgYwUAADgSKrwAAMBRSTStp8JrULvkkkt0zz33aMyYMTr++OP1/vvv64EHHtA111wjSbJYLLr55pt19913a/LkyRo/frxuv/12VVVV6bLLLjN38AAAAB0QeAEAgKNCD6/c8NBDD+n222/XDTfcoLq6OlVVVenb3/627rjjjsQxP/jBD+Tz+XT99deroaFBZ599tl566SW53W4TRw4AANAZgRcAADgq3KUxNxQWFmrx4sVavHhxt8dYLBbddddduuuuuwZuYAAAAH3AlSkAADgq7RVeXFYAAAAgO3BlCgAAjkqQwAsAAABZhitTAABwVOjhBQAAgGxD4AUAAI5KoocXd2kEAABAluDKFAAAHJVAiCWNAAAAyC5cmQIAgKMSjMQCLyeBFwAAALIEV6YAAOCotFd40cMLAAAA2YHACwAAHJVEDy8qvAAAAJAluDIFAABHhbs0AgAAINsQeAEAgKPSFnjRwwsAAADZgitTAABwVIJh7tIIAACA7MKVKQAAOCqJHl4OLisAAACQHbgyBQAAfRaJGgpFDEn08AIAAED2IPACAAB91racUWJJIwAAALIHV6YAAKDPkgMvmtYDAAAgW3BlCgAA+qytf5fVItmtFpNHAwAAAMQQeAEAgD4LJO7QaJPFQuAFAACA7EDgBQAA+ow7NAIAACAbcXUKAAD6zB+KVXg5bVxSAAAAIHtwdQoAAPosGIkvaaTCCwAAAFmEq1MAANBngVB7Dy8AAAAgWxB4AQCAPkv08LJzSQEAAIDswdUpAADos/a7NHJJAQAAgOzB1SkAAOizYDzwchJ4AQAAIItwdQoAAPqsvcKLHl4AAADIHgReAACgz+jhBQAAgGzE1SkAAOizxF0aHVR4AQAAIHsQeAEAgD4LRuI9vGxcUgAAACB7cHUKAAD6rL3Ci0sKAAAAZA+uTgEAQJ/RwwsAAADZiKtTAADQZ9ylEQAAANmIwAsAAPQZFV4AAADIRlydAgCAPgvGK7ycBF4AAADIIlydAgCAPmtf0sglBQAAALIHV6cAAKDP2u/SSA8vAAAAZA8CLwAA0Gf08AIAAEA24uoUAAD0WTDCkkYAAABkH65OAQBAnyWWNBJ4AQAAIItwdQoAAPqsvWk9PbwAAACQPQi8AABAn9HDCwAAANmIq1MAANBnwbYKLweXFAAAAMgeXJ0CAIA+a1vS6LSxpBEAAADZg8ALAAD0WYAKLwAAAGQhrk4BAECfBUL08AIAAED24eoUAAD0GXdpBAAAQDYi8AIAAH0SiRoKRw1JkpMKLwAAAGQRrk4BAECftN2hUWJJIwAAALILV6cAAKBPAuFI4jGBFwAAALIJV6cAAKBP2vp32awW2W1cUgAAACB7cHUKAAD6JJhoWM/lBAAAALILV6gAAKBP2pY00rAeAAAA2YYrVAAA0Cf+EBVeAAAAyE69vkLdu3evvva1r6msrEwej0fTpk3Te++9l9hvGIbuuOMOVVZWyuPxaPbs2dq8eXNGBw0AAMwXSCxptJk8EgAAACBVrwKvw4cPa9asWXI4HPr73/+u9evX6+c//7mGDRuWOOb+++/Xgw8+qCVLlmjFihXKz8/X3Llz5ff7Mz54AABgnrYljVR4AQAAINvYe3Pwfffdp+rqaj3++OOJbePHj088NgxDixcv1o9+9CNdeumlkqQnnnhC5eXleu6553TllVdmaNgAAMBsbU3r6eEFAACAbNOrwOuFF17Q3Llz9cUvflGvv/66Ro0apRtuuEHXXXedJGn79u2qqanR7NmzE68pLi7WjBkztHz58i4Dr0AgoEAgkHju9XolSaFQSKFQqE9vKp22c/bHuXMB85Me85Me85Me85Me85NeNs5Piz82FqfNYvq4zP7+AAAAyC69Cry2bdumRx55RAsWLNAPf/hDrVq1St/73vfkdDo1b9481dTUSJLKy8tTXldeXp7Y19GiRYt05513dtr+8ssvKy8vrzfD65WlS5f227lzAfOTHvOTHvOTHvOTHvOTXjbNz5qDFkk2NTce1osvvmj2cAAAAICEXgVe0WhU06dP17333itJOuWUU7Ru3TotWbJE8+bN69MAFi5cqAULFiSee71eVVdXa86cOSoqKurTOdMJhUJaunSpLrjgAjkcjoyff7BjftJjftJjftJjftJjftLLxvlpXbNX2vyxqspH6rOfPdXUsYRCIT3//POmjgEAAADZo1eBV2VlpY477riUbccee6z+9Kc/SZIqKiokSbW1taqsrEwcU1tbq5NPPrnLc7pcLrlcrk7bHQ5Hv17Q9/f5BzvmJz3mJz3mJz3mJz3mJ71smp+ILJIkl8OWNWMCAAAApF7epXHWrFnauHFjyrZNmzZp7NixkmIN7CsqKrRs2bLEfq/XqxUrVmjmzJkZGC4AAMgWgVCsab3LYTN5JMiUvXv36mtf+5rKysrk8Xg0bdo0vffee4n9hmHojjvuUGVlpTwej2bPnq3NmzebOGIAAICu9SrwuuWWW/Tuu+/q3nvv1ZYtW/TUU0/p0Ucf1fz58yVJFotFN998s+6++2698MIL+uijj/SNb3xDVVVVuuyyy/pj/AAAwCSB+F0aXdylMSccPnxYs2bNksPh0N///netX79eP//5zzVs2LDEMffff78efPBBLVmyRCtWrFB+fr7mzp0rv99v4sgBAAA669WSxtNPP11//vOftXDhQt11110aP368Fi9erKuuuipxzA9+8AP5fD5df/31amho0Nlnn62XXnpJbrc744MHAADmCYQjkgi8csV9992n6upqPf7444lt48ePTzw2DEOLFy/Wj370I1166aWSpCeeeELl5eV67rnnenU3bgAAgP7W6yvUz33uc/roo4/k9/u1YcMGXXfddSn7LRaL7rrrLtXU1Mjv9+uVV17RlClTMjZgAACQHYKJCi+WNOaCF154QdOnT9cXv/hFjRw5Uqeccop+85vfJPZv375dNTU1mj17dmJbcXGxZsyYoeXLl3d5zkWLFqm4uDjxp7q6ut/fBwAAgNSHwAsAAEBqX9LopMIrJ2zbtk2PPPKIJk+erH/84x/613/9V33ve9/Tf//3f0uSampqJEnl5eUprysvL0/s62jhwoVqbGxM/Nm9e3f/vgkAAIC4Xi1pBAAAaMOSxtwSjUY1ffp03XvvvZKkU045RevWrdOSJUs0b968Pp2zu7txAwAA9DeuUAEAQJ+036WRy4lcUFlZqeOOOy5l27HHHqtdu3ZJkioqKiRJtbW1KcfU1tYm9gEAAGQLrlABAECfBOjhlVNmzZqljRs3pmzbtGmTxo4dKynWwL6iokLLli1L7Pd6vVqxYoVmzpw5oGMFAAA4EpY0AgCAPgnSwyun3HLLLTrrrLN077336ktf+pJWrlypRx99VI8++qik2I2Jbr75Zt19992aPHmyxo8fr9tvv11VVVW67LLLzB08AABABwReAACgT+jhlVtOP/10/fnPf9bChQt11113afz48Vq8eLGuuuqqxDE/+MEP5PP5dP3116uhoUFnn322XnrpJbndbhNHDgAA0BmBFwAA6JP2JY0EXrnic5/7nD73uc91u99iseiuu+7SXXfdNYCjAgAA6D2uUAEAQJ/QwwsAAADZigovAADQI3VNfv3u7R1atqFOtU1+NbSEJFHhBQAAgOxD4AUAADoxDENrdh3W3ga/DvuCWre3Uc+v3adgJJpyXGWxWyeMKjZplAAAAEDXCLwAAEAnv1+1Wwuf/ajT9uljh+mbs8bpmPJCjSh0qdjjkMViMWGEAAAAQPcIvAAAQIpo1NCjb2yTJB1fVaSxZXkaUeDS50+u0mljS00eHQAAAHBkBF4AACDFqxvrtP2gT0Vuu/747ZnKd3G5AAAAgMGFLrMAACDFY29tlyR95YwxhF0AAAAYlAi8AABAwob9Xr2z9ZBsVou+cdY4s4cDAAAA9AmBFwAASHj87Vh114XHV2hUicfk0QAAAAB9wzoFAACGsG0HmvWHVbu1t6FVwXBUr208IEm65uxx5g4MAAAAOAoEXgAADBF1Xr/W7GpQMBKVPxjR39ft16vxgCvZqWNKdOqYYSaMEAAAAMgMAi8AAIaIL/16uXYcaknZZrFInzlmpM6aNFxuh1Vuu03nTBkhi8Vi0igBAACAo0fgBQDAEFDvCybCrpkTyuS0WzVpZIG+fuZYjRueb/LoAAAAgMwi8AIAYAjYdqBZklRV7NbT159p8mgAAACA/sVdGgEAGAK2HfBJkiaMKDB5JAAAAED/I/ACAGAI2HowVuE1YQTLFwEAAJD7CLwAABgCEhVe9OsCAADAEEDgBQDAENDWw2viSJY0AgAAIPcReAEAkOPCkah21cfu0EgPLwAAAAwFBF4AAOS43YdbFYoYcjusqixymz0cAAAAoN8ReAEAkOO21sWWM44fXiCr1WLyaAAAAID+R+AFAECO28YdGgEAADDEEHgBAJDj2u7QOJE7NAIAAGCIIPACACDHtQVeNKwHAADAUEHgBQBAjmNJIwAAAIYaAi8AAHJYY2tIB5uDkqTxLGnEEGGYPQAAAGA6Ai8AAHLYtgOx6q7yIpcK3Q6TRwP0r+EFLknSG5sOyDCIvQAAGMoIvAAAyGGJ/l3D6d+F3PfVGWPkdli1dneDXt1YZ/ZwAACAiQi8AADIYfTvwlAystCteTPHSZIeWLqJKi8AAIYwAi8AAHKIYRj6eF+jnli+Q//3zW16beMBSdyhEUPHtz89UflOm9bt9eofH9eaPRwAAGASu9kDAAAAR6c1GNGGGq9WbKvXc+/v1cbapk7HTB5J4IWhoTTfqWvOHq+H/rlFv1i6SXOOK5fVajF7WAAAYIAReAEAkCUMw9Cu+hbtbWhVsz8sXzCs5kBEzf6wvK0BfbzdqtefXaeWYFS+YFhN/rAaW0PaecinaNLKLafdqlkTy1Tscchmtaq61KNZk4ab98aAAfatsyfod+/s0MbaJv31o/36/ElVZg8JAAAMMAIvAABMZBiGnlm9R8+u2aOP93nV5A+nOdoq1ezrcs+IQpeOryrS3OMr9NlplSr2cEdGDF3FeQ5d96kJemDpJi16cYM+3N2gyhKPqordqizxqLLYreEFLtmo/AIAIGcReAEAYJKGlqBu+9NHeunjmsQ2py1WkVXodqjAZVeBy658l115Dotq9uzUSccdo+I8p/Lj+wrdDk0cma+RhW4T3wmQfa6eNU6/e2eH9jf69X/f2t5pv91qUXmRW5XxEKyq2K2KYrcqiz2qKol9Lct3shwSAIBBisALAAATrNxer5t//772NfrlsFn0vc9M1vnHlmtyeYEcts73lAmFQnrxxe367Dnj5XBQvQUcSaHboT9++0y9vumg9je0an+jX/sbY19rvX6Fo4b2NrRqb0OrtPNwl+dw2qwqL3apsjhWFdYWhlUUuVUVrxQrzXfKYiEUAwAg2xB4AQAwgMKRqH756hY9uGyzooY0rixPD33lVE0bXWz20ICcM2lkoSaNLOy0PRyJ6kBzQPsaYiFYTaM/8Xhfo181ja2qawooGIlqd32rdte3dvs9nHZrPAxztwdjJR5VFrlVWeJWVbFHJXkOQjEAAAYYgRcAAP1s5yGfth/0qabRr2fX7NXKHfWSpCtOHa07Lz1eBS7+OQYGkt1mjYdTHknDujwmFImq1uuPhWGN/k5VYvsa/DrYHFAwHNXOQy3aeail2+/ndlhTqsRioVgsDKsojn0t8tgJxQAAyCCusAEA6EcvfrRfNzy5JmVbvtOme/5lmi47ZZRJowJwJA6bVaOH5Wn0sLxujwmGY6HYvoZW1Xjbq8QSwViDX4d8QflDUW0/GAu+u5PntHWuEotXjlWVxIKxIjfLmQEA6CkCLwAA+tGLH+2XJI0q8WhKeYHGlObpmrPHa2xZvskjA3C0nHarqkvzVF3afSjmD0XioVjnMKzt8eGWkFqCEW094NPWA92HYgUuuyrjzfWrij2qLHGn9BarLPYon4pRAAAkEXgBANBvDMPQu9tiyxcf+NJJmjGhzOQRARhobodNY8vy04bcrcFIey+xtuWT3uRllH41tobUHAhrc12zNtc1d3uuQrc9JQy79ORROpO/ewAAQxCBFwAA/WTrgWYdbA7IZbfq5DElZg8HQJbyOG2aMKJAE0YUdLnfH4po2wGf1u5u0Jpdh/X+rsPdVoI1+cPa6G/SxtomSdLa3Y36+02f6rexAwCQrQi8AADoJ8vj1V2njR0ml91m8mgAZBvDMFTvC6rG6483yA+opjHWD6zGG1Bto1813lh1V09YLdKIQpcqitwqL4otfbzkpKp+fhcAAGQnAi8AAPrJu9sOSRLLiYAhKBCOqM4biIVXjf7EHR/3e/2JIKvOG1AwEu3R+TyOWFP7tiCrvMitiiJX4nFlsUfDC5yy26z9/M4AABgcCLwAAOgHhmFoRTzwmjmRwAvIFYZhqLE11CHICqjGG+vBVeMNqNbrV70v2ONzDi9wxkOrtiDLrfLi2Ne2QKvIbZfFYunHdwYAQG4h8AIAoB9sqWvWweag3A6rThxdbPZwAPRAKBJVXVMgpSKr1htrGt+27LDW65c/1LOqLKfdGgut4gFWcqBVUexSeZFbIwvdctqpygIAINMIvAAA6AfL49Vd08eW0r8LMJlhGGoKhBNLCROBVvxxTbxK65AvIMPo2TmH5TkSywuTe2a1Pa8ocqskz0FVFgAAJiHwAgCgH7T37yo1eSTA0OELhPXCB/u0q74lJdyq8frVEoz06BwOm0UjCzsGWa5EZVZlsUcji1xyOwiyAQDIZgReAABkWDRq6N34HRrp3wUMnN++tV0/X7qpT68tcts1YUSBqkvzVOi2K99pU57TrnxX7KthSI2tIYWjhg75Aipw2ZXnaj+OZYkAAGQXAi8AADJs3b5G1fuC8jhsmjaqxOzhAEPGZ0+s1Oa6Zh1uCcoXCKslGFFz/KsvEFYg3H3vLa8/rLW7G7R2d0OfvrfDZlG+y658p115TpvyXHYVxMOy/PjzfKet/RiXTflOe/x5+/48l10F8f0O7rgIAECfEXgBAJBBG/Z79a3/fk+SNGvScKo+gAE0cUSBHvzKKd3uD0ei8gUjagmG5Qu0f/UFwvIF24OxxL5gWC2BiHxtx6U8D8sXjCgYD9FCEUMNLSE1tIQy9n6cNmtSMJZacVbgigVr+fGvsedJx7WFa4nnhGgAgKGFwAsAgAx5d9shXffEe2ryh3VMeaHuvuwEs4cEIIndZlWxx6pijyNj5wxFoomgrKXLYCySst2XVHHWXoGWelwwEgvRgpGogi3RzIZodmvKcs3kqrS2gCz2PDVk61SVlhSs2QnRAABZiMALAICj1NAS1EP/3KInlu9QKGLojHGl+s03pqs4L3M/VAPITo5+CNGC4ahag/HgLBhWcyCilnhFWex5e6DWXpUWTqleSwRr8fAtEaKFowqGozqcwRDNZbcmKs3agrGCDs87BmtdVaXF+qLZlOcgRAMAHD0CLwAAeiEaNbT9kE87DvpU6w1oV32LnlqxU15/WJJ08bRK/fxLJ3EHNwB95rRb5bRbMxqaB8PR+DLN9vCsLShLDsYSFWfx45oTyztjz9uq0nyBsMJRQ5IUCEcVCAdV78vYcBMhWnvFWef+Z931QevuOJvVkrkBAgCyHoEXAABpBMIRvb+rQcu3HtJ7O+v14Z5GNcXDrWRTKwq18LPH6tNTRpgwSgBILxaiOVWSl7lztoVozR2Wabb1ROvcJ61j/7TOAVx/hmhuh7VTUNa5Ki3p5gFJSzqPqyxWRbE7c4MBAPQ7Ai8AACT5AmGt2lGv5VsP6YM9DfK2xn6Iq/X6O93ZzWW3anJ5gSqK3BpZ5Nb0scN06cmjqB4AMKS0h2jOjJzPMAwFI9GU5Zqpyzc7LO/sIljrGMD5ghFF4iGaPxSVPxTUoT6EaIUuu967fbZcdqp3AWCwIPACAAxJ7+86rNv+9JH2NbQqEI4m+tt0ZXiBUzMnDteM8aU6ZUyJppQXcqczAMgwi8Uil90ml92mYfldh2iGYaglGFG9L6jDLcGkryEd9gVV3xJUQ9t2X0iHfEEd8gVkGH0bU6HLrmH5Tp01sUxO/t4HgEHlqAKvn/70p1q4cKFuuukmLV68WJLk9/v1b//2b/r973+vQCCguXPn6le/+pXKy8szMV4AAI7ah3sa9I3fruy0NHFUiUezJpXp9HGlGlHoUqHbrtJ8l8aV5clioXoLADLNH4qFVykBli+o+pb2AOuwL6jDSc+D4e5/QZFOntOmYXlOleY7NSzfqdI8R/yrUyXxr8PyHSpt25bnlNNOyAUAg1WfA69Vq1bp17/+tU488cSU7bfccov+9re/6ZlnnlFxcbFuvPFGXX755Xr77bePerAAABytj/c16uuPxcKu08cN06LLT1Se06Y8py1jy3IAYCgKhCM67AvpcEswJayqj2/rHGoF5Q/1Lbxy2q0qy3d2HWDlx8Kq5ABrWJ6Tm4kAwBDTp8CrublZV111lX7zm9/o7rvvTmxvbGzUY489pqeeekqf+cxnJEmPP/64jj32WL377rs688wzMzNqAAD6oK7Jr68/tlKNrSGdOqZEj199hgpcrO4HgI5CkWg8uAqp3hdfJniEAMsXjPTpezlslvbgKhFgOeKBVfv2YUkBlsdho/IWAJBWn67y58+fr4svvlizZ89OCbxWr16tUCik2bNnJ7ZNnTpVY8aM0fLly7sMvAKBgAKBQOK51+uVJIVCIYVCob4ML622c/bHuXMB85Me85Me85Me85PeQMzP797arnpfUMeUF+j/fv0UuazGoPnvwecnPeYl82hdkTsiUUMNLe29rpKDqoaWrgOsru5G2xM2q0XD8hzxcMqZFFo5OlRjtVViOVTgshNeAQAyrteB1+9//3utWbNGq1at6rSvpqZGTqdTJSUlKdvLy8tVU1PT5fkWLVqkO++8s9P2l19+WXl5GbxvcgdLly7tt3PnAuYnPeYnPeYnPeYnvf6an1BUemK1TZJFZxU36s1/Ds7/Dnx+MBBoXZHdWoJh1TT6OzVrP5zohZW6rLCxNdSnpu0WizQsLxZKlfYgwBqW71Shyy4rd6wFAGSBXgVeu3fv1k033aSlS5fK7XZnZAALFy7UggULEs+9Xq+qq6s1Z84cFRUVZeR7JAuFQlq6dKkuuOACORyOjJ9/sGN+0mN+0mN+0mN+0uvv+fl/a/aqecXHqip26wdfPVv2QXa3LT4/6YVCIT3//PNmDyMnZLJ1RXeV/EjPH4pob0Ordte3aM/hVu0+HPu6J/78kC/Yp/MWe9r6WSVVYCWWEXYOsIo8DtkIrwAAg1SvAq/Vq1errq5Op556amJbJBLRG2+8oV/+8pf6xz/+oWAwqIaGhpQqr9raWlVUVHR5TpfLJZfL1Wm7w+Ho1wv6/j7/YMf8pMf8pMf8pMf8pNcf82MYhp54d7ck6RtnjZPH3fnfncGCzw/6WyZbV3RXyT/UBcNR7WtoTQqzWrS7vlV74sFWXVPgiOcocNlTm7UfIcAq8TgGXdAPAMDR6FXgdf755+ujjz5K2Xb11Vdr6tSpuvXWW1VdXS2Hw6Fly5bpiiuukCRt3LhRu3bt0syZMzM3agAAemHF9npt2O+V22HVladXmz0cIGtlunVFd5X8uS4ciWp/o7+9MiupOmv34RbVeP1HXGKY77SpujRPo4d5NHpY+9fq0tjXYg/BNwAA6fQq8CosLNQJJ5yQsi0/P19lZWWJ7ddee60WLFig0tJSFRUV6bvf/a5mzpzJHRoBAKb53ds7JEmXnzpaJXlOcwcDZKn+aF3RXSX/YBeJGqpr8ieqshJf4wHX/ka/ItH0iZbbYU0EWdVtX0vbn5fkOWjkDgDAUcj4vdh/8YtfyGq16oorrki5ew8AAGb4n3d36h/rY9Un3zxrnLmDAbJYf7SuGKwMw9CBpoB2H25fZtjWT2vP4RbtbWhVKJI+0HLarBo1zNOpMqst0Bpe4CTQAgCgHx114PXaa6+lPHe73Xr44Yf18MMPH+2pAQDoM8Mw9J8vb9TDr26VFAu7ppQXmjwqIHsNpdYVhmGo3hdMaQif3CB+7+FWBcLRtOewWy2qKmkLtOJVWqVt1Vp5Glno4m6FAACYKOMVXgAAZIM7nv9Y//PuTknSLbOn6HvnTzJ5REB2y6XWFYZhyNsa7tQQPrliqyUYSXsOq0WqLPZoVFdLDkvzVF7oogk8AABZjMALAJBzdhz06X/e3SmLRVr0L9N05RljzB4SkBMGQ+uK1mBElz38tjbWNh3x2PIiV6cwa/SwPFUPy1NFsVtOO4EWAACDFYEXACDnPLtmjyTpU5NHEHYBR2Ewtq74YE9DIuwaXuBs75tVmtogvqrEI7fDZvJoAQBAfyHwAgDklGjU0J/W7JUkXXHqKJNHA2CgbT3QLEk675gRevzqM0weDQAAMAt12gCAnLJie732NrSq0GXX3ONz685xAI5sa51PkjRxRIHJIwEAAGYi8AIA5JQ/xZczXnxiJcuVgCGorcJr4kgCLwAAhjICLwBAzmgJhvX3j/ZLkq44bbTJowFghkTgRYUXAABDGoEXACBnvLSuRr5gRGPL8jR97DCzhwNggLUGI9rb0CpJmjgi3+TRAAAAM9G0HgCQtXbXt+hXr21RnTegQDiqYDiqQDiS9LjtTyTxXJIuP2W0LBaLyaMHMNC2H/TJMKSSPIdK851mDwcAAJiIwAsAkJVe/rhG//7MB/L6w7163bA8h750OssZgaGo1uuXJI0e5iH0BgBgiCPwAgBkjeaQ9Mbmg3rlk4N6euUuSdJJ1SX6yunVcjtsctqtctmtctmTHjusctqscjlsctmtKnI75LSzYh8YilpDEUlSnoNLXAAAhjquBgAApvMFwrr68VVaucMuvbcmsf26T43X9+dOJcAC0CP+eODlcvB3BgAAQx2BFwDAdA8s3aSVOw5LkiYMz9OJo0t0+amjdc6UESaPDMBg0lbh5XbYTB4JAAAwG4EXAMBUH+1p1ONvb5ckfWdqRP921dlyOBwmjwrAYOQPxW5c4SHwAgBgyKPeGwBgmnAkqoV//lBRQ/rctAodO8wwe0gABjF/osKLS1wAAIY6KrwAAAMiGjXU0BrSweaAGlpC8gXDemfLQa3b61WR267/+OwxWvnGHrOHCWAQC7CkEQAAxBF4AQAyzjAMef1h7Tjo0+ubDui1jXX6cE+jwtGuK7gWfvZYDS9wDfAoAeSaQDi2pNHFjS4AABjyCLwAABmzZtdhLfjDWu1taFUo0nW4VZLn0LA8p/JdNuU57DpxdLG+PL1akUh4gEcLINdEjdjfO1aLxeSRAAAAsxF4AQAy5n/f3akdh1oSz4vcds2cWKbzjhmpsyYOV0WxW85uKi8ikYEaJYBc1VZEaiHwAgBgyCPwAgBkhGEYemfLIUnSr646VZ+ZOpI+OgAGVLzAS1byLgAAhjwaHAAAMmLrAZ9qvH457VbCLgCmYEkjAABoQ+AFAMiId7YelCRNHzuMsAuAKYxE4GXyQAAAgOlY0ggAyIi3t8QCr1mThps8EgBDVVsPr/d2Htbv3t6uIo9DxR6HijwOFbkdKvLYVeR2KM9po88XAAA5jsALAHDUIlFDy7fG+nedNbHM5NEAGKryXLHq0ne2HtI78b+TumKzWlTktncKwlIeezo8JjADAGBQIfACABy1dXsb5fWHVei2a9qoYrOHA2CIuvbs8XLZrDrQHJTXH5K3Nf7HH45/DSkUMRSJGjrcEtLhllCfvk/HwKy4UzjWRZiW9NzjIDADAKC/EXgBAI7a2/H+XWdOKJPdRntIAOYYWejWgjnHdLvfMAz5Q9H2MMwfkrc1nPQ8nHZ7Y2tI4ejRB2Z2q+UIwVjn7clLM90OK4EZAABHQOAFADhq72yJLR2axXJGAFnMYrHI47TJ47SpvMjd69cbhqHWUKRDGNbz0KyxNaRI1FA4aqjeF1S9L9in9+GwWXpQTdZ9lRmBGQBgKCDwAgAcFX8oolU76iVJZ0+mYT2A3GWxWJTntCvPaVdF8dEHZo2tHUKzjkFZF6FZJGooFDF0yBfUof4IzBKPHR2WbbaHZi47gRkAIPsReAEA0tp1qEUf7m1QJL6MZ+ehFn28r1Hr93nVFIj98BUIRzWy0KWJIwrMHi4AZK1MBGYtwUinIKyxi15l3YVmUUNHHZg5bdZE1VhhD6vMkgMzt8PWp+8LAEBvEHgBALoVDEd12a/e7tGym8tOGcVv/AGgH1ksFuW77Mp32VXZh/uDGIYhXzDSTVVZaq+yrkKzJn8sMAtGojrYHNTB5j4GZnbrEZdgjirx6KITKuW00xcSANA3BF4AgG59uKdB9b6gPA6bThlTIpvVohGFLk0bVazjq4o1vMApm9Uip92qij70wwEADByLxaICl10FLruq5On166NRQ75gOCUYO9AU0J7DrdpzuCXlayAc7fY8wXBUB5sDOtgcSP8Nr5QuPXlUr8cJAIBE4AUASGPF9lhvrnOPGaFHvnaayaMBAGSaP9Re8dXYw+b7TUnbg5Hug62eslktnSq8Kos9msmNUAAAR4HACwDQrbbAa8b4UpNHAgDoij/UuadXurCq4/ZgmkqsnrJa1LmHV9oli7F9xfHneU4bS+IBABlH4AUA6FI4EtXq+N0XZ0zgt+wA0B8C4cgRQ6l0TekzFVgVpg2r4k3nOwRWbY/zCawAAFmIwAsA0KV1+7zyBSMq9jh0THmh2cMBgKzUk8Cq8/b25+l6XfWUxSIVupLviNg5sOqy2ir+ON9pl9VKYAUAyC0EXgCALq3YdkiSdPq4Un4QApCzguFoD8Kq7sMrfyjzgVWPwqq27R6HCgisAADohMALANCllfH+XWdOoH8XgMEhGjX03s7DOtAU6BRWNXZRXZWpwEqSCt1pQqluthfHq6wIrAAAyDwCLwBAJ5GooZVt/bvG078LwODw6ze26b6XPunTa9sqrAq7CauK0wRZBS67bARWAABkFQIvAEAnG/Z71eQPq8Bl17GV9O8CMDjsbWiRJI0e5tHUiqL0SwHbel25HSpwE1gBAJBrCLwAAJ2siC9nnD5umOw2q8mjAYDeueLU0brlgilmDwMAAJiIwAsAICm2jHH9Pq+2HmjWc+/vlcRyRgAAAACDE4EXAECSdPMf1uovH+xL2XbWRAIvAAAAAIMPgRcAQJtrm/SXD/bJYpFOH1eqiSPyNX1sqU6qLjF7aAAAAADQawReAAD9+o1tkqQ5x5Xr11+fbvJoAAAAAODo0IkYAIa4fQ2tiZ5d3/n0RJNHAwAAAABHj8ALAIa4x97arnDU0JkTSnXKmGFmDwcAAAAAjhqBFwAMYQ0tQT29cpck6V/PnWTyaAAAAAAgM+jhBQA5rrElpE11TdpS1yxfIKxAOCpfIKwtdc36eJ9XLcGIjqss0jmTh5s9VAAAAADICAIvAMhRdU1+XfWbFdpc15z2OKtFuuWCKbJYLAM0MgAAAADoXwReAJCjHn19WyLsqip2a3J5oYblOeS0W+Vx2DRueL6OKS/UMRWFKitwmTxaAAAAAMgcAi8AyEENLUE9Fe/N9dtvTtdnppabPCIAAAAAGDg0rQeAHPTE8p1qCUZ0bGWRzjtmpNnDAYABtftwiwzDMHsYAADARAReAJBjWoJhPf72dknSv547kd5cAIaMU6qHSZKeXbNX3/6f1WpsDZk8IgAAYBYCLwDIMX9YtVuHW0IaU5qnz55QYfZwAGDAXH7qKN192Qly2qx6eX2tPv/Lt/TxvkazhwUAAExADy8AGOQMw9ALH+zTqh31OtgU1LvbD0mSvv3pCbLb+L0GgKHDYrHoa2eO1Ymji/Wv/7tGOw+16F9+9Y7u+vzx+vLp1VS8AgAwhBB4AcAgVuf16wd/+lCvbTyQsr2iyK0rTh1t0qgAwFwnji7R3753tv7tjx9o2Sd1uu3Zj7Rqx2HdfdkJ8jhtZg8PAAAMAAIvABhEDMPQnsOt2lTbpE9qmvSbN7epoSUkp92qr80Yq7FleRpe4NJpY4fJ7eCHOgBDV0meU7/5xnQteWOr/vMfG/WnNXv08b5G/eqqUzVhRIHZwwMAAP2MwAsABoHd9S16ds1ePfv+Hu081JKy7/iqIv3iyydrSnmhSaMDgOxktVp0w7mTdEr1MH336ff1SU2TPv/Lt3XfFSfq4hMrzR4eAADoRwReAJCFttQ16bG3tmtrnU87632q9QYS+5w2qyaMyNfk8kJNHztMXzljjJx2enUBQHdmTizTi987W999+n2t2F6v+U+tUa33OF1z9nizhwYAAPoJgRcAZBHDMPTE8p2698UNCoSjie0WizRzQpm+cNpoXXhChfKc/PUNAL0xssitr84YoxXb6yVJq3bUE3gBAJDD+IkJALLEoeaA/v2ZD/RqvAH9pyYP1xdOG62xZfkaV5ankjynySMEgMEpFInqp3//RI+9tV1S7BcIP7nsBJNHBQAA+hOBFwBkgY/2NOo7/7taexta5bRbtfCiqZo3c5ysVovZQwOAQa2uya8bn3xfK3fEKru+8+mJ+vc5U2S3sRQcAIBcxr/0AGCy/7d6j65Y8o72NrRq/PB8vXDjLF09azxhF4ABtWjRIp1++ukqLCzUyJEjddlll2njxo0px/j9fs2fP19lZWUqKCjQFVdcodraWpNGfGSrdtTr4gff0sod9Sp02fXrr5+m2y6aStgFAMAQ0Kt/7XPxQggAzPTm5gP692c+UDAc1flTR+q5+bM0taLI7GEBGIJef/11zZ8/X++++66WLl2qUCikOXPmyOfzJY655ZZb9Je//EXPPPOMXn/9de3bt0+XX365iaPummEYeuyt7frKo+/qQFNAx5QX6vkbZ2nu8RVmDw0AAAyQXi1pbLsQOv300xUOh/XDH/5Qc+bM0fr165Wfny8pdiH0t7/9Tc8884yKi4t144036vLLL9fbb7/dL28AAAYrwzB030ufSJK+PL1aiy6fRlUXANO89NJLKc9/97vfaeTIkVq9erXOOeccNTY26rHHHtNTTz2lz3zmM5Kkxx9/XMcee6zeffddnXnmmZ3OGQgEFAi032XW6/X275uQ5AuEdeufPtRfP9wvSbr05CotunwaN/sAAGCI6dW//P1xIQQAQ9VL62q0bq9X+U6bbr1oKmEXgKzS2NgoSSotLZUkrV69WqFQSLNnz04cM3XqVI0ZM0bLly/v8jpv0aJFuvPOOwdmwJIaWoL6wpLl2lLXLLvVots/d5y+MXOsLBb+fgUAYKg5ql91ZeJCqLvf/IVCIYVCoaMZXpfaztkf584FzE96zE96zE96yfMTiRr62T9iS8KvPmusCp2WIT9vfH7SY37SY14yKxqN6uabb9asWbN0wgmxuxnW1NTI6XSqpKQk5djy8nLV1NR0eZ6FCxdqwYIFieder1fV1dX9Nu6/fbRfW+qaNbzApV9//VSdNra0374XAADIbn0OvDJ1IdTdb/5efvll5eXl9XV4R7R06dJ+O3cuYH7SY37SY37SW7p0qVbWWbTtoE15dkOjfZv04oubzB5W1uDzkx7zg4Ewf/58rVu3Tm+99dZRncflcsnlcmVoVEfWEohIks6ZPJywCwCAIa7PgVemLoS6+83fnDlzVFSU+cbNoVBIS5cu1QUXXCCHw5Hx8w92zE96zE96zE/XolFDh1uC2nvYp+dfXSH78HH6R22tpKDmf2aKrvjUeLOHmBX4/KTH/KQXCoX0/PPPmz2MnHDjjTfqr3/9q9544w2NHj06sb2iokLBYFANDQ0pv9ysra1VRUV2NIMPRqKSJAd3YQQAYMjrU+CVyQuh7n7z53A4+vWCvr/PP9gxP+kxP+kN9fnZ29Cq597fqw37vVq/36tdh1oUjhrxvTZp825J0qgSj645e6IcDpt5g81CQ/3zcyTMD/qLYRj67ne/qz//+c967bXXNH58ahh/2mmnyeFwaNmyZbriiiskSRs3btSuXbs0c+ZMM4bcSagt8LLTswsAgKGuV4FXLlwIAUB/MgxD3/mf1fpob2OnfWX5ThVaAjrruGpNGz1Mc4+vkMdJ2AUgO8yfP19PPfWUnn/+eRUWFibaURQXF8vj8ai4uFjXXnutFixYoNLSUhUVFem73/2uZs6cmTU3JgpR4QUAAOJ6FXjlwoUQAPSndXu9+mhvo5w2q26+YLKOqyzSpJEFKi9yS9GIXnzxRX32s8dRoQMg6zzyyCOSpHPPPTdl++OPP65vfvObkqRf/OIXslqtuuKKKxQIBDR37lz96le/GuCRdi8UiVXTOgm8AAAY8noVeOXChRAA9Kffr9olSbrwhArdcO6klH2haMSMIQFAjxiGccRj3G63Hn74YT388MMDMKLeC4ZjFV5rdh3WI69tVZ7TJo/Tprz4H4/D3v7YaVOeM/bcZbfKYmEZJAAAuaTXSxqPJNsvhACgv7QGI3ph7T5J0pWnV5s8GgAYetzxnoirdhzWqh2He/w6q0XyOGzyONsDsbx4INZVYJayzWlXnqNziOZJvMZGmAYAgAn6fJdGAECqFz/ar6ZAWGNK83TmhDKzhwMAQ841Z4+Ty27V4ZagWoIRtQYjagmGY49DkU7bAvGKsKgh+YIR+YL9U4nrSQnEYiFZfofArFOIlhSYdQre2oI1h01WK2EaAABdIfACgAz5w6rY3Re/fHo1P4AAgAlGFrp1ywVTenx8JGqoJRiOh2DxQCwUbn+c2B4/JtQhRGvbH4qoNRiWL9AWrIXlD0UT36c1FNsuX+bfs8tu7VSN5nHYlO+KP3d0EaIdYYln2+vs9EIDAAxiBF4AkAFbDzRr5Y56WS3SF04bbfZwAAA9YLNaVOh2qNCd+RuJRKNGalVZKNx1iBZsD8lS93fYFko9vq3TSCAcVSAc1eGWUMbfg9NmTapKi1eYOezKc3UXmB15iWeeIxbEOe2EaQCA/kXgBQBHYeuBZv39o/169v29kqTPTB0ZuyMjAGBIs1otynfZle/K/OW2YRjyh6KdlmumhGhtgVkoNWRr6bi/i+We0XiYFoxEFWyNqrE182Ga3WpJWqZp77TsMz+5Yq2LJZ7d9VLzcBMCAEAcgRcA9MF7O+r1X8s2683NBxPbPA6b/rXDnRkBAMg0iyUWFnmcNmW6Y6RhGAqEo0lLOLtZ4hmKxJZwtu3vZrlnx+q1cDxNC0cNNfnDavKHJQUy+h6sFnVa4pm87DO/qyWeHXqlpVSnJb3e7SBMA4DBgsALAHrIH4rolQ21evLdXVq+7ZCk2HKYT00erotOqNAFx1WoNN9p8igBAOg7i8Uit8Mmt8OmYf1w/mAiTOt6iacv2CFk62K5Z1evbw1GFIy034SgORBWcyCc8fFbLB1uQuDoeEOB+JJOV2r12qljh+nk6pKMjwcA0D0CLwA4gt31Lfr1G1v1wtp98vpjF88Om0VfOK1aN5w7UdWleSaPEACAwcFpt8ppt6pYme+bFopEYzcI6LCEszkQVn1zUId8AR1qDupg0uNDzQEd9AUVDEeP/A0kGYYSYVxv5Dlt+vDHc7gRAAAMIAIvAOhGXZNfD/9zi55auUuhSGwJRlWxW/9y6ih9dcZYjSrxmDxCAAByQzgS7dRvrKu7Z6Zs62I5ZVdN/wM9DLOOltthTVR05btSG/jPnFhG2AUAA4zACwA6CEei+t07O/TA0k2J3+DOmlSmG86dpJkTymS10rsDADD0HGk5Ynv41EU4Feq6YX7H5Yj9qeNyxJTG+F3ccbJzM/2u7zjZ1ueL6wMAyC4EXgCQ5MM9DbrtTx9p/X6vJOmk6hLdOvcYnTVpuMkjAwAgvd40nO9YJdVdYOXrouF8f+rYcD6v0x0aaTgPAOgZAi8AiPv7R/t10+/XKhiJqtjj0MKLpupL06v5jS0AIGMMw5A/FG1vxB7qXCWVqIAKdb2Er2Ng5Yu/pjUUUWQAQim71dJ1IJW0hM/j7CKccnRRJZX0Oo/TJpedUAoAkBkEXgAg6Q+rdmnhsx8pakizjy3XT6+YpuEFLrOHBQAwQTRqJEKlIy3h6+ough17Tvk6HG/0fyYlp82aUvGU7+xwN8EulvB1Dqy6XsLntNOLCgCQ/Qi8AOS8SNTQgaaA9jW26mBTQOGooUjUUGswohqvX9sP+vTn9/dKkq48vVr3/Ms02ajqAoBBJxo1tHrXYdX7gp3DqVBq76iOS/h8gfbwyh8amCbnLru1iyqppOV53S3h6zKwSgqnHDYapAMAhjwCLwCDzt6GVm2ta5YhKWoY2nWoRR/va9TGmqZEk/moEQu0mgOxpR49WeLxnU9P1K0XHsNSCgAYpH779nbd/bcNA/b9Clx2Dct3qDTfpdI8h4blO1Wa51RpQexrSZ5TBa4OlVVJzdD55QoAAP2HwAvAoLLjoE8XP/imfPFgq6dsVovKC10aUeiS026VzWqR22FTeaFb5cVunVJdovOmjuynUQMABsJpY4fppNHFamgNyR+KyB+Kyh+KKBDun4qt5kBYzYGwdte3dnuM02aV22GV22GL/0l+bJPbHnvuSdrnij/2JL/GbpPbaYt9TTpHyuvofwUAQAKBF4BBIxyJ6uY/rJUvGNGIQpdGFLhksUjlRW4dX1Wk4yqLVJznkCRZZFGe06YCt12FLrvKClz8Jh0ActwpY4bp+RvP7rQ9Go3dvdAfisgfji1n9Iei8ocjsUAsFFVrKJISkvnDEfmDEfnbXheKqDXU/jgQf31rMH5sW7gWiioYaQ/YgpHYc68/PCBz4LLHend1DMYSj+PbYw3iU/elhme2+HlSz+FJCuScNgI2AED2IvACMGj86rWtWru7QYVuu56bP0ujSjxmDwkAMAhY43cV9DhtA/L9IlEjEYz5w9F4wBZRICkYa00O1+JVaG3HtQVoraGIAh1CuLawLpB0TPKy/UA4Gq9oC/X7+7RYJHeHYMzlsMnTMWjrsjotOTyLvd7jtKWEch3DOgd9yQAAvUDgBWBQ+GB3g/5r2WZJ0k8uPYGwCwCQtWxWi/JdduW7BuZSOxSJplanhVIr2PxJFWyBDs9Tj+viHEnBXVtQ13aXScOQWuPbBoLNaklUobnsHSvT4tvjgZnHGQ/aulpGGg/hYgGbNaXSzZN0HJXhADC4EXgBMJ1hxO6iWOsNqNbrV70vqIhhyDCkmsZWvb31kD7Y3aBI1NDFJ1bq0pOrzB4yAABZw2GzymGzqtDd/9/LMAwFI9FYlVmn8Cx1CWhyaJZc0ZaoTkssB00f1rWJRI1437T+f5+S5LBZOleqxQM2V3KIllKdZo1XuXUM25KP7xzWuexWWQnYACCjCLwADKgP9zRoU22z6pr8qmn0a2NNkz6paVJj65GXXpxcXaJ7LjuBfiEAAJjEYrHIZY/1+JLH0e/fzzCS+q+lhGdJwViHPmqpy0Hbj2tN6skW6GJZqT8cVTDpBgehiKFQJKymAeq/5rRbU/qopevD5nHGQ7dECJd6XMfQzeOwafSwPKrWAAwpBF4ABoRhGLrvpY1a8vrWLvdbLdKIQpfKi9wqzXfKbrXKapEK3HbNGF+qsyYOV3Vp3gCPGgAAmMliaa+yGgiRqJHSay25Oi2Q6KOWdGODUNcVbW292wIdlpV6W0NqaAml3NigTTAeuDV2f9PPo/LpKSP039ec0T8nB4AsROAFoN+FI1H9x5/X6Q/v7ZYkzZxQpsoSt8qL3Jo0okDHVhZp4sj82G+LAQAABoBhGApFjPR364wHXoGkZZatwdQll4Guqs66CsLC7f3PBprDZtEYfnEIYIgh8ALQb+p9Qb277ZD+sGq3Xt90QFaLtOjyafry6WPMHhoAAMhC4Ui02z5gnZY0huNhUxdLGtuXKaa+LtAhyIqaFEDZrJZE7y8a8ANA/yDwApBx9b6g/v1Pa/Tm5oOJbU67VQ9eeYouPKHCxJEBAIDeiEaNTqFRcmVTT5f4JZYFdnHnyOQqqbBJCZTFovbgyB4LlVxd9MZy221yOzv217LGG9l31cS+6z5cDpvVlPcJAEMJgReAjKoPSF/5vyu17WCLJOmY8kLNnFimL5w2WieMKjZ5dAAADG4dm7i3dgyXetjEPblKqtsm7qFol72mBoqrLThK07zd7UiukupYKRUPrjoEWZ1DKqucNis3xQGAHEPgBeCo+EMR7W/0yxcIq6bBp8XrbGoMtqiq2K3fXXOGppQXmj1EAAD6jWEYCkaiHXo5pd75r2Pzcn+H4wJtdxAMpfaG6hRkxR+bxWmzdh08JVU2dQqbOi2vS16el3qO1JCKAAoAcHQIvAD0yGFfUM+s3q1DvqCa/WEdaApoc12zdhzydWjAatHEEfn632/NUGWxx6zhAgCGsFCki15OKU3HUyubOi6v67pKKqLWtmV8SZVRrSHzGpHbrZaUKqWOvZvaekN5uguoHKm9oTqeK6WflJ0+UACAwYXAC0CP3P+PT/T0yt1d7st32lTgtivfaVOx0axff+t0jSTsAgAMsNc21um7T7+vJn/YlO9vtSilUbirQ0Pxriub2o9xJVVJdax46qqayk4fKAAAukXgBeCIAuGI/vbhfknSF08brVHDPCrxODRpZKGmVBRoRIFLFotFoVBIL774ooblOU0eMQBgKFq/39sp7PJ0CIu66gvVqUqq7c548ceuLiqe2l/XfpzDZmEZHgAAWYLAC8ARvbHpoLz+sMqLXPrpFSeypAEAkJU+NWmE7tdGuR1WrfyP2Sp02QmgAAAYoqiDBnBEL3ywT5L0uROrCLsAAFnrhFFFKi9yyR+Kas3Ow4RdAAAMYQReANJqCYb1yvpaSdLnT6oyeTQAAHTPYrHoM1PLJUnLNtSZPBoAAGAmAi8AaS1dX6vWUERjy/J04uhis4cDAEBaFxw3UpK0bEOtDLNunwgAAExH4AUgrb/ElzN+/qQqloYAALLeWROHy+2wal+jXxv2N5k9HAAAYBICLwDdqvP69fqmA5KkS1jOCAAYBNwOm86eNEKS9Nhb2+UPRUweEQAAMAN3aQSGKMMw5AtGtO1Asz7c06h1ext1sDmg5kBYzYGw9jX4Ve8LSpKmVhRqSnmhySMGAKBn/uWUUXplQ63+tGaP3ttZrx9fclyitxcAABgaCLyAIeDplbv0p9V71OSPhVlN/pCaA2FFe9DaZHiBU987f3L/DxIAgAy5+MRKhaMn694XN2jnoRZd87v3dP7UkbrjkuM0tizf7OEBAIABQOAF5Lhf/nOz/vPlTd3uL/Y4dOLoYp00ukSjh3mU57Ir32lTRbFbY0rzVOh2DOBoAQDIjEtPHqXzjy3XQ//crMfe3K5ln9Tpzc0Hdf05E3TDeROV5+QyGACAXMa/9EAOSw67vvPpiTp70nAVuO0qdNtV6LKrwG2Xx2GjGT0AICcVuOxaeNGx+uJp1brzLx/rzc0H9ctXt+jZNXv0o88dp4tOqODfQAAAchSBFzCINflDWr/Pq09qmnSwOaB6X1CHW4Kq9wV1qDmozXXNkqTvzz1G88+bZPJoAQAwx6SRBXrimjP0j49r9ZO/rtfehlbd8OQazZpUpv9zyfGaTJ9KAAByDoEXMAhEooY21jRp7e4Gba5r0o6DPm0/6NOOQy1HfO0PLjxGN5xL2AUAGNosFosuPKFCn54yQo+8vlVLXt+qt7cc0kX/9aaunjVO3zt/Msv4AQDIIQReQJYKRaJatqFOz7y3W+9uOyRfsOvbqlcVu3VcVZEqiz0alu9UaZ5DpQUuleY5NaY0T2PK8gZ45AAAZC+P06YFF0zRF04drbv+ul6vbKjVb97crufW7tMjV52q6eNKzR4iAADIAAIvwGSGYehAU0BbDjRr16EW7W/0a19Dq17fdEB1TYHEcQUuu06qLtZxlUUaNzxf48vyNaWiUMMLXCaOHgCAwWlMWZ6uP2eCNuz3am9Dqw40BfTs+3sJvAAAyBEEXsAAavKH9NrGA1q+7ZD2NbSqptGvvYdb1RQId3n88AKnvnBatT5/UpWOqSiUzUpjXQAAjtZ7O+r1i1c26e0thyRJDptFX5perX+fc4zJIwMAAJlC4AX0E8MwtGZXgzbs92rnIZ8+qWnSim31CkainY61WqQxpXkaNzxflcUeVRW7NaWiUOcdM1JOu9WE0QMAkHtW7zysxa9s0pubD0qS7FaLvji9WvPPm6jRw2gBAABALiHwAvpBMBzVHc+v0+9X7e60b8KIfF1wbLkmjMhXRTzcGlOWJ5fdZsJIAQDIfe/vOqxfvLJZb2w6IKkt6BqtG86dpOpSgi4AAHIRgReQYfW+oL7zv6u1cnu9rBbp01NGaPzwAo0bnqezJg7XpJEFZg8RAIAhYe3uBi1+ZZNe2xgLumxWi75w6mjd+BmCLgAAch2BF9ALB5oCWn3QovoVu9QciKqxNdTpz76GVnn9YRW67Hrwq6fovGNGmj1sAACGlA/3NOgXSzfp1aSg64pTR+nG8yZz92IAAIYIAi+gF771P2u0fr9N2vxJ2uPGlObpsXnTNbm8cIBGBgAADjUH9IP/96GWfVKX2Fbscej6cybouMoi7WtsVUNrUHlOu/KcNuU5bfI4bXLarLJYuDEMAAC5hMAL6KEdB31av79JVhm64LhyDct3qcjjULHHkfja9ufYykJ6cgEAMMB+v2p3StglSY2tIf3sHxvTvs5utcgTD8BSwzC78hw25bna93kcNuW72vclHjtt8X32RJCW57DJbuPmMwAAmIHAC+ihVzbUSpImFhn65VdOlsPhMHlEAAAg2ZdPr1YoElW9L6iWYEQtwXD8a0StwYh8wbBa489bgmGFIoYkKRw11OQPq8kflhTI6JicdmssLHPEQrB8Vyw0y3PalOeKB2rxcC2/LShLCt3ynPZEGJef9NjjsMlqpSoNAIDuEHgBPbRsQ+w3xieUGiaPBAAAdGV4gUs3z57S4+NDkWi3YVh7UBZ77Et63HZ82+Pk17Tti8YvF4LhqILhqBoUyvj79Tjal2V2rE5re9xd5Vp34Vqe0yaXnSWeAIDBj8AL6IHG1pBW7aiXJJ0wjMALAIBc4LBZVeyxqtiT2aptwzAUCEdjYVgoopZAUoAWCssXaAvKwvH97ftagpHY/vjjlkBELaHkMC6S+D6toYhaQxHJl9Hhy2pRSmVZSlgWX7bZtmQzUakW35fntMeWgDraq9PyXTblOWKPnXaWeAIABgaBF9ADr286oHDU0MQR+RrubjR7OAAAIItZLBa5HTa5HTYNy/C5o1FD/nDkiJVmbdVpvmBSuJZ0TGuiaq19XyAcjX0PQ2oOhNUcCGd49JLDZmkPxpzx/miObsK1Dss4O1eupVaq2VjiCQBIQuAF9MCyeP+uzxwzQooQeAEAAHNYrZZ46JP5y/hI1EiEYS2dlnmmLtvsFK4lV66FwomqtbZjwvE1nqGIoVAkLK8/82Gaq61fWltlWVIY1vY8UbnmaK9M61S51kW/NJZ4AsDgQ+AFHEEoEtWr8Ts+fWbqCNV9vMXkEQEAAGSezWpRoduhQnfmb8wTTCzxTFqqGV/S2RqMyBcIqzWUFK4F0uxLCd7a+6UFwlEFwlEdbsl8v7TkqrM8hz1x506PIxaa5bvaH7eFay6HTQ6bVQ6bJf71yI+ddqvsVoscdquc8X1UrgFA3xB4AUfw3o7D8vrDGpbn0CnVJfrHx2aPCAAAYHBx2mNhTrH6p19ax6WaLd3ehCDdDQpS97WG2vuldeyfNpAslli/OWc8GLMnPXbYrPHn7eGZ3WZJhGUOu1UOa3yf/cjBmzP++sTz+Gvs1vbHnV/f+XxUxAHIBgReQAfhSFRL19fqn5/U6XBLSFsPNEuSzps6kt+wAQAAZJHkfmml+c6MnjsaNRKVZW3Vack3HGgNxW4wkAjXkm9QEIooEIoqHI0qFIkqFDYUjMSfhw2FItHY80j741D8edvyzzaG0X63z8HC3hay2SzxqrX2wK1TqJb0PDnUcyQFeSmhXrwKzmnvPnCjig6AROCFISoQjqjJ3/4bvHpfUAebA9pc16w/rtqtGq+/02s+d2KlCSMFACD7PPzww/rZz36mmpoanXTSSXrooYd0xhlnmD0sIKOsVovyXXbluwb2R6Zo1FAoGlUoYiicCMMMhcKxwCwYD8xCbdsTj1OfB+Ovb9seDLcFcLHHoW4Ct2DSa8Lx84S62h+OKhQ1ugziwlFD4WhErZlfXdovzK6is1vTV9SlfD+q6IAeI/BCzjIMQyu312vDfq+a/GE1BcLaecinTbXN2nnIpw6/PEtRlu/UFaeN1tiyPBW6HRpV4tGpY0oUDme+wSoAAIPJH/7wBy1YsEBLlizRjBkztHjxYs2dO1cbN27UyJEjzR4eMOhZrRa5rDYNcM7WZ4ZhKBI1YmFbNB6ERTqHcp0q2sJRhaPx5/HXxAK99td0DNzSBXyhcFtQ2F5FF+qmoi7XqugSQRhVdECKQfLXKNBzB5oCeuGDfXpyxU5tO+BLe6zHEWssWpLn0PB8l0YUuTTnuHJdeEKFXHbbAI0YAIDB44EHHtB1112nq6++WpK0ZMkS/e1vf9Nvf/tb3XbbbSnHBgIBBQKBxHOv1zugYwXQ/ywWi+w2i+w2yaPBcf3cVkWXWuHWXhHXVkXXm4q6jo8TFXVHWMIaSgr1wm2v7VBRRxUdVXToGwIvZC3DiPUwCEfa/0EKR2Kl023/GPhDEdU0+rWvsVWbapv07rZ6balrTpwj32nTrEnDNSzPqUK3XZUlHh1TXqgp5QUaXuCSld80AADQY8FgUKtXr9bChQsT26xWq2bPnq3ly5d3On7RokW68847B3KIAHBEg7WKLhw1EtVxVNH1v3RVdA5r9+Fb5wq7rivq2nrcUUXXfwbJ/+IYarbUNelb//2edhxq6dPrp40q1pdPr9Zlp4xSwWD5lwwAgCx38OBBRSIRlZeXp2wvLy/XJ5980un4hQsXasGCBYnnXq9X1dXV/T5OAMglyVV0bgdVdP1RRReKRGV0aHlDFV3vquicNptOGVMy4H0P0+m3kdDMFH3lC4T1nf9d02XYZbVIdlsszbZbLXI5bKosdquy2K3qYXk6fXypzhhXqmEZvksPAADoPZfLJZfLZfYwAAADLFer6LqrqGuroutUUZd0E4lQJCp//M6vvkA4cRdYX9vdXYMRtQbDaglFOoVvPXsP5lfRnTWxTE9dd6Zp37+jfvn40cwUfWUYhhY++5G21DVrZKFL/+87Z6mswBlLj61WliACAGCi4cOHy2azqba2NmV7bW2tKioqTBoVAGAoMgwjsUwzFE6+w+gR7mCaVF2W/LrUu5R2sZyzh1VkXd7RNN6mpy9BVjaxWS2xSi9rvCos5WYIVl10QnZdC/RL4NWbZqZAG8Mw9N/v7NALH+yTzWrRL796qsaU5Zk9LAAAEOd0OnXaaadp2bJluuyyyyRJ0WhUy5Yt04033mju4AAAR6VtGWIoYsRCmuTHXYY9RrxPWDzsiR8XiiY9jgdB4aQQKdhxaWH82HDUiFdGHbm6qu3xYJfcJ6ytn5e9Q4jU8Q6bnZYYxvt4tfX06q5vWLfLGO3t53Mm9xdLOq5tnIOtT1jGA6/eNjPt7u49oVBIoVDmF8u2nbM/zp0LBnp+fIGw1u3z6rVNB/Xy+lrtqm+VJH1/zmSdMrow6/478flJj/lJj/lJj/lJj/lJj3kZOAsWLNC8efM0ffp0nXHGGVq8eLF8Pl/iF50AgJhItItlcN2EO+HkEKnLqqSul8h1bjp/5EbzXVUphSKxJX2DXSIUssfvxGhrD3Ts1uQm8V00gu/UlyreC8va/rjr3lYd+ltxJ8iskfHAq7fNTLu7e8/LL7+svLz+q+5ZunRpv507Fxzt/AQjUnNY8oUkX9giX9Lj5pDUHJLq/Bbtb5EMtf9P7rAYOqvcUEXDer344vqjfRv9hs9PesxPesxPesxPeswPzPblL39ZBw4c0B133KGamhqdfPLJeumllzpd+wFAJrX1WOrYcDwltGkLd5L6L6Vfgtb13Qa77M/UMWzqMrhKXco22POjjk3Quwx77PFqImv745S7EHYTCjns8WVxtuSqpM7VTSnN17u5Q2Lb6+xWAiSkMr2FXHd375kzZ46Kiooy/v1CoZCWLl2qCy64QA6HI+PnH+yOND+GYWjLAZ8ONAVS7o4RikTlC0b08T6vVu9q0NYDvh5/z8pit6aPLdEFx47UOZOHZ9VdHTri85Me85Me85Me85Me85NeKBTS888/b/Ywhowbb7yRJYzAINfW/6jr3kWdex6Fo6n9j47Y86ibu+h1fp3R+bXh+LK4pP5HwYh5jbgzJdH/qJtwJ+V5N1VCne+W1750LWVZXEo1U8fqoi6qm+ydjxtsy9eAjjKeLPS2mWl3d+9xOBz9ekHf3+cf7JLnxx+KaOn6Wi3bUKu3thzSwebAEV4d47RbNSzPoWF5TpXEvw7Ld2p4vlOl+U5Vlnh0cnWJyovc/flW+gWfn/SYn/SYn/SYn/SYHwAYPLz+kOqbgyk9j1KWrnXT8yjUobKobelax55Hqa/r2V3kEk2+c6D/UVtY09twJ9G7KKn/kcPWdc+jjtVMya/rrudRpz5J3IALMEXGAy+amQ5+UUPa29CqfY2NWvZJnf60Zo8aWtp7o7gdVo0pzUv8o+KM/0bAZbdq0shCnTZ2mE4dU6LSfCclpQAAABiSttQ16bMPvqVgePBUJiX3P4otHUtuaJ1cTdR9lZLDfuSeR517JXVXzZRa2ZT6OpavAUivX9aO0cw0+xmGoe0HfdpxyKcdB1u0q75FOw/54n9sirz7ZsrxlcVuXXbKKJ0zeYROHVsil91m0sgBAACA7JfvsquiyK3dh1tkmFhMlee0qdjjUJHbEfvqsavI40jZVuxxqNBtV57TLo/TKo/DrjynTXlOmzxOm/Kcdpa3ARh0+iXwoplpdjrQFNDbWw7qjU0H9OaWgzrQ1N3SxNja8urSPB1bUaQvnDZa50wZwT9yAAAAQA9VFnv0xg/Ok2EY8oeiagmG1RKMyB+KqCUY+9MaCrc/TmyLqDV+bNu2lvi21lCHY4ORI/a2ajv//kb/Ub0fp90qjyM5BLMpz2FPPPY42rd7nEmBmcOWNkjLc9rkslup1gKQcf3WHZxmpgMrEI5od32rGltDisbvYHLYF9T+Rr/2HG7Vu9sOaf1+b8pr3A6rxg8v0LiyPI0py9PY0nyNLnFq6wcr9JVLL5Lb5TTp3QAAAAC5wWKxyBMPeMr64fzhSFQtoYj8wZ4HaYnt3QRpscdhtYQiieq0YDjWA6yxNZR+QH1gsSg1TOsmSIsFZzblOdoDs3RBmjt+TofNmvExA8h+2Xs7PBxRY2tID7y8Ua9sqNO+xtYelUqfMKpIn5o8Qp+aPFynjR3WaWliKBTS4U9ENRcAAAAwCNhtVhXZrCpyZ/6GJoZhKBCOpoRjRwzSQpGuq9mCse1t52oJRhL9zQyjvRKtPzhslkSlWVso1l2Q1lah5nFY27d3CNKSX+u222hGD2QpAq9B6h8f1+j259apLmlZYr7TprICl2xWi6wWqcjjUFWxR5XFbp0wqlhnTx6u4QWd74gJAAAAAB1ZLBa5HbFKqWH9cP5I1EgEZKnVZZFO25ODtPb97UFaS9I2fzxUi0RjFQGxu1aG5fWH++FdpFanpS777CpI67AsNGkJaEoIF69yc9qpTgP6isBrkPjjqt169M1tavKH1BKMqCn+l/WE4fn6j4uP1YmjSzS8gLsiAgAAABgcbFaLClx2Fbgy/2OpYRgKRqJHDNJaOvRMSxekJb/eH2rvndYaim2XL+NvQ3arpUOQ1rE/2hGCtHiI1l3fNarTkMsIvLKcYRh66J9b9MDSTSnb7VaLrj9ngr53/mS5HdwxEQAAAADaWCwWuew2uew2leRl/vzRRHVaWz+0sLytYXlbQ2psDcnrj39tDSeee1tD8vpjx3hbQ2oKHLniLBw11OQPJwoeMm3G+FI9fd2ZBF/ISQReWaixNaQmf0itwYieXLFLv3tnhyTphnMn6rPTKuVx2jS8wKViT+bX6QMAAACAGYz4zbfC0Vh1VjhiKBSJxv8YCkeiHbbHvoajUQXDhsLRaOr2jsdFogrGzxOKRBWKGgqFo0nfr/3YLl8fNRQMR+PfJ/W4cLQHDZWzUK3Xr6hhyCoCL+QeAi8TRKKGGltDOtwSVG2jX3sbWrXncKs27Pdq3d5G7evilsE/vuQ4XT1rvAmjBQAAADCYRKOGQtGuAp94sBONKhSOHZMa+HQOmHoa+HQXRCW+T7rj2h5Hoz26Eddg4bBZ5LBZZbda5LRbZbda5bBb5LBaY9vj+xPH2axy2izx46xyWC3pj7NZU7cnf5/Eazq/Pvk8IwpcsnMXS+QoAq8BtO1As/79mQ+0ZlfDEY912a3yOG0qzXPqptmTdenJo/p/gAAAAAAkxaqNwlEjKaTpEPIkwp8jhEQdXx8PgcLRbqqVjvh9ujkuEWYZiWbtucBmtcSCnJTgJhbadB/4HOE4W/x8KQGUJR4kdR0wpQRIab5P2+vtVgv9lQGTEXgNkBc+2KeFf/pQvqRb7Ra67BpR5NKoEo9GlXg0aWSBpo0q1nFVRSrsh9sKAwAAAAMtGm0PYkLhaLzap3Mw1HkZWvy4aLTbwCfYTcCUulztSMclBVAdvk8u6TrIiT+PBz926xECn+TjUgKfI1cexSqWjhQstQVb7d+P3lIA+orAqx81+UN6beMB/eWDfXp5fa2kWFPA+79woqpKPHJQOgoAAIAeaKs2ShcMJVcExYKlIwc+wW6Cn7QVSt1UHnV3XA4VG8lmtSQFP0cKfFIrglKXqHUTDCUqj5IConjAlDYY6qryKClgslFtBGAIIvDKEK8/pMff2qGlG2rUEowoGI6qzhtQMBK7Xa3FIs0/d5Junj2ZNdIAAAAmiXSx/Ky3TbDT9jrqSRPsboKobs8dza1qI4tF8aoiSzw06iYk6lj9Y7XKmTb86bC8rIs+SB1Doq6CpW4roag2AoBBhcDrKO2ub9Fz7+/Vb97cJm8Xt4qdMDxfc46v0OdOrNQJo4pNGCEAAMDQsHpnvX78wsfytobbl7vFm3KH4kFXLjXEtncKcnoa+HTTd6iLKqSOAVSXYdARg6jU89gIjQAAA4DAqw8aW0P6+csb9drGA9pV35LYPmlkgb7z6YmqHuaR025Vab5TY0rzKB8GAAAYAM+v3ad1e729ek1btVFKkJOoPOq+IXVPA5+e3oGty+VpHQOmlOVtLFEDACAdAq8+eHbNHj2xfKek2G/WThlToq+dOVafO7GK31gBAACYJBov3/rS9NH62pljU4OoDgET1UYAAOQ2Aq8+2FTbLEn68vRq3X7JcSpwMY0AAADZorLYoxNHl5g9DAAAYCK6p/fBtgOxwGvGhFLCLgAAAAAAgCxD4NUHWw/4JEkTRxSYPBIAAAAAAAB0RODVS42tIR1sDkiSJozIN3k0AAAAAAAA6IjAq5faljOOLHSp0O0weTQAAAAAAADoiMCrl7axnBEAACCrGWYPAAAAmI7Aq5e2xiu8WM4IAACQXWwWiyQpGiXyAgBgqCPw6iUqvAAAALKT1RoLvCIGgRcAAEMdgVcvUeEFAACQnajwAgAAbQi8eiEciWrnoRZJVHgBAABkG1tbhReBFwAAQx6BVy/sOdyqYCQql92qUSUes4cDAACAJCxpBAAAbexmD2Aw2XYwtpxx/PD8xAUVAAAAskPbksaP93r1p9V7NLLIpZGFbo0sdKkkzyGLhes3AACGCgKvXthaR8N6AACAbFXojl3artxRr5U76lP2OW1WjSh0aUShSyMLXSlhWPLjsgJXYmkkAAAYvAi8eqGtwmsiDesBAACyzpWnj1FrKKJd9S060BRQnTeguia/DreEFIxEtbehVXsbWtOew2qRygrioVgiIHPHQzGXRiSFZC67bYDeGQAA6C0Cr15oq/CaQIUXAABA1inOc+jm2VM6bQ+EIzrYHFSd16+6poDqmgI6kPS4rsmvOm9AB5sDihrSgaaADjQF9PGRvp/H0alCbEShSyOL3InAbGSRWwUuLrkBABho/OvbC1sPtFV4EXgBAAAMFi67TaNKPEe86VAkauiQL1YZdiApCEuEYk3t+4KRqBpbQ2psDWlzXXPa8+Y5bfEAzK0RRa7E4+SwbEShS8PoMwYAQMYQePVQQ0tQh3xBSdIEljQCAADkHJvVEg+i3GmPMwxDja2hRABW1+SPB2TxP972582BsFqCEe041KIdh1rSntdhs2hEgUsjkivEkpZTtj0uy3fKbuNm6wAApEPglUYkauhn/9ioD3Y3JPo9VBS5lU9ZOgAAwJBlsVhUkudUSZ5TU8oL0x7bEgynVoklPe7YZywUMbSv0a99jf6057RapNL85Ob7XfcZG1HokttBnzEAwNBEcpPGXz/cpyWvb03Z9ukpI0waDQAAAAabPKdd44bbNW54+hUCwXBUB5oDPe4zdrA59nj9/vTfv7s+YxdNqzziEk8AAAYzAq9uRKOGfvVqLOy68vRqff7kKpUXuTXhCBcrAAAAQG857dZe9xnbc7hV6/c16uN9Xn28z6sab+fKsO76jL2x+aCeuOaMjL4HAACyCYFXN5Z9UqeNtU0qcNm18KJjVZznMHtIAAAAGAKaA+GUSq/knmDJyyIbW0M9PqfVIpUVuBK9wb5x1rj+ewMAAGQBAq8uGIahX766RZL09ZljCbsAAABwVAzD0OGWUIfeXe13fzzgbX/cEoz0+LxOm1Uj4ssUOy5dTH5cSqN7AMAQQ+DVhbe3HNIHuxvkdlh17dnjzR4OAAAAslQ4EtUhXzDRfD75zo3JvbgONAcUihg9Pm++06aRRe72IKuLuzWOLHSp2OOQxWLpx3cIAMDgROAVt36fV+9sPagDTQG9sqFWknTl6WM0vMBl8sgAAAAw0PyhSGIZ4YGugqx4hdYhX0BGz3MsDctzJAKrEfHwqj3UcmlkUawii7uCAwBwdIb8v6SBcET/9cpmLXl9q6JJFytOm1XXnzPBvIEBAAAgowzDiPXHSgqvDiT1yUrumeX1h3t8XpvVouEFzkSA1RZejShypwRZwwucctlt/fgOAQBAmyEXeG2ua9aGBouMj2rUHIzqf9/dqU9qmiRJ50wZockjCzSi0KUzxpeqils1AwAADBq761u0/aAvpbn7gabUqqzWUC/6Y9mt7YFVhyWFI5Iel+Y7ZbOyrBAAgGwy5AKvn728Sa9utEkbPkxsK8136p7LTtBF0ypNHBkAAACOxn8t26z/t3rPEY8rdNlTAqvkBu/JPbOKPHb6YwEAMEgNucBrXFm+RuUd0OjyUhV7nBpTmqfvnDuRXl0AAACD3Pjh+ZpaUdi+tLCLJu8jC93yOFlWCABArhtygdcPLzpGJxtb9dnPni6Hw2H2cAAAAJAh88+bpPnnTTJ7GAAAIAtYzR4AAAAAAAAAkEkEXgAAAAAAAMgpBF4AAAAAAADIKQReAAAAAAAAyCkEXgAAAAAAAMgpBF4AAAAAAADIKQReAAAAAAAAyCkEXgAAAAAAAMgpBF4AAAAAAADIKQReAAAAAAAAyCkEXgAAAAAAAMgpBF4AAAAAAADIKQReAAAAQ9yOHTt07bXXavz48fJ4PJo4caJ+/OMfKxgMphz34Ycf6lOf+pTcbreqq6t1//33mzRiAACA9OxmDwAAAADm+uSTTxSNRvXrX/9akyZN0rp163TdddfJ5/PpP//zPyVJXq9Xc+bM0ezZs7VkyRJ99NFHuuaaa1RSUqLrr7/e5HcAAACQKusCL8MwJMUuqvpDKBRSS0uLvF6vHA5Hv3yPwYz5SY/5SY/5SY/5SY/5SY/5Sa9tfqT2awn03IUXXqgLL7ww8XzChAnauHGjHnnkkUTg9eSTTyoYDOq3v/2tnE6njj/+eK1du1YPPPBAt4FXIBBQIBBIPG9sbJTUf9d5AAAgN7VdO/TmOi/rAq+mpiZJUnV1tckjAQAAg1FTU5OKi4vNHsag19jYqNLS0sTz5cuX65xzzpHT6Uxsmzt3ru677z4dPnxYw4YN63SORYsW6c477+y0nes8AADQF725zsu6wKuqqkq7d+9WYWGhLBZLxs/v9XpVXV2t3bt3q6ioKOPnH+yYn/SYn/SYn/SYn/SYn/SYn/Ta5mf9+vWqqqoyeziD3pYtW/TQQw8lqrskqaamRuPHj085rry8PLGvq8Br4cKFWrBgQeJ5NBpVfX29ysrKuM7LcsxlZjGfmcNcZg5zmTnMZeZ0N5eGYaipqalX13lZF3hZrVaNHj26379PUVERH8Q0mJ/0mJ/0mJ/0mJ/0mJ/0mJ/0Ro0aJauVe/K0ue2223TfffelPWbDhg2aOnVq4vnevXt14YUX6otf/KKuu+66o/r+LpdLLpcrZVtJSclRnbMn+P8kc5jLzGI+M4e5zBzmMnOYy8zpai57W8GfdYEXAAAAMuPf/u3f9M1vfjPtMRMmTEg83rdvn8477zydddZZevTRR1OOq6ioUG1tbcq2tucVFRWZGTAAAECGEHgBAADkqBEjRmjEiBE9Onbv3r0677zzdNppp+nxxx/vVCk3c+ZM/cd//IdCoVDixglLly7VMccc0+VyRgAAADMNuZp/l8ulH//4x53K6xHD/KTH/KTH/KTH/KTH/KTH/KTH/BydvXv36txzz9WYMWP0n//5nzpw4IBqampUU1OTOOarX/2qnE6nrr32Wn388cf6wx/+oP/6r/9K6dFlNj4HmcNcZhbzmTnMZeYwl5nDXGZOJufSYnDvbgAAgCHtd7/7na6++uou9yVfKn744YeaP3++Vq1apeHDh+u73/2ubr311oEaJgAAQI8ReAEAAAAAACCnDLkljQAAAAAAAMhtBF4AAAAAAADIKQReAAAAAAAAyCkEXgAAAAAAAMgpQy7wevjhhzVu3Di53W7NmDFDK1euNHtIpli0aJFOP/10FRYWauTIkbrsssu0cePGlGPOPfdcWSyWlD/f+c53TBrxwPo//+f/dHrvU6dOTez3+/2aP3++ysrKVFBQoCuuuEK1tbUmjnjgjBs3rtPcWCwWzZ8/X9LQ+9y88cYbuuSSS1RVVSWLxaLnnnsuZb9hGLrjjjtUWVkpj8ej2bNna/PmzSnH1NfX66qrrlJRUZFKSkp07bXXqrm5eQDfRf9JNz+hUEi33nqrpk2bpvz8fFVVVekb3/iG9u3bl3KOrj5zP/3pTwf4nfSPI31+vvnNb3Z67xdeeGHKMUP18yOpy7+LLBaLfvaznyWOyeXPDzrjOu/o9eQaEX3z05/+VBaLRTfffLPZQxmU9u7dq6997WsqKyuTx+PRtGnT9N5775k9rEEnEono9ttv1/jx4+XxeDRx4kT95Cc/EfeyO7JMXPcjJhM/I/TEkAq8/vCHP2jBggX68Y9/rDVr1uikk07S3LlzVVdXZ/bQBtzrr7+u+fPn691339XSpUsVCoU0Z84c+Xy+lOOuu+467d+/P/Hn/vvvN2nEA+/4449Pee9vvfVWYt8tt9yiv/zlL3rmmWf0+uuva9++fbr88stNHO3AWbVqVcq8LF26VJL0xS9+MXHMUPrc+Hw+nXTSSXr44Ye73H///ffrwQcf1JIlS7RixQrl5+dr7ty58vv9iWOuuuoqffzxx1q6dKn++te/6o033tD1118/UG+hX6Wbn5aWFq1Zs0a333671qxZo2effVYbN27U5z//+U7H3nXXXSmfqe9+97sDMfx+d6TPjyRdeOGFKe/96aefTtk/VD8/klLmZf/+/frtb38ri8WiK664IuW4XP38IBXXeZnR02tE9M6qVav061//WieeeKLZQxmUDh8+rFmzZsnhcOjvf/+71q9fr5///OcaNmyY2UMbdO677z498sgj+uUvf6kNGzbovvvu0/3336+HHnrI7KFlvUxc9yMmUz8jHJExhJxxxhnG/PnzE88jkYhRVVVlLFq0yMRRZYe6ujpDkvH6668ntn360582brrpJvMGZaIf//jHxkknndTlvoaGBsPhcBjPPPNMYtuGDRsMScby5csHaITZ46abbjImTpxoRKNRwzCG9udGkvHnP/858TwajRoVFRXGz372s8S2hoYGw+VyGU8//bRhGIaxfv16Q5KxatWqxDF///vfDYvFYuzdu3fAxj4QOs5PV1auXGlIMnbu3JnYNnbsWOMXv/hF/w4uC3Q1P/PmzTMuvfTSbl/D5yfVpZdeanzmM59J2TZUPj/gOq+/dHWNiN5pamoyJk+ebCxdunRIXycdjVtvvdU4++yzzR5GTrj44ouNa665JmXb5Zdfblx11VUmjWhw6st1P7rW158RemLIVHgFg0GtXr1as2fPTmyzWq2aPXu2li9fbuLIskNjY6MkqbS0NGX7k08+qeHDh+uEE07QwoUL1dLSYsbwTLF582ZVVVVpwoQJuuqqq7Rr1y5J0urVqxUKhVI+S1OnTtWYMWOG3GcpGAzqf//3f3XNNdfIYrEktg/lz02y7du3q6amJuWzUlxcrBkzZiQ+K8uXL1dJSYmmT5+eOGb27NmyWq1asWLFgI/ZbI2NjbJYLCopKUnZ/tOf/lRlZWU65ZRT9LOf/UzhcNicAZrgtdde08iRI3XMMcfoX//1X3Xo0KHEPj4/7Wpra/W3v/1N1157bad9Q/nzM1Rwndd/urtGRM/Nnz9fF198ccrnE73zwgsvaPr06friF7+okSNH6pRTTtFvfvMbs4c1KJ111llatmyZNm3aJEn64IMP9NZbb+miiy4yeWSDW0+u+9F33f2McCT2/hlO9jl48KAikYjKy8tTtpeXl+uTTz4xaVTZIRqN6uabb9asWbN0wgknJLZ/9atf1dixY1VVVaUPP/xQt956qzZu3Khnn33WxNEOjBkzZuh3v/udjjnmGO3fv1933nmnPvWpT2ndunWqqamR0+ns9D9beXm5ampqzBmwSZ577jk1NDTom9/8ZmLbUP7cdNT2eejq7522fTU1NRo5cmTKfrvdrtLS0iH3efL7/br11lv1la98RUVFRYnt3/ve93TqqaeqtLRU77zzjhYuXKj9+/frgQceMHG0A+PCCy/U5ZdfrvHjx2vr1q364Q9/qIsuukjLly+XzWbj85Pkv//7v1VYWNhpeflQ/vwMJVzn9Y/urhHRc7///e+1Zs0arVq1yuyhDGrbtm3TI488ogULFuiHP/yhVq1ape9973tyOp2aN2+e2cMbVG677TZ5vV5NnTpVNptNkUhE99xzj6666iqzhzao9eS6H33T3c8IPTFkAi90b/78+Vq3bl1KjypJKT1gpk2bpsrKSp1//vnaunWrJk6cONDDHFDJv+E48cQTNWPGDI0dO1Z//OMf5fF4TBxZdnnsscd00UUXqaqqKrFtKH9u0HehUEhf+tKXZBiGHnnkkZR9CxYsSDw+8cQT5XQ69e1vf1uLFi2Sy+Ua6KEOqCuvvDLxeNq0aTrxxBM1ceJEvfbaazr//PNNHFn2+e1vf6urrrpKbrc7ZftQ/vwAR6u7a0T0zO7du3XTTTdp6dKlnf5uQu9Eo1FNnz5d9957ryTplFNO0bp167RkyRICr1764x//qCeffFJPPfWUjj/+eK1du1Y333yzqqqqmEtknXQ/I/TEkFnSOHz4cNlstk530qutrVVFRYVJozLfjTfeqL/+9a969dVXNXr06LTHzpgxQ5K0ZcuWgRhaVikpKdGUKVO0ZcsWVVRUKBgMqqGhIeWYofZZ2rlzp1555RV961vfSnvcUP7ctH0e0v29U1FR0amhcjgcVn19/ZD5PLX9Q7Zz504tXbr0iL+5mTFjhsLhsHbs2DEwA8wiEyZM0PDhwxP/P/H5iXnzzTe1cePGI/59JA3tz08u4zov83pzjYiurV69WnV1dTr11FNlt9tlt9v1+uuv68EHH5TdblckEjF7iINGZWWljjvuuJRtxx57bKLlCHru+9//vm677TZdeeWVmjZtmr7+9a/rlltu0aJFi8we2qDWk+t+9E5vf0boypAJvJxOp0477TQtW7YssS0ajWrZsmWaOXOmiSMzh2EYuvHGG/XnP/9Z//znPzV+/Pgjvmbt2rWSYv/gDDXNzc3aunWrKisrddppp8nhcKR8ljZu3Khdu3YNqc/S448/rpEjR+riiy9Oe9xQ/tyMHz9eFRUVKZ8Vr9erFStWJD4rM2fOVENDg1avXp045p///Kei0WgiLMxlbf+Qbd68Wa+88orKysqO+Jq1a9fKarV2Wso3FOzZs0eHDh1K/P801D8/bR577DGddtppOumkk4547FD+/OQyrvMypy/XiOja+eefr48++khr165N/Jk+fbquuuoqrV27VjabzewhDhqzZs3Sxo0bU7Zt2rRJY8eONWlEg1dLS4us1tQYwGazKRqNmjSi3NCT6370XF9+RujKkFrSuGDBAs2bN0/Tp0/XGWecocWLF8vn8+nqq682e2gDbv78+Xrqqaf0/PPPq7CwMLGuuLi4WB6PR1u3btVTTz2lz372syorK9OHH36oW265Reecc86QuJ3yv//7v+uSSy7R2LFjtW/fPv34xz+WzWbTV77yFRUXF+vaa6/VggULVFpaqqKiIn33u9/VzJkzdeaZZ5o99AERjUb1+OOPa968ebLb2/8aGYqfm+bm5pTqte3bt2vt2rUqLS3VmDFjdPPNN+vuu+/W5MmTNX78eN1+++2qqqrSZZddJin228kLL7xQ1113nZYsWaJQKKQbb7xRV155ZcpS0cEq3fxUVlbqC1/4gtasWaO//vWvikQiib+LSktL5XQ6tXz5cq1YsULnnXeeCgsLtXz5ct1yyy362te+lhO3Ik83P6Wlpbrzzjt1xRVXqKKiQlu3btUPfvADTZo0SXPnzpU0tD8/Y8aMkRS7mHzmmWf085//vNPrc/3zg1Rc52XGka4R0XOFhYWdep/l5+errKyMnmi9dMstt+iss87Svffeqy996UtauXKlHn30UT366KNmD23QueSSS3TPPfdozJgxOv744/X+++/rgQce0DXXXGP20LLe0V73o93R/ozQY328c+Sg9dBDDxljxowxnE6nccYZZxjvvvuu2UMyhaQu/zz++OOGYRjGrl27jHPOOccoLS01XC6XMWnSJOP73/++0djYaO7AB8iXv/xlo7Ky0nA6ncaoUaOML3/5y8aWLVsS+1tbW40bbrjBGDZsmJGXl2f8y7/8i7F//34TRzyw/vGPfxiSjI0bN6ZsH4qfm1dffbXL/5fmzZtnGEbsFsW33367UV5ebrhcLuP888/vNG+HDh0yvvKVrxgFBQVGUVGRcfXVVxtNTU0mvJvMSzc/27dv7/bvoldffdUwDMNYvXq1MWPGDKO4uNhwu93Gsccea9x7772G3+83941lSLr5aWlpMebMmWOMGDHCcDgcxtixY43rrrvOqKmpSTnHUP38tPn1r39teDweo6GhodPrc/3zg864zjt6R7pGxNH59Kc/bdx0001mD2NQ+stf/mKccMIJhsvlMqZOnWo8+uijZg9pUPJ6vcZNN91kjBkzxnC73caECROM//iP/zACgYDZQ8t6mbjuR8zR/ozQUxbDMIyex2MAAAAAAABAdhsyPbwAAAAAAAAwNBB4AQAAAAAAIKcQeAEAAAAAACCnEHgBAAAAAAAgpxB4AQAAAAAAIKcQeAEAAAAAACCnEHgBAAAAAAAgpxB4AQAAAAAAIKcQeAEAAAAAACCnEHgBAAAAAAAgpxB4AQAAAAAAIKf8fwslh2tfDcjsAAAAAElFTkSuQmCC",
       "text/plain": [
-       "<Figure size 640x480 with 1 Axes>"
+       "<Figure size 1500x500 with 2 Axes>"
       ]
      },
      "metadata": {},
@@ -243,13 +236,16 @@
    "source": [
     "import matplotlib.pyplot as plt\n",
     "\n",
+    "plt.figure(figsize=(15, 5))\n",
+    "plt.subplot(121)\n",
     "plt.plot(np.sort(dh.eig()))\n",
-    "plt.ylim(None, 20)\n",
+    "plt.grid()\n",
+    "plt.subplot(122)\n",
+    "DOS = sisl.physics.electron.DOS(np.linspace(-15, 85, 50), dh.eig())\n",
+    "plt.plot(DOS, np.linspace(-15, 85, 50))\n",
     "print(np.real(dh.eig()).min())\n",
     "print(np.imag(dh.eig()).min())\n",
     "\n",
-    "import matplotlib.pyplot as plt\n",
-    "\n",
     "coords = dh.xyz[-3:]\n",
     "\n",
     "\n",
@@ -271,14 +267,14 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 4,
+   "execution_count": 44,
    "metadata": {},
    "outputs": [
     {
      "name": "stdout",
      "output_type": "stream",
      "text": [
-      "Hamiltonian and exchange field rotated. Elapsed time: 2.613376875 s\n",
+      "Hamiltonian and exchange field rotated. Elapsed time: 5448.302713666 s\n",
       "================================================================================================================================================================\n"
      ]
     }
@@ -380,14 +376,14 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 5,
+   "execution_count": 45,
    "metadata": {},
    "outputs": [
     {
      "name": "stdout",
      "output_type": "stream",
      "text": [
-      "Site and pair dictionaries created. Elapsed time: 2.628747208 s\n",
+      "Site and pair dictionaries created. Elapsed time: 5448.330635 s\n",
       "================================================================================================================================================================\n"
      ]
     }
@@ -495,21 +491,21 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 6,
+   "execution_count": 46,
    "metadata": {},
    "outputs": [
     {
      "name": "stderr",
      "output_type": "stream",
      "text": [
-      "k loop:   0%|          | 0/400 [00:00<?, ?it/s]"
+      "k loop:   0%|          | 0/1000 [00:00<?, ?it/s]"
      ]
     },
     {
      "name": "stdout",
      "output_type": "stream",
      "text": [
-      "k set created. Elapsed time: 2.648811208 s\n",
+      "k set created. Elapsed time: 5448.339522166 s\n",
       "================================================================================================================================================================\n"
      ]
     }
@@ -530,14 +526,14 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 7,
+   "execution_count": 47,
    "metadata": {},
    "outputs": [
     {
      "name": "stdout",
      "output_type": "stream",
      "text": [
-      "Rotations done perpendicular to quantization axis. Elapsed time: 2.901047458 s\n",
+      "Rotations done perpendicular to quantization axis. Elapsed time: 5448.55170175 s\n",
       "================================================================================================================================================================\n"
      ]
     }
@@ -593,7 +589,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 8,
+   "execution_count": 48,
    "metadata": {},
    "outputs": [
     {
@@ -601,15 +597,15 @@
      "output_type": "stream",
      "text": [
       "Starting matrix inversions\n",
-      "Total number of k points: 400\n",
-      "Number of energy samples per k point: 100\n",
+      "Total number of k points: 1000\n",
+      "Number of energy samples per k point: 500\n",
       "Total number of directions: 3\n",
-      "Total number of matrix inversions: 120000\n",
+      "Total number of matrix inversions: 1500000\n",
       "The shape of the Hamiltonian and the Greens function is 84x84=7056\n",
       "Memory taken by a single Hamiltonian is: 0.015625 KB\n",
       "Expected memory usage per matrix inversion: 0.5 KB\n",
-      "Expected memory usage per k point for parallel inversion: 150.0 KB\n",
-      "Expected memory usage on root node: 58.59375 MB\n",
+      "Expected memory usage per k point for parallel inversion: 750.0 KB\n",
+      "Expected memory usage on root node: 732.421875 MB\n",
       "================================================================================================================================================================\n"
      ]
     },
@@ -617,14 +613,14 @@
      "name": "stderr",
      "output_type": "stream",
      "text": [
-      "k loop: 100%|██████████| 400/400 [03:18<00:00,  2.01it/s]"
+      "k loop: 100%|██████████| 1000/1000 [39:28<00:00,  2.37s/it]"
      ]
     },
     {
      "name": "stdout",
      "output_type": "stream",
      "text": [
-      "Calculated Greens functions. Elapsed time: 201.17534975 s\n",
+      "Calculated Greens functions. Elapsed time: 7816.996438041 s\n",
       "================================================================================================================================================================\n"
      ]
     },
@@ -728,7 +724,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 9,
+   "execution_count": 49,
    "metadata": {},
    "outputs": [
     {
@@ -759,11 +755,11 @@
       "[0 0 1]  --»  [array([1, 0, 0]), array([0, 1, 0])]\n",
       "================================================================================================================================================================\n",
       "Parameters for the contour integral:\n",
-      "Number of k points:  20\n",
-      "k point directions:  xy\n",
+      "Number of k points:  1000\n",
+      "k point directions:  x\n",
       "Ebot:  -13\n",
-      "Eset:  100\n",
-      "Esetp:  100\n",
+      "Eset:  500\n",
+      "Esetp:  10000\n",
       "================================================================================================================================================================\n",
       "Atomic information: \n",
       "----------------------------------------------------------------------------------------------------------------------------------------------------------------\n",
@@ -781,65 +777,42 @@
       "Magnetic entity1        Magnetic entity2       [i  j  k]       d [Ang]\n",
       "----------------------------------------------------------------------------------------------------------------------------------------------------------------\n",
       "[3]Fe(2)        [4]Fe(2)        [0 0 0]        d [Ang] Not yet.\n",
-      "Isotropic:  -63.6758032703783\n",
-      "DMI:  [-9.33104051e-01 -1.01015977e-03 -2.89340046e-06]\n",
-      "Symmetric-anisotropy:  [-3.45194760e+00  1.37702774e+00  1.29440230e-03 -8.59287160e-04\n",
-      "  1.03571867e-04]\n",
+      "Isotropic:  -29.367092812942317\n",
+      "DMI:  [-8.89792781e-01  4.16580540e+00 -3.59948517e-06]\n",
+      "Symmetric-anisotropy:  [-4.56670448  1.17667365  0.33235705  4.15578143 -0.03089926]\n",
       "Energies for debugging: \n",
-      "array([[-6.22987755e-02, -9.33207623e-04,  9.33000479e-04,\n",
-      "        -6.16008834e-02],\n",
-      "       [-6.77166238e-02,  8.59287160e-07, -1.16103239e-06,\n",
-      "        -6.71277509e-02],\n",
-      "       [-4.95869253e-02, -1.29729570e-06, -1.29150890e-06,\n",
-      "        -4.95886050e-02]])\n",
+      "array([[-0.02819042, -0.00085889,  0.00092069, -0.02597706],\n",
+      "       [-0.03690207, -0.00415578,  0.00417583, -0.0339338 ],\n",
+      "       [-0.03160427, -0.00033236, -0.00033235, -0.03122101]])\n",
       "J_ii for debugging: (check if this is the same as in calculate_exchange_tensor)\n",
-      "array([-0.06712775, -0.04958693, -0.06229878])\n",
+      "array([-0.0339338 , -0.03160427, -0.02819042])\n",
       "Test J_xx = E(y,z) = E(z,y)\n",
-      "-0.06712775087066476 -0.04958860499691915\n",
+      "-0.03393379729617979 -0.03122101431812066\n",
       "\n",
       "[4]Fe(2)        [5]Fe(2)        [0 0 0]        d [Ang] Not yet.\n",
-      "Isotropic:  -60.98280300207994\n",
-      "DMI:  [-3.79822521e+00  6.14860175e+00  5.42292960e-03]\n",
-      "Symmetric-anisotropy:  [0.06335389 0.13420785 0.07070523 6.23730266 0.03599778]\n",
+      "Isotropic:  -74.98944973858333\n",
+      "DMI:  [-6.82212886e+00  1.12891691e+01  5.77778279e-04]\n",
+      "Symmetric-anisotropy:  [-1.03411109  0.63716653 -0.16031601 10.67952816 -0.32649278]\n",
       "Energies for debugging: \n",
-      "array([[-6.08485951e-02, -3.83422299e-03,  3.76222743e-03,\n",
-      "        -6.11803648e-02],\n",
-      "       [-6.06181097e-02, -6.23730266e-03,  6.05990085e-03,\n",
-      "        -6.09194491e-02],\n",
-      "       [-5.99054512e-02, -6.52823043e-05, -7.61281635e-05,\n",
-      "        -5.98235524e-02]])\n",
+      "array([[-0.07435228, -0.00649564,  0.00714862, -0.07459251],\n",
+      "       [-0.07593353, -0.01067953,  0.01189881, -0.07602356],\n",
+      "       [-0.0795701 ,  0.00016089,  0.00015974, -0.07975458]])\n",
       "J_ii for debugging: (check if this is the same as in calculate_exchange_tensor)\n",
-      "array([-0.06091945, -0.05990545, -0.0608486 ])\n",
+      "array([-0.07602356, -0.0795701 , -0.07435228])\n",
       "Test J_xx = E(y,z) = E(z,y)\n",
-      "-0.06091944910723846 -0.059823552399443855\n",
-      "\n",
-      "[3]Fe(2)        [5]Fe(2)        [0 0 0]        d [Ang] Not yet.\n",
-      "Isotropic:  -60.99404713886926\n",
-      "DMI:  [ 3.78607092e+00 -6.13231029e+00  5.42245509e-03]\n",
-      "Symmetric-anisotropy:  [ 0.05074819  0.13690999  0.07070521 -6.22528695 -0.04239585]\n",
-      "Energies for debugging: \n",
-      "array([[-6.08571371e-02,  3.82846677e-03, -3.74367506e-03,\n",
-      "        -6.11817053e-02],\n",
-      "       [-6.06264662e-02,  6.22528695e-03, -6.03933362e-03,\n",
-      "        -6.09432989e-02],\n",
-      "       [-5.99054277e-02, -6.52827577e-05, -7.61276679e-05,\n",
-      "        -5.98235222e-02]])\n",
-      "J_ii for debugging: (check if this is the same as in calculate_exchange_tensor)\n",
-      "array([-0.0609433 , -0.05990543, -0.06085714])\n",
-      "Test J_xx = E(y,z) = E(z,y)\n",
-      "-0.060943298947041494 -0.059823522197803146\n",
+      "-0.07602356083253432 -0.07975458355080561\n",
       "\n",
       "================================================================================================================================================================\n",
       "Runtime information: \n",
-      "Total runtime: 199.56628429199998 s\n",
+      "Total runtime: 2369.9501148750005 s\n",
       "----------------------------------------------------------------------------------------------------------------------------------------------------------------\n",
-      "Initial setup: 0.150397667 s\n",
-      "Hamiltonian conversion and XC field extraction: 0.802 s\n",
-      "Pair and site datastructure creatrions: 0.015 s\n",
-      "k set cration and distribution: 0.020 s\n",
-      "Rotating XC potential: 0.252 s\n",
-      "Greens function inversion: 198.274 s\n",
-      "Calculate energies and magnetic components: 0.052 s\n"
+      "Initial setup: 0.17246508400057792 s\n",
+      "Hamiltonian conversion and XC field extraction: 0.904 s\n",
+      "Pair and site datastructure creatrions: 0.028 s\n",
+      "k set cration and distribution: 0.009 s\n",
+      "Rotating XC potential: 0.212 s\n",
+      "Greens function inversion: 2368.445 s\n",
+      "Calculate energies and magnetic components: 0.180 s\n"
      ]
     }
    ],
@@ -925,7 +898,596 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 10,
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "Magnetic entities integrated.\n",
+    "Pairs integrated.\n",
+    "Magnetic parameters calculated.\n",
+    "##################################################################### GROGU OUTPUT #############################################################################\n",
+    "================================================================================================================================================================\n",
+    "Input file: \n",
+    "Not yet specified.\n",
+    "Output file: \n",
+    "test.pickle\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",
+    "Parameters for the contour integral:\n",
+    "Number of k points:  1000\n",
+    "k point directions:  x\n",
+    "Ebot:  -13\n",
+    "Eset:  500\n",
+    "Esetp:  10000\n",
+    "================================================================================================================================================================\n",
+    "Atomic information: \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",
+    "================================================================================================================================================================\n",
+    "Exchange [meV]\n",
+    "----------------------------------------------------------------------------------------------------------------------------------------------------------------\n",
+    "Magnetic entity1        Magnetic entity2       [i  j  k]       d [Ang]\n",
+    "----------------------------------------------------------------------------------------------------------------------------------------------------------------\n",
+    "[3]Fe(2)        [4]Fe(2)        [0 0 0]        d [Ang] Not yet.\n",
+    "Isotropic:  -29.367092812942317\n",
+    "DMI:  [-8.89792781e-01  4.16580540e+00 -3.59948517e-06]\n",
+    "Symmetric-anisotropy:  [-4.56670448  1.17667365  0.33235705  4.15578143 -0.03089926]\n",
+    "Energies for debugging: \n",
+    "array([[-0.02819042, -0.00085889,  0.00092069, -0.02597706],\n",
+    "       [-0.03690207, -0.00415578,  0.00417583, -0.0339338 ],\n",
+    "       [-0.03160427, -0.00033236, -0.00033235, -0.03122101]])\n",
+    "J_ii for debugging: (check if this is the same as in calculate_exchange_tensor)\n",
+    "array([-0.0339338 , -0.03160427, -0.02819042])\n",
+    "Test J_xx = E(y,z) = E(z,y)\n",
+    "-0.03393379729617979 -0.03122101431812066\n",
+    "\n",
+    "[4]Fe(2)        [5]Fe(2)        [0 0 0]        d [Ang] Not yet.\n",
+    "Isotropic:  -74.98944973858333\n",
+    "DMI:  [-6.82212886e+00  1.12891691e+01  5.77778279e-04]\n",
+    "Symmetric-anisotropy:  [-1.03411109  0.63716653 -0.16031601 10.67952816 -0.32649278]\n",
+    "Energies for debugging: \n",
+    "array([[-0.07435228, -0.00649564,  0.00714862, -0.07459251],\n",
+    "       [-0.07593353, -0.01067953,  0.01189881, -0.07602356],\n",
+    "       [-0.0795701 ,  0.00016089,  0.00015974, -0.07975458]])\n",
+    "J_ii for debugging: (check if this is the same as in calculate_exchange_tensor)\n",
+    "array([-0.07602356, -0.0795701 , -0.07435228])\n",
+    "Test J_xx = E(y,z) = E(z,y)\n",
+    "-0.07602356083253432 -0.07975458355080561\n",
+    "\n",
+    "================================================================================================================================================================\n",
+    "Runtime information: \n",
+    "Total runtime: 2369.9501148750005 s\n",
+    "----------------------------------------------------------------------------------------------------------------------------------------------------------------\n",
+    "Initial setup: 0.17246508400057792 s\n",
+    "Hamiltonian conversion and XC field extraction: 0.904 s\n",
+    "Pair and site datastructure creatrions: 0.028 s\n",
+    "k set cration and distribution: 0.009 s\n",
+    "Rotating XC potential: 0.212 s\n",
+    "Greens function inversion: 2368.445 s\n",
+    "Calculate energies and magnetic components: 0.180 s"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 50,
+   "metadata": {},
+   "outputs": [
+    {
+     "ename": "SyntaxError",
+     "evalue": "invalid syntax (1614865498.py, line 1)",
+     "output_type": "error",
+     "traceback": [
+      "\u001b[0;36m  Cell \u001b[0;32mIn[50], line 1\u001b[0;36m\u001b[0m\n\u001b[0;31m    Magnetic entities integrated.\u001b[0m\n\u001b[0m             ^\u001b[0m\n\u001b[0;31mSyntaxError\u001b[0m\u001b[0;31m:\u001b[0m invalid syntax\n"
+     ]
+    }
+   ],
+   "source": [
+    "Magnetic entities integrated.\n",
+    "Pairs integrated.\n",
+    "Magnetic parameters calculated.\n",
+    "##################################################################### GROGU OUTPUT #############################################################################\n",
+    "================================================================================================================================================================\n",
+    "Input file: \n",
+    "Not yet specified.\n",
+    "Output file: \n",
+    "test.pickle\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",
+    "Parameters for the contour integral:\n",
+    "Number of k points:  1000\n",
+    "k point directions:  x\n",
+    "Ebot:  -13\n",
+    "Eset:  100\n",
+    "Esetp:  10000\n",
+    "================================================================================================================================================================\n",
+    "Atomic information: \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",
+    "================================================================================================================================================================\n",
+    "Exchange [meV]\n",
+    "----------------------------------------------------------------------------------------------------------------------------------------------------------------\n",
+    "Magnetic entity1        Magnetic entity2       [i  j  k]       d [Ang]\n",
+    "----------------------------------------------------------------------------------------------------------------------------------------------------------------\n",
+    "[3]Fe(2)        [4]Fe(2)        [0 0 0]        d [Ang] Not yet.\n",
+    "Isotropic:  -29.36445059067917\n",
+    "DMI:  [-8.91784089e-01  4.16579134e+00 -3.59901207e-06]\n",
+    "Symmetric-anisotropy:  [-4.5685362   1.18301202  0.3323615   4.15576769 -0.02915828]\n",
+    "Energies for debugging: \n",
+    "array([[-0.02818144, -0.00086263,  0.00092094, -0.02597893],\n",
+    "       [-0.03690122, -0.00415577,  0.00417581, -0.03393299],\n",
+    "       [-0.03160383, -0.00033237, -0.00033236, -0.03122058]])\n",
+    "J_ii for debugging: (check if this is the same as in calculate_exchange_tensor)\n",
+    "array([-0.03393299, -0.03160383, -0.02818144])\n",
+    "Test J_xx = E(y,z) = E(z,y)\n",
+    "-0.03393298678958339 -0.031220578858898954\n",
+    "\n",
+    "[4]Fe(2)        [5]Fe(2)        [0 0 0]        d [Ang] Not yet.\n",
+    "Isotropic:  -74.99083333178646\n",
+    "DMI:  [-6.82222254e+00  1.12891474e+01  5.95795483e-04]\n",
+    "Symmetric-anisotropy:  [-1.03214008  0.64046775 -0.16031723 10.67951504 -0.32709547]\n",
+    "Energies for debugging: \n",
+    "array([[-0.07435037, -0.00649513,  0.00714932, -0.07459916],\n",
+    "       [-0.07593293, -0.01067952,  0.01189878, -0.07602297],\n",
+    "       [-0.0795695 ,  0.00016091,  0.00015972, -0.07975399]])\n",
+    "J_ii for debugging: (check if this is the same as in calculate_exchange_tensor)\n",
+    "array([-0.07602297, -0.0795695 , -0.07435037])\n",
+    "Test J_xx = E(y,z) = E(z,y)\n",
+    "-0.07602297341278706 -0.07975398680085131\n",
+    "\n",
+    "================================================================================================================================================================\n",
+    "Runtime information: \n",
+    "Total runtime: 501.2763705420002 s\n",
+    "----------------------------------------------------------------------------------------------------------------------------------------------------------------\n",
+    "Initial setup: 0.12531970899999578 s\n",
+    "Hamiltonian conversion and XC field extraction: 0.888 s\n",
+    "Pair and site datastructure creatrions: 0.062 s\n",
+    "k set cration and distribution: 0.024 s\n",
+    "Rotating XC potential: 0.245 s\n",
+    "Greens function inversion: 499.875 s\n",
+    "Calculate energies and magnetic components: 0.056 s"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [
+    {
+     "ename": "SyntaxError",
+     "evalue": "invalid syntax (1097899960.py, line 1)",
+     "output_type": "error",
+     "traceback": [
+      "\u001b[0;36m  Cell \u001b[0;32mIn[40], line 1\u001b[0;36m\u001b[0m\n\u001b[0;31m    Magnetic entities integrated.\u001b[0m\n\u001b[0m             ^\u001b[0m\n\u001b[0;31mSyntaxError\u001b[0m\u001b[0;31m:\u001b[0m invalid syntax\n"
+     ]
+    }
+   ],
+   "source": [
+    "Magnetic entities integrated.\n",
+    "Pairs integrated.\n",
+    "Magnetic parameters calculated.\n",
+    "##################################################################### GROGU OUTPUT #############################################################################\n",
+    "================================================================================================================================================================\n",
+    "Input file: \n",
+    "Not yet specified.\n",
+    "Output file: \n",
+    "test.pickle\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",
+    "Parameters for the contour integral:\n",
+    "Number of k points:  1000\n",
+    "k point directions:  x\n",
+    "Ebot:  -13\n",
+    "Eset:  50\n",
+    "Esetp:  10000\n",
+    "================================================================================================================================================================\n",
+    "Atomic information: \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",
+    "================================================================================================================================================================\n",
+    "Exchange [meV]\n",
+    "----------------------------------------------------------------------------------------------------------------------------------------------------------------\n",
+    "Magnetic entity1        Magnetic entity2       [i  j  k]       d [Ang]\n",
+    "----------------------------------------------------------------------------------------------------------------------------------------------------------------\n",
+    "[3]Fe(2)        [4]Fe(2)        [0 0 0]        d [Ang] Not yet.\n",
+    "Isotropic:  -29.359042196157713\n",
+    "DMI:  [-8.93813379e-01  4.16574774e+00 -3.59758916e-06]\n",
+    "Symmetric-anisotropy:  [-4.57147523  1.18770348  0.33237646  4.15572505 -0.02850636]\n",
+    "Energies for debugging: \n",
+    "array([[-0.02817134, -0.00086531,  0.00092232, -0.02597527],\n",
+    "       [-0.03689861, -0.00415573,  0.00417577, -0.03393052],\n",
+    "       [-0.03160252, -0.00033238, -0.00033237, -0.03121927]])\n",
+    "J_ii for debugging: (check if this is the same as in calculate_exchange_tensor)\n",
+    "array([-0.03393052, -0.03160252, -0.02817134])\n",
+    "Test J_xx = E(y,z) = E(z,y)\n",
+    "-0.03393051742760642 -0.031219274782673133\n",
+    "\n",
+    "[4]Fe(2)        [5]Fe(2)        [0 0 0]        d [Ang] Not yet.\n",
+    "Isotropic:  -74.98947011735306\n",
+    "DMI:  [-6.82126017e+00  1.12890803e+01  6.51055432e-04]\n",
+    "Symmetric-anisotropy:  [-1.03170277  0.64211178 -0.16032097 10.67947464 -0.32714492]\n",
+    "Energies for debugging: \n",
+    "array([[-0.07434736, -0.00649412,  0.00714841, -0.07459988],\n",
+    "       [-0.07593108, -0.01067947,  0.01189869, -0.07602117],\n",
+    "       [-0.07956766,  0.00016097,  0.00015967, -0.07975216]])\n",
+    "J_ii for debugging: (check if this is the same as in calculate_exchange_tensor)\n",
+    "array([-0.07602117, -0.07956766, -0.07434736])\n",
+    "Test J_xx = E(y,z) = E(z,y)\n",
+    "-0.07602117288317486 -0.07975215648287909\n",
+    "\n",
+    "[3]Fe(2)        [5]Fe(2)        [0 0 0]        d [Ang] Not yet.\n",
+    "Isotropic:  -75.04936117965987\n",
+    "DMI:  [ 6.85701206e+00 -1.13051160e+01  6.39732773e-04]\n",
+    "Symmetric-anisotropy:  [ -1.05318057   0.68103129  -0.16032    -10.70364791   0.31624841]\n",
+    "Energies for debugging: \n",
+    "array([[-0.07436833,  0.00654076, -0.00717326, -0.07467721],\n",
+    "       [-0.07601601,  0.01070365, -0.01190658, -0.07610254],\n",
+    "       [-0.07956796,  0.00016096,  0.00015968, -0.07975245]])\n",
+    "J_ii for debugging: (check if this is the same as in calculate_exchange_tensor)\n",
+    "array([-0.07610254, -0.07956796, -0.07436833])\n",
+    "Test J_xx = E(y,z) = E(z,y)\n",
+    "-0.07610254174486064 -0.07975245001820208\n",
+    "\n",
+    "================================================================================================================================================================\n",
+    "Runtime information: \n",
+    "Total runtime: 261.890385875 s\n",
+    "----------------------------------------------------------------------------------------------------------------------------------------------------------------\n",
+    "Initial setup: 0.1216103749998183 s\n",
+    "Hamiltonian conversion and XC field extraction: 1.045 s\n",
+    "Pair and site datastructure creatrions: 0.011 s\n",
+    "k set cration and distribution: 0.006 s\n",
+    "Rotating XC potential: 0.224 s\n",
+    "Greens function inversion: 260.450 s\n",
+    "Calculate energies and magnetic components: 0.033 s"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [
+    {
+     "ename": "SyntaxError",
+     "evalue": "invalid syntax (1463926646.py, line 1)",
+     "output_type": "error",
+     "traceback": [
+      "\u001b[0;36m  Cell \u001b[0;32mIn[30], line 1\u001b[0;36m\u001b[0m\n\u001b[0;31m    Magnetic entities integrated.\u001b[0m\n\u001b[0m             ^\u001b[0m\n\u001b[0;31mSyntaxError\u001b[0m\u001b[0;31m:\u001b[0m invalid syntax\n"
+     ]
+    }
+   ],
+   "source": [
+    "Magnetic entities integrated.\n",
+    "Pairs integrated.\n",
+    "Magnetic parameters calculated.\n",
+    "##################################################################### GROGU OUTPUT #############################################################################\n",
+    "================================================================================================================================================================\n",
+    "Input file: \n",
+    "Not yet specified.\n",
+    "Output file: \n",
+    "test.pickle\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",
+    "Parameters for the contour integral:\n",
+    "Number of k points:  3000\n",
+    "k point directions:  x\n",
+    "Ebot:  -13\n",
+    "Eset:  50\n",
+    "Esetp:  10000\n",
+    "================================================================================================================================================================\n",
+    "Atomic information: \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",
+    "================================================================================================================================================================\n",
+    "Exchange [meV]\n",
+    "----------------------------------------------------------------------------------------------------------------------------------------------------------------\n",
+    "Magnetic entity1        Magnetic entity2       [i  j  k]       d [Ang]\n",
+    "----------------------------------------------------------------------------------------------------------------------------------------------------------------\n",
+    "[3]Fe(2)        [4]Fe(2)        [0 0 0]        d [Ang] Not yet.\n",
+    "Isotropic:  -29.437239263932703\n",
+    "DMI:  [-8.91096571e-01  4.16386016e+00 -3.68076797e-06]\n",
+    "Symmetric-anisotropy:  [-4.68204753  1.22980695  0.33079224  4.15200152 -0.03180686]\n",
+    "Energies for debugging: \n",
+    "array([[-0.02820743, -0.00085929,  0.0009229 , -0.025985  ],\n",
+    "       [-0.03707978, -0.004152  ,  0.00417572, -0.03411929],\n",
+    "       [-0.03167345, -0.0003308 , -0.00033079, -0.03129215]])\n",
+    "J_ii for debugging: (check if this is the same as in calculate_exchange_tensor)\n",
+    "array([-0.03411929, -0.03167345, -0.02820743])\n",
+    "Test J_xx = E(y,z) = E(z,y)\n",
+    "-0.03411928679618277 -0.03129215156448653\n",
+    "\n",
+    "[4]Fe(2)        [5]Fe(2)        [0 0 0]        d [Ang] Not yet.\n",
+    "Isotropic:  -74.98086703681366\n",
+    "DMI:  [-6.81158085e+00  1.13265508e+01  5.68487147e-04]\n",
+    "Symmetric-anisotropy:  [-1.03582622  0.63764125 -0.15959679 10.72110977 -0.32326625]\n",
+    "Energies for debugging: \n",
+    "array([[-0.07434323, -0.00648831,  0.00713485, -0.07458268],\n",
+    "       [-0.07592777, -0.01072111,  0.01193199, -0.07601669],\n",
+    "       [-0.07959289,  0.00016017,  0.00015903, -0.07977653]])\n",
+    "J_ii for debugging: (check if this is the same as in calculate_exchange_tensor)\n",
+    "array([-0.07601669, -0.07959289, -0.07434323])\n",
+    "Test J_xx = E(y,z) = E(z,y)\n",
+    "-0.07601669325607426 -0.07977653333847946\n",
+    "\n",
+    "================================================================================================================================================================\n",
+    "Runtime information: \n",
+    "Total runtime: 808.34989875 s\n",
+    "----------------------------------------------------------------------------------------------------------------------------------------------------------------\n",
+    "Initial setup: 0.1594020829999998 s\n",
+    "Hamiltonian conversion and XC field extraction: 0.995 s\n",
+    "Pair and site datastructure creatrions: 0.021 s\n",
+    "k set cration and distribution: 0.021 s\n",
+    "Rotating XC potential: 0.232 s\n",
+    "Greens function inversion: 806.884 s\n",
+    "Calculate energies and magnetic components: 0.037 s"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "Magnetic entities integrated.\n",
+    "Pairs integrated.\n",
+    "Magnetic parameters calculated.\n",
+    "##################################################################### GROGU OUTPUT #############################################################################\n",
+    "================================================================================================================================================================\n",
+    "Input file: \n",
+    "Not yet specified.\n",
+    "Output file: \n",
+    "test.pickle\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",
+    "Parameters for the contour integral:\n",
+    "Number of k points:  2000\n",
+    "k point directions:  x\n",
+    "Ebot:  -13\n",
+    "Eset:  50\n",
+    "Esetp:  10000\n",
+    "================================================================================================================================================================\n",
+    "Atomic information: \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",
+    "================================================================================================================================================================\n",
+    "Exchange [meV]\n",
+    "----------------------------------------------------------------------------------------------------------------------------------------------------------------\n",
+    "Magnetic entity1        Magnetic entity2       [i  j  k]       d [Ang]\n",
+    "----------------------------------------------------------------------------------------------------------------------------------------------------------------\n",
+    "[3]Fe(2)        [4]Fe(2)        [0 0 0]        d [Ang] Not yet.\n",
+    "Isotropic:  -29.313466612791018\n",
+    "DMI:  [-8.85106686e-01  4.16843252e+00 -3.50602459e-06]\n",
+    "Symmetric-anisotropy:  [-4.7530289   1.27630805  0.34469396  4.15760367 -0.03005581]\n",
+    "Energies for debugging: \n",
+    "array([[-0.02803716, -0.00085505,  0.00091516, -0.02583675],\n",
+    "       [-0.03702114, -0.0041576 ,  0.00417926, -0.0340665 ],\n",
+    "       [-0.03178903, -0.0003447 , -0.00034469, -0.03139132]])\n",
+    "J_ii for debugging: (check if this is the same as in calculate_exchange_tensor)\n",
+    "array([-0.0340665 , -0.03178903, -0.02803716])\n",
+    "Test J_xx = E(y,z) = E(z,y)\n",
+    "-0.03406649551707385 -0.031391324924774755\n",
+    "\n",
+    "[4]Fe(2)        [5]Fe(2)        [0 0 0]        d [Ang] Not yet.\n",
+    "Isotropic:  -74.97640334609967\n",
+    "DMI:  [-6.84069197e+00  1.12809486e+01  8.75098125e-04]\n",
+    "Symmetric-anisotropy:  [-1.03102826  0.63307456 -0.16123618 10.67285449 -0.32447457]\n",
+    "Energies for debugging: \n",
+    "array([[-0.07434333, -0.00651622,  0.00716517, -0.07457845],\n",
+    "       [-0.07591065, -0.01067285,  0.01188904, -0.07600743],\n",
+    "       [-0.07959609,  0.00016211,  0.00016036, -0.07978159]])\n",
+    "J_ii for debugging: (check if this is the same as in calculate_exchange_tensor)\n",
+    "array([-0.07600743, -0.07959609, -0.07434333])\n",
+    "Test J_xx = E(y,z) = E(z,y)\n",
+    "-0.07600743160328334 -0.07978158717380361\n",
+    "\n",
+    "================================================================================================================================================================\n",
+    "Runtime information: \n",
+    "Total runtime: 534.5927785000001 s\n",
+    "----------------------------------------------------------------------------------------------------------------------------------------------------------------\n",
+    "Initial setup: 0.14542879100008577 s\n",
+    "Hamiltonian conversion and XC field extraction: 0.903 s\n",
+    "Pair and site datastructure creatrions: 0.016 s\n",
+    "k set cration and distribution: 0.014 s\n",
+    "Rotating XC potential: 0.242 s\n",
+    "Greens function inversion: 533.223 s\n",
+    "Calculate energies and magnetic components: 0.050 s"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [
+    {
+     "ename": "SyntaxError",
+     "evalue": "invalid syntax (577185453.py, line 1)",
+     "output_type": "error",
+     "traceback": [
+      "\u001b[0;36m  Cell \u001b[0;32mIn[10], line 1\u001b[0;36m\u001b[0m\n\u001b[0;31m    Magnetic entities integrated.\u001b[0m\n\u001b[0m             ^\u001b[0m\n\u001b[0;31mSyntaxError\u001b[0m\u001b[0;31m:\u001b[0m invalid syntax\n"
+     ]
+    }
+   ],
+   "source": [
+    "Magnetic entities integrated.\n",
+    "Pairs integrated.\n",
+    "Magnetic parameters calculated.\n",
+    "##################################################################### GROGU OUTPUT #############################################################################\n",
+    "================================================================================================================================================================\n",
+    "Input file: \n",
+    "Not yet specified.\n",
+    "Output file: \n",
+    "test.pickle\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",
+    "Parameters for the contour integral:\n",
+    "Number of k points:  1000\n",
+    "k point directions:  x\n",
+    "Ebot:  -13\n",
+    "Eset:  50\n",
+    "Esetp:  10000\n",
+    "================================================================================================================================================================\n",
+    "Atomic information: \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",
+    "================================================================================================================================================================\n",
+    "Exchange [meV]\n",
+    "----------------------------------------------------------------------------------------------------------------------------------------------------------------\n",
+    "Magnetic entity1        Magnetic entity2       [i  j  k]       d [Ang]\n",
+    "----------------------------------------------------------------------------------------------------------------------------------------------------------------\n",
+    "[3]Fe(2)        [4]Fe(2)        [0 0 0]        d [Ang] Not yet.\n",
+    "Isotropic:  -29.359042196157713\n",
+    "DMI:  [-8.93813379e-01  4.16574774e+00 -3.59758916e-06]\n",
+    "Symmetric-anisotropy:  [-4.57147523  1.18770348  0.33237646  4.15572505 -0.02850636]\n",
+    "Energies for debugging: \n",
+    "array([[-0.02817134, -0.00086531,  0.00092232, -0.02597527],\n",
+    "       [-0.03689861, -0.00415573,  0.00417577, -0.03393052],\n",
+    "       [-0.03160252, -0.00033238, -0.00033237, -0.03121927]])\n",
+    "J_ii for debugging: (check if this is the same as in calculate_exchange_tensor)\n",
+    "array([-0.03393052, -0.03160252, -0.02817134])\n",
+    "Test J_xx = E(y,z) = E(z,y)\n",
+    "-0.03393051742760642 -0.031219274782673133\n",
+    "\n",
+    "[4]Fe(2)        [5]Fe(2)        [0 0 0]        d [Ang] Not yet.\n",
+    "Isotropic:  -74.98947011735306\n",
+    "DMI:  [-6.82126017e+00  1.12890803e+01  6.51055432e-04]\n",
+    "Symmetric-anisotropy:  [-1.03170277  0.64211178 -0.16032097 10.67947464 -0.32714492]\n",
+    "Energies for debugging: \n",
+    "array([[-0.07434736, -0.00649412,  0.00714841, -0.07459988],\n",
+    "       [-0.07593108, -0.01067947,  0.01189869, -0.07602117],\n",
+    "       [-0.07956766,  0.00016097,  0.00015967, -0.07975216]])\n",
+    "J_ii for debugging: (check if this is the same as in calculate_exchange_tensor)\n",
+    "array([-0.07602117, -0.07956766, -0.07434736])\n",
+    "Test J_xx = E(y,z) = E(z,y)\n",
+    "-0.07602117288317486 -0.07975215648287909\n",
+    "\n",
+    "================================================================================================================================================================\n",
+    "Runtime information: \n",
+    "Total runtime: 272.760337666 s\n",
+    "----------------------------------------------------------------------------------------------------------------------------------------------------------------\n",
+    "Initial setup: 0.13726933300000033 s\n",
+    "Hamiltonian conversion and XC field extraction: 0.970 s\n",
+    "Pair and site datastructure creatrions: 0.013 s\n",
+    "k set cration and distribution: 0.017 s\n",
+    "Rotating XC potential: 0.217 s\n",
+    "Greens function inversion: 271.373 s\n",
+    "Calculate energies and magnetic components: 0.034 s"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
    "metadata": {},
    "outputs": [
     {
diff --git a/test.pickle b/test.pickle
index bbb3a28..292d728 100644
Binary files a/test.pickle and b/test.pickle differ