{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:26.831051Z",
     "start_time": "2024-03-30T23:39:23.276549Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "LightSimBackend import error: No module named 'grid2op'\n",
      "PhysicalLawChecker import error: No module named 'grid2op'\n",
      "TimeSerie import error: cannot import name 'TimeSerie' from 'lightsim2grid.timeSerie' (C:\\ProgramData\\miniconda3\\envs\\py3.11\\Lib\\site-packages\\lightsim2grid\\timeSerie.py)\n",
      "ContingencyAnalysis import error: cannot import name 'ContingencyAnalysis' from 'lightsim2grid.contingencyAnalysis' (C:\\ProgramData\\miniconda3\\envs\\py3.11\\Lib\\site-packages\\lightsim2grid\\contingencyAnalysis.py)\n",
      "rewards import error: No module named 'grid2op'\n"
     ]
    }
   ],
   "source": [
    "import pandas as pd\n",
    "pd.options.display.float_format = '{:,.5f}'.format\n",
    "import igraph\n",
    "import pandapower as pp\n",
    "import pandapower.plotting\n",
    "from pandapower.test.shortcircuit.test_1ph import *"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:26.837253Z",
     "start_time": "2024-03-30T23:39:26.834035Z"
    }
   },
   "outputs": [],
   "source": [
    "# beauftragt: Y0y0d5, YN0y0d5, Y0yn0d5, YN0yn0d5, Y0y0y0, Y0d5d5, YN0d5d5, Y0d5y0, Y0y0d11 und D0d0d0 \n",
    "# implementiert zusätzlich: YNdy, Ydyn, YNdyn, YNyy"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:26.896786Z",
     "start_time": "2024-03-30T23:39:26.838329Z"
    }
   },
   "outputs": [],
   "source": [
    "vg = \"yyy\"\n",
    "net = single_3w_trafo_grid(vg)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:27.057796Z",
     "start_time": "2024-03-30T23:39:26.898800Z"
    }
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "1 in collections is of unknown type. Skipping\n"
     ]
    },
    {
     "data": {
      "text/plain": "<Figure size 1000x800 with 1 Axes>",
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA+gAAAMDCAYAAAAxID+lAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAACWvUlEQVR4nOzdd3hUZd7G8XuSSSGQ0EJJCL0EqaGDIAiooEAorosVRVAU+6u79rq6iu6urg1FBXWtq6KEoqggIkgRlEAE6T0hSCiBhNSZ94/ZOZlJzsQAIWdm8v1c117OTIbJLzCr3PM8z31sTqfTKQAAAAAAYKkQqwcAAAAAAAAEdAAAAAAA/AIBHQAAAAAAP0BABwAAAADADxDQAQAAAADwAwR0AAAAAAD8AAEdAAAAAAA/QEAHAAAAAMAPENABAAAAAPADBHQAAAAAAPwAAR0AAAAAAD9AQAcAAAAAwA8Q0AEAAAAA8AMEdAAAAAAA/AABHQAAAAAAP0BABwAAAADADxDQAQAAAADwAwR0AAAAAAD8AAEdAAAAAAA/QEAHAAAAAMAPENABAAAAAPADdqsHAAAAAPzZwYOfaNeuR1RUdNzqUQKe3R6tFi3+poYN/2T1KIBfsjmdTqfVQwAAAAD+avXqc5Sb+5vVYwSNqKj26t17k9VjAH6JFXQAAACgHCUr567ToTZbiMLCGkgKUXHxcRUXZys0NEahodEev8qhwsLf5XQ6FBbWQDabXQ5HroqKjigkpKbs9joez3WqsPCQnM5ChYXFymYLl8ORp6KiLIWERMpur19qniw5HHmy2+srJCRSTmeBCgsPyWYLU1hYrCSbx3OPyuHIkd1eVyEhUXI6i1RY+LvXz+BW9mcp+zMYP12Zn6Xsz1DyXNfPUvb3E0BpBHQAAACggqKiEtW162JFRDQ2Htu160nt2vWwmjb9i1q0eEj5+QeUmjpEkk1JSd8pKqqd8dyMjJnavHmyYmOvUdu2r6i4OEfr1w9XcfEJde36g2JiehvPzcqar7S0cYqJ6aUOHT6WJG3cOF5ZWQvUufM81a8/wnhudvZqpaZeqMjIFurS5SuFhtbU1q23KD39dSUmvqW4uOuN5+bmbtG6dYNlt9c2fhb3z9Cixd/UosVDxnPdP0tR0TElJS1WVFQ742eIj79Jbdu+IpvNFfKLio5r/frhyslJU9eu8xQT09v4GWJjxyg7+ycVFOyv9D8TIJiwxR0AAAAox48/JqigYL9sNrv69t3rFc7d3AG3SZNbdeTIov8FWu9w7uYOuA0bXqm8vB3KyflVXbt+4xXO3dwBt27dCyVJR458o06dZnuFczd3SK9Zs6MiI1vp4MEPlJj4plc4d/MM6XXrDtX+/S+XCeduniG9ceMJ2rNnmuLjp3iFczfPkJ6QcKf27HlG9etfog4dPtbKla1UULBf4eFNdO65+3z/hgPVGAEdAAAAKIc7oIeHx+ncc9N9Pm/r1tu0f//LkqTevTebhnO3ffte1LZtd0iSundfqZiYPj6f+/vvX+jXX8dKkjp2/FwNGozx+dzs7FX6+ee+kqQ2bf6thITbfT43N3eLVq9OlCQ1aXKr2rZ9yedz8/MPaMWKOElSnTpD1LXrN2XCuVtR0XEtWxYjSQoPj1ffvjsVEhLu8ftIQAd84TJrAAAAQAU4HHk+v5aff0BHjiwy7h88+F+fzy0qOq6DBz827h848LacToeP71mgzMx3jPuZme/I4Sgwfa7T6dCBA297zPBxuee9PWc8cmSR8vMP+Hzu4cMLjNu5uZt08uQ2n889dmypcbuo6IhOnFjn87kAvBHQAQAAgAooKjqijIyZZR733ALeu/dmtWjxN+3a9bB27XrS5DVKtoB3775KiYlvKT39dW3dekuZkO5wFHidOe/ceZ6yshZo48bxZUK60+nwOnPevfsq5eSkaf364aYh3fPMee/em1VUdEypqUNMQ7rnmfN+/dJlt9fRunWDlZu7pcxzPc+c9+9/SLVqdVNq6oXKzl79h7+/AAjoAAAAQIWEhNTU5s2TvUK6d4ma68x5ixYPmYZ07xI115nzuLjrlZj4ZpmQ7hnO3WfO69cfoU6dZpcJ6d7h3HXmPCamt7p2/cY0pJcuhIuKaqekpO9MQ3pJOHedOY+IiFPXrotlt9cuE9Ld4dx95jwsrL66dPlKNWt2UmrqhXI6zVf+AZQgoAMAAAAVYLfXUXz8FCOkm4Vzt9Ih3Sycu5UO6Q5Hfplw7lY2pOeXCeduZiHdV1u7WUgvHc7dZ84jIhqXCemlw3lISPj/fs+ijZBeWHjorP3ZAMGCkjgAAACgHJ7lZv367flfIH5NkqsEzVdbu1SyWi1JoaExPtvaJfdq9SRJks0Wpk6dPjdta5fcq9Vj5XQWSlKZS6l5cre7FxdnS5LPtnappN29oMBVhlf6Umqe3B9Q5OZukiTFxo7xCueeioqOa/nyWDmdBZTEAeVgBR0AAACoIJstRM2bP2Lcb9x4Qrlt7QkJdxi369Q532c4l6RGja42boeHN1a9ehf7fG69ehcrPLzkcm+ev7a0mJjeqlPnfNOZSouKaqfGjScY95s3f8RnW3tERGM1bXqPcb9lyydNw7nkWkm322N8fl8ALgR0AAAAoIJcq8ZDFR4epzp1hmjPnmmmxXFSyZnz0NAY1a+frKysFNPiOKnkzLnNFq7Y2EuVn7/PtDhOKjlznp+/T7Gxl8pmCzctjnPbtetJZWWlqH79ZIWGxvgsjpNcq/h79kxTnTpDFB4ep9TUoT7b3bOy5mvLlpsVE9NPkZGtlZp6kWlxnORaxS8szDL9GoASBHQAAACgQhweZ86XqGvXb7zOpHsqfea8c+c5PtvdSxfCder0qWlxnFS2EK5Tp09Ni+PcPM+cd+48x2dxnORdCNe16zdKSlris93d88x5UtISdeu2zLQ4TirZYm+zhZ327zxQXRDQAQAAgAooLPzdqxDOZgtR27avlAnpvgrhzNrdzdrapbLFcU6nw7StXTIrjnOFdLNCOF/t7maFcL7a3c0K4cyK46SScF6zZieFhcVWxR8TENAoiQMAAADK4S6Jk0LUu/emMmfOPYNzmzYv6ODBj03b2t3cwbl580eUk7PetK3dzR2c4+Ju/N/9GWXa2t08g3PNml20e/cTPgvhPINzw4bjtW3bnWXa2t3cxXF2e201bXqPtmy5uUxbu5tns32rVtO0destqlmzk7p0+UqrV59jlO1REgeYI6ADAAAA5XAH9LCwRurf3/w8ttPp0KZNE3Tw4PuSpO7dV5VbCLdz56PavfsJSVLnzvN8trVLUkbGW9q8ebIk/S+cT/L53Kys+dqwYaQkV8Fby5aP+3xudvZq/fxzH0lSw4ZX6Zxz3vVZCJebu0WrVydKkmJi+ikpaYnPQrj8/ANasSJOkhQW1kB9+myX3R7t1YZPQAfMscUdAAAAqACns8jn14qLc5SXt8O4n5OT5vO5DkeBcnLWG/ePH/+lnO/p0PHjazyeu8a0OM7stXJy1vssjis9Y17eDhUX5/h87smTW43bBQUHVVh42Odz8/P3GLcdjjwVFGT4fC4AbwR0AAAAoAKKirKUlTXf5HH3mfNf1b37SsXH32RaHCd5nznv3Hmez+I4qXQh3FtKTHzLtDjOzbsQbp7P4jjJ88z5TerefaVycn712e7u3jofGztGvXqlyeE4aVocJ5VsnY+JOVe9e29RRESCaXEcAHN2qwcAAAAAAkFISKTS0sZ5nRc3K4SLju4lSca2dPd5cbNCOPfr7Nr1sCQZ58V9FcJ5vq7neXGzQrhOnWYrLW2cNm4c73Ve3KwQrmvXb5SaeqHWrx+uLl2+kt0eLcm8EC4p6TutWzdYqalD1LXrYkVEuK7H7nmu3f0aXbsuVmrqEK1bN1hOZ/HZ+GMBggor6AAAAEAF2O31Vb/+JUpLG6esrPk+29rN2t19tbVLZdvdywvnZu3uZuFcMm93Nwvnknm7u1k4l2Ta7m4WziV5tbsXFv5+1v+MgEBHSRwAAABQDs9ys759d2jjxvE6dOgLSVJoaIzPtvaSoP2aJMlmC/fZ1i6VrIK72Hy2tUslq+CS66/yvtrapZJVcKfTtdU9Pv4m07Z2qWQVvLg4W5IUGzvGtK1dKml3LyhIlyTFxJzrFc495ecf0MqVTeV0FlESB5SDFXQAAACggkJCwpWY+KZxPyHhTp9t7TZbiNq0ecG4X7/+qHLb2ps3f8C4HRnZ0mc4l1wr6ZGRLU1/bWmurfSjjPtt2rzgs609Jqa3EhLuNO4nJr7ps609KqqdWrWaZtxv3/5t03AuuVbS7fa6PmcE4EJABwAAACqoqOi4NmxIVkhIDYWHx2vPnmdMi+Mk95nzy2WzhSkioqkOHZptWhwnlay2SzZFRrZSXt4O0+I4t127nlRe3g5FRraSZPNZHCe5VtsPHZqtiIimstnCtHHj5T7b3bOy5mvPnmcUHh6vkJAa2rAh2bQ4TnKttm/deovCwhooNDRaaWmjTYvjJNdqe2Fhls+fB4ALAR0AAACoEKdx5jwpaYn69t3pdSbdk/eZ88/Vt+8urzPpXq9a6sx5377by2139zxz3rfv9jJn0j15njnv23eXOnX63Ge7u+eZ8759dyopaYnXmXRPnmfO+/TZrh491nidSffk3grva9UeQAn+XwIAAABUQGHhIa9CuJCQcHXo8HGZkG5WCGdWHCf5bmsvXRznZlYIZ1YcJ5m3tZsVx0nmbe1mxXGSeVu7WXGcVBLO7fbaCgtrUAV/SkBgoyQOAAAAKIe7JE6yqXv3lWXOnHsG8g4dPlZm5jumbe1S6UD+ho4fX2Pa1u7mGchd98u2tbt5BvLo6J7avPmGMm3tbp6BvFGja7Vx4/gybe1unoG8ZcsnlJY2rkxbu5tnIG/X7nVt3Hi57Pba6tp1sdau7WmU7VESB5gjoAMAAADlcAf0sLAG6t//oOlzHI4CpaWN0+HDrlX0zp3n+SyEczod2rJlqjIyXpckJSa+VW4hnGe7e3lt7ZI7pE+SJMXFTVG7dq/63FqelTVfGzaMlCTVq+daWfdVCJedvVo//9xHklSzZhd167bMZyFcbu4WrV6dKEmy2+upV69fFRHR2KsNn4AOmGOLOwAAAAAAfoCADgAAAFRAYeEhZWevLvO4e4v7kSPfqGPHzxUbO8a0OE4q2eKekTFDiYlvKj7+JtPiODfvM+e+i+Mkzy3uNykx8U1lZMzw2e7u3uIeGztGHTt+riNHvjEtjpNKtrjHxJyrrl2/VV7eLtPiOKlki3tU1DlKSlqqkJBI0+I4AObsVg8AAAAABAKbLUypqRcaJXGSeSFc/fqXaOPG8UpLG+d1Dt2sEK5x44mSpM2bJ0uS11Z3s0I41+Pu7e4lj5kVwkk243U9z6GbFcJ16jRbaWnjtHHjeK9z6GaFcF27fqPU1Au1fv1wr3PonufPu3ZdrIiIxkpK+k7r1g1WauoQSeaXgQNQghV0AAAAoALCwmJVs2YnpaZeqOzs1abhXJJpu7uvtnZf7e6+wrlZu7t5ODdvdzcL55JM293Nwrkk03Z3s3AuyavdvbDw9yr4UwICGyVxAAAAQDk8y816996k9euH68SJX2S311Vh4SHTtnbJc3V9vsLDGys/f5/PtnbPAB8Z2VJ5eTvKLYRzB/jIyFbKy9vps61dKgnwEREJKig4oPr1R5i2tUslq+thYbEqKjqiWrW6mba1SyWr6yEhEXI48hQRkeAVzj25iuPOkeSgJA4oByvoAAAAQAXZ7dHq3DlFDsdJFRSkq1mz+3y2tbtW0j+S01mo/Py9io0d57Ot3b2SLjmVl7dDkZGtym1rb9Hiof+F8x2SnD7DueRaSY+NHaf8/L1yOgvVocNHPtva69cfoWbN7lNBQbocjpPq3DnFZ1t7TExvtW37igoLf1dx8XF16jTHNJxLrpX0sLD6Pn8eAC4EdAAAAKCCHI4C41y3JO3b94JpcZzkWhXftu1O435W1lzT4ji33bv/btzOy9vpszhOcq2K5+XtNP21pWVlzVdW1lzj/rZtd5oWx0muVfF9+14w7m/ePNm0OE5yrYrv2HGvcf+3364zLY6TpPz8AyoqOuJzRgAuBHQAAACggtxnzjt3nqcBA7K9zqR78j5z/pYGDsz3OpNemueZ80GDisucSffkeeZ80KDictvdPc+cDxyYr8TEt7zOpHvyPHM+YEC2Onee53Um3ZPnmfN+/TLUvfsqrzPpnvLzDyg1dYjPDwUAlCCgAwAAABVQVJTlVQhnt0erS5evyoR0s0I4s+I4t9KFcL6K4yTzQjiz4jjJvK3drDhOMm9rNyuOk8zb2s2K46SScF5UdExhYQ3O+p8REOgoiQMAAADK4S6Jk6TOneeVOXNeVHRc69cPV05Omrp2/VoHDrxdpq3drXTz+/Hjv5i2tUtlg74k07Z2N8+gHx3dzbSt3c0z6DdufJ1SUy8q09bu5hn0W7Z8UqmpF5Vpa3fzDPrt27+ttLTRKio69r/LrQ0xyvYoiQPMEdABAACAcrgDut1eXwMGHDJ9jiukD1N29gpJUmLiWz4L4dwh/dChLySp3Lb2kpD+miQpPv6mcgvh3CFdkmJjx/hsa5fcIX2SJCkmpp+6dFnosxAuK2u+NmwYKUmKjGytbt2W+SyEy85erZ9/7iNJCg2NVo8eaxQV1c6rDZ+ADphjizsAAABQATab3efXQkNrKjKylXG/Zs1OPp8bEhKumjW7GPejo7uV8z1DFB3d0+O5PX2G89KvVbNmF5/hvPSMkZGtFBpa0+dza9Roa9wOD2+osLB6Pp8bEdHMuB0SEqnw8DifzwXgjYAOAAAAVEBh4e/Kzd1S5nH3KvfBgx+oTZt/KybmXNPiOLddu57U7t1PqHnzRxQbO8ZncZzkXuW+QXFxUxQXN0WbN9/gs93dvRU9NnaMmjd/RLt3P2FaHCeVbEWPiTlXbdr8WwcPfmBaHCeVnDmPijpHiYlv6fjxtabFcVLJmfPw8Hi1b/8fORz5psVxAMz5/hgQAAAAgMFmC9G6dYOVlPSdoqLaSTIvhGvceKLWrx+u1NQL1bXrN4qJ6W28RulCOPd297S0cUb5nFvpQjjXDDbjMm+eW+jNCuFstjBju7vnFnqzQrjQ0FrG63puoTcrhAsPb6S0tHHauHG81xZ6z0I49+9RVFQ7paZeqPXrh0viZC3wR1hBBwAAACogLKyB7PbaWrdusHJzt5iGc0k+291Lh3NJPtvdzdrafbW7m4VzSabt7mbhXJJpu7tZOJdk2u5uFs4lebW7Fxaan98HUIKSOAAAAKAcnuVmPXqs+V8QPaqoqHN09Oh3pm3tkne7e5065ysrK8VnIZxnu3v9+qN06NBsn23tnh8MxMaOU1bWXJ9t7VLJBwP16yfr6NElPtvapZIPBurUGazc3E2y2+uYtrVLJR8MREf3UEHBQTkcJ73CuSdXcVxfSU5K4oByENABAACAcpRuH8/Pz9CKFfGSpGbN7lOrVk/7/LVFRce1bFmMJKl+/WR17jzH53MdjgItXRohSYqIaKq+fXf5LIRzOh1aubKF8vP3SpIGDswvtxBuw4bRyspKkSQNGJDts61dknbsuF979jwjSerXL10REb5L3jyb4Hv1SlPNmh19Pnf58gYqLDxEQAfKwRZ3AAAAoIKcTod2737CuH/gwLumxXFu+/b927h99OgSn8VxkpSZ+Z5xu6DggA4f/tLncw8f/lIFBQdMf21p2dmrdfToEtOZSsvN3aIDB9417u/e/YRpcZzkOnO+d+8/jPs7dz5kWhwnuT6oKCrK9vl9AbgQ0AEAAIAKKjlz/pb69cvwOpNemueZ8wEDssucSfdUcub8Jg0cmKf69Uf4bHcvOXM+QgMH5ik+/iavM+mePM+cDxiQXeZMuifPM+f9+mUoMfEtrzPpnjzPnPfuvVmdO8/zOpPuyb3V3+ks/MPfX6C6I6ADAAAAFVBUdNSrEC4iorG6dl1sGtJLF8L5Ko6TyhbChYREmBbHSWaFcBGmxXGSeSGcWXGcZN7WblYcJ5m3tZsVx7l+z0rO4YeFxZ61PxsgWHAGHQAAACiH+wy6JCUmvlWmEK50YD148L9l2trdPAOru928dFu7m2dxXKdOsyXJtK1dKnu5N/cHAb4K4Tw/QGjY8M+mbe1unh8gNG/+iFJTh5Zpa3fz/AAhMfFNbdiQbPysaWnjvM7yAyiLgA4AAACUwx3Q7fa6GjDgsOlz3CE9N3eTJPlsa5dKQnp29o+SpPj4m0zb2qWSkH7o0BeSpNjYMT7b2ktC+muSpJiYc322tUslIV2SoqLO8dnWLnmXwYWHxykpaYlpW7vkCukbNoyUJIWE1FBS0hLFxPQuU7YHoCy2uAMAAAAVEBIS6fNrERGNVbfuUON+w4Z/9vlcuz1aDRuON+43bnydz7b2kJBwNWp0rXG/UaNrfba122whatz4Oo8Zxpfb1u45Y926Q32Gc0mqV+8S43ZU1DmqUaONz+fWrj3QuG2311WtWkk+nwvAGwEdAAAAqIDCwt+Vn3/A9Gu7dj2p/ftfVpMmtyoq6hyfxXGSazV627Y71bDhVYqJ6afU1It8trtnZc3Xxo3jVa/eCNWrN+J/W97LFsdJ7jPnFykmpp8aNrxK27bdaVocJ5WcOY+KOkdNmtyq/ftfNi2Ok0p2B4SHx6tZs/t09Oh3psVxUsnugNDQGDVv/ogKCw+ZFscBMGe3egAAAAAgEDidDqWmDimzFbx0IZw70K5bN7jMOe3ShXDFxTlav364UlMvVNeu3ygmprfx3NKFcJK0ceN4paWNU6dOs1W//gjjuaUL4UJDa8puj9bmzZMlyevcvFkhXFhYI2O7u+fWfLNCuBo12hqv67k1v/T5+piY3oqJ6a20tHHauLFkxwAA31hBBwAAACogLKyBioqOKTV1iLGSXjqcS/LZ7l46nNtsIT7b3cu2tYcrJCTctN3drK3dZgsxbXc3C+eSTNvdzcK5JNN2d7NwLsmr3b2oKKsq/piAgEZJHAAAAFAOz3KzpKTFRsCtW3eo9u9/2WchnGfAbdx4gvbsmWba1i55rz4nJNypPXueMW1rl7zb3Zs1u0/79r3gs63ds929WbN7deDAuz7b2qWSDxyaNLlVR44s8tnWLpV84NCw4ZXKy9uhnJxfy+wCcPMsjqMkDvCNgA4AAACUo3T7eG7uFq1enShJatLkVrVt+5LPX5uff0ArVsRJkurUGaKuXb/xWQhXVHRcy5bFSJLCw+PVt+9On4VwDkeBVq5sqYKCdEnSgAHZPgvhXFvzL9TRo4slSf36ZZRbCLd1623av/9lSVLv3pt9trVL0r59L2rbtjskSd27r1RMTB+fz122LFZFRVkEdKAcbHEHAAAATsHBg/81bh85sshncZwkHT68wLidm7tJJ09u8/ncY8eWGreLio7oxIl1Pp974sQ6FRUdMf21pZ08uc24/FvpmUrLzz+gI0cWGfc9f9bSioqO6+DBj437Bw68bVocJ7k+UHA4cny+FgAXAjoAAABQQZ5nznv33lzmTLqnkjPnN6lfv3TZ7XV8tru7z5zHxo5R//6HVKtWN68z6Z7cZ85r1eqm/v0PKTZ2jNeZdE8lZ87rqF+/dMXH3+R1Jt2T55b83r03lzmT7slzS3737quUmPiW15l0T+4t+Q5HXrm/twAI6AAAAECFFBcf9yqEi4pqp6Sk70xDeulCuIiIONPiOKlsIVxYWH3T4jipbCFcWFh90+I4yawQLs60OE4yL4QzK46TzNvazYrjJO/z8nZ7/bPy5wIEE86gAwAAAOVwn0GXZFoIVzoIHz68oExbu1vpIHzy5NYybe1upYOwpDJt7W6eQbhTp9mqUaOtaVu75F0cl5j4purVu8S0rd3Nc9dAQsIdpm3tbp4fTLRp84I2brzcmGnz5ileZ/kBlEVABwAAAMrhDuihoTE677xjps9xh3R3aVt8/E2mbe1SSUh3nwuPjR1j2tYulYT07OwfJUkxMeeatrVLJSH90KEvJElRUef4bGsvCemvSXKV0vlqa5dKQrokhYbG+Gxrl9whfZIkyWYLU6dOn6t+/RFlyvYAlMUWdwAAAKACQkNr+fxaVFQ7NW48wbjfvPkjPtvaIyIaq2nTe4z7LVs+6bOt3W6PVsuWT3g89wmfbe0hIeFq2bJkK3rTpvf4bGu32ULUvPkjxv3GjSeU29aekHCHcbtOnfN9hnNJatToauN2eHhj1at3sc/nAvBGQAcAAAAqoLDwkIqKjpt+LSNjpvbsmaY6dYYoPDxOqalDfba7Z2XN15YtNysmpp8iI1srNfUi0+I4yXXmPC1tnGrW7KKaNbsoLW2caXGc5FrFT029SJGRrRUT009bttxsWhwnuVfxhyo8PE516gzRnj3TTIvjpJJV/NDQGNWvn6ysrBTT4jipZBXfZgtXbOylys/fZ1ocB8AcAR0AAACoAKezUOvXDy8T0j3PXXft+o2Skpb4bHf3LIRLSlqibt2WmRbHSd6FcN26LVO3bstMi+Mk73Pw3botU1LSEtPiOKn0Ofgl6tr1G9PiOKnsOfjOnef4bHcvfQ6+U6dPvYrjAPwxAjoAAABQAWFhscrJSfMK6aXb2m22EJ/t7qXb2kNCwhUR0di03b10W7vdHi27Pdq03b1sW3tjhYSEm7a7m7W122whpu3uZm3tkkzb3UuH8/r1R0iSV7t7UdHRqvmDAgIYJXEAAABAOTzLzTp1mm0E54YNx2vbtjtN29ol7+DctOk92rLlZtO2dsk7OLdqNU1bt95i2tYueQfntm1f0Y4d95q2tUvewbldu+nau/cfPtvaPdvd27R5QQcPfuyzrV0qKY5r3vwR5eSsLxPOPXkWx1ESB/hGQAcAAADKUbp9PDt7tX7+uY8kqWHDq3TOOe/6LITLzd2i1asTJUkxMf2UlLTEZyFcfv4BrVgRJ0kKC2ugPn22+yyEKyo6rlWrWquw8HdJUr9+GT4L4RyOAq1bd76ys1dIknr33uyzEM7pdGjTpgk6ePB9SVL37qvKLYTbufNR7d7tKrHr3HmeaTh3W7asnoqKjhDQgXKwxR0AAAA4BTk5acbtvLwdKi7O8fnckye3GrcLCg6qsPCwz+fm5+8xbjsceSooyPD53IKCDDkceaa/trTCwsMqKDhoOlNpxcU5ysvbYdz3/FlLczgKlJOz3rh//PgvPp/rdDrkcBT4/DoAFwI6AAAAUEElZ85vUvfuK5WT86tpcZxUcuY8NnaMevVKk8Nx0rQ4Tio5cx4Tc656996iiIgE0+I4qWTrfEREgnr33qKYmHNNi+Okkq3zDsdJ9eqVptjYMabFcZLn1vlf1b37SsXH32RaHCd5b53v3Hmez+I4qWTrvMPh+4MMAC4EdAAAAKACHI5cr0K4mJg+6tr1mzLFcVLZQriaNTuaFsdJZQvhoqLamhbHSWUL4aKi2poWx0llC+Fq1uxoWhwnmRXC9TEtjnP9PpQthDMrjpO8z7Xb7XUr/c8ECDYEdAAAAKACioqOlCmEi4npXSakm7W1SzJtdzdra5dk2u5u1tYuybTd3aytXZJpu7uvtnazdndfbe1S2XZ3z3CemPimQkKiquzPCghUlMQBAAAA5XCXxIWE1NR552WbFsK5g3ZxcbYkKTZ2jGlbu1SyCl5QkC5Jiok517StXSpZBc/N3SRJioo6x7StXSpZBc/O/lGSFB4eb9rWLpWsgh869IUkKTQ0xmdbe0nQfk2SZLOF+2xrl0ra3V1sSkx8U3Fx15cp2wNQFivoAAAAQAXY7XV8trXHxPRWQsKdxn3XirF5W3tUVDu1ajXNuN++/ds+29ojIhqrXbvXjfvt2r3us63dbo9W+/ZvG/dbtZrms609JCRciYlvGvcTEu702dZus4WoTZsXjPv1648qt629efMHjNuRkS0VF3e9z+cC8EZABwAAACqgqCjLZxN5VtZ87dnzjMLD4xUSUkMbNiSbFsdJrtX2rVtvUVhYA4WGRistbbRpcZzkWm3fuPFy2e31ZLfX08aNl5sWx0mu1fa0tNEKDY1WWFgDbd16i2lxnOtnOa4NG5IVElJD4eHx2rPnGdPiOMm92n65bLYwRUQ01aFDs02L46SS1XbJpsjIVsrL22FaHAfAHAEdAAAAqACHI08bN44vE9I9z5z37btTSUlLTIvjJO9CuD59tqtHjzWmxXGSdyFcr16/qlevX02L4yTvQrgePdaoT5/tpsVxknchXFLSEvXtu9O0OM71M3ueOf9cffvuMi2Ok1TmzHnfvtvLbXcHUBYBHQAAAKgAu72+srIWeIV0s0I4s+I4qWxbu90ebVocJ5Vta4+IaGxaHCeVbWuPimpnWhwnmbW19zYtjpPM29rNiuOksuHcva3dsziuuNh8RwGAEpTEAQAAAOXwLDdLTHzdCOSNGl2rjRvHl2lrd/MM5C1bPqG0tHFl2trdPAN5u3av/29be23TQjjPQN6hw0fasmVKmbZ2N89A3qnTbO3c+UiZtnY3z0DeocPHysx8x7StXSodyN/Q8eNryoRzT57FcZTEAb4R0AEAAIBylG4fz8qarw0bRkqS6tUboU6dZvsshMvOXq2ff+4jSapZs4u6dVvmsxAuN3eLVq9OlCTZ7fXUq9evPgvh8vMP6KefOqqo6LAkqXfvzT4L4YqKjuuXXwYoJ2e9JKl791U+C+EcjgKlpY3T4cOuVfTOnef5LIRzOh3asmWqMjJcJXaJiW+VWwj3ww+1VVycTUAHysEWdwAAAAAA/AABHQAAAKgg95nz2Ngx6tjxcx058o1pcZxUssU9JuZcde36rfLydpkWx0klW9yjos5RUtJShYREmhbHSSVb3ENCIpWUtFRRUeeYFsdJJVvc8/J2qWvXbxUTc65pcZxUssX9yJFv1LHj54qNHWNaHCeVbHHPyJihxMQ3FR9/k2lxnNuuXU8a14gH4BsBHQAAAKgAhyPPqxCuQYMx6tRpdpniOKlsIVzdukNNi+OksoVwdeqcZ1ocJ5UthKtT5zzT4jipbCFc3bpDTYvjXD+bdyFcgwZjTIvjJLNCuEmmxXFu7vPnoaExlfrnAQQjAjoAAABQAUVFWWUK4erXH1EmpJu1tUsybXc3a2uXZNrubtbWLsm03d2srV2Sabu7WVu7JNN2d19t7b7a3d3hvEWLvyk01PzsPYASlMQBAAAA5XCXxIWERGrAgGOmhXDure9hYbEqKjqiWrW6mba1SyWr6yEhEXI48hQRkWDa1i6VrK47HHmS9L9t7WXb2qWS1fX8/H0KCYmUw5Fv2tYulayunzjxi+z2uiosPGTa1i55rq7PV3h4Y+Xn7/PZ1u4Z4CMjWyovb4datPibWrR4qEzZHoCyCOgAAABAOdzBUgpReHicz+cVF2cb1/p2Pc/3ZlWHI1dFRUckSWFhjWSz2X0+1+ksUGHh7/97bgPZbGU/IFi8+KQ++OC4rrwySgMGHJMk2e11FRISVc5P5lBBQYYkKTQ0+g+3oLt+D6SQkBqy2+tV6LlSqMLDG//vsQxJDgI6UA4COgAAAFCO1avPUW7ub1aPUa5rr5X27JGaNZPeecfqacoXFdVevXtvsnoMwC/5/qgOAAAAgFq0+Jt27XrYtH3dX5w86VqdPnmy/FV+q9nt0WrR4m9WjwH4LVbQAQAAgACXkJCg/fv3q0mTJtq3j+3jQKCixR0AAAAAAD9AQAcAAAAAwA8Q0AEAAAAA8AMEdAAAAAAA/AABHQAAAAAAP0BABwAAAADADxDQAQAAAADwAwR0AAAAAAD8AAEdAAAAAAA/QEAHAAAAAMAPENABAAAAAPADBHQAAAAAAPwAAR0AAAAAAD9AQAcAAAAAwA8Q0AEAAAAA8AMEdAAAAAAA/AABHQAAAAAAP0BABwAAAADADxDQAQAAAADwAwR0AAAAAAD8AAEdAAAAAAA/QEAHAAAAAMAPENABAAAAAPADBHQAAAAAAPwAAR0AAAAAAD9AQAcAAAAAwA8Q0AEAAAAA8AMEdAAAAAAA/AABHQAAAAAAP0BABwAAAADADxDQAQAAAADwAwR0AAAAAAD8AAEdAAAAAAA/QEAHAAAAAMAPENABAAAAAPADBHQAAAAAAPwAAR0AAAAAAD9AQAcAAAAAwA8Q0AEAAAAA8AMEdAAAAAAA/AABHQAAAAAAP0BABwAAAADADxDQAQAAAADwAwR0AAAAAAD8AAEdAAAAAAA/QEAHAAAAAMAPENABAAAAAPADBHQAAAAAAPwAAR0AAAAAAD9AQAcAAAAAwA8Q0AEAAAAA8AMEdAAAAAAA/AABHQAAAAAAP0BABwAAAADADxDQAQAAAADwAwR0AAAAAAD8AAEdAAAAAAA/QEAHAAAAAMAPENABAAAAAPADBHQAAAAAAPwAAR0AAAAAAD9AQAcAAAAAwA8Q0AEAAAAA8AMEdAAAAAAA/AABHQAAAAAAP0BABwAAAADADxDQAQAAAADwAwR0AAAAAAD8AAEdAAAAAAA/QEAHAAAAAMAPENABAAAAAPADBHQAAAAAAPwAAR0AAAAAAD9gt3oAIJAUFEj790sOhxQXJ0VFWT0RAACo7vLzpeJi122n09pZAJwZVtCBCti3T7rvPqlRI6lVK6lNG6l+femmm6SNG62eDgAAVEd79kj33CM1bCgdOOB6LD1duuUWafNma2cDcHpsTiefswHl+ekn6aKLpOPHSz6ddrPbJZtN+u9/pTFjLBkPAABUQ8uXSxdfLOXmuv9+kiBpv6Qmstv3KTRUmj1buuQSa+cEcGoI6EA50tOlTp2k7Oyy4dzNZpNCQ6UVK6SePat2PgAAUP3s3i117izl5LiO3bmUBHRpn2w2KSzMtdDQpYtlowI4RWxxB8rxyivlh3Op5KzXM89UzUwAAKB6e/FF18p5STgvy+l0ff3ZZ6tuLgBnjhV0wIeiIteZ88OHK/b8kBDXWfW4uLM7FwAAqL7y811nzrOzS3/FewXdzW53nU+vX7/qZgRw+lhBB3zIzKx4OJdcn1L/9tvZmwcAAGDfPrNw7ltRkbRly9mbB0DlIqADPpS3bcyX8rbCAwAAnCnffz+JLvXPivwaAP6GgA740KiRVKvWqfwKp/btWyxOjQAAgLOlSRMpMtLsK3+T1Pt//yxhs0mtW1fBYAAqBQEd8CE8XJo0ydXQ/seKJC3UxIlD1adPHy1YsICgDgAAKl1UlHTttVJISOll8T9JWvW/f7rY7VJystS4cVVOCOBMENCBctx2myuo22x/9MwQSdMkST/99JNGjBhBUAcAAGdF06az5XAUSip/77rDIf31r1UzE4DKQUAHytG6tfTFF1JEhOtT6NLsdld4nz5dmj37dnXt2tX4GkEdAABUtlmzZunhh/8k10p5oWy2sgU4drvr6jKzZknnnlvlIwI4AwR04A9cdJH000/S5Zd7h3SbTRo2TFqyRLrpphCNHTtWP//8s2bPnk1QBwAAlW7WrFmaNGnS//4uMU9XXvmixo8PKfP3k0sukX74QZowwbJRAZwmroMOnIKsLNel1BwOqVUrV1GLGYfDoTlz5ujxxx9Xamqq19d69eqlxx57TBdffLFsf7x3HgAAoFQ4l+644w49//zzstls+v1316XUHA6pTRspLs7iYQGcNgI6cBYR1AEAwJkqL5wDCC4EdKAKENQBAMDpIJwD1QsBHahCBHUAAFBRhHOg+iGgAxYgqAMAgPIQzoHqiYAOWIigDgAASiOcA9UXAR3wAwR1AAAgEc6B6o6ADvgRgjoAANUX4RwAAR3wQwR1AACqF8I5AImADvg1gjoAAMGPcA7AjYAOBACCOgAAwYlwDsATAR0IIAR1AACCB+EcQGkEdCAAEdQBAAhshHMAZgjoQAAjqAMAEHgI5wB8IaADQYCgDgBAYCCcAygPAR0IIgR1AAD8F+EcwB8hoANBiKAOAIB/IZwDqAgCOhDECOoAAFiPcA6gogjoQDVAUAcAwBqEcwCngoAOVCMEdQAAqg7hHMCpIqAD1RBBHQCAs4twDuB0ENCBaoygDgBA5SOcAzhdBHQABHUAACoJ4RzAmSCgAzAQ1AEAOH2EcwBnioAOoAyCOgAAp4ZwDqAyENAB+ERQBwDgjxHOAVQWAjqAP0RQBwDAHOEcQGUioAOoMII6AAAlCOcAKhsBHcApI6gDAKo7wjmAs4GADuC0EdQBANUR4RzA2UJAB3DGCOoAgOqCcA7gbCKgA6g0BHUAQDAjnAM42wjoACodQR0AEGwI5wCqAgEdwFlDUAcABAPCOYCqQkAHcNYR1AEAgYpwDqAqEdABVBmCOgAgkBDOAVQ1AjqAKkdQBwD4O8I5ACsQ0AFYhqAOAPBHhHMAViGgA7AcQR0A4C8I5wCsREAH4DcI6gAAKxHOAViNgA7A7xDUAQBVjXAOwB8Q0AH4LYI6AKAqEM4B+AsCOgC/R1AHAJwthHMA/oSADiBgENQBAJWJcA7A3xDQAQQcgjoA4EwRzgH4IwI6gIBFUAcAnA7COQB/RUAHEPAI6gCAiiKcA/BnBHQAQYOgDgAoD+EcgL8joAMIOgR1AEBphHMAgYCADiBoEdQBWMnhcCgzM1MZGRlKT0/3+qf7dmZmpvLy8lRUVKSioiIVFhbKZrPJbrfLbrcrLCxMNWrUUOPGjRUXF6f4+HjTfzZo0IB/l5WDcA4gUBDQAQQ9gjqAs83hcGjLli1as2aN1q5dq7Vr1+qXX37RiRMnquT7161bV927d1ePHj3Us2dP9ejRQy1btuTfayKcAwgsBHQA1QZBHUBlycrK0tdff63Vq1efdhgPCQlRgwYNVLNmTWPF3G63S5Kxml5UVKQTJ07o999/P+UZPUN73759dcEFFyg6OvqUXyeQEc4BBBoCOoBqh6AO4HRs3bpVKSkpSklJ0bJly+RwOMp9frNmzdS5c2clJCSYbktv2LChQkNDK/S9CwsLfW6X37t3r1JTU3XgwIFyXyM8PFxDhgxRcnKyRo0apYSEhAr/7IGIcA4gEBHQAVRbBHUA5SkuLtaKFSs0d+5cpaSk6LfffvP53GbNmhlby93/i42NrcJppfT0dGN7/dq1a7VmzZpyQ3u3bt2UnJys5ORkdevWLaj+XUc4BxCoCOgAqj2COgBP27Zt02uvvaZ3333X59byxMREjRo1SkOHDlXPnj2rPIxXVHp6utasWaOvv/5aKSkp2rt3r+nzmjZtquuvv1433HCDmjRpUsVTVi7COYBARkAHgP8hqAPVV3FxsebPn69XX31VCxcuLPP1kJAQ9e/f39genpiYaMGUZ8bpdCo1NdXYEbBmzZoyzwkNDdWYMWM0depUDR48OOD+fUc4BxDoCOgAUApBHag+MjMz9dZbb+n111/Xnj17vL4WHh6ukSNHavTo0brkkkv8dpX8dO3fv1/z5s3TnDlztHDhwjJn6tu3b6+pU6dqwoQJql27tkVTVhzhHEAwIKADgA8EdSB4/fbbb/rb3/6mTz75RIWFhV5fa9GihW6++WZNnDhRDRo0sGjCqrV371698cYbmjFjhjIzM72+FhUVpQkTJuiBBx5Q06ZNLZqwfIRzAMGCgA4Af4CgDgSPffv26bHHHtOsWbO8VoxtNpsuueQSTZ06VcOGDatwu3qwKSgo0Oeff65XX31VS5cu9fpaRESEbr31Vt1///2qX7++RROWRTgHEEwI6ABQQQR1IHAdPnxYzzzzjF566SXl5eUZj9evX1+TJ0/WlClT1LJlSwsn9D9paWmaPn263n33Xa9rvMfExOjee+/VHXfcoZo1a1o4IeEcQPAhoAPAKSKoA4EjJydHL774oqZNm6Zjx44Zj/tTyPR3WVlZmjZtWpkPNxo3bqxHHnlEkydPVlhYWJXPRTgHEIwI6ABwmgjqgP9yOp165513dP/993tdCzwiIkK33Xab7rvvPr/aph0I9u3bp8cff1wzZ870Oh7QunVrPf/88xo1alSVzUI4BxCsCOgAcIYI6oB/2bdvn2644QZ99dVXxmMhISGaOHGiHn30Ub8tOgsUv/32mx566CF99tlnXo9fffXV+ve//6169eqd1e9POAcQzAjoAFBJCOqAtZxOp2bNmqW77rpL2dnZxuNjx47VU089pXPOOcfC6YLP6tWrdd999+m7774zHmvcuLFmzJhx1lbTCecAgh0BHQAqGUEdqHpmq+bx8fGaMWOGRowYYeFkwc3pdOrdd9/VHXfc4XXG/2ysphPOAVQHBHQAOEsI6sDZ52vV/Nprr9Xzzz+vunXrWjhd9bF//35NmTJF8+fPNx6rzNV0wjmA6oKADgBnGUEdODuOHTuma665RnPnzjUeY9XcOr5W0ydPnqyXX35ZERERp/W6hHMA1QkBHQCqCEEdqDxbt25VcnKyfvvtN+MxVs39w/79+3XjjTdqwYIFxmP9+vXT7Nmz1bhx41N6LcI5gOqGgA4AVYygDpyZr7/+WuPHj9fRo0clSfXq1dM777yjkSNHWjsYDE6nU2+//bZuueUWnTx5UpKUkJCgL774Qj169KjQaxDOAVRHBHQAsAhBHTg1TqdT//73v3X33Xcb1+Hu2LGjUlJS1KpVK4ung5mff/5Zo0eP1r59+yRJkZGRmjVrli6//PJyfx3hHEB1FWL1AABQXYWEhGjs2LH6+eefNXv2bHXt2tX42k8//aQRI0aoT58+WrBggfgsFdVdfn6+rr/+et11111GOE9OTtaKFSsI536se/fuWrNmjc4991xJUl5enq644go9+OCDxp9jaYRzANUZK+gA4CdYUQfMZWZmauzYsVqxYoXx2IMPPqgnnnhCISGsNQSC/Px8TZ06VTNnzjQeGzVqlN5//31FR0cbjxHOAVR3BHQA8DMEdaDEvn37NHToUG3ZskWSVKNGDc2aNUvjx4+3eDKcKqfTqZdeeslrF0S/fv305Zdfqnbt2oRzABABHQD8FkEd1d2uXbs0ZMgQ7dy5U5KrZGzOnDnq3r27xZPhTHz77be67LLLjJK/Hj166Nprr9Udd9xBOAdQ7RHQAcDPEdRxtjmdThUVFamwsFBOp1N2u112u12hoaGWzbRt2zYNGTJEe/fulSS1bt1aixYtUvPmzS2bCZUnNTVVF154oX7//fcyXyOcA6jOCOgAECAI6jgVhYWF2rlzp9LT05Wenq6MjAxlZGR43c7IyNCJEyd8lnXZbDaFh4erQYMGio+PV1xcnOLi4srcbtGiRaVee3znzp0aOHCg0fzdvn17ffvtt2rSpEmlfQ9Yb+PGjRoyZIgyMzONx6ZMmaLp06fz7zEA1RYBHQACDEEdpRUWFmrjxo1as2aN1q5dq7Vr1yo1NVX5+flVNkPLli3Vs2dP9ejRQz169FD37t1Vr169U36dffv2aeDAgca29k6dOunbb79Vo0aNKntk+IEtW7aoU6dOKiwslORqfV+0aJHq1Klj7WAAYBECOgAEKIJ69ZWVlaWvvvpKy5cvP+0wHh0drbi4ONWuXdvY0m6322Wz2VRcXGxsec/Ly1NmZqYOHjx4ypf7c4f2nj17avjw4ercuXO578cDBw5o0KBBRiFc+/bt9f3336thw4an9H0RWDZt2qRBgwYZ29379u2rr7/+2qvdHQCqCwI6AAQ4gnr1sHXrVqWkpCglJUXLli3zuS3drW3bturWrZuaNWtmuj29Vq1ap/T9i4qKlJmZWWab/P79+7Vx40atW7dOubm55b5G8+bNlZycrFGjRmnQoEEKDw83vpabm6uBAwdq7dq1klxnzpcuXar4+PhTmhOBaePGjRo0aJAOHTokSbrkkkuUkpJiaQ8CAFiBgA4AQYKgHlyKi4u1YsUKI5Rv3rzZ53Pbtm1rbC3v2bOnunXrptq1a1fhtK4A/9tvvxlb7NeuXVtuaI+JidHw4cOVnJys4cOH69Zbb9VHH30kSWrWrJl++OEHNWvWrCp/BFgsNTVVgwcP1pEjRyRJf/3rXzVt2jSLpwKAqkVAB4AgQ1APbHv27NGMGTP05ptvepVneWrXrp0RbHv27FnlYbyi3KF9yZIlmjt3rr777jvjrLGnkJAQY0dAzZo1tXLlSnXq1Kmqx4UfWLx4sS666CIVFxdLkv7zn//o6quvtngqAKg6BHQACFIE9cDhcDj0zTff6NVXX9W8efPKbF8PCQlR//79je3hiYmJFk16ZrKzs7Vw4UKlpKRo/vz5xkqpp/j4eN1zzz267rrrKrUZHoHjlVde0a233ipJioiI0NKlS9W7d2+LpwKAqkFAB4AgR1D3X4cPH9bbb7+t6dOna9u2bV5fs9vtGjlypMaOHatLLrlEsbGxFk15dhQVFem9997TjTfeaLqqXqNGDV1xxRWaOnWqevToYcGEsIrT6dRNN92kGTNmSHJ9aPPTTz/RRwCgWiCgA0A1QVD3H3v37tXf/vY3/ec//1FeXp7X15o0aaIpU6Zo8uTJiouLs2jCsy8rK0u9e/fWjh07JEn9+/dXZGSkFi1aVOa5vXv31sMPP6wRI0bw/qwmCgoKdMEFF+iHH36Q5HoPfP/994qMjLR4MgA4uwjoAFDNENStk5WVpaefflovv/xymcuiXXDBBZo6dapGjRolu91u0YRVo7i4WBdddJEWL14sSerWrZuWLVumqKgo/fbbb5o+fbrefvttZWdne/26AQMG6JlnnlH//v2tGBtV7ODBg+rVq5f27NkjSZowYYLeeecdi6cCgLOLgA4A1RRBverk5OTohRde0LPPPusVOmNiYnT99dfrpptuCthz5afjH//4h/7yl79Ikho2bKg1a9aoadOmXs/JycnRhx9+qJdffrnM+3PUqFH6+9//TpFcNbBu3Tr179/fuBrAf//7X1122WUWTwUAZw8BHQCqOYL62VNYWKg33nhDTzzxhFcje2RkpG6//Xbdd9991a4I7bffflNSUpLy8/Nls9m0ZMkSDRw40OfznU6nPv/8cz3wwANel5qz2Wy65ppr9MQTT6h58+ZVMTos8sEHH+iqq66SJMXGxurXX39Vw4YNLZ4KAM6OEKsHAABYKyQkRGPHjtXPP/+s2bNnq2vXrsbXfvrpJ40YMUJ9+vTRggULxGe6Fffll1/qnHPO0S233GKE89DQUN1www3aunWrpk2bVu3CeXFxsa6//npje/+dd95ZbjiXXEF83LhxSktL05tvvqkmTZpIcgX3d999V+3atdNf/vIXn9dbR+C74oorNG7cOEnSoUOHjIZ3AAhGrKADALywon5mjh49qrvuuktvv/221+OXXnqpnnzySbVv396awfyA59b2tm3bat26dYqKijql1zh58qReeeUV/f3vf/e6TFvbtm01a9YszqcHqczMTHXs2FFZWVmS2OoOIHgR0AEApgjqp27BggW68cYbtX//fuOx888/X88++6x69epl4WTWK721/YcffjijMH306FE999xz+te//mU04dtsNt1555168sknTzn4w/999NFHuuKKKySx1R1A8GKLOwDAFFvfK+7o0aOaOHGiRowYYYTzmJgYvfnmm1q8eHG1D+dmW9vPdKW7Tp06euqpp7Ru3Tr169dPkmvb+/PPP6+kpCQtX778jOeGfxk/frzXVvdbbrnF4okAoPKxgg4AqBBW1M2ZrZoPGzZMb7zxRplm8urqn//8p+655x5Jp7+1vTzFxcV64YUX9NBDD7GaHuRKb3X/+OOP9ec//9niqQCg8hDQAQCnhKDuUlhYqP/7v//Tyy+/bDwWExOjf/3rX7r++uuD/uevqMzMTLVu3Vo5OTmVsrW9PJs3b9Z1112nlStXGo8lJiZqzpw51eoydsHOc6t7QkKCtmzZoho1alg8FQBUDra4AwBOCVvfXdtrhw0b5hXOhw0bprS0NE2aNIlw7uGpp55STk6OJGnKlClntcQtMTFRy5Yt03PPPaeIiAhJrtDep08fffnll2ft+6JqjR8/Xpdccokkad++fXr11VctnggAKg8r6ACAM1LdVtTT0tKUnJysnTt3SpLCw8P14osv6sYbbwyan7Gy7NixQ+3bt1dhYaGioqK0fft2NW7cuEq+92+//abx48dr/fr1klwfLD377LP6v//7P/6cgsCGDRvUtWtXOZ1O1atXTzt27FDt2rWtHgsAzhgr6ACAM1KdVtTnzJmjfv36GeG8UaNGWrJkiaZMmULoM/Hoo4+qsLBQkvR///d/VRbOJal9+/Zavny5xo4dK8n1QdI999yja6+91jinjsDVuXNnXX311ZKkw4cP69lnn7V4IgCoHKygAwAqVTCuqDudTj311FN6+OGHjcd69OihL774QgkJCRZO5r/Wr1+vpKQkOZ1O1a9fX9u3b7dkhdPhcOiJJ57Q448/bjzWu3dvff7554qPj6/yeVB5du3apXbt2qmwsFA1atTQ9u3bFRcXZ/VYAHBGWEEHAFSqYFtRz8vL0+WXX+4Vzi+//HItXbqUcF6O+++/3/jzfeCBByzbfhwSEqLHHntMn3zyidHmvnr1avXq1Utr1qyxZCZUjhYtWujmm2+WJJ08eVJPPPGExRMBwJljBR0AcFYF8op6Tk6ORo8erUWLFklyXbrr73//u+69916/nNdfLF26VIMGDZIkNW3aVFu2bFFkZKTFU0nr1q3T6NGjtWfPHklSdHS0FixYoAEDBlg8GU7XwYMH1bp1a504cUKhoaHatGmT2rZta/VYAHDaWEEHAJxVgbqinp2drYsvvtgI5zVr1tScOXN03333Ec7/wEMPPWTcfuKJJ/winEtSUlKSfvrpJyOQHz9+XMOGDdPixYstngynq2HDhrrnnnskScXFxV5HGQAgELGCDgCoUoGwon706FENHz5cq1atkuS6vvlXX32lfv36WTZToPjll1/UvXt3Sa6itrS0NIWGhlo8lbeTJ09q7NixWrhwoSQpMjJSn3/+uYYPH27xZDgdx48fV8uWLZWVlaWwsDDt2bOnSgsJAaAysYIOAKhS/r6ifvz4cV188cVGOK9Xr54WL15MOK+g6dOnG7fvuOMOvwvnklSjRg3NmTNHo0aNkuTqGRgzZoy+/fZbiyfD6YiOjtYNN9wgSSosLNRbb71l8UQAcPpYQQcAWMqfVtRzc3N18cUXa+nSpZKkBg0a6Ntvv1WXLl3O+vcOBkePHlWTJk2Um5ur6Oho7d+/X9HR0VaP5VNBQYGuuuoqffrpp5Jcwf2rr77SwIEDLZ4Mp2rXrl1q1aqVnE6nmjZtqh07dshut1s9FgCcMlbQAQCW8pcV9by8PI0ePdoI5/Xq1SOcn6J3331Xubm5kqQJEyb4dTiXpPDwcH3wwQcaPXq0JNfW9xEjRmjlypUWT4ZT1aJFC40YMUKStHfvXs2fP9/iiQDg9BDQAQB+wcqg7nQ6NXnyZGOLc0xMjL7++mvC+SlwOp169dVXjfvuy1/5u7CwMH388cfG+fMTJ05o5MiR2rlzp8WT4VRNnTrVuO35XgSAQEJABwD4FSuC+j/+8Q+9//77kqSoqCh9+eWX6tGjR6W8dnXx3XffafPmzZKkQYMGqWPHjhZPVHERERGaPXu2Bg8eLEnKysrS6NGjdeLECYsnw6kYNmyYWrZsKUn6+uuvtXXrVosnAoBTR0AHAPilqgrqCxYs0L333mvc/89//qNzzz33jGavjjxXLD1XMgNFjRo19NlnnxnX0N6wYYMmTJggh8Nh8WSoqJCQEK+dG6+99pqF0wDA6aEkDgAQEM5GmdymTZvUt29fZWdnS5Iee+wxPfroo5U6d3Vw8OBBxcfHq7i4WI0bN9bu3bsVHh5u9Vin5bffflOfPn2M98Sjjz6qxx57zNqhUGGHDh1SQkKC8vPzVbduXR04cCBg34sAqidW0AEAAaGyV9SPHDmi5ORkI4hdeumlevjhh8/a/MFs3rx5Ki4uluQqhwvkQNS+fXt99NFHxgc9jz/+uNHyDv8XGxurcePGSXL9f/yHH36weCIAODUEdABAQKmMoF5UVKTx48dr27ZtkqSuXbvqnXfeUUgI/1k8HSkpKcbtsWPHWjhJ5bj44ov17LPPGvevvfZarVu3zrqBcErcrfyS93sTAAIBW9wBAAHtdLa+/+Uvf9E//vEPSa5rnf/0009q3rx5lc4dLE6ePKnY2Fjl5uaqYcOGysjICIoPOpxOp6699lr95z//kSQ1a9ZMP//8s+rXr2/xZPgjx44dU2xsrIqKitSyZUtt3779lI6+AICVAv+/oACAau1UV9SXLFlihHO73a7PPvuMcH4GFi9ebFz7fOTIkUERziXJZrNpxowZ6t27tyRpz549uv322y2eChVRu3ZtnX/++ZKknTt36tdff7V2IAA4BcHxX1EAQLVXkaA+btw4TZw40Xj82Wef1XnnnWfFuEHDcwtxcnKyhZNUvsjISM2ePVt16tSRJH3wwQf64osvLJ0JFeP5XmSbO4BAwhZ3AEBQKm/ruySdd955WrJkSdCs+FrB4XAoISFBGRkZioyM1KFDh1SzZk2rx6p0//nPfzRhwgRJUqNGjfTrr7+y1d3P7dq1y7gmep8+fbRy5UqLJwKAiuFvJQCAoFR6Rd295VVyXfN65syZhPMztHbtWmVkZEiSLrjggqAM55J09dVXa9SoUZKkzMxMtroHgBYtWqhz586SpFWrVunAgQMWTwQAFcPfTAAAQS0kJEQXXnihdu/ebTz29NNPq02bNhZOFRzmzZtn3A627e2ebDabXn/9dba6BxjP9+T8+fMtnAQAKo6ADgAIevfff7927twpybW1/bbbbrN4ouCwYsUK4/awYcMsnOTsi4uL04svvmjcv+mmm5SVlWXhRPgjw4cPN26zxR1AoOAMOgAgqH3//ffG9vYaNWpo/fr1rJ5XAqfTqfr16+vIkSNq2LChDhw4EPSXsnI6nRo9erTmzp0rSbryyiv1/vvvWzwVfMnJyVFMTIwcDoeSkpL0yy+/WD0SAPwhVtABAEHL4XDojjvuMO6ztb3y7Nq1S0eOHJEk9ezZM+jDuWS+1d1zFwH8S82aNXXOOedIktLS0pSXl2fxRADwxwjoAICg9dFHHxkN7j169GBreyVau3atcbtHjx4WTlK14uLi9Mwzzxj377vvPrEZ0X+535tFRUXasGGDxdMAwB8joAMAglJBQYEefvhh4/4zzzxDa3slqq4BXZImTZqkdu3aSZKWLl2qL7/80uKJ4Ivne9PzPQsA/oq/qQAAgtIbb7yhHTt2SHJdAuyCCy6weKLgUp0Dut1u11NPPWXcv//+++VwOCycCL4Q0AEEGkriAABB58SJE2rdurUOHjwoSfrpp5/Us2dPi6cKHtWxIK40p9Op3r17a82aNZKk9957T1dddZXFU6E0iuIABBpW0AEAQeeFF14wwvlll11GOK9k1bEgrjSbzeZ1Fv3hhx9WQUGBhRPBDEVxAAINAR0AEFQOHTqkZ599VpIUGhqqJ5980uKJgs+mTZuM20lJSdYNYrGhQ4fqwgsvlCTt3LlTr7/+usUTwYz7PVpUVKRt27ZZOwwA/AECOgAgqPzjH//Q8ePHJUmTJ082yrxQedLT043bzZo1s3AS6z399NPG7SeffJIVWj/k+R7NyMiwcBIA+GMEdABA0Dh58qTeeOMNSVJ4eLgeeeQRiycKTp4hJy4uzsJJrNejRw/96U9/kiQdPHhQn376qcUToTTP9ygBHYC/I6ADAILGf//7Xx0+fFiS9Oc//1nx8fEWTxScPFfQq3tAl6Q777zTuP3qq69aNwhMeb5HPd+7AOCPCOgAgKDhGY6mTp1q4STBzXMVkg9BpHPPPVddunSRJK1YsYKmcD/j+R5lBR2AvyOgAwCCwpo1a7R69WpJrlKovn37WjxR8HKvQtpsNjVq1Mjiaaxns9m8PhCaPn26hdOgNFbQAQQSAjoAICh4hqKpU6dWy0t/VRX3KmSDBg1kt9stnsY/XHXVVYqOjpYkvf/++zp69Ki1A8HAGXQAgYSADgAIeEeOHNEHH3wgSYqJidGVV15p8UTBy+Fw6MCBA5LY3u6pVq1auvbaayVJubm5evfddy2eCG6RkZGqW7euJAI6AP9HQAcABLy3337buLzVddddp5o1a1o8UfDKyspSUVGRJAriSrv55puN26+++qqcTqeF08CT+73KFncA/o6ADgAIeG+++aZx2zMkofLl5OQYt2NiYiycxP906NBB559/viRp8+bNWrZsmbUDweB+r+bl5fHBCQC/RkAHAAS0LVu2aOPGjZJcbdrt27e3eKLg5l49l8T5cxOTJk0ybs+ZM8fCSeDJ873q+R4GAH9DQAcABLS5c+cat8eMGWPdINUEAb18l1xyiUJDQyW5Ajqrtf6BgA4gUBDQAQABLSUlxbidnJxs4STVg2e4CQsLs3AS/1SvXj0NGDBAkrRt2zZt3rzZ4okgeb9XCegA/BkBHQAQsLKysoxzvu3atVNiYqLFEwU/VoT/mOcHRZ4fIME6npdddDgcFk4CAOUjoAMAAtaCBQuMv2yPGjXK4mmqB8+twsXFxRZO4r8834sEdP/Azg8AgYKADgAIWJ7nz9neXjU8A3phYaGFk/ivtm3b6pxzzpEkrVixQr///rvFE8HzvUp3AgB/RkAHAASk/Px8ffXVV5Jc537PPfdciyeqHijbqhj3KrrD4dCCBQssngaUGwIIFAR0AEBA+uGHH3T8+HFJ0ogRI/hLdxUJDw83bp88edLCSfyb546OefPmWTgJJNf1zyUpNDRUISH89ReA/+LfUACAgLRq1Srj9kUXXWThJNVLgwYNjNuZmZkWTuLf+vbtq+joaEne71VY48CBA5KkRo0aWTwJAJSPgA4ACEhr1qwxbvfs2dPCSaqX8PBwxcbGSpLS09MtnsZ/hYaGqnv37pKkvXv36uDBgxZPVH0VFxcbHybFxcVZPA0AlI+ADgAISGvXrpUk1apVS+3atbN4muolPj5ekmtVksuu+dajRw/jtvv9iqp38OBB42oP7vcuAPgrDuwBAALO77//rr1790qSunfvzpnSKhYXF6f169eroKBAhw8fVv369a0eqVKdOHFCBw4cUF5engoLC1VUVKSioiI5nU6FhYXJbrfLbrerRo0aatSokbGVvbTSAf3iiy+uqh8BHjIyMozbrKAD8HcEdABAwPFcjfQMQaganiEnPT09oAJ6fn6+NmzYoC1btigjI0Pp6ele/8zIyDDKByuqVq1aiouLU1xcnOLj441/eu4uYAXdOp5HMQjoAPwdAR0AEHAI6Nby3CackZGhzp07WziNb+4wvmbNGq1du1Zr165VWlpapV+//cSJE9q6dau2bt3q8znz58/XpEmT1KNHD/Xs2VNdunRRZGRkpc4Bc54r6GxxB+DvCOgAgIDjWRBHQK96pVfQ/cXJkye1aNEizZ8/X6tWrTrlMB4dHW2sgMfFxalmzZpeW9olGdvdi4qKlJOTY6y6p6enKzs72+drFxYWaubMmZo5c6Yk17W4O3bsqD59+mjEiBG64IILFBUVdWa/ATDFCjqAQEJABwAEHArirOW5Cml1QM/MzNT8+fOVkpKir7/+utxrs4eEhKh9+/bq0aOHunbtqoSEBK9t6TVr1jyjWXJzc722y+/bt09vvfWWNm7cWOa5RUVFSk1NVWpqqmbMmKHIyEhdcMEFSk5O1siRIwmSlcjzPcoKOgB/R0AHAASU3NxcoyCuS5cuFMRZoFWrVsbttLS0Kv/+mzZt0pw5c5SSkqKVK1eaNsl7hvGePXuqR48eSkpKOuMQXp6oqCi1bt1arVu3Nh6LjY3VtddeK0m6/fbb1bZtW2PL/caNG4128by8PM2bN0/z5s2TJPXu3VvJyclKTk722yMEgcLzPdqyZUsLJwGAP0ZABwAEFM/zpAkJCRZOUn116NBBERERys/Pr7Lys7y8PH3yySd69dVXtXLlStPnNGrUSCNHjlRycrKGDBmiWrVqVcls5fF8j0ZEROjWW2817ufk5GjJkiWaO3euUlJSvN7bq1ev1urVq/XQQw+pZ8+emjp1qsaPH882+FNUXFysdevWSXKF83r16lk7EAD8AZYdAAABhcIn64WHh6tLly6SpC1btpR79vpM7dy5U/fee68SEhI0YcKEMuG8U6dOeuCBB7Ry5Uqlp6frzTffVHJysl+Ec6lsoZ6nmjVrasSIEXrttde0b98+/fTTT3r44YfVtWtXr+etWbNG119/vRISEnT33XeXW0YHb7/99ptyc3Ml0VcBIDAQ0AEAAYXCJ//gGXZ++eWXSn3t4uJizZ8/XyNGjFDr1q317LPPKisry/h6586d9fzzz2v79u3asGGDnnrqKfXp08cvjztUtFAvJCREPXv21BNPPKF169Zp165devHFF9WtWzfjOUeOHNG//vUvtWvXTsOGDdOcOXNUVFR0VucPdFzxAUCg8b//kgEAUA5W0P2DZ9iprG3uxcXFeu+995SYmKiRI0dqwYIFxvnysLAwXXHFFVq2bJlSU1N15513ep2F91cxMTHGtvTSK+jlad68uW677TatXbtWK1eu1IQJExQREWF8/euvv9aYMWPUpk0bzZw5k6DuAwEdQKAhoAMAAgor6P7BM+x4XvbudDidTs2fP1/dunXTNddco+3btxtfa9q0qZ566int3btXH3zwgfr37y+bzXZG368q2Ww24316Oo33NptNffr00TvvvKN9+/Zp2rRpXkVnu3fv1qRJk9SlSxd9/vnnpoV51Znne7N79+4WTgIAFUNABwAEFM9VSAK6dTp27Gis6J7JCvry5cs1cOBAjRw5Uhs2bDAeHzJkiObMmaOdO3fqgQceUKNGjc54Zqu436fHjh0zzkOfjtjYWP31r3/V1q1bNX/+fA0bNsz42qZNmzRu3Dj169dPS5YsOdORg0Lpgrj69etbOxAAVAABHQAQULimsX8406K4tLQ0JScna8CAAVq2bJnxeK9evbRo0SItWrRIycnJCg0NrdS5rVBeUdzpCA0N1SWXXKKvvvpKS5cu1bnnnmt8bdWqVRo8eLAuvvhiI5xWVxTEAQhEBHQAQEA5cOCAJCkyMlK1a9e2eJrqrVevXsbt77//vkK/5vjx45o6daq6dOmiuXPnGo8nJibq008/1apVqzRkyJBKn9VKnjs93O/fynLeeedp2bJlmjNnjjp27Gg8/tVXX6lbt26aNGmSjh49WqnfM1B89913xu2ePXtaOAkAVBwBHQAQUE6ePClJqlWrVkCdRQ5Gw4cPN26npKT84fMXLVqkzp07a/r06cZZ6SZNmuiNN95QWlqaLr300qD8M/W85Jv7/VuZbDabkpOTlZqaqrffflvNmjUzvjZz5kx16tRJX375ZaV/X3/n+Z68+OKLLZwEACqOgA4ACCjutmq73W7xJBg6dKhq1KghSZo7d64cDofp844fP66bb75ZF1xwgXbv3i3JdQ3wadOmaevWrZo8eXJQ/3l6/mxns209NDRU1157rbZs2aLnn39e0dHRkqT9+/frkksu0fXXX19tVtOzs7ONs/jNmzdX586drR0IACqIgA4ACCjugBMWFmbxJIiKitKFF14oScrMzNRPP/1U5jnuVfPXXnvNeOz888/X+vXr9de//tUI+MHM871aFZdDi4iI0J133qm0tDRddNFFxuOzZs2qNqvpCxcuVGFhoSQpOTk5KHdmAAhOBHQAQEBhBd2/JCcnG7c9txT7WjV/+eWXtWjRooC4hnllqaoV9NKaNWumr776Sm+88Ua1W033fC96vkcBwN8R0AEAAYXrPPuXkSNHGquT7lD022+/qUePHqar5rfccotCQqrXXz88V299HQM4m9978uTJpqvp3bp187q0XbAoKirS/PnzJUkxMTEaOHCgxRMBQMVVr/9CAgACnns1sipXIuFbo0aN1Lt3b0muS6fNmjVLffr00datWyVV31VzT+6t1pJ1RzPMVtN37dqlfv366YsvvrBkprNl+fLlOnLkiCRXkWF4eLjFEwFAxRHQAQABhYDufzy3EE+aNMm4JnqXLl2UmppaLVfNPXm+V608muG5mu6+LnhOTo7Gjh2rJ598Mmh2p3hevo/t7QACTfX9ryUAICC5VyAJ6P5j2LBhxm13yBs3bpyWL1+u1q1bWzWW3/B8r/pDuWGzZs30ww8/6IorrjAee/jhhzV+/Hjl5ORYONmZKyws1AcffCDJ1WrP5dUABBoCOgAgoEREREiScnNzLZ4EkpSenq6pU6d6PXbDDTfok08+8br+d3Xm+V71l+3WNWrU0Pvvv69nnnnGOCP/ySefaMCAAdqzZ4/F052+OXPmKCMjQ5I0evRo1atXz+KJAODUENABAAGlUaNGklxbc48fP27xNNXbunXr1KtXL61evdrr8YKCgmq9pb20AwcOGLcbN25s4STebDab7r33XqWkpBjn0tetW6eePXuW+TMNFK+++qpxu/QHRwAQCPivJwAgoMTFxRm33StlqHqrV6/W4MGDlZ6eLklq2rSpEfI++ugjZWVlWTmeX3H/Hkne719/MXLkSK1cudI4jvD777/rggsu0LJlyyye7NRs2rRJ3333nSQpMTFRQ4YMsXgiADh1BHQAQECJj483bhPQrbFs2TJdcMEFxnW0+/XrpzVr1mjy5MmSpPz8fM2aNcvCCf2L+31aq1Yt40MMf9OhQwetXr1a559/viTXdeyHDRumRYsWWTvYKZg+fbpx++abb/a6vB0ABAoCOgAgoLCCbq0lS5Zo2LBhxvGCwYMH6+uvv1bDhg110003Gc977bXXqvya3/7K/T71x9VzT/Xq1dOCBQs0fPhwSa6z8yNGjNDChQstnuyPnThxQu+8844k1/n6a6+91uKJAOD0ENABAAHFM+R4bh3G2bd8+XKNHDnSKD0bNmyY5s+fb5TBtWvXThdeeKEkafv27frmm28sm9Vf5OTkGJed8/eALrnC7RdffGFcniw/P19jxozR4sWLLZ6sfB988IHx+3zVVVepTp061g4EAKeJgA4ACChscbfG6tWrdfHFFxuX4Ro5cqTmzJmjGjVqeD3Ps5jrpZdeqtIZ/ZHne9TzvevPIiIi9Omnn+rSSy+VJOXl5WnUqFF+eybd4XDo5ZdfNu7ffPPNFk4DAGeGgA4ACCisoFe91NRUr23tw4YN06effmpc8s7TyJEj1bRpU0nS/PnztWrVqiqd1d/4e0GcL2FhYfrwww+NlfTc3Fxdcskl+umnnyyerKzPPvtMGzZskCT17dtX3bt3t3giADh9BHQAQEDxXIXcv3+/hZNUD5mZmRo1apRRCHf++edr9uzZpuFckux2ux588EHj/n333Sen01kVo/olz4AeKCvobmFhYfrvf/9rnEk/fvy4kpOTtW/fPosnK1FYWOj1fnvkkUcsnAYAzhwBHQAQUKKjo9WwYUNJ0oYNG6p1+Dvb8vPzNW7cOO3du1eS1KdPH82dO1dRUVHl/rrrr79ebdu2leQqlQuEkrGzZf369cZt92XMAklERIRmz56tgQMHSnJd033MmDE6efKkxZO5zJw5U1u3bpUkDRw40PgwAQACFQEdABBwevToIUk6fPiwdu/ebfE0wcnpdGrq1Kn68ccfJUkJCQn64osvjEK48oSFhenJJ5807t9///3VttF97dq1xu2ePXtaOMnpq1Gjhj777DO1bNlSkutnmjRpkuUfjuXm5urxxx837k+bNo1LqwEIeAR0AEDAcQd0yTsAofK89NJLmjlzpiQpMjJSX3zxhRo3blzhX/+nP/3JOAu8bt06ffzxx2dlTn/mdDqN92eDBg2UkJBg8USnLzY2VnPmzFHNmjUlSR9++KGmTZtm6UwvvviiUcI3ZswY9e3b19J5AKAyENABAAGHgH52ffPNN7rrrruM+zNnzvT6Pa+IkJAQPfPMM8b9hx56SAUFBZU2YyDYs2ePsrKyJLnes4G+utu5c2e99957xv0HHnhAc+fOtWSWI0eOGB8QhISE6KmnnrJkDgCobAR0AEDA8QyLa9assXCS4LN161aNHz/e2JJ+//3364orrjit17rwwgs1dOhQSdKOHTv05ptvVtqcgcDzw6NT/YDDX40ZM0Z/+9vfJLl2CFx55ZX69ddfq3yOZ555xiguvPbaa9WhQ4cqnwEAzgYCOgAg4CQkJBhFcWvXrrX8LGywOHnypMaMGaMjR45IkkaNGuV1lvx0PP3008bthx9+WAcOHDij1wsknh8eBer5czMPPvigLrvsMknSiRMnNHr0aOMSfFXh119/1QsvvCDJVWLneQ4dAAIdAR0AEHBsNhtFcWfBww8/rI0bN0qSOnTooPfee08hIWf2V4VevXrpyiuvlOT6s7r55purzQcqwbiCLrn+/zdr1iwlJSVJkrZv36577723Sr53UVGRrrvuOuO4xF133aWmTZtWyfcGgKpAQAcABCS2uVeuH3/8Uf/6178kuVYlP/30U8XExFTKa7/wwgtq0KCBJOmLL77Qhx9+WCmv68+CqSDOTM2aNfXZZ58Zl9ybPn26Fi9efNa/73PPPWf8/719+/Z69NFHz/r3BICqREAHAASk3r17G7erIhgEs5MnT2rixInGyvYTTzyhc845p9Jev0GDBpo+fbpx/7bbbgv6re7r1683CuJ69eoV8AVxZlq1auXV5D5p0iSdOHHirH2/tLQ0PfbYY5JcxXCzZs1SZGTkWft+AGAFAjoAICANHjxY4eHhkqSUlJRqs236bHj44Ye1ZcsWSVKfPn109913V/r3uPTSSzV+/HhJrq3uN910U1D/maWkpBi3L7nkEgsnObumTp2qQYMGSZJ27dqlv/71r2fl+xQVFWnixInG1va7776by6oBCEo2ZzD/1xEAENQuueQSffnll5Jc533d191Gxf34448aMGCAnE6nIiIi9Msvv1Tq6rmn33//XR07dtTvv/8uSXr//feN8+nBplevXsZW7N27d6tZs2YWT3T27NixQ507d1Zubq4kadGiRRoyZEilfo+nn35aDzzwgCTX1vZffvmF1XMAQYkVdABAwEpOTjZue65YomLO9tb20qrLVvf09HQjnCclJQV1OJfKbnW//vrrK7XV/ddff2VrO4Bqg4AOAAhYI0eONG4T0E/do48+eta3tpdWeqv7+PHjjW3LwWLevHnGbc8PkYKZ51b33bt36/7776+U1z127Jguu+wytrYDqDbY4g4ACGg9e/Y02rL37NnDJZcqaNeuXWrXrp0KCwvP+tb20g4dOqQuXbooIyNDknTTTTd5rawHupEjR2r+/PmSXFcYCKZLrJXHc6t7SEiIfv31V7Vv3/60X6+4uFijR482fi87duyoNWvWsHoOIKixgg4ACGijRo0ybnuuXKJ8jzzyiAoLCyW5ViWrKpxLUmxsrD7//HNFRERIkl577bWgCeg5OTn69ttvJUnx8fHVqhehVatWxsq5w+HQQw89dEav99BDDxnhvG7dupozZw7hHEDQI6ADAAKa5xbiOXPmWDhJ4Fi/fr3ee+89SVK9evXOWvN2efr06aMZM2YY92+//XYtWbKkyueobN98843y8/MluT48CsbLq5XnrrvuUqNGjSRJn332mVavXn1ar/PBBx/omWeekSSFhobqk08+UevWrSttTgDwVwR0AEBAS0pKMra1f/vtt9q3b5/FE/m/Bx980CiGu//++1W7dm1L5pgwYYJx7r2oqEh/+tOftHPnTktmqSzvvPOOcbu6nD/3VLNmTT3yyCPG/fvuu++UL6e3Zs0aTZo0ybj//PPPa+jQoZU2IwD4MwI6ACCg2Ww2TZw4UZLrzOobb7xh8UT+bdmyZcZRgISEBN1yyy2WzjNt2jQNHz5ckpSVlaXRo0frxIkTls50uvbu3WuUFcbHx+uiiy6yeCJrTJ48Wa1atZIkfffdd/rmm28q/GszMjI0ZswY5eXlGa916623npU5AcAfEdABAAHvhhtuUGhoqCRpxowZxtlqeHM6nbrvvvuM+4899phq1Khh4USu7csffvih2rVrJ0nasGGDrrrqKhUVFVk61+mYMWOGHA6HJGnKlCmy2+0WT2SN8PBwPfnkk8b9+++/3/h9Kc+JEyc0duxY7d+/X5LUv39/vfLKK9XumACA6o2ADgAIeAkJCRo9erQk6cCBA/riiy+sHchPzZs3T8uXL5cktW/fXtdee63FE7nUqVNHKSkpxlb7lJQUTZgwQcXFxRZPVnEFBQXG7g273a7JkydbPJG1xo8fr6SkJEnSzz//rE8++aTc5+fm5mrUqFFatWqVJKlp06b67LPPFB4efrZHBQC/QkAHAASFqVOnGrdfffVVCyfxT06nU48//rhx/6mnnvKrFd7ExER99tlnRrP7hx9+qMmTJ1do5dUffP7558rMzJQkjR07VvHx8RZPZK2QkBA9/fTTxv3HH3/c51n0vLw8jR071igJrFOnjubOnWuUzQFAdUJABwAEhSFDhigxMVGStGTJEm3cuNHiifzL6tWrjevFd+vWTWPHjrV4orKGDh2qTz/91Pjg4O2339bNN98cECH9lVdeMW57flhUnQ0bNkz9+/eXJG3atMm0pT8/P19//vOf9fXXX0uSoqOjtXDhQnXt2rUqRwUAv0FABwAEBZvNpptvvtm4HyzX1a4snrsKbrvtNr891zty5Eh99NFHXp0CkydP9uvt7hs2bNAPP/wgSTrnnHM0aNAgiyfyDzabTbfddptxv/TOlpMnT2rMmDGaO3euJCkqKkoLFixQ7969q3ROAPAnNuepXvsCAAA/dfToUcXHx+vkyZOKjo7W3r17LbuEmD85dOiQEhISlJ+fr7p162rfvn2Kioqyeqxyffjhh7rmmmuMYH7llVfqnXfe8att+W5Tpkwxrun+0ksv0TruoaCgQM2aNVNmZqZCQ0O1Z88excfHKycnR8nJyVq8eLEkVzifN2+eBg8ebPHEAGAtVtABAEGjTp06uvrqqyVJx48f13PPPWfxRP5h1qxZys/PlyRNnDjR78O5JF1xxRX6+OOPjUD+wQcf6LLLLlNOTo7Fk3nbvn27Zs6cKUmqVauWrrnmGosn8i/h4eG64YYbJJVcBvHQoUMaPny4Ec6jo6P11VdfEc4BQKygAwCCzM6dO5WYmKjCwkJFRUVp27ZtiouLs3osyzgcDrVp00Y7d+6UJG3ZskVt27a1eKqKmzt3rv70pz+poKBAkpSUlKQ5c+aoWbNmFk/mcuWVV+rDDz+UJD388MN64oknLJ7I/+zdu1ctWrSQw+FQw4YNVaNGDe3evVuSVLt2bS1cuFB9+vSxeEoA8A+soAMAgkrLli2Ns+i5uble12OujhYuXGiE84suuiigwrkkjRo1SnPnzlV0dLQkad26derZs6eWLVtm8WSuWdzhPDY2Vvfcc4/FE/mnpk2bKjk5WZJ08OBBI5w3btxYixcvJpwDgAcCOgAg6Dz44IOqVauWJFfJ2LZt2yyeyDqexVyB2i5+0UUXaeXKlWrdurUk6ffff9eQIUP05ptvWjrX/fffb9x+8MEHFRMTY+E0/svpdJb5venZs6fWrFmj7t27WzQVAPgntrgDAILSY489Zlz3+/LLLzdWOquTzMxMxcXFyel0qmnTptqxY4dflqxV1OHDhzV+/Hh9++23xmO33Xab/vWvf1X5z7VkyRLjzHSzZs20ZcsW4xruKJGTk6OJEyfqk08+8Xp88+bNateunUVTAYD/YgUdABCU7r77bsXGxkqSPvroI/3yyy8WT1T15s+fL/fn8FdddVVAh3NJqlevnr788kvdcccdxmMvvfSShg0bpkOHDlXZHE6nU/fdd59x/4knniCcm9i1a5cGDBhQJpxL0jfffGPBRADg/wjoAICgFB0drYceesi477kdubpISUkxbo8ePdrCSSqP3W7XCy+8oLfeekthYWGSpMWLF6tjx46aPXt2lczwxRdfaNWqVZKkTp06GVcOgIvT6dSMGTPUpUsXrVu3TpLr/48vvvii8RzP9yYAoARb3AEAQSs/P1/t27fXrl27JEn//e9/ddlll1k7VBU5efKkYmNjlZubq4YNGyojI0MhIcH1ufzy5cs1btw4HTx40Hjs8ssv10svvWTsnqhs2dnZ6ty5s/bs2SPJFTRHjRp1Vr5XINq9e7cmT57sdQyhTZs2SklJUfv27dWyZUvt3r1bYWFhOnToEOf2AaCU4PovNQAAHiIiIjRt2jTj/tSpU/X7779bOFHVWbx4sXJzcyW5mtCDLZxLUv/+/ZWamuq1O+Cjjz46q6vpf/3rX41wPmTIEI0cOfKsfJ9A414179y5s1c4nzRpktasWaNzzjlHNpvNaHMvLCzUwoULrRoXAPxW8P3XGgAAD5dddpnGjRsnSTp06JBuueUWiyeqGp5biN2hKBg1btxYn3/+ud577z3VrVtXkutSXpdeeqmuuOKKSj2b/u233+r111+XJNWsWVNvvvmmbDZbpb1+oNq9e7cuuugiTZkyRcePH5ckJSQk6Msvv9Sbb76p2rVrG8/1fC/OnTu3ymcFAH/HFncAQNDLzMxUx44dlZWVJSn4t7o7HA4lJCQoIyNDkZGRysrKUlRUlNVjnXUHDhzQTTfdpDlz5hiPNWzYUNOmTdM111yj0NDQ037t0lvbX331Vd18881nPHMgy8/P1/Tp0/XII48YwVxyrZr/85//9ArmbgUFBWrQoIGys7NVr149ZWZmBnx5IQBUJlbQAQBBr1GjRnr55ZeN+1OnTvU6txxsfv75Z2VkZEiSLrjggmoRziXfq+kTJ05U165dlZKSotNdl/jLX/5ihPPBgwdrypQplTZ3oCkuLta7776rxMRE3XXXXX+4au4pPDxcw4cPl+S6bN6PP/5YZXMDQCAgoAMAqoXx48d7bXW/9dZbLZ7o7FmwYIFxO5i3t5ux2Wy66qqrtHHjRo0dO9Z4/Ndff9Xo0aM1YMAA/fDDD6f0mt9++61mzJghybW1/a233grKM/1/xOl0au7cuUpKStK1116r3bt3G1+bNGmS0tLSjPBdHs/3pOd7FQBAQAcAVBM2m02vvvqq6tevL0n65JNP9N///tfiqc6O1atXG7eHDBli4STWady4sWbPnq0lS5aob9++xuM//vijBg4cqJEjR2r9+vV/+DrZ2dmaNGmScf+5555Ty5Ytz8rM/mzZsmU677zzlJycrLS0NOPx4cOH6+effy531bw0z/ek53sVAMAZdABANfPRRx/piiuukCTFxMRo1apVat++vcVTVa64uDgdOHBAdevWVVZWVrUvMnM6nZozZ44eeOABbdq0yXjcZrPpz3/+s26//Xb169evzO+Tw+HQuHHjjDPtgwcP1rffflttVs8dDoe++eYb/fvf/9aXX37p9bU+ffromWee0fnnn39ar92kSROlp6erdu3aOnLkSLV/jwKAW/X4LwwAAP8zfvx4jR8/XpJrdTQ5OVlHjhyxeKrKk56ergMHDkiSunfvTvCRK4iPGTNGGzZs0MyZM9W0aVNJruD+8ccfq3///urWrZtmzJihEydOGL/u0UcfNcJ5nTp1qs3W9sOHD+uf//yn2rVrp+HDh3uF8/bt22v27NlasWLFaYdzSerRo4ck6dixY9q+ffuZjgwAQSP4/ysDAIAHm82mt956S127dpUkbd26VZdffrmKioosnqxyrF271rjtDkFwCQ0N1cSJE7Vlyxb985//VGxsrPG11NRUTZkyRU2aNNHtt9+uf/3rX3ryySclSSEhIfr444+Dfmv7mjVrNHHiRDVp0kT33HOPV3Bu2rSp3nzzTW3YsEFjx4494w9+PN+bnu9ZAKjuCOgAgGqnZs2amjNnjhHQvv76a917770WT1U5COh/LDIyUv/3f/+nPXv26J133lGfPn2Mr2VnZ+ull17S3XffbTz2+OOP66KLLrJi1LMuPT1dr732mnr37q1evXrp7bffVl5envH1oUOH6rPPPtOOHTs0adKkSrskGgEdAMxxBh0AUG0tXbpUQ4cONVbPZ82apeuuu87aoc7QqFGjNG/ePEnS9u3b1apVK4snCgxr167V9OnT9f7773sFVMm162LAgAFKTk5WcnKy2rVrZ9GUZ87pdGr9+vVKSUlRSkqK1qxZU+Y5tWvX1nXXXaebbrrprPUzZGRkKD4+XpKrNG7RokVn5fsAQKAhoAMAqrUZM2YY17QODw/XkiVL1K9fP4unOn3x8fHKyMigIO405Ofna9CgQVq1alW5z0tMTFRycrJGjRqlfv36Vdqq8tmSn5+v77//XikpKZo7d65xPffSkpKSdMstt+iKK65QzZo1z/pc7vdqnTp1dPjwYd6rACACOgAAuvXWW/XKK69Ikho1aqTly5erdevWFk916jxXJYcOHapvv/3W4okCh8Ph0MSJE/Xuu+9KcrWMv/zyy1qxYoVSUlL022+/mf66GjVqKCkpST179lSPHj3Uo0cPtW/f3rLQXlBQoLS0NK1du1Zr167VmjVrtGHDBhUUFJg+v1u3bsbOgG7dulVpSPbc7bFt27aA/P8cAFQ2AjoAoNorLCzUsGHD9N1330mSmjVrpqVLl6p58+YWT3ZqfvjhBw0cOFCSdMcdd+iFF16wdqAA4XQ6NXXqVL322muSXGfUf/jhB/Xs2dN4zpYtWzR37lylpKRo2bJlcjgcPl/PM7R37txZTZo0UVxcnOLj49WgQYMzboIvKirSwYMHlZ6eroyMDO3bt0+pqal/GMYl1y6RwYMHGzsA3I32VnjggQf09NNPS5IWLlwYtOf8AeBU+PeeLAAAqkBYWJg++eQTnX/++UpLS9OePXs0dOhQff/992rSpInV41VYRkaGcTuQ5raS0+nUXXfdZYTz0NBQvf/++17hXJLatWunu+++W3fffbeysrK0YMECLViwQKtXr9aOHTu8nnvy5EmtWLFCK1asKPP9QkND1bhxY8XFxRn/i4qKkt1uN/4nuUJ4YWGhioqKdOLECWVkZBj/y8zMVEXXV2w2mxITE9WnTx+NGjVKF110kaKjo0/nt6rSuXd7SN7vXQCozgjoAABIql+/vr799lsNHDhQW7Zs0fbt2zV06FAtXrzYK0j4s/T0dON2oMxsJafTqXvvvVf//ve/JbnC7Lvvvqtx48aV++vq16+va665Rtdcc40k6ciRI/r555+NLeVr164tE9rdiouLtX//fu3fv79yf5j/zd++fXtjq32PHj2UlJTkN4G8NM/3qOd7FwCqMwI6AAD/06hRIy1evFgDBw7Ujh07tHnzZg0cOFCLFy9Ws2bNrB7vD3muQsbFxVk4if9zOp2688479eKLLxqPvfXWW7ryyitP+bXq1q2roUOHaujQocZj7tC+bds2ZWRkGNvR3f/MzMwsd5u8Gbvdbqy+x8fHe/0zMTHRr8O4Gc/3KCvoAOBCQAcAwEOTJk20ePFiDR48WDt37tT27duNkO7vlywjoFeMw+HQ1KlT9frrrxuPvfbaa5o4cWKlfQ+z0O6puLhYBw8e1IEDB5Sfn6+ioiJjW7vNZjO2u4eFhSkyMlJxcXGKjY094/Pr/oQVdAAoi4AOAEApzZs31/fff6+hQ4dq69at2r17twYOHKh58+YpKSnJ6vF8Yov7H8vLy9PkyZP1/vvvS3JtC585c6auu+66Kp0jNDTUOINeXTVu3Ni4zQo6ALgEz8ewAABUoqZNm+r7779Xhw4dJEn79+9X//799cknn1g8mW/ukFOjRg3FxMRYPI3/SU9P16BBg4xw7i6Eq+pwDpeIiAjVr19fEivoAOBGQAcAwIe4uDgtWbJEvXv3liTl5ubqz3/+sx599NFTPj9cFdwBPT4+vkqvZx0IfvrpJ/Xq1UurV6+WJEVFRenTTz/VFVdcYfFk1Zt7B0FGRkaFm+kBIJgR0AEAKEeDBg30/fffG43dkvTEE0/oT3/6k06cOGHhZN4KCgp05MgRSa6yO5R47733dN555xmrtM2aNdPy5cs1ZswYaweDEdDz8/N17Ngxi6cBAOsR0AEA+AORkZF655139I9//MMo6fr888917rnnaufOnRZP51JQUGDcjoqKsnAS/1FcXKx7771X11xzjfLz8yVJAwYM0E8//eTXXQLVSY0aNYzbnu9hAKiuCOgAAFSAzWbT3XffrXnz5ql27dqSpA0bNqhXr15atGiRxdNJRUVFxm27nQ7Yw4cPKzk5Wc8++6zx2A033KBFixapYcOGFk4GT57vVc/3MABUVwR0AABOwcUXX6xVq1apXbt2kqSsrCxdcMEFuu222yzd8l5YWGjcru4Bfc6cOerQoYMWLFggyVUG9/LLL+v1119XeHi4xdPBk+d71fM9DADVFQEdAIBTlJiYqFWrVmn48OHGYy+//LK6dOmiJUuWWDITK+iuVfOrr75aY8aMUWZmpiSpXr16+vrrr3XLLbdQnOeHwsLCjNsEdAAgoAMAcFrq1Kmj+fPn64UXXjDO0e7cuVODBw+2ZDW9uodP96q5+xJqkjRixAht2LBBQ4YMsXAylMezud3d7wAA1Rn/JgQA4DSFhITojjvu0Pr163XeeecZj1uxml5dz/IePnxY11xzjdeqeZ06dfTOO+9o7ty5io+Pt3hClIedHwDgjYAOAMAZatOmjZYsWWK6mn7jjTcal/c6m6rbVuHi4mK9/fbb6tixo9577z3j8REjRigtLU0TJkyo9rsKAoHne9XzPQwA1RUBHQCASuBeTU9NTdWAAQOMx9944w21adNG999/v44ePXrWvn91WUF3Op2aM2eOunTpookTJ+rAgQOSpNq1a+vtt9/W3Llz1aRJE4unREWxgg4A3gjoAABUorZt2+r777/XCy+8oOjoaEnSyZMn9cwzz6hVq1Z67rnndPLkyUr/vhEREcaK8fHjxyv99f3B0qVL1b9/f40ZM0YbN240Hk9OTtavv/6qa6+9llXzAOP5XvW8JjoAVFcEdAAAKpl7NX379u268847jUt7HTlyRH/961/Vtm1bvfnmm5W60m2329WgQQNJUkZGRqW9rj9Yv369RowYoUGDBmnFihXG4/369dP333+vOXPmsGoeoNzv1Vq1aqlWrVoWTwMA1iOgAwBwljRo0EDPP/+8Nm/e7HUmev/+/brhhhvUsWNHvfrqq8rOzq6U7xcXFydJOnDggBwOR6W8plWcTqe+//57XXbZZUpKSjKuaS5JHTp00BdffKHly5dr4MCBFk6JM+XuZ6DMDwBcCOgAAJxlLVq00DvvvKP169dr1KhRxuNbtmzRLbfcoiZNmmjq1KlKS0s7o+/jDjmFhYXKyso6o9eySnZ2tl555RV16tRJ559/vj799FPjUlxNmzbVrFmztH79eo0ePZrt7AHuxIkTxhZ394dLAFDdEdABAKginTp1UkpKipYtW6ZBgwYZj584cULTp09X586dNXDgQH300UcqKCg45df3DDmBts19w4YNuvnmmxUfH69bb73V64x5w4YN9c9//lNbtmzRddddp9DQUAsnRWXxfI+ygg4ALgR0AACqWP/+/bVkyRKtW7dOU6ZMUc2aNY2v/fDDD7riiivUrFkz3XvvvVq5cmWFt6t7BvSquLTbmUpPT9eMGTN03nnnqUuXLnrttdeUk5NjfP28887Thx9+qL179+r//u//FBkZaeG0qGyeAZ0VdABwIaADAGCRrl276rXXXtP+/fv10ksvqUOHDsbXMjMz9eyzz6pfv36Ki4vT5MmTlZKSotzcXJ+v57kK6Y8r6E6nU6mpqfrb3/6mXr16qUmTJpoyZYqWLVtmPKdWrVq6+eabtX79ei1dulSXX365UbKH4OL5IRIr6ADgwgUnAQCwWO3atXXrrbfqlltu0dKlS/Xqq69q9uzZRsv7wYMH9dZbb+mtt95SZGSkLrjgAiUnJ+viiy9WQkKC8Tr+uMU9NzdXy5cvV0pKilJSUrRnzx7T53Xo0EG33HKLrr76asXExFTxlLACK+gAUBYBHQAAP2Gz2TRo0CANGjRIBw8e1Pz585WSkqKvv/7aWDnPy8vTvHnzNG/ePElS48aN1aNHD/Xo0cPrMlV79+6t8vlzc3OVmpqqNWvWaO3atVq7dq02bdqk4uJi0+d369ZNo0aN0qhRo9SjRw9K36oZz/coAR0AXGxOdzUqAADwSydPntTixYuVkpKiuXPnVmh1vE6dOrrlllvUrFkzxcfHKy4uTnFxcWrYsKHs9tP7fN7pdCo7O1vp6enKyMhQRkaG0tPT9euvv2rt2rXauHFjueflw8LCNGTIECUnJ2vkyJFq1qzZac2B4DBkyBB99913klzb3QnpAEBABwAgoDgcDv3888+aO3eufvzxR61du1ZHjhyp8K8PCQlRw4YNjdBeu3Zt2e12438hISEqKipSUVGRCgsLlZeXp8zMTCOMnzx5ssLfy263q2PHjurZs6eGDx+uiy66iO3rkOT6sKdu3bo6duyY4uPjtX//fqtHAgC/QEAHACCAOZ1O7dq1y9hW/vbbbyszM7PK5/AM4+4t9126dKF5Haa2bdumtm3bSpJGjRqllJQUiycCAP/AGXQAAAKYzWZTy5Yt1bJlS1122WWqX7++/vrXv0qSbr/9dnXs2NFrO7r7n5mZmT7PhpupU6eO4uLijJV3z3+2aNGCMI5TsnbtWuN2jx49LJwEAPwLAR0AgCDiGXaKi4t14403mj6vuLhYv//+u3Jycry2tDscDoWFhclutyssLExhYWFq2LChatSoUVU/AqoBAjoAmCOgAwAQRLp3727c9gxBpYWGhqpx48ZVMRJQBgEdAMxxBh0AgCDTpk0bbd++XTVq1FB2dvZpt7YDZ4NnQVxcXJzS09OtHgkA/EaI1QMAAIDK5V6RPHnypH799VeLpwG8bdu2TceOHZPE6jkAlEZABwAgyPTr18+4vXDhQgsnAcryfE96vlcBAAR0AACCzsiRI43bc+fOtXASoCzP9+SoUaMsnAQA/A9n0AEACEIdOnTQpk2bFBISogMHDqhBgwZWjwQoOztbsbGxKiwsVPPmzbVz507ZbDarxwIAv8EKOgAAQSg5OVmS5HA4tGDBAounAVwWLlyowsJCSa73KOEcALwR0AEACELugC5JKSkpFk4ClPB8L3q+RwEALmxxBwAgCBUXF6tx48Y6dOiQatasqUOHDikyMtLqsVCNFRUVqWHDhjpy5IhiYmL0+++/Kzw83OqxAMCvsIIOAEAQCg0NNcricnJy9N1331k8Eaq75cuX68iRI5Kkiy++mHAOACYI6AAABCnPLcRffPGFdYMAkubMmWPcZns7AJhjizsAAEHqxIkTatiwoU6ePKmYmBjt379ftWrVsnosVEN5eXlq2rSpDh06pLCwMGVmZqpu3bpWjwUAfocVdAAAglStWrV0xRVXSHJd3uqDDz6weCJUV59++qkOHTokSbr00ksJ5wDgAwEdAIAgNnXqVOP2K6+8IjbOwQqvvvqqcdvzPQkA8MYWdwAAglyfPn20evVqSa6irnPPPdfiiVCd/PLLL+revbskqVOnTlq/fj3XPwcAH1hBBwAgyHmuWHquZAJVYfr06cbtqVOnEs4BoBysoAMAEOROnjyphIQEHT58WOHh4dq7d68aNmxo9VioBo4ePaomTZooNzdXtWrVUnp6uqKjo60eCwD8FivoAAAEuRo1auj666+XJBUUFGjmzJkWT4Tq4t1331Vubq4kacKECYRzAPgDrKADAFANbN++XW3atJEkNW3aVFu3blVERITFUyGYFRUV6ZxzztG2bdskSRs2bFCnTp0sngoA/Bsr6AAAVAOtW7fWiBEjJEl79+7Va6+9ZvFECHazZs0ywvngwYMJ5wBQAaygAwBQTXi2acfGxmrHjh1sOcZZkZubq7Zt2yo9PV0SVw8AgIpiBR0AgGqiW7duuvzyyyVJhw4d0j//+U+LJ0Kwevnll41wnpycTDgHgApiBR0AgGpk27ZtOuecc1RUVKRatWpp+/btNLqjUh05ckStWrXS0aNHZbPZtGHDBnXs2NHqsQAgILCCDgBANdKmTRvdcMMNkqQTJ07oqaeesngiBJtp06bp6NGjklzN7YRzAKg4VtABAKhmMjIy1KZNG+Xm5iosLExbtmxRixYtrB4LQWD//v1q06aN8vLyFB4eri1btqh58+ZWjwUAAYMVdAAAqpm4uDjdeeedkqTCwkI9/PDD1g6EoPH4448rLy9PkjR16lTCOQCcIlbQAQCoho4dO6ZWrVrp8OHDkqSFCxfqoosusngqBLKlS5dq0KBBkqTo6Gjt2LFDsbGxFk8FAIGFFXQAAKqh2rVr6+9//7txf/LkyTp27JiFEyGQ5eTk6PrrrzfuP/bYY4RzADgNBHQAAKqpG2+8UUOHDpUk7d27V3/5y18sngiB6oEHHtD27dslSf3799cdd9xh8UQAEJjY4g4AQDW2a9cude7cWSdOnJDEVnecOs+t7ZGRkUpNTVW7du0sngoAAhMr6AAAVGMtWrTQc889Z9xnqztORemt7X//+98J5wBwBgjoAABUc1OmTGGrO07Lgw8+6LW1/fbbb7d4IgAIbGxxBwAAZba6f/nllxo+fLjFU8Gfff/99xo8eLCcTidb2wGgkrCCDgAAymx1v+qqq4yVUaC0PXv26M9//rPc6zxsbQeAysEKOgAAkCQ5nU4lJydr3rx5kqQOHTpoxYoViomJsXgy+JPc3FwNGDBAv/zyiyRp6NChWrhwoUJDQy2eDAACHyvoAABAkmSz2fTee++pffv2kqSNGzfq6quvlsPhsHgy+Aun06nrr7/eCOetWrXSxx9/TDgHgEpCQAcAAIbatWsrJSVFderUkSTNnTtXjzzyiLVDwW88/fTT+vjjjyVJtWrVUkpKiurXr2/xVAAQPNjiDgAAyvjmm280fPhwY/X8o48+0vjx4y2eClaaM2eOxowZI8m122LOnDkaNWqUtUMBQJBhBR0AAJRx4YUX6l//+pdxf+LEiVq7dq2FE8FKaWlpuvrqq437Tz31FOEcAM4CVtABAIApp9OpyZMna+bMmZKk+Ph4/fDDD2rVqpXFk6Eq7d27V+edd552794tSbr88sv1wQcfyGazWTwZAAQfAjoAAPApPz9fQ4YM0Y8//ihJat68uZYuXapmzZpZPBmqQkZGhgYOHKht27ZJknr06KGlS5cqKirK4skAIDixxR0AAPgUERGhOXPmqGPHjpKk3bt3a8iQIUpPT7d4MpxtBw8e1NChQ41w3rZtW82dO5dwDgBnEQEdAACUKzY2Vt9++63atWsnSdq+fbuGDBmi/fv3WzwZzpbMzEwNHTpUmzZtkiS1aNFCixYtUlxcnMWTAUBwY4s7AACokH379mngwIHauXOnJKl169ZavHgx292DTHp6uoYOHarffvtNkpSQkKClS5eqZcuWFk8GAMGPFXQAAFAhCQkJ+u6774ySuO3bt2vgwIHavn27xZOhsuzZs0cDBw40wnnTpk313XffEc4BoIoQ0AEAQIW5S+Lc2913796t8847T6tXr7Z4MpypdevWacCAAcYHLi1bttTSpUvVpk0biycDgOqDgA4AAE5JkyZN9P333xvFce6m7/fee8/iyXC6Pv30U/Xv31979+6V5CqEW7p0qVq0aGHtYABQzRDQAQDAKWvcuLGWLFmi8847T5LrcmzXXHON7r33XhUXF1s8HSrK4XDo0Ucf1WWXXabc3FxJUp8+fbR06VIlJCRYPB0AVD+UxAEAgNNWUFCg2267TTNmzDAeu+SSS/TBBx+odu3aFk6GP3LixAlNmDBBn3/+ufHYhAkT9PrrrysyMtLCyQCg+iKgAwCAM+J0OjV9+nTdfvvtxup5+/btlZKSorZt21o8Hczs3LlTo0eP1oYNGyRJISEheu6553TXXXfJZrNZPB0AVF8EdAAAUCkWL16syy67TIcPH5Yk1alTR++8846Sk5MtngyevvrqK1199dXKysqSJNWuXVsfffSRhg8fbvFkAADOoAMAgEoxZMgQ/fTTT0Z53NGjRzV69GhNmDBBR44csXg6HDt2TJMnT9bFF19shPN27dpp1apVhHMA8BMEdAAAUGlatWqlFStWaOzYscZj//nPf9SxY0fNnTvXwsmqt4ULF6pTp0566623jMcuvvhirVq1SomJiRZOBgDwREAHAACVKjo6Wp999plmzZplFMVlZGQoOTmZ1fQqduzYMd1www0aPny49u3bJ0mqVauWXn/9dc2fP1916tSxdkAAgBfOoAMAgLNm//79uvHGG7VgwQLjsbi4OL3++usaNWqUhZMFv4ULF2ry5MlGMJck+/+3d38hVtf5H8ffMx5dWysDdxZNQzfLC00vnPBCLLyYDIlGt98SsSDkj/6gN/GT6qIuyi6i2DvDMjPKZW9aLLaZoshouwhJSMICTTeCUBsJArMVdWea+V24Tk7OOGfOnPP9fr7f7+MBYsyM+r05L3rzPI61Wtx5553R09MT7e06DUBqLDMA0DJz586Nd955Z9Safs8998RXX32V8xOWz9dffx333XffZdV89erVMTAwEO+++27s2bMn56cEYDQKOgCQidFqent7e2zcuDGefvrpmDdvXo5PV3x9fX3xzDPPxK5du2JgYGD4411dXbFr1644evRorFmzJiIiFi9eHF9++aWKDpAYqwwAZOJiTX/99ddj9uzZERExODgYr776atx0003x2GOPDf8TbdTv1KlT8cQTT8TChQtjx44dw8d5R0dH7Ny5Mz744IOYP39+dHV1xcqVKyMi4tChQyo6QIIUdAAgc2fOnIlt27bF888/Hz/++OPwx2fOnBmPP/54PPLIIzFjxowcnzB9Z8+eje3bt8ezzz474hvvXX311fHoo4/Gli1b4pprrhnxa/bu3auiAyTMgQ4A5OaHH36I5557Ll544YU4f/788Mc7OjrigQceiIcffjjmz5+f4xOm5/jx47Fz58545ZVX4uTJk8Mfnzp1amzevDmefPLJ6OjoGPXXDg0NxapVq2Lfvn0REfHGG2/Evffem8lzAzA+BzoAkLtjx47F1q1b47XXXovBwcHhj7e3t8ddd90VmzdvjjVr1lS29g4NDcVHH30U27dvj56envj555+HP9fW1hYbNmyIrVu3xoIFC8b9vVR0gHQ50AGAZBw+fDi2bt0ab7755ohvdBYRsXDhwti0aVPcf//9MWvWrJyeMFunTp2K3bt3x0svvRRHjhwZ8bkpU6bE+vXr46mnnoqlS5fW/Xuq6ADpcqADAMnp6+uLXbt2xcsvvxwnTpwY8bnp06fHunXrYt26dbF27dq47rrr8nnIFjl9+nS8//770dPTE2+99VacPXt2xOdnz54dDz30UDz44IMNf+d7FR0gTQ50ACBZAwMD0dvbGy+++GJ8+OGHl32+VqvF7bffHt3d3XH33XfHjTfemMNTTt63334bvb290dPTEx9//HH09/df9jWrV6+OzZs3x/r162Pq1KmT+vNUdIA0OdABgEI4cuRI7NixI3bv3j3iu5ZfasmSJdHd3R133HFHdHZ2xrXXXpvxU9bnp59+is8//zz27t0bvb29cfDgwVG/bubMmbFhw4bYtGlTLF68uKnPoKIDpMeBDgAUSn9/f3zyySfR29sbb7/9dnzzzTdjfu2iRYuis7Nz+Mfy5cszP9ovHuMHDhwY/nHkyJEY63/BFixYEN3d3dHd3R233XZbTJs2rSXPpaIDpMeBDgAU1tDQUBw+fDh6enqip6cnPv300zEP34sWLVoUy5Yti3nz5sWcOXPi+uuvH/HzzJkzo62tre4///Tp09HX1xfffffdiJ9PnDgRX3zxxRWP8YtWrFgxfJTfcsstdf/5k6WiA6TFgQ4AlMb3338f7733Xuzfvz8OHDgQBw8eHPHvq9dj+vTpMWfOnJgxY0ZMnTo1arVa1Gq1iLjwd+Iv/jhz5kz09fVd9k3cxjNt2rRYtmxZdHZ2xooVK2Lt2rUxZ86cCf0ezaKiA6TFgQ4AlFZ/f38cOnQoPvvss+G3lzdytDfq0mO8s7Mzbr311liyZEnL3rbeCBUdIB0OdACgUvr7++PYsWOjvi394n+fPHkyzp07F/39/TEwMBCDg4MREdHe3j5c1KdPnx6zZ88efmv8aG+Xv+GGG5I6xkejogOkw4EOADCOSw/0MlLRAdJgeQEAxtHe3l7qg7WrqytWrlwZERGHDh2KPXv25PxEANWkoAMAoKIDJMDqAgCgogMkQEEHACAiVHSAvFlcAAAiQkUHyJuCDgDAMBUdID/WFgCAYSo6QH4UdAAARlDRAfJhaQEAGEFFB8iHgg4AwGVUdIDsWVkAAC6jogNkT0EHAGBUKjpAtiwsAACjUtEBsqWgAwAwJhUdIDvWFQCAManoANlR0AEAuCIVHSAblhUAgCtS0QGyoaADADAuFR2g9awqAADjUtEBWk9BBwCgLio6QGtZVAAA6qKiA7SWgg4AQN1UdIDWsaYAANRNRQdoHQUdAIAJUdEBWsOSAgAwISo6QGso6AAATJiKDtB8VhQAgAlT0QGaT0EHAKAhKjpAc1lQAAAaoqIDNJeCDgBAw1R0gOaxngAANExFB2geBR0AgElR0QGaw3ICADApKjpAcyjoAABMmooOMHlWEwCASVPRASZPQQcAoClUdIDJsZgAADSFig4wOQo6AABNo6IDNM5aAgDQNCo6QOMUdAAAmkpFB2iMpQQAoKlUdIDGKOgAADSdig4wcVYSAICmU9EBJk5BBwCgJVR0gImxkAAAtISKDjAxCjoAAC2jogPUzzoCANAyKjpA/RR0AABaSkUHqI9lBACgpVR0gPoo6AAAtJyKDjA+qwgAQMup6ADjU9ABAMiEig5wZRYRAIBMqOgAV6agAwCQGRUdYGzWEACAzKjoAGNT0AEAyJSKDjA6SwgAQKZUdIDRKegAAGRORQe4nBUEACBzKjrA5RR0AAByoaIDjGQBAQDIhYoOMJKCDgBAblR0gF9YPwAAcqOiA/xCQQcAIFcqOsAFlg8AgFyp6AAXKOgAAORORQdQ0AEASICKDqCgAwCQCBUdqDqLBwBAElR0oOoUdAAAkqGiA1Vm7QAASIaKDlSZgg4AQFJUdKCqLB0AAElR0YGqUtABAEiOig5UkZUDACA5KjpQRQo6AABJUtGBqrFwAAAkSUUHqkZBBwAgWSo6UCXWDQCAZKnoQJUo6AAAJE1FB6rCsgEAkDQVHagKBR0AgOSp6EAVWDUAAJKnogNVoKADAFAIKjpQdhYNAIBCUNGBslPQAQAoDBUdKDNrBgBAYajoQJkp6AAAFIqKDpSVJQMAoFBUdKCsFHQAAApHRQfKyIoBAFA4KjpQRgo6AACFpKIDZWPBAAAoJBUdKBsFHQCAwlLRgTKxXgAAFJaKDpSJgg4AQKGp6EBZWC4AAApNRQfKQkEHAKDwVHSgDKwWAACFp6IDZaCgAwBQCio6UHQWCwCAUlDRgaJT0AEAKA0VHSgyawUAQGmo6ECRKegAAJSKig4UlaUCAKBUVHSgqBR0AABKR0UHishKAQBQOio6UEQKOgAApaSiA0VjoQAAKCUVHSgaBR0AgNJS0YEisU4AAJSWig4UiYIOAECpqehAUVgmAABKTUUHikJBBwCg9FR0oAisEgAApaeiA0WgoAMAUAkqOpA6iwQAQCWo6EDqFHQAACpDRQdSZo0AAKgMFR1ImYIOAEClqOhAqiwRAACVoqIDqVLQAQCoHBUdSJEVAgCgclR0IEUKOgAAlaSiA6mxQAAAVJKKDqRGQQcAoLJUdCAl1gcAgMpS0YGUKOgAAFSaig6kwvIAAFBpKjqQCgUdAIDKU9GBFFgdAAAqT0UHUqCgAwBAqOhA/iwOAACEig7kT0EHAID/UtGBPFkbAAD4LxUdyJOCDgAAl1DRgbxYGgAAuISKDuRFQQcAgF9R0YE8WBkAAPgVFR3Ig4IOAACjUNGBrFkYAAAYhYoOZE1BBwCAMajoQJasCwAAjEFFB7KkoAMAwBWo6EBWLAsAAFyBig5kRUEHAIBxqOhAFqwKAACMQ0UHsqCgAwBAHVR0oNUsCgAA1EFFB1pNQQcAgDqp6EArWRMAAKiTig60koIOAAAToKIDrWJJAABgAlR0oFUUdAAAmCAVHWgFKwIAABOkogOtoKADAEADVHSg2SwIAAA0QEUHmk1BBwCABqnoQDNZDwAAaJCKDjSTgg4AAJOgogPNYjkAAGASVHSgWRR0AACYJBUdaAarAQAAk6SiA82goAMAQBOo6MBkWQwAAGgCFR2YLAUdAACaREUHJsNaAABAk6jowGQo6AAA0EQqOtAoSwEAAE2kogONUtABAKDJVHSgEVYCAACaTEUHGqGgAwBAC6jowERZCAAAaAEVHZgoBR0AAFpERQcmwjoAAECLqOjARCjoAADQQio6UC/LAAAALaSiA/VS0AEAoMVUdKAeVgEAAFpMRQfqoaADAEAGVHRgPBYBAAAyoKID41HQAQAgIyo6cCXWAAAAMqKiA1eioAMAQIZUdGAslgAAADKkogNjUdABACBjKjowGisAAAAZU9GB0SjoAACQAxUd+DULAAAAOVDRgV9T0AEAICcqOnApr34AAMiJig5cSkEHAIAcqejARV75AACQIxUduEhBBwCAnKnoQISCDgAAuVPRgQgFHQAAkqCiA17xAACQABUdUNABACARKjpUm1c7AAAkQkWHalPQAQAgISo6VJdXOgAAJERFh+pS0AEAIDEqOlSTVzkAACRGRYdqUtABACBBKjpUj1c4AAAkSEWH6lHQAQAgUSo6VItXNwAAJEpFh2pR0AEAIGEqOlSHVzYAACRMRYfqUNABACBxKjpUg1c1AAAkTkWHalDQAQCgAFR0KD+vaAAAKAAVHcpPQQcAgIJQ0aHcvJoBAKAgVHQoNwUdAAAKREWH8nKgAwBAgQwNDcWqVati3759ERHx17/uic7O/4nBwYj58yOuuSbnBwQa5kAHAICCuVDRN0XE/0Vb2//G0NBVERExbVrEn/8csWVLxNKl+T4jMHEOdAAAKJh//nMourrOxeBgLSKmjvhcrRbR1hbx979HrF+fy+MBDXKgAwBAgXzzzYU6fvbsUAwNtY36NW1tEVOmROzfH7F8ecYPCDTMd5MAAIAC2bYt4j//iTGP84iIiwnuL3/J6KGAplDQAQCgIM6di+joiPj3v+v7+lotoq8v4ne/a+1zAc2hoAMAQEEcP17/cR4RMTAQcfRo654HaC4HOgAAlJj3y0JxONABAKAg5s6N+O1v6//6KVMibr65dc8DNJcDHQAACuKqqyI2brzwd8vHU6tF/PGPEb//feufC2gO3yQOAAAK5F//ili2LOL8+Su/fX3KlIh9+yJWrMju2YDJUdABAKBAbr454h//iPjNb0Yv6bXaheP8b39znEPRONABAKBg7rwz4sCBiA0bIqZN++XjtVrEn/4UsX9/xH335fd8QGO8xR0AAArsxx8jvv76wtvd//CHiFmz8n4ioFEOdAAAAEiAt7gDAABAAhzoAAAAkAAHOgAAACTAgQ4AAAAJcKADAABAAhzoAAAAkAAHOgAAACTAgQ4AAAAJcKADAABAAhzoAAAAkAAHOgAAACTAgQ4AAAAJcKADAABAAhzoAAAAkAAHOgAAACTAgQ4AAAAJcKADAABAAhzoAAAAkAAHOgAAACTAgQ4AAAAJcKADAABAAhzoAAAAkID/BxTd7UDzLniRAAAAAElFTkSuQmCC"
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/plain": "<Axes: >"
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pp.plotting.simple_plot(net, ext_grid_size=5)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:28.728039Z",
     "start_time": "2024-03-30T23:39:27.058805Z"
    }
   },
   "outputs": [],
   "source": [
    "%%capture\n",
    "net = single_3w_trafo_grid(\"Yyy\")\n",
    "sc.calc_sc(net, fault=\"1ph\", case=\"max\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:28.740814Z",
     "start_time": "2024-03-30T23:39:28.730053Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": "   ikss_ka  rk0_ohm   xk0_ohm   rk_ohm    xk_ohm\n0  1.51934 15.80517 158.05171 15.80517 158.05171\n1  0.00000      inf       inf  1.52174  22.24258\n2  0.00000      inf       inf  0.11730   2.66335",
      "text/html": "<div>\n<style scoped>\n    .dataframe tbody tr th:only-of-type {\n        vertical-align: middle;\n    }\n\n    .dataframe tbody tr th {\n        vertical-align: top;\n    }\n\n    .dataframe thead th {\n        text-align: right;\n    }\n</style>\n<table border=\"1\" class=\"dataframe\">\n  <thead>\n    <tr style=\"text-align: right;\">\n      <th></th>\n      <th>ikss_ka</th>\n      <th>rk0_ohm</th>\n      <th>xk0_ohm</th>\n      <th>rk_ohm</th>\n      <th>xk_ohm</th>\n    </tr>\n  </thead>\n  <tbody>\n    <tr>\n      <th>0</th>\n      <td>1.51934</td>\n      <td>15.80517</td>\n      <td>158.05171</td>\n      <td>15.80517</td>\n      <td>158.05171</td>\n    </tr>\n    <tr>\n      <th>1</th>\n      <td>0.00000</td>\n      <td>inf</td>\n      <td>inf</td>\n      <td>1.52174</td>\n      <td>22.24258</td>\n    </tr>\n    <tr>\n      <th>2</th>\n      <td>0.00000</td>\n      <td>inf</td>\n      <td>inf</td>\n      <td>0.11730</td>\n      <td>2.66335</td>\n    </tr>\n  </tbody>\n</table>\n</div>"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(net.res_bus_sc)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:28.821936Z",
     "start_time": "2024-03-30T23:39:28.742606Z"
    }
   },
   "outputs": [],
   "source": [
    "%%capture\n",
    "net = single_3w_trafo_grid(\"Ydd\")\n",
    "sc.calc_sc(net, fault=\"1ph\", case=\"max\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:28.832169Z",
     "start_time": "2024-03-30T23:39:28.823073Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": "   ikss_ka  rk0_ohm   xk0_ohm   rk_ohm    xk_ohm\n0  1.51934 15.80517 158.05171 15.80517 158.05171\n1  0.00000      inf       inf  1.52174  22.24258\n2  0.00000      inf       inf  0.11730   2.66335",
      "text/html": "<div>\n<style scoped>\n    .dataframe tbody tr th:only-of-type {\n        vertical-align: middle;\n    }\n\n    .dataframe tbody tr th {\n        vertical-align: top;\n    }\n\n    .dataframe thead th {\n        text-align: right;\n    }\n</style>\n<table border=\"1\" class=\"dataframe\">\n  <thead>\n    <tr style=\"text-align: right;\">\n      <th></th>\n      <th>ikss_ka</th>\n      <th>rk0_ohm</th>\n      <th>xk0_ohm</th>\n      <th>rk_ohm</th>\n      <th>xk_ohm</th>\n    </tr>\n  </thead>\n  <tbody>\n    <tr>\n      <th>0</th>\n      <td>1.51934</td>\n      <td>15.80517</td>\n      <td>158.05171</td>\n      <td>15.80517</td>\n      <td>158.05171</td>\n    </tr>\n    <tr>\n      <th>1</th>\n      <td>0.00000</td>\n      <td>inf</td>\n      <td>inf</td>\n      <td>1.52174</td>\n      <td>22.24258</td>\n    </tr>\n    <tr>\n      <th>2</th>\n      <td>0.00000</td>\n      <td>inf</td>\n      <td>inf</td>\n      <td>0.11730</td>\n      <td>2.66335</td>\n    </tr>\n  </tbody>\n</table>\n</div>"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(net.res_bus_sc)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:28.919184Z",
     "start_time": "2024-03-30T23:39:28.833158Z"
    }
   },
   "outputs": [],
   "source": [
    "%%capture\n",
    "net = single_3w_trafo_grid(\"Ddd\")\n",
    "sc.calc_sc(net, fault=\"1ph\", case=\"max\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:28.933415Z",
     "start_time": "2024-03-30T23:39:28.923169Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": "   ikss_ka  rk0_ohm   xk0_ohm   rk_ohm    xk_ohm\n0  1.51934 15.80517 158.05171 15.80517 158.05171\n1  0.00000      inf       inf  1.52174  22.24258\n2  0.00000      inf       inf  0.11730   2.66335",
      "text/html": "<div>\n<style scoped>\n    .dataframe tbody tr th:only-of-type {\n        vertical-align: middle;\n    }\n\n    .dataframe tbody tr th {\n        vertical-align: top;\n    }\n\n    .dataframe thead th {\n        text-align: right;\n    }\n</style>\n<table border=\"1\" class=\"dataframe\">\n  <thead>\n    <tr style=\"text-align: right;\">\n      <th></th>\n      <th>ikss_ka</th>\n      <th>rk0_ohm</th>\n      <th>xk0_ohm</th>\n      <th>rk_ohm</th>\n      <th>xk_ohm</th>\n    </tr>\n  </thead>\n  <tbody>\n    <tr>\n      <th>0</th>\n      <td>1.51934</td>\n      <td>15.80517</td>\n      <td>158.05171</td>\n      <td>15.80517</td>\n      <td>158.05171</td>\n    </tr>\n    <tr>\n      <th>1</th>\n      <td>0.00000</td>\n      <td>inf</td>\n      <td>inf</td>\n      <td>1.52174</td>\n      <td>22.24258</td>\n    </tr>\n    <tr>\n      <th>2</th>\n      <td>0.00000</td>\n      <td>inf</td>\n      <td>inf</td>\n      <td>0.11730</td>\n      <td>2.66335</td>\n    </tr>\n  </tbody>\n</table>\n</div>"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(net.res_bus_sc)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:29.026189Z",
     "start_time": "2024-03-30T23:39:28.935408Z"
    }
   },
   "outputs": [],
   "source": [
    "%%capture\n",
    "net = single_3w_trafo_grid(\"Yynd\")\n",
    "sc.calc_sc(net, fault=\"1ph\", case=\"max\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:29.037081Z",
     "start_time": "2024-03-30T23:39:29.028198Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": "   ikss_ka  rk0_ohm   xk0_ohm   rk_ohm    xk_ohm\n0  1.51934 15.80517 158.05171 15.80517 158.05171\n1  3.33940  0.10727  18.19496  1.52174  22.24258\n2  0.00000      inf       inf  0.11730   2.66335",
      "text/html": "<div>\n<style scoped>\n    .dataframe tbody tr th:only-of-type {\n        vertical-align: middle;\n    }\n\n    .dataframe tbody tr th {\n        vertical-align: top;\n    }\n\n    .dataframe thead th {\n        text-align: right;\n    }\n</style>\n<table border=\"1\" class=\"dataframe\">\n  <thead>\n    <tr style=\"text-align: right;\">\n      <th></th>\n      <th>ikss_ka</th>\n      <th>rk0_ohm</th>\n      <th>xk0_ohm</th>\n      <th>rk_ohm</th>\n      <th>xk_ohm</th>\n    </tr>\n  </thead>\n  <tbody>\n    <tr>\n      <th>0</th>\n      <td>1.51934</td>\n      <td>15.80517</td>\n      <td>158.05171</td>\n      <td>15.80517</td>\n      <td>158.05171</td>\n    </tr>\n    <tr>\n      <th>1</th>\n      <td>3.33940</td>\n      <td>0.10727</td>\n      <td>18.19496</td>\n      <td>1.52174</td>\n      <td>22.24258</td>\n    </tr>\n    <tr>\n      <th>2</th>\n      <td>0.00000</td>\n      <td>inf</td>\n      <td>inf</td>\n      <td>0.11730</td>\n      <td>2.66335</td>\n    </tr>\n  </tbody>\n</table>\n</div>"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(net.res_bus_sc)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:29.127456Z",
     "start_time": "2024-03-30T23:39:29.039080Z"
    }
   },
   "outputs": [],
   "source": [
    "%%capture\n",
    "net = single_3w_trafo_grid(\"Ydyn\")  # Extra vector group\n",
    "sc.calc_sc(net, fault=\"1ph\", case=\"max\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:29.138548Z",
     "start_time": "2024-03-30T23:39:29.129469Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": "   ikss_ka  rk0_ohm   xk0_ohm   rk_ohm    xk_ohm\n0  1.51934 15.80517 158.05171 15.80517 158.05171\n1  0.00000      inf       inf  1.52174  22.24258\n2  8.83645  0.00670   1.13718  0.11730   2.66335",
      "text/html": "<div>\n<style scoped>\n    .dataframe tbody tr th:only-of-type {\n        vertical-align: middle;\n    }\n\n    .dataframe tbody tr th {\n        vertical-align: top;\n    }\n\n    .dataframe thead th {\n        text-align: right;\n    }\n</style>\n<table border=\"1\" class=\"dataframe\">\n  <thead>\n    <tr style=\"text-align: right;\">\n      <th></th>\n      <th>ikss_ka</th>\n      <th>rk0_ohm</th>\n      <th>xk0_ohm</th>\n      <th>rk_ohm</th>\n      <th>xk_ohm</th>\n    </tr>\n  </thead>\n  <tbody>\n    <tr>\n      <th>0</th>\n      <td>1.51934</td>\n      <td>15.80517</td>\n      <td>158.05171</td>\n      <td>15.80517</td>\n      <td>158.05171</td>\n    </tr>\n    <tr>\n      <th>1</th>\n      <td>0.00000</td>\n      <td>inf</td>\n      <td>inf</td>\n      <td>1.52174</td>\n      <td>22.24258</td>\n    </tr>\n    <tr>\n      <th>2</th>\n      <td>8.83645</td>\n      <td>0.00670</td>\n      <td>1.13718</td>\n      <td>0.11730</td>\n      <td>2.66335</td>\n    </tr>\n  </tbody>\n</table>\n</div>"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(net.res_bus_sc)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:29.242692Z",
     "start_time": "2024-03-30T23:39:29.140543Z"
    }
   },
   "outputs": [],
   "source": [
    "%%capture\n",
    "net = single_3w_trafo_grid(\"YNynd\")\n",
    "sc.calc_sc(net, fault=\"1ph\", case=\"max\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:29.257054Z",
     "start_time": "2024-03-30T23:39:29.244810Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": "   ikss_ka  rk0_ohm  xk0_ohm   rk_ohm    xk_ohm\n0  1.78326  5.12452 88.22838 15.80517 158.05171\n1  3.49934  0.21021 15.31723  1.52174  22.24258\n2  0.00000      inf      inf  0.11730   2.66335",
      "text/html": "<div>\n<style scoped>\n    .dataframe tbody tr th:only-of-type {\n        vertical-align: middle;\n    }\n\n    .dataframe tbody tr th {\n        vertical-align: top;\n    }\n\n    .dataframe thead th {\n        text-align: right;\n    }\n</style>\n<table border=\"1\" class=\"dataframe\">\n  <thead>\n    <tr style=\"text-align: right;\">\n      <th></th>\n      <th>ikss_ka</th>\n      <th>rk0_ohm</th>\n      <th>xk0_ohm</th>\n      <th>rk_ohm</th>\n      <th>xk_ohm</th>\n    </tr>\n  </thead>\n  <tbody>\n    <tr>\n      <th>0</th>\n      <td>1.78326</td>\n      <td>5.12452</td>\n      <td>88.22838</td>\n      <td>15.80517</td>\n      <td>158.05171</td>\n    </tr>\n    <tr>\n      <th>1</th>\n      <td>3.49934</td>\n      <td>0.21021</td>\n      <td>15.31723</td>\n      <td>1.52174</td>\n      <td>22.24258</td>\n    </tr>\n    <tr>\n      <th>2</th>\n      <td>0.00000</td>\n      <td>inf</td>\n      <td>inf</td>\n      <td>0.11730</td>\n      <td>2.66335</td>\n    </tr>\n  </tbody>\n</table>\n</div>"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(net.res_bus_sc)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:29.368422Z",
     "start_time": "2024-03-30T23:39:29.258160Z"
    }
   },
   "outputs": [],
   "source": [
    "%%capture\n",
    "net = single_3w_trafo_grid(\"YNdyn\")  # Extra vector group\n",
    "sc.calc_sc(net, fault=\"1ph\", case=\"max\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:29.380999Z",
     "start_time": "2024-03-30T23:39:29.370437Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": "   ikss_ka  rk0_ohm  xk0_ohm   rk_ohm    xk_ohm\n0  1.79376  4.86639 85.86380 15.80517 158.05171\n1  0.00000      inf      inf  1.52174  22.24258\n2  9.04239  0.01219  0.98955  0.11730   2.66335",
      "text/html": "<div>\n<style scoped>\n    .dataframe tbody tr th:only-of-type {\n        vertical-align: middle;\n    }\n\n    .dataframe tbody tr th {\n        vertical-align: top;\n    }\n\n    .dataframe thead th {\n        text-align: right;\n    }\n</style>\n<table border=\"1\" class=\"dataframe\">\n  <thead>\n    <tr style=\"text-align: right;\">\n      <th></th>\n      <th>ikss_ka</th>\n      <th>rk0_ohm</th>\n      <th>xk0_ohm</th>\n      <th>rk_ohm</th>\n      <th>xk_ohm</th>\n    </tr>\n  </thead>\n  <tbody>\n    <tr>\n      <th>0</th>\n      <td>1.79376</td>\n      <td>4.86639</td>\n      <td>85.86380</td>\n      <td>15.80517</td>\n      <td>158.05171</td>\n    </tr>\n    <tr>\n      <th>1</th>\n      <td>0.00000</td>\n      <td>inf</td>\n      <td>inf</td>\n      <td>1.52174</td>\n      <td>22.24258</td>\n    </tr>\n    <tr>\n      <th>2</th>\n      <td>9.04239</td>\n      <td>0.01219</td>\n      <td>0.98955</td>\n      <td>0.11730</td>\n      <td>2.66335</td>\n    </tr>\n  </tbody>\n</table>\n</div>"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(net.res_bus_sc)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:29.483326Z",
     "start_time": "2024-03-30T23:39:29.383016Z"
    }
   },
   "outputs": [],
   "source": [
    "%%capture\n",
    "net = single_3w_trafo_grid(\"YNyd\")\n",
    "sc.calc_sc(net, fault=\"1ph\", case=\"max\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:29.495096Z",
     "start_time": "2024-03-30T23:39:29.484338Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": "   ikss_ka  rk0_ohm  xk0_ohm   rk_ohm    xk_ohm\n0  1.78326  5.12452 88.22838 15.80517 158.05171\n1  0.00000      inf      inf  1.52174  22.24258\n2  0.00000      inf      inf  0.11730   2.66335",
      "text/html": "<div>\n<style scoped>\n    .dataframe tbody tr th:only-of-type {\n        vertical-align: middle;\n    }\n\n    .dataframe tbody tr th {\n        vertical-align: top;\n    }\n\n    .dataframe thead th {\n        text-align: right;\n    }\n</style>\n<table border=\"1\" class=\"dataframe\">\n  <thead>\n    <tr style=\"text-align: right;\">\n      <th></th>\n      <th>ikss_ka</th>\n      <th>rk0_ohm</th>\n      <th>xk0_ohm</th>\n      <th>rk_ohm</th>\n      <th>xk_ohm</th>\n    </tr>\n  </thead>\n  <tbody>\n    <tr>\n      <th>0</th>\n      <td>1.78326</td>\n      <td>5.12452</td>\n      <td>88.22838</td>\n      <td>15.80517</td>\n      <td>158.05171</td>\n    </tr>\n    <tr>\n      <th>1</th>\n      <td>0.00000</td>\n      <td>inf</td>\n      <td>inf</td>\n      <td>1.52174</td>\n      <td>22.24258</td>\n    </tr>\n    <tr>\n      <th>2</th>\n      <td>0.00000</td>\n      <td>inf</td>\n      <td>inf</td>\n      <td>0.11730</td>\n      <td>2.66335</td>\n    </tr>\n  </tbody>\n</table>\n</div>"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(net.res_bus_sc)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:29.590956Z",
     "start_time": "2024-03-30T23:39:29.496106Z"
    }
   },
   "outputs": [],
   "source": [
    "%%capture\n",
    "net = single_3w_trafo_grid(\"YNdy\")  # Extra vector group\n",
    "sc.calc_sc(net, fault=\"1ph\", case=\"max\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:29.602393Z",
     "start_time": "2024-03-30T23:39:29.592968Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": "   ikss_ka  rk0_ohm  xk0_ohm   rk_ohm    xk_ohm\n0  1.79376  4.86639 85.86380 15.80517 158.05171\n1  0.00000      inf      inf  1.52174  22.24258\n2  0.00000      inf      inf  0.11730   2.66335",
      "text/html": "<div>\n<style scoped>\n    .dataframe tbody tr th:only-of-type {\n        vertical-align: middle;\n    }\n\n    .dataframe tbody tr th {\n        vertical-align: top;\n    }\n\n    .dataframe thead th {\n        text-align: right;\n    }\n</style>\n<table border=\"1\" class=\"dataframe\">\n  <thead>\n    <tr style=\"text-align: right;\">\n      <th></th>\n      <th>ikss_ka</th>\n      <th>rk0_ohm</th>\n      <th>xk0_ohm</th>\n      <th>rk_ohm</th>\n      <th>xk_ohm</th>\n    </tr>\n  </thead>\n  <tbody>\n    <tr>\n      <th>0</th>\n      <td>1.79376</td>\n      <td>4.86639</td>\n      <td>85.86380</td>\n      <td>15.80517</td>\n      <td>158.05171</td>\n    </tr>\n    <tr>\n      <th>1</th>\n      <td>0.00000</td>\n      <td>inf</td>\n      <td>inf</td>\n      <td>1.52174</td>\n      <td>22.24258</td>\n    </tr>\n    <tr>\n      <th>2</th>\n      <td>0.00000</td>\n      <td>inf</td>\n      <td>inf</td>\n      <td>0.11730</td>\n      <td>2.66335</td>\n    </tr>\n  </tbody>\n</table>\n</div>"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(net.res_bus_sc)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:29.708673Z",
     "start_time": "2024-03-30T23:39:29.603391Z"
    }
   },
   "outputs": [],
   "source": [
    "%%capture\n",
    "net = single_3w_trafo_grid(\"YNdd\")\n",
    "sc.calc_sc(net, fault=\"1ph\", case=\"max\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:29.721244Z",
     "start_time": "2024-03-30T23:39:29.710690Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": "   ikss_ka  rk0_ohm  xk0_ohm   rk_ohm    xk_ohm\n0  1.84354  3.76841 75.01993 15.80517 158.05171\n1  0.00000      inf      inf  1.52174  22.24258\n2  0.00000      inf      inf  0.11730   2.66335",
      "text/html": "<div>\n<style scoped>\n    .dataframe tbody tr th:only-of-type {\n        vertical-align: middle;\n    }\n\n    .dataframe tbody tr th {\n        vertical-align: top;\n    }\n\n    .dataframe thead th {\n        text-align: right;\n    }\n</style>\n<table border=\"1\" class=\"dataframe\">\n  <thead>\n    <tr style=\"text-align: right;\">\n      <th></th>\n      <th>ikss_ka</th>\n      <th>rk0_ohm</th>\n      <th>xk0_ohm</th>\n      <th>rk_ohm</th>\n      <th>xk_ohm</th>\n    </tr>\n  </thead>\n  <tbody>\n    <tr>\n      <th>0</th>\n      <td>1.84354</td>\n      <td>3.76841</td>\n      <td>75.01993</td>\n      <td>15.80517</td>\n      <td>158.05171</td>\n    </tr>\n    <tr>\n      <th>1</th>\n      <td>0.00000</td>\n      <td>inf</td>\n      <td>inf</td>\n      <td>1.52174</td>\n      <td>22.24258</td>\n    </tr>\n    <tr>\n      <th>2</th>\n      <td>0.00000</td>\n      <td>inf</td>\n      <td>inf</td>\n      <td>0.11730</td>\n      <td>2.66335</td>\n    </tr>\n  </tbody>\n</table>\n</div>"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(net.res_bus_sc)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:29.840707Z",
     "start_time": "2024-03-30T23:39:29.723258Z"
    }
   },
   "outputs": [],
   "source": [
    "%%capture\n",
    "net = single_3w_trafo_grid(\"YNyy\")  # Extra vector group\n",
    "sc.calc_sc(net, fault=\"1ph\", case=\"max\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:29.854031Z",
     "start_time": "2024-03-30T23:39:29.842704Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": "   ikss_ka  rk0_ohm   xk0_ohm   rk_ohm    xk_ohm\n0  1.51934 15.80517 158.05171 15.80517 158.05171\n1  0.00000      inf       inf  1.52174  22.24258\n2  0.00000      inf       inf  0.11730   2.66335",
      "text/html": "<div>\n<style scoped>\n    .dataframe tbody tr th:only-of-type {\n        vertical-align: middle;\n    }\n\n    .dataframe tbody tr th {\n        vertical-align: top;\n    }\n\n    .dataframe thead th {\n        text-align: right;\n    }\n</style>\n<table border=\"1\" class=\"dataframe\">\n  <thead>\n    <tr style=\"text-align: right;\">\n      <th></th>\n      <th>ikss_ka</th>\n      <th>rk0_ohm</th>\n      <th>xk0_ohm</th>\n      <th>rk_ohm</th>\n      <th>xk_ohm</th>\n    </tr>\n  </thead>\n  <tbody>\n    <tr>\n      <th>0</th>\n      <td>1.51934</td>\n      <td>15.80517</td>\n      <td>158.05171</td>\n      <td>15.80517</td>\n      <td>158.05171</td>\n    </tr>\n    <tr>\n      <th>1</th>\n      <td>0.00000</td>\n      <td>inf</td>\n      <td>inf</td>\n      <td>1.52174</td>\n      <td>22.24258</td>\n    </tr>\n    <tr>\n      <th>2</th>\n      <td>0.00000</td>\n      <td>inf</td>\n      <td>inf</td>\n      <td>0.11730</td>\n      <td>2.66335</td>\n    </tr>\n  </tbody>\n</table>\n</div>"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(net.res_bus_sc)  # warum dasselbe Ergebnis wie Yyy?"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Single-bus fault"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "all-bus fault vs single-bus fault vs LU factorization"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Ynynd"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:30.077108Z",
     "start_time": "2024-03-30T23:39:29.855926Z"
    }
   },
   "outputs": [],
   "source": [
    "%%capture\n",
    "net = single_3w_trafo_grid(\"Ynynd\")  \n",
    "sc.calc_sc(net, fault=\"1ph\", case=\"max\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:30.096432Z",
     "start_time": "2024-03-30T23:39:30.085182Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": "   ikss_ka  rk0_ohm  xk0_ohm   rk_ohm    xk_ohm\n0  1.78326  5.12452 88.22838 15.80517 158.05171\n1  3.49934  0.21021 15.31723  1.52174  22.24258\n2  0.00000      inf      inf  0.11730   2.66335",
      "text/html": "<div>\n<style scoped>\n    .dataframe tbody tr th:only-of-type {\n        vertical-align: middle;\n    }\n\n    .dataframe tbody tr th {\n        vertical-align: top;\n    }\n\n    .dataframe thead th {\n        text-align: right;\n    }\n</style>\n<table border=\"1\" class=\"dataframe\">\n  <thead>\n    <tr style=\"text-align: right;\">\n      <th></th>\n      <th>ikss_ka</th>\n      <th>rk0_ohm</th>\n      <th>xk0_ohm</th>\n      <th>rk_ohm</th>\n      <th>xk_ohm</th>\n    </tr>\n  </thead>\n  <tbody>\n    <tr>\n      <th>0</th>\n      <td>1.78326</td>\n      <td>5.12452</td>\n      <td>88.22838</td>\n      <td>15.80517</td>\n      <td>158.05171</td>\n    </tr>\n    <tr>\n      <th>1</th>\n      <td>3.49934</td>\n      <td>0.21021</td>\n      <td>15.31723</td>\n      <td>1.52174</td>\n      <td>22.24258</td>\n    </tr>\n    <tr>\n      <th>2</th>\n      <td>0.00000</td>\n      <td>inf</td>\n      <td>inf</td>\n      <td>0.11730</td>\n      <td>2.66335</td>\n    </tr>\n  </tbody>\n</table>\n</div>"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(net.res_bus_sc)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:30.204778Z",
     "start_time": "2024-03-30T23:39:30.098446Z"
    }
   },
   "outputs": [],
   "source": [
    "%%capture\n",
    "net = single_3w_trafo_grid(\"Ynynd\")  \n",
    "sc.calc_sc(net, fault=\"1ph\", case=\"max\", bus=0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:30.219950Z",
     "start_time": "2024-03-30T23:39:30.206744Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": "   ikss_ka  rk0_ohm  xk0_ohm   rk_ohm    xk_ohm\n0  1.78326  5.12452 88.22838 15.80517 158.05171",
      "text/html": "<div>\n<style scoped>\n    .dataframe tbody tr th:only-of-type {\n        vertical-align: middle;\n    }\n\n    .dataframe tbody tr th {\n        vertical-align: top;\n    }\n\n    .dataframe thead th {\n        text-align: right;\n    }\n</style>\n<table border=\"1\" class=\"dataframe\">\n  <thead>\n    <tr style=\"text-align: right;\">\n      <th></th>\n      <th>ikss_ka</th>\n      <th>rk0_ohm</th>\n      <th>xk0_ohm</th>\n      <th>rk_ohm</th>\n      <th>xk_ohm</th>\n    </tr>\n  </thead>\n  <tbody>\n    <tr>\n      <th>0</th>\n      <td>1.78326</td>\n      <td>5.12452</td>\n      <td>88.22838</td>\n      <td>15.80517</td>\n      <td>158.05171</td>\n    </tr>\n  </tbody>\n</table>\n</div>"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(net.res_bus_sc)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:30.340774Z",
     "start_time": "2024-03-30T23:39:30.222984Z"
    }
   },
   "outputs": [],
   "source": [
    "%%capture\n",
    "# now with LU factorization instead of Y inversion\n",
    "net = single_3w_trafo_grid(\"Ynynd\")  \n",
    "sc.calc_sc(net, fault=\"1ph\", case=\"max\", bus=0, inverse_y=False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:30.355680Z",
     "start_time": "2024-03-30T23:39:30.342811Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": "   ikss_ka  rk0_ohm  xk0_ohm   rk_ohm    xk_ohm\n0  1.78326  5.12452 88.22838 15.80517 158.05171",
      "text/html": "<div>\n<style scoped>\n    .dataframe tbody tr th:only-of-type {\n        vertical-align: middle;\n    }\n\n    .dataframe tbody tr th {\n        vertical-align: top;\n    }\n\n    .dataframe thead th {\n        text-align: right;\n    }\n</style>\n<table border=\"1\" class=\"dataframe\">\n  <thead>\n    <tr style=\"text-align: right;\">\n      <th></th>\n      <th>ikss_ka</th>\n      <th>rk0_ohm</th>\n      <th>xk0_ohm</th>\n      <th>rk_ohm</th>\n      <th>xk_ohm</th>\n    </tr>\n  </thead>\n  <tbody>\n    <tr>\n      <th>0</th>\n      <td>1.78326</td>\n      <td>5.12452</td>\n      <td>88.22838</td>\n      <td>15.80517</td>\n      <td>158.05171</td>\n    </tr>\n  </tbody>\n</table>\n</div>"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(net.res_bus_sc)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:30.464411Z",
     "start_time": "2024-03-30T23:39:30.357692Z"
    }
   },
   "outputs": [],
   "source": [
    "%%capture\n",
    "# now with LU factorization instead of Y inversion\n",
    "net = single_3w_trafo_grid(\"Ynynd\")  \n",
    "sc.calc_sc(net, fault=\"1ph\", case=\"max\", bus=1, inverse_y=False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:30.476434Z",
     "start_time": "2024-03-30T23:39:30.466426Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": "   ikss_ka  rk0_ohm  xk0_ohm  rk_ohm   xk_ohm\n1  3.49934  0.21021 15.31723 1.52174 22.24258",
      "text/html": "<div>\n<style scoped>\n    .dataframe tbody tr th:only-of-type {\n        vertical-align: middle;\n    }\n\n    .dataframe tbody tr th {\n        vertical-align: top;\n    }\n\n    .dataframe thead th {\n        text-align: right;\n    }\n</style>\n<table border=\"1\" class=\"dataframe\">\n  <thead>\n    <tr style=\"text-align: right;\">\n      <th></th>\n      <th>ikss_ka</th>\n      <th>rk0_ohm</th>\n      <th>xk0_ohm</th>\n      <th>rk_ohm</th>\n      <th>xk_ohm</th>\n    </tr>\n  </thead>\n  <tbody>\n    <tr>\n      <th>1</th>\n      <td>3.49934</td>\n      <td>0.21021</td>\n      <td>15.31723</td>\n      <td>1.52174</td>\n      <td>22.24258</td>\n    </tr>\n  </tbody>\n</table>\n</div>"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(net.res_bus_sc)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 35,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:30.583181Z",
     "start_time": "2024-03-30T23:39:30.478447Z"
    }
   },
   "outputs": [],
   "source": [
    "%%capture\n",
    "# now with LU factorization instead of Y inversion\n",
    "net = single_3w_trafo_grid(\"Ynynd\")  \n",
    "sc.calc_sc(net, fault=\"1ph\", case=\"max\", bus=2, inverse_y=False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:30.595832Z",
     "start_time": "2024-03-30T23:39:30.585193Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": "   ikss_ka  rk0_ohm  xk0_ohm  rk_ohm  xk_ohm\n2  0.00000      inf      inf 0.11730 2.66335",
      "text/html": "<div>\n<style scoped>\n    .dataframe tbody tr th:only-of-type {\n        vertical-align: middle;\n    }\n\n    .dataframe tbody tr th {\n        vertical-align: top;\n    }\n\n    .dataframe thead th {\n        text-align: right;\n    }\n</style>\n<table border=\"1\" class=\"dataframe\">\n  <thead>\n    <tr style=\"text-align: right;\">\n      <th></th>\n      <th>ikss_ka</th>\n      <th>rk0_ohm</th>\n      <th>xk0_ohm</th>\n      <th>rk_ohm</th>\n      <th>xk_ohm</th>\n    </tr>\n  </thead>\n  <tbody>\n    <tr>\n      <th>2</th>\n      <td>0.00000</td>\n      <td>inf</td>\n      <td>inf</td>\n      <td>0.11730</td>\n      <td>2.66335</td>\n    </tr>\n  </tbody>\n</table>\n</div>"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(net.res_bus_sc)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Ydd"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 37,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:30.706810Z",
     "start_time": "2024-03-30T23:39:30.596844Z"
    }
   },
   "outputs": [],
   "source": [
    "%%capture\n",
    "net = single_3w_trafo_grid(\"Ydd\")  \n",
    "sc.calc_sc(net, fault=\"1ph\", case=\"max\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:30.722321Z",
     "start_time": "2024-03-30T23:39:30.708824Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": "   ikss_ka  rk0_ohm   xk0_ohm   rk_ohm    xk_ohm\n0  1.51934 15.80517 158.05171 15.80517 158.05171\n1  0.00000      inf       inf  1.52174  22.24258\n2  0.00000      inf       inf  0.11730   2.66335",
      "text/html": "<div>\n<style scoped>\n    .dataframe tbody tr th:only-of-type {\n        vertical-align: middle;\n    }\n\n    .dataframe tbody tr th {\n        vertical-align: top;\n    }\n\n    .dataframe thead th {\n        text-align: right;\n    }\n</style>\n<table border=\"1\" class=\"dataframe\">\n  <thead>\n    <tr style=\"text-align: right;\">\n      <th></th>\n      <th>ikss_ka</th>\n      <th>rk0_ohm</th>\n      <th>xk0_ohm</th>\n      <th>rk_ohm</th>\n      <th>xk_ohm</th>\n    </tr>\n  </thead>\n  <tbody>\n    <tr>\n      <th>0</th>\n      <td>1.51934</td>\n      <td>15.80517</td>\n      <td>158.05171</td>\n      <td>15.80517</td>\n      <td>158.05171</td>\n    </tr>\n    <tr>\n      <th>1</th>\n      <td>0.00000</td>\n      <td>inf</td>\n      <td>inf</td>\n      <td>1.52174</td>\n      <td>22.24258</td>\n    </tr>\n    <tr>\n      <th>2</th>\n      <td>0.00000</td>\n      <td>inf</td>\n      <td>inf</td>\n      <td>0.11730</td>\n      <td>2.66335</td>\n    </tr>\n  </tbody>\n</table>\n</div>"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(net.res_bus_sc)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 39,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:30.848396Z",
     "start_time": "2024-03-30T23:39:30.724283Z"
    }
   },
   "outputs": [],
   "source": [
    "%%capture\n",
    "net = single_3w_trafo_grid(\"Ydd\")  \n",
    "sc.calc_sc(net, fault=\"1ph\", case=\"max\", bus=0, inverse_y=False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 40,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:30.862671Z",
     "start_time": "2024-03-30T23:39:30.850391Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": "   ikss_ka  rk0_ohm   xk0_ohm   rk_ohm    xk_ohm\n0  1.51934 15.80517 158.05171 15.80517 158.05171",
      "text/html": "<div>\n<style scoped>\n    .dataframe tbody tr th:only-of-type {\n        vertical-align: middle;\n    }\n\n    .dataframe tbody tr th {\n        vertical-align: top;\n    }\n\n    .dataframe thead th {\n        text-align: right;\n    }\n</style>\n<table border=\"1\" class=\"dataframe\">\n  <thead>\n    <tr style=\"text-align: right;\">\n      <th></th>\n      <th>ikss_ka</th>\n      <th>rk0_ohm</th>\n      <th>xk0_ohm</th>\n      <th>rk_ohm</th>\n      <th>xk_ohm</th>\n    </tr>\n  </thead>\n  <tbody>\n    <tr>\n      <th>0</th>\n      <td>1.51934</td>\n      <td>15.80517</td>\n      <td>158.05171</td>\n      <td>15.80517</td>\n      <td>158.05171</td>\n    </tr>\n  </tbody>\n</table>\n</div>"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(net.res_bus_sc)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 41,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:30.972483Z",
     "start_time": "2024-03-30T23:39:30.864684Z"
    }
   },
   "outputs": [],
   "source": [
    "%%capture\n",
    "net = single_3w_trafo_grid(\"Ydd\")  \n",
    "sc.calc_sc(net, fault=\"1ph\", case=\"max\", bus=1, inverse_y=False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 42,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:30.985449Z",
     "start_time": "2024-03-30T23:39:30.974495Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": "   ikss_ka  rk0_ohm  xk0_ohm  rk_ohm   xk_ohm\n1  0.00000      inf      inf 1.52174 22.24258",
      "text/html": "<div>\n<style scoped>\n    .dataframe tbody tr th:only-of-type {\n        vertical-align: middle;\n    }\n\n    .dataframe tbody tr th {\n        vertical-align: top;\n    }\n\n    .dataframe thead th {\n        text-align: right;\n    }\n</style>\n<table border=\"1\" class=\"dataframe\">\n  <thead>\n    <tr style=\"text-align: right;\">\n      <th></th>\n      <th>ikss_ka</th>\n      <th>rk0_ohm</th>\n      <th>xk0_ohm</th>\n      <th>rk_ohm</th>\n      <th>xk_ohm</th>\n    </tr>\n  </thead>\n  <tbody>\n    <tr>\n      <th>1</th>\n      <td>0.00000</td>\n      <td>inf</td>\n      <td>inf</td>\n      <td>1.52174</td>\n      <td>22.24258</td>\n    </tr>\n  </tbody>\n</table>\n</div>"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(net.res_bus_sc)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 43,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:31.121024Z",
     "start_time": "2024-03-30T23:39:30.986463Z"
    }
   },
   "outputs": [],
   "source": [
    "%%capture\n",
    "net = single_3w_trafo_grid(\"Ydd\")  \n",
    "sc.calc_sc(net, fault=\"1ph\", case=\"max\", bus=2, inverse_y=False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 44,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:31.137500Z",
     "start_time": "2024-03-30T23:39:31.123064Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": "   ikss_ka  rk0_ohm  xk0_ohm  rk_ohm  xk_ohm\n2  0.00000      inf      inf 0.11730 2.66335",
      "text/html": "<div>\n<style scoped>\n    .dataframe tbody tr th:only-of-type {\n        vertical-align: middle;\n    }\n\n    .dataframe tbody tr th {\n        vertical-align: top;\n    }\n\n    .dataframe thead th {\n        text-align: right;\n    }\n</style>\n<table border=\"1\" class=\"dataframe\">\n  <thead>\n    <tr style=\"text-align: right;\">\n      <th></th>\n      <th>ikss_ka</th>\n      <th>rk0_ohm</th>\n      <th>xk0_ohm</th>\n      <th>rk_ohm</th>\n      <th>xk_ohm</th>\n    </tr>\n  </thead>\n  <tbody>\n    <tr>\n      <th>2</th>\n      <td>0.00000</td>\n      <td>inf</td>\n      <td>inf</td>\n      <td>0.11730</td>\n      <td>2.66335</td>\n    </tr>\n  </tbody>\n</table>\n</div>"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(net.res_bus_sc)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## YNyd"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 45,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:31.269076Z",
     "start_time": "2024-03-30T23:39:31.140340Z"
    }
   },
   "outputs": [],
   "source": [
    "%%capture\n",
    "net = single_3w_trafo_grid(\"YNyd\")  \n",
    "sc.calc_sc(net, fault=\"1ph\", case=\"max\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 46,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:31.282373Z",
     "start_time": "2024-03-30T23:39:31.271090Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": "   ikss_ka  rk0_ohm  xk0_ohm   rk_ohm    xk_ohm\n0  1.78326  5.12452 88.22838 15.80517 158.05171\n1  0.00000      inf      inf  1.52174  22.24258\n2  0.00000      inf      inf  0.11730   2.66335",
      "text/html": "<div>\n<style scoped>\n    .dataframe tbody tr th:only-of-type {\n        vertical-align: middle;\n    }\n\n    .dataframe tbody tr th {\n        vertical-align: top;\n    }\n\n    .dataframe thead th {\n        text-align: right;\n    }\n</style>\n<table border=\"1\" class=\"dataframe\">\n  <thead>\n    <tr style=\"text-align: right;\">\n      <th></th>\n      <th>ikss_ka</th>\n      <th>rk0_ohm</th>\n      <th>xk0_ohm</th>\n      <th>rk_ohm</th>\n      <th>xk_ohm</th>\n    </tr>\n  </thead>\n  <tbody>\n    <tr>\n      <th>0</th>\n      <td>1.78326</td>\n      <td>5.12452</td>\n      <td>88.22838</td>\n      <td>15.80517</td>\n      <td>158.05171</td>\n    </tr>\n    <tr>\n      <th>1</th>\n      <td>0.00000</td>\n      <td>inf</td>\n      <td>inf</td>\n      <td>1.52174</td>\n      <td>22.24258</td>\n    </tr>\n    <tr>\n      <th>2</th>\n      <td>0.00000</td>\n      <td>inf</td>\n      <td>inf</td>\n      <td>0.11730</td>\n      <td>2.66335</td>\n    </tr>\n  </tbody>\n</table>\n</div>"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(net.res_bus_sc)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 47,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:31.412232Z",
     "start_time": "2024-03-30T23:39:31.284383Z"
    }
   },
   "outputs": [],
   "source": [
    "%%capture\n",
    "net = single_3w_trafo_grid(\"YNyd\")  \n",
    "sc.calc_sc(net, fault=\"1ph\", case=\"max\", bus=0, inverse_y=False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 48,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:31.429355Z",
     "start_time": "2024-03-30T23:39:31.414245Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": "   ikss_ka  rk0_ohm  xk0_ohm   rk_ohm    xk_ohm\n0  1.78326  5.12452 88.22838 15.80517 158.05171",
      "text/html": "<div>\n<style scoped>\n    .dataframe tbody tr th:only-of-type {\n        vertical-align: middle;\n    }\n\n    .dataframe tbody tr th {\n        vertical-align: top;\n    }\n\n    .dataframe thead th {\n        text-align: right;\n    }\n</style>\n<table border=\"1\" class=\"dataframe\">\n  <thead>\n    <tr style=\"text-align: right;\">\n      <th></th>\n      <th>ikss_ka</th>\n      <th>rk0_ohm</th>\n      <th>xk0_ohm</th>\n      <th>rk_ohm</th>\n      <th>xk_ohm</th>\n    </tr>\n  </thead>\n  <tbody>\n    <tr>\n      <th>0</th>\n      <td>1.78326</td>\n      <td>5.12452</td>\n      <td>88.22838</td>\n      <td>15.80517</td>\n      <td>158.05171</td>\n    </tr>\n  </tbody>\n</table>\n</div>"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(net.res_bus_sc)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 49,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:31.541083Z",
     "start_time": "2024-03-30T23:39:31.431371Z"
    }
   },
   "outputs": [],
   "source": [
    "%%capture\n",
    "net = single_3w_trafo_grid(\"YNyd\")  \n",
    "sc.calc_sc(net, fault=\"1ph\", case=\"max\", bus=1, inverse_y=False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 50,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:31.554596Z",
     "start_time": "2024-03-30T23:39:31.542078Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": "   ikss_ka  rk0_ohm  xk0_ohm  rk_ohm   xk_ohm\n1  0.00000      inf      inf 1.52174 22.24258",
      "text/html": "<div>\n<style scoped>\n    .dataframe tbody tr th:only-of-type {\n        vertical-align: middle;\n    }\n\n    .dataframe tbody tr th {\n        vertical-align: top;\n    }\n\n    .dataframe thead th {\n        text-align: right;\n    }\n</style>\n<table border=\"1\" class=\"dataframe\">\n  <thead>\n    <tr style=\"text-align: right;\">\n      <th></th>\n      <th>ikss_ka</th>\n      <th>rk0_ohm</th>\n      <th>xk0_ohm</th>\n      <th>rk_ohm</th>\n      <th>xk_ohm</th>\n    </tr>\n  </thead>\n  <tbody>\n    <tr>\n      <th>1</th>\n      <td>0.00000</td>\n      <td>inf</td>\n      <td>inf</td>\n      <td>1.52174</td>\n      <td>22.24258</td>\n    </tr>\n  </tbody>\n</table>\n</div>"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(net.res_bus_sc)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 51,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:31.667552Z",
     "start_time": "2024-03-30T23:39:31.556612Z"
    }
   },
   "outputs": [],
   "source": [
    "%%capture\n",
    "net = single_3w_trafo_grid(\"YNyd\")  \n",
    "sc.calc_sc(net, fault=\"1ph\", case=\"max\", bus=2, inverse_y=False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 52,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:31.680527Z",
     "start_time": "2024-03-30T23:39:31.669953Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": "   ikss_ka  rk0_ohm  xk0_ohm  rk_ohm  xk_ohm\n2  0.00000      inf      inf 0.11730 2.66335",
      "text/html": "<div>\n<style scoped>\n    .dataframe tbody tr th:only-of-type {\n        vertical-align: middle;\n    }\n\n    .dataframe tbody tr th {\n        vertical-align: top;\n    }\n\n    .dataframe thead th {\n        text-align: right;\n    }\n</style>\n<table border=\"1\" class=\"dataframe\">\n  <thead>\n    <tr style=\"text-align: right;\">\n      <th></th>\n      <th>ikss_ka</th>\n      <th>rk0_ohm</th>\n      <th>xk0_ohm</th>\n      <th>rk_ohm</th>\n      <th>xk_ohm</th>\n    </tr>\n  </thead>\n  <tbody>\n    <tr>\n      <th>2</th>\n      <td>0.00000</td>\n      <td>inf</td>\n      <td>inf</td>\n      <td>0.11730</td>\n      <td>2.66335</td>\n    </tr>\n  </tbody>\n</table>\n</div>"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(net.res_bus_sc)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## YNdd"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 53,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:31.792388Z",
     "start_time": "2024-03-30T23:39:31.682538Z"
    }
   },
   "outputs": [],
   "source": [
    "%%capture\n",
    "net = single_3w_trafo_grid(\"YNdd\")  \n",
    "sc.calc_sc(net, fault=\"1ph\", case=\"max\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 54,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:31.807739Z",
     "start_time": "2024-03-30T23:39:31.794399Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": "   ikss_ka  rk0_ohm  xk0_ohm   rk_ohm    xk_ohm\n0  1.84354  3.76841 75.01993 15.80517 158.05171\n1  0.00000      inf      inf  1.52174  22.24258\n2  0.00000      inf      inf  0.11730   2.66335",
      "text/html": "<div>\n<style scoped>\n    .dataframe tbody tr th:only-of-type {\n        vertical-align: middle;\n    }\n\n    .dataframe tbody tr th {\n        vertical-align: top;\n    }\n\n    .dataframe thead th {\n        text-align: right;\n    }\n</style>\n<table border=\"1\" class=\"dataframe\">\n  <thead>\n    <tr style=\"text-align: right;\">\n      <th></th>\n      <th>ikss_ka</th>\n      <th>rk0_ohm</th>\n      <th>xk0_ohm</th>\n      <th>rk_ohm</th>\n      <th>xk_ohm</th>\n    </tr>\n  </thead>\n  <tbody>\n    <tr>\n      <th>0</th>\n      <td>1.84354</td>\n      <td>3.76841</td>\n      <td>75.01993</td>\n      <td>15.80517</td>\n      <td>158.05171</td>\n    </tr>\n    <tr>\n      <th>1</th>\n      <td>0.00000</td>\n      <td>inf</td>\n      <td>inf</td>\n      <td>1.52174</td>\n      <td>22.24258</td>\n    </tr>\n    <tr>\n      <th>2</th>\n      <td>0.00000</td>\n      <td>inf</td>\n      <td>inf</td>\n      <td>0.11730</td>\n      <td>2.66335</td>\n    </tr>\n  </tbody>\n</table>\n</div>"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(net.res_bus_sc)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 55,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:31.949699Z",
     "start_time": "2024-03-30T23:39:31.809752Z"
    }
   },
   "outputs": [],
   "source": [
    "%%capture\n",
    "net = single_3w_trafo_grid(\"YNdd\")  \n",
    "sc.calc_sc(net, fault=\"1ph\", case=\"max\", bus=0, inverse_y=False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 56,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:31.986344Z",
     "start_time": "2024-03-30T23:39:31.951695Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": "   ikss_ka  rk0_ohm  xk0_ohm   rk_ohm    xk_ohm\n0  1.84354  3.76841 75.01993 15.80517 158.05171",
      "text/html": "<div>\n<style scoped>\n    .dataframe tbody tr th:only-of-type {\n        vertical-align: middle;\n    }\n\n    .dataframe tbody tr th {\n        vertical-align: top;\n    }\n\n    .dataframe thead th {\n        text-align: right;\n    }\n</style>\n<table border=\"1\" class=\"dataframe\">\n  <thead>\n    <tr style=\"text-align: right;\">\n      <th></th>\n      <th>ikss_ka</th>\n      <th>rk0_ohm</th>\n      <th>xk0_ohm</th>\n      <th>rk_ohm</th>\n      <th>xk_ohm</th>\n    </tr>\n  </thead>\n  <tbody>\n    <tr>\n      <th>0</th>\n      <td>1.84354</td>\n      <td>3.76841</td>\n      <td>75.01993</td>\n      <td>15.80517</td>\n      <td>158.05171</td>\n    </tr>\n  </tbody>\n</table>\n</div>"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(net.res_bus_sc)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 57,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:32.119098Z",
     "start_time": "2024-03-30T23:39:31.988392Z"
    }
   },
   "outputs": [],
   "source": [
    "%%capture\n",
    "net = single_3w_trafo_grid(\"YNdd\")  \n",
    "sc.calc_sc(net, fault=\"1ph\", case=\"max\", bus=1, inverse_y=False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 58,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:32.130783Z",
     "start_time": "2024-03-30T23:39:32.121109Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": "   ikss_ka  rk0_ohm  xk0_ohm  rk_ohm   xk_ohm\n1  0.00000      inf      inf 1.52174 22.24258",
      "text/html": "<div>\n<style scoped>\n    .dataframe tbody tr th:only-of-type {\n        vertical-align: middle;\n    }\n\n    .dataframe tbody tr th {\n        vertical-align: top;\n    }\n\n    .dataframe thead th {\n        text-align: right;\n    }\n</style>\n<table border=\"1\" class=\"dataframe\">\n  <thead>\n    <tr style=\"text-align: right;\">\n      <th></th>\n      <th>ikss_ka</th>\n      <th>rk0_ohm</th>\n      <th>xk0_ohm</th>\n      <th>rk_ohm</th>\n      <th>xk_ohm</th>\n    </tr>\n  </thead>\n  <tbody>\n    <tr>\n      <th>1</th>\n      <td>0.00000</td>\n      <td>inf</td>\n      <td>inf</td>\n      <td>1.52174</td>\n      <td>22.24258</td>\n    </tr>\n  </tbody>\n</table>\n</div>"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(net.res_bus_sc)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 59,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:32.238871Z",
     "start_time": "2024-03-30T23:39:32.132794Z"
    }
   },
   "outputs": [],
   "source": [
    "%%capture\n",
    "net = single_3w_trafo_grid(\"YNdd\")  \n",
    "sc.calc_sc(net, fault=\"1ph\", case=\"max\", bus=2, inverse_y=False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 60,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2024-03-30T23:39:32.250262Z",
     "start_time": "2024-03-30T23:39:32.239867Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": "   ikss_ka  rk0_ohm  xk0_ohm  rk_ohm  xk_ohm\n2  0.00000      inf      inf 0.11730 2.66335",
      "text/html": "<div>\n<style scoped>\n    .dataframe tbody tr th:only-of-type {\n        vertical-align: middle;\n    }\n\n    .dataframe tbody tr th {\n        vertical-align: top;\n    }\n\n    .dataframe thead th {\n        text-align: right;\n    }\n</style>\n<table border=\"1\" class=\"dataframe\">\n  <thead>\n    <tr style=\"text-align: right;\">\n      <th></th>\n      <th>ikss_ka</th>\n      <th>rk0_ohm</th>\n      <th>xk0_ohm</th>\n      <th>rk_ohm</th>\n      <th>xk_ohm</th>\n    </tr>\n  </thead>\n  <tbody>\n    <tr>\n      <th>2</th>\n      <td>0.00000</td>\n      <td>inf</td>\n      <td>inf</td>\n      <td>0.11730</td>\n      <td>2.66335</td>\n    </tr>\n  </tbody>\n</table>\n</div>"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(net.res_bus_sc)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "For further vector groups, we included the check in the test, all the vector groups passed."
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.9.13"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}