diff --git a/Project.toml b/Project.toml
index 86e1560..262c349 100644
--- a/Project.toml
+++ b/Project.toml
@@ -1,7 +1,7 @@
name = "ElectronGas"
uuid = "8c98e03e-8755-48b1-a5ef-1e2e7812a7c9"
authors = ["Xiansheng Cai, Pengcheng Hou, Kun Chen"]
-version = "0.2.5"
+version = "0.2.6"
[deps]
CompositeGrids = "b5136c89-beeb-4521-9139-60d2cac8be56"
diff --git a/example/G0W0_2D.ipynb b/example/G0W0_2D.ipynb
index cd7a26e..f93ccfb 100644
--- a/example/G0W0_2D.ipynb
+++ b/example/G0W0_2D.ipynb
@@ -26,7 +26,6 @@
"metadata": {},
"outputs": [],
"source": [
- "using Revise\n",
"using ElectronGas\n",
"using GreenFunc\n",
"using LaTeXStrings\n",
@@ -37,15 +36,7 @@
"cell_type": "code",
"execution_count": 2,
"metadata": {},
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}()\n"
- ]
- }
- ],
+ "outputs": [],
"source": [
"beta = 1000.0\n",
"rs = 1.0\n",
@@ -55,17 +46,16 @@
},
{
"cell_type": "code",
- "execution_count": 3,
+ "execution_count": 4,
"metadata": {},
"outputs": [],
"source": [
- "# sigma = SelfEnergy.G0W0(para, minK = 1e-6, Nk=12, order=4);\n",
- "# Nk, order, minK, rtol = 11, 8, 1e-8, 1e-10\n",
"Nk, order, minK, rtol = 8, 8, 1e-7, 1e-10\n",
- "sigma = SelfEnergy.G0W0(para, 100*para.EF, rtol, Nk, 10 * para.kF, minK * para.kF, order, :rpa)\n",
- "sigma_wn = GreenFunc.toMatFreq(sigma);\n",
- "dlr = sigma_wn.dlrGrid\n",
- "kgrid = sigma_wn.spaceGrid\n",
+ "sigma_dyn, sigma_ins = SelfEnergy.G0W0(para, 100*para.EF, rtol, Nk, 10 * para.kF, minK * para.kF, order, :rpa)\n",
+ "sigma_ωn = sigma_dyn |> to_dlr |> to_imfreq;\n",
+ "wgrid, kgrid = sigma_ωn.mesh\n",
+ "dlr = wgrid.representation\n",
+ "ωn_grid = wgrid.representation.ωn;\n",
"kFidx = searchsortedfirst(kgrid.grid, para.kF);"
]
},
@@ -78,250 +68,324 @@
},
{
"cell_type": "code",
- "execution_count": 8,
+ "execution_count": 5,
"metadata": {},
"outputs": [
{
"data": {
+ "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlgAAAGQCAIAAAD9V4nPAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nOzdd1QUVxsG8HdmC7sLS0dEFAEbghVEUcHeS2JP7Mau+dSYxBqNmmJiojFqYqJGY2yxG3uwYEVj772CVKnCNrbMfH+sQYMoirCzsM/v5Hh2h8vsCxP22Xtn7h2G53kCAACwVazQBQAAAAgJQQgAADYNQQgAADYNQQgAADYNQQgAADYNQQgAADYNQQgAADYNQQgAADYNQQgAADYNQQgAADatZAehXq+fPHmy0FUAEZHJZBK6BCDCgbAaHMcJXQIQEb3OMqIlOwizsrKWL18udBVARKTRaIQuAYhwIKwGDoQ14DjOaDQW2KxkByEAAMBbQhACAIBNQxACAIBNQxACAIBNQxACAIBNQxACAIBNQxACAIBNKz1BmJGRcfHiRZ1OJ3QhAABQkpSSIIzcfyCgYas2ExYG1AtPTEwUuhwAACgxxEIX8LZyar074rhp2+Q5KcO2kZMXc3bT2O+XT576mbc94yknRujyAACKSVpa2sOHD4WuwnIYhqlTpw7LFn3/rcQHIfskIcSd2SMVk1FPRGTQXX4iGn7clKDm03PIS8GUt6fy9kw5Bfk4MOUU5G3PVLAnLwUjKSWdYQCwUV9++eXWrVvLlCkjdCEWcuvWrf3794eFhRX5nkt8EEpizgwPYGvOn9FtUE/etaJTTsrJ/TtdXcVEpOcoVccnaihBwydq6H42fz6VErXc/Sx6pOYdJOQlZ8rZk7+S8VJQOcXTfys7Mk5SoX8qAICCmEymiRMn/u9//xO6EAsJDw8vpjXlS3wQmjUMC3t46Z/k5OTy5cvndpylLJVTMOUUFJLfEGlGzrOATNDw51L5+9l8gpri1Lyeo3IKxt/xaVI+n5deCgbDrQAApUkpCUIisrOz8/Hxef32LnbkYscEuRC9EJNaIyVq+ftZz5IyOpnuZ3MJGj4j51nf0fxvbl5WdGBECEkAgJKm9ARhEZKLyV/J+CvpxYzUmShBwyeoKVH79N9N9ylBwyVqKEbFO0qfDbQ+P+Lqq2TsX+83nZOTk52d7e7uXuQ/FAAA5AtB+GZkopdmpIGjlP+ekryWQYlaLkFNMSpeKsr/lGQlR8b531OSa9Zv+vjzrxilR60Kbn9vWScSiSz80wEA2CAEYZGRvPyUJE+UrKV4NR+v5h+pKUHDH0/i49R8opZiVbyEpfL2THl7OvbZN5pPjpJUcWrLJwcOHGjbtq0gPwgAgE1BEFoCQ1RWTmXlTIh7PmcRM/UUr+ZjVfQPy5NYSkQqkXLgAXV7uamZF9PMi6nogHOPAFD04uPjr127lpGRERgYGBAQIJFIhK5IGAhC4TlLyVnKBLnQpJEf/LC4A5Wp5KW6t2HitDOZTHQy//k5zsDx4Z5sK2+msScT5IJQBIC3defOnZEjR0ZFReVu8fLymj179qBBgwq3w+joaJPJ1KRJk6Kpz7IQhFZkyidj+/R49/Hjx8HBwSKRqLoHDahCRHQ/mz+exEcn819fRCgCwNu6fv16REREeno6SeXkU5fkTvToUmJiwgcffBAXFzdt2rRC7PPYsWMGgwFBCEWgYsWKFStWzLPRX8n4K5k8oTj7IqdHKALAmxs6dGh6ejoFNKfhq0lZhoiI5yhqMW34ZObMme+8806tWrXedJ+TJ0/Os0Wj0UilUrFYTEQGg4Hnean0P4uVZGdnK5XKPN+l0+nEYrH5uywG64yVMP5KZkAVdkm46OH74uOdxZ18mHOpfMdIU7l1hl4HTUtvctcyeKFrBADrde3atZMnT5LckUatf5qCRMSw1PJ/1GiAyWRasWJFIXb79ddfz5w5k4gOHz4cERExaNCggICAMmXK7Nq16+uvv65ataqHh8eCBQvMjefPn+/j4xMSElK2bNlFixaZN6rV6p49e/r5+QUGBn7++efVqlV7+x/2NaFHWIK9rKeYY+IjyqKnCGBzbj/hY1UFtDkYdZ6IqGoEKVzyfq1uF4r+I+qfCwfiC/g8LWapqdd/ltnSaDR6vZ6IDAZDdHT05MmTV65cuW7dun79+k2YMOHBgwfXr18PDQ0dMmSIg4NDp06dxowZIxaLExISQkND27VrV6VKle+//16lUsXExIhEot69ez958uQNf/rCQxCWEvmG4jeXOJ3xWSgGumB9OIDSbON9/kgS9+o28XfURESyvGOSRERyJRHFpKnmXC5gSU8xQ3XcxM4vWZa5UqVKHTt2JKIWLVo8efJkzJgxRBQYGOjs7BwTExMUFFSxYsWtW7dev35dp9PJ5fILFy5UqVLl77//njJlinn4dPTo0UePHi3gpy06CMJS6PlQTNDw0cn8gfhnodjYkwkvywS7IxQBSptpddlpBZ3w2ksVOywiiruaz9firhBRq1q+W9q/VTQ4OjqaH5hT7fmnOTk5RNSlSxcHB4euXbuKxeIDBw5kZ2cTUVZWVm5LJyentyngTSEIS7lyCqanH9PTj+i5UFxwjdMiFAFsUtOmTZVKZXb8VfpnHYX1efaF7BSK/IGIOnXqVKwFaDSayMhItVotk8k4jpsyZYp5e0BAwNmzZ5s3b05EZ86cKdYa8kAQ2pCXhaLGyDdBKALYBoVCMXPmzE8++YR+H0L3/qHgLiRzoAdnKXIupcfVqlWrX79+xVqAXC4vX778/PnzmzdvvmLFiszMTPP2SZMmde7cmeM4iUSydu1axoJvRQhCG4VQBLBZ48ePT05O/v777/nDv9LhX3O3165de+fOnYVbX6Zx48bmmwX6+vrmzsqXyWSTJk3KbTNy5EhPT0+GYXbv3j137twzZ8707NmzSZMm5gtEGzRocPDgwe3bt8tksq+//jq3p2gBDM+X4KvtU1NTq1evnpKSInQhpUduKB5P4tNy8obiN/MWrFy/uUb1gCU/fJPnFhn5TgkCy8OBsBIqlcrBwaFYX2LMmDHVqlUr9I15z5w589tvv509e1ar1VaqVKlbt259+/bNM9XPkrRa7ePHjytWrKjT6fr16+fr6zt37tznG4SHh8+ZM6dx48avv0+O40wmU4HRjh4h/MfzPcWH2fyRJP5w4tNzipVj9p7de1Y3aOfd6/s1oz7au2mN0MUCQOGFhoaGhoYKXcUzWq22Y8eOKpXKYDC0bNly+vTpFntpBCG8lK+S8VUyA6sQEcWo+LHTr+pqdiY7e67Ou6fnfZepp5ddPA0AJcLdu3evXLmSlZVVvXr1mjVryuVyAYtxdXW9ejW/a1mLnzAry2zZsqVdu3bt2rXbvHnzK5rt2rVrxIgRjx49slhh8DIVHZhpfVq7Hv+Jbh6W7ZjuWbeZ/wZDr4OmA/EleWwdwFZdu3atYcOGVapU6dat26BBgxo0aODj47N48WLLnyy7devW6tWrLfyieQjQIzxy5MiwYcP++OMPlmUHDBjg7u7erFmzF5slJSVNnTr1xo0bI0eOrFChgsXLhLxCQ0P/Wjz7j43bgztXGz54kIoTb7zPfXrKpDXR4Kpsz3IMzkwBlAiXLl1q2rTpkydPHKTiumWdHO3EVx9nx6SmfvjhhwkJCV999ZUli7lx48by5cv79+9vyRfNQ4Ag/Omnn8aOHdu5c2ci+uijjxYtWpRvEH744YezZs3q3bu3peuDl4sID48IDzc/diYaHsAOD2DPpfJLb3L19kpbeZv6V2E6VGBFuNgUwIoNHTr0yZMnzX3df2pX29FOTEQ80fqrcVOirn/zzTfdunULDg4u3J55nk9KSnJ3dzdfn/LkyRO5XP78BTgcx6WkpLi6ulrVvQ8FGBq9cOFCWFiY+XHDhg0vXLjwYpsNGzawLNu1a1fLlgaFEeLOLAkXXe2sb+XNzDrP+a43Tj5jilFhxBTAGl2+fPns2bOOdpKf/01BImKIetco/16QN8dxv//+eyF2+9VXXw0cODA0NDQ8PPzIkSNXrlypV69ew4YN/f39P/roI/OI68qVK/39/du1a1euXLmPP/64KH+qtyNAjzA5OdnF5elir66urklJSXkapKWlTZ8+/dChQwXuymAwZGZmPv/hpVevXoW+mBjeBpuj6lOe+pSnG0+YPx+yIdtEtVz4DyqZOntzYtzjxILUarUlZyLDy6jV6uJ+CYPBkHdLwgPj47hXf9eZ3ZFE1MDbRWmX9/2/tX+ZdVfjLpw4qr147NU7YSRSWfVQYp/9bet0ul27dp09e9bPz0+v19esWdPcudRqtW3atNm4ceN7773XqlWrfv36icVilUrVsGHDQ4cOmdeReU08z2u1WpWqoGXFn2Oenm+N0yeUSqVWqzU/VqvVL64pN27cuIkTJ3p7exe4K4lEolQqly1blrulYsWKxT13B/LF87z5Nx/qQKHeNDuMdsZyS2+KJl+g/lWY4QGssyFjyhff3n8UP25Iv04d2gtdb6mVeyBAcMV9IF58f9fdvqB/eOPV35V99xoRKSSiF79kLxERkTo9TXOxgAWvGZaVVqrByuyf39i5c2c/Pz8iunTpUnp6uqOj44EDBziO8/f3j4qKeu+99zw9PTdt2nT58uX09HSTyXThwoU3CkKGYeRy+Rv9Vs3zCAtsJkAQVqxY8d69e+bzgvfu3fPx8cnTIDo6eu3atcOGDTM/DQ4OXrRo0cv6eSKRKCQkpDjrhcKQiainH9vTj72Zya+8wzXcYTT8POxJjS5c7X4Xpow9VKF8zZo1ha4RoLRRNutWYJvAcvto5V/XU7Jf/NLVlGwiqlyvodugzwrx6m5ubuYHqampRHTgwAHzUy8vrxo1ahDR0KFDMzMzBwwY4OTklJCQ8EZ9u2IlQBD27t17+fLlAwYMIKLly5fnXg7z448/du3atWLFig8ePMhtLJPJTp48WbduXcvXCUUiwJn5NlQ0M1hU7osYrv77RJReq+fRE6cQhACCiIiIcHZ2vpOeue5qXJ8a5XO3J6tzlpx7QETvvPNO4facOyZfvXp1jUYzYcKE3Gg027lz59mzZ/39/YnIkvPlCyTA2Zthw4a5uLj4+/tXqlRJqVSOGDHCvH3atGm3b9+2fD1gATIR1anqJ/pnDcVflZzbOCOz3pxLnNYodFkAtkcul5snSEw+eO3jfVf23X98Ii598dkH7dadSFbn1KtX7+2v1TcvN9q5c+dt27bt37//u+++27ZtGxHVqFFj7ty5hw8f/uijj+7fv18EP0wREaBHKJfLd+/eHR8fz/N8+fLPPo9kZ2e/eJJfp9NZtjooLltXL/vsyzl3r0SNn/+Zd2jtLy9wvhsMH9cQjQliFVjgCMCCPvzww7S0tC+++GLTjYRNNxJytzdo0OCvv/4SiwvzBxkREWE0Pvts+/PPP2/atCkyMlKj0VSvXt28ltv69evnzp27bNmytm3btmjRwt7enoiqVasm7CRCwqLbUFQKsdbz1Qz+i/PckSQOcViEsOi2lbD+RbevXr26YsWKs2fP5uTk+Pv7d+nSpUePHiJRPhfRWAksug2lUA0XZmNL0dUM9rtLnP8Gw3jEIYAF1ahR44cffhC6CquAGV4gsBouzKpmogMdxNcyeP8NhjmXOA3OHQKABSEIwSq8Ig5zZ50CABQHBCFYkTxxOO1gQtWQcN+wNlWDG+EmJABQTBCEYHXMcbivvXjtL3PvNPrk8Ziou80/nzjrW6HrAoDSCUEIVqqWKxNsrybHMkTEO3rGpDwRuiIAKJ1wiR5Yryljhx/vO0wX2I6uRt7uM3/yGdOsYJGd9V7dDWBpFy5c2LRpk9BVWIh55bbigCAE61UvJOTqschz587VrTuOHD0/PGEK+cu4oomovgfurgBALVq0+PPPP20nCOvXr29enq3IYUI9FA3LzOPe9IAbd9LU04/9JlSEGYf5woR6K2GBCfVQoNecUI9zhFCS9PRjL3aVZORQ7a3GI4kl+DMcAFgPBCGUMGXktKqZaH6YqN9h04jjJjVm3wPA20EQQonUyYe50l1MRLW2GA+hawgAbwFBCCWVs5SWhIsWNRINOmIacdykMghdEACUTAhCKNk6VGAudxMTUa2txgPx6BoCwBtDEEKJ5ySlJeGiXxqLhh4zjThuykLXEADeBIIQSom25ZnL3cUyEQVuNu6I4YQuBwBKDAQhlB6OElrQUPRnc9Gnp7leB03pOUIXBAAlAYIQSpuIsszFrmJ/R6q5xfgXuoYAUBAEIZRCCjF9Gyra2FI0+TTX66ApVSd0QQBgxRCEUGo19mQudBX7O1LtrcYtD9A1BID8IQihNJOL6dtQ0eZWomnnuF4HTdv3H23TvV//kWMTEhKELg0ArAWCEEq/hmWYC13F7prYbqMm7q/16Vq7Fu169hO6KACwFghCsAkyEbWnq7Labck7iK/ZMTlTbTKZhC4KAKwCghBsRd26dR1u/k2xF5gL23USZQ6PO/wCABGCEGxH+fLlty6b3/HOLx+w/3Setar5bmMKriYFANyhHmxK40aNdjVqREQ80azzpoY7jH+3E1V2xP3uAWwaeoRgixiimcGiKbXZ5rtN51OxVDeATUMQgu0aUo1d1Iht97fx7zhkIYDtQhCCTetSkd3RRjzkqGntXcy4B7BRCEKwdWFlmAMdRJ+d5Waex4QKAFuEIASg6s7MyXfE22P4sSdNHEZJAWwMghCAiMhLQcc6iW9l8j0OmrRGoasBAAtCEAI85SChnW3FMhG13GNMw70MAWwGghDgGSlLa5uLIsoyTXYaY1UYJAWwCQhCgP9giObUFw0PYBvtNF1MQxYClH4IQoB8jKvBLghj2/9tPJKILAQo5RCEAPnr7seuayHuFWVcfw9TDAFKMwQhwEs192IOdhBPOsN9f5kjIpVKJXRFAFD0EIQAr1LDhTnRWfTHxTT32k0rhXesWKPejRs3hC4KAIoSghCgAN72TOu7S9NDBj7+34HY7r9+OHmm0BUBQFFCEAIUTK9R80p3IiKlR1Z2ttDlAEBRwv0IAQr20YgPtnbupXt4TH3tqNfomUKXAwBFCT1CgIJVqVLl+smozWPbXYnaluDf+ssLuI4UoPRAjxDgtbi4uLRs2ZKI/i5DETuNSgl9VAOfIwFKA/wlA7wZDxnt7yD68Sq38jb6hQClAXqEAG+sgj2zr72oxR6Th5zpWIERuhwAeCvoEQIURlUnZlsr0eCjxiOJ3M2bN5OTk4WuCAAKCT1CgEIK9WBWR/Ct23V0cHQUZTyaMmrgxx+OELooAHhj6BECFB53PUpSvnrGwDWpYw58v+hXocsBgMJAEAIUHsuyEjIREfE8z+M+FQAlEoIQoPBatGhRm48ts7y7fF5zt45jjbiMFKAEwjlCgMITi8VH9mx7+PCh0sl54GmHiadNP4SJhC4KAN4MeoQAb8vX19fNxXldC/G+eH7JTfQKAUoYBCFA0XCU0I42olnnTYdwU3uAEgVBCFBk/JXM+hbiPlHGE7fiDx8+nJmZKXRFAFAwBCFAUWpSlump3t+kY7eu83YEhjW/d++e0BUBQAEEC8KUlJTo6OhXrMfx4MGDkydPJiYmWrIqgLd3fN1C06hNme/OSWwzY8GS34UuBwAKIEwQrl+/PiAgYMaMGUFBQStWrMjzVZ7nAwICIiIiJk6cGBgYOH78eEGKBCgcpYMDqdKIiFWlOirthS4HAArAWH4WcE5OToUKFdauXdu6detTp061adMmPj7ewcHh+Tb79+9v1aoVwzAxMTFBQUH79u1r1KjRi7tKTU2tXr16SkqKpWqHl8rOzlYqlUJXYRWuXbvWrmd/nUSZaWD27virVWVnS746DoSVUKlUed7WwPI4jjOZTBKJ5NXNBOgRHjlyRCaTtW7dmogaNGjg7e0dGRmZp03r1q0ZhiGiihUrenh4pKamWr5OgMIJCgp6dP38/eO7NuyJGnneIcsgdEEA8EoCTKh/9OiRr69v7lM/P7/Y2NiXNd6zZ49Go4mIiHhZA4PBcODAgdynlStXfn7nAEJRKpXdlLQ/nh92zLShBWbZA1gvAYJQo9FIpdLcp3Z2dhqNJt+W165dGzx48PLly11cXPJtYDAYtFrt7Nmzc7f06tWrX79+RVswvA61Wm3uxMPzvqxBLfZLl1zR9/UzWeYVcSCshFqtFroEII7jJBJJgUOjAgRh2bJl09LScp+mpaV5eXm92Oz27dtt27adO3dup06dXrYriUTi6OgYFRVVLIXCm+B5HmdEXuRAtL4V32wX07SCLMDZEvmEA2E9cCAEZz5HWGAzAc4R1qtX78aNGxkZGUSk0WjOnTtXv379PG1iYmLatm37+eefo3sHJV2gM/NVPVGvKJPWKHQpAJAfAYLQz8+vU6dO/fr127t374ABAxo3blyjRg0iWrp0qfkKGp1O16RJEzc3t4yMjDlz5syZM+f8+fOWrxOgqAwPYGu4MJPPWGh0FADeiDB3n/jjjz/mzp27fPny6tWrT5gwwbwxKChIr9cTEc/zvXv3JiJzr5GIcnJyBKkToKgsbiwK3mZsUY57tyKWcwKwLgLMIyxCmEdoPTB9rUCnU/jO+4xn3hX7OBTjyUIcCCuBeYTWwHrnEQLYpvoezLggUf/DJlMJ/vAJUAohCAEsZ3JtViqi2Rdxz0IAK4IgBLAclqE1zcQ/n4zvPHzCe0NGX716VeiKAABBCGBZnnKSLH1/lyRso3vP1r0G5l4RBgBCQRACWFR2draekVLdd6hKeE6lJpcvXxa6IgBbhyAEsCilUmnPaejuCUq8yd0+HhgYKHRFALYOQQhgafu2rO0Us7Ju9BfSPvPEjh5ClwNg6zCPEIoGpq8VwpgTJq2JfosoyntT4EBYCcwjtAaYRwhg7b6tLzqUwC85eHnfvn1arVbocgBslDBLrAEAEdmLqcG570bvO+HgU91j4ucXju1HZw7A8tAjBBDS4Z2buQ+3ZL0zO6FKx8jISKHLAbBFCEIAIUnEYtJmEZHoSYKjo6PQ5QDYIgyNAgjpl7lfDRnTSs3KJBVqtGzVWuhyAGwRghBASB3atkm83UaXk9Nyn+iPu/zgqpa4iz0APA9DowDCk9nZ/dpYNPWMKVUndCkAtgc9QgCrUNOV6eVL732zto7xdr8eXerWrSt0RQC2Aj1CAGuh3zbrUPTJH7LqtOk3CmuQAlgMeoQA1uLvffv5MUeJYVON+i0799aqVUvoigBsAnqEANaikp8vcyOKOKP8ZmTdoAChywGwFQhCAGux5tcfWz5Y5bWwqbiMf5tO7wpdDoCtwNAogLXw8vLav+1PInovyjT/KvdZHXxOBbAE/KUBWJ1vQ9kfr5oSNULXAWAbEIQAVsdPyQyqws48bxK6EACbgCAEsEaf1RVtj+GupJfg24UClBQIQgBr5CylqXVEk86gUwhQ7BCEAFZqVHX2XhatO/vo7t27QtcCUJq96qrRnJycpKSk+Ph4IipXrlzZsmVlMpmlCgOwdRKWAo59OXDGMRdnx3reDrs3rmYYLMkNUPTyCcKsrKxly5atX78+PT29TJkyzs7ORJSRkfH48WM/P7/u3bv369cPN04DKG5arfbk/j3Gj6NTGObM2g8uXbpUp04doYsCKIXyBuG2bdtWr17dvXv3HTt2eHl55flqfHz8wYMH+/btO2DAgJ49e1qqSABbxDAMz3Pmx7zRwLI4kQFQLP7zp7Vw4UKTybR169a+ffu+mIJE5O3tPWDAgJ07d/I8//3331uqSABbJJPJRvZ/z2Nhc+mizuXdlFh6FKCYMDz/9PrszMzM5OTkatWqveZ33rp1y93d3c3NrdhqK1hqamr16tVTUlIErAHMsrOzlUql0FWUQqmpqReTtH0jVV84/xNWL7h27dqvbo8DYSVUKpWDg4PQVdg6juNMJpNEInl1s2c9Qmdn59dPQSKqVq2asCkIYAvc3d3tMmMzfx00+lB2y4Efbdr6l9AVAZQ2r3vWYejQoUuXLr148WKxVgMAL/p55Z/6nvO4ZqPS+iz/cdkqocsBKG1ed9FtHx+f4cOHF2spAJAvv/JekoeXDf4NmJjz3p4eQpcDUNq8bo/Q19eXiB48eLB169bExMRirAgA/uuzT8eFpx9y/LIWv+GTg2cud3pvAMdxQhcFUHq8bhCaL9328/NTKpXz5s27efNmcVYFAM84ODhE7diktFfQrEvp449GZzsePHhQ6KIASo/XHRqNjIx0cnKKiIioUqWKXC4PCAggoidPnjg5ORVneQDwlInjSCIjIoPMRa1WC10OQOnxukHI8/zvv/8+dOhQNze30NDQzMzMiIiIHTt29O/fv1jrAwCz4QP7Ll76TiYvNyReX+f0ODQ01NvbW+iiAEqDZ/MIX+3SpUu1a9fmef7q1auHDh06cuRIdHS0WCyOi4sr7hJfAfMIrQemr1nA4cOH3xn+afaAP5jU+zXOzL8cHfViGxwIK4F5hNbgNecRvm6P0DyNl2GYmjVr1qxZc+zYsTzPf/fdd29bJgC8tqysLL52JypblS9bNXn/lxzHYd01gLdX+L8ihmF69OhRhKUAwKsFBwcrru2hh2eZM5vKOCqELgeglHirj5OVKlUqqjoAoEDly5ff/vuiVld/ZP+aHquXV6rdIDY2VuiiAEq8lwahXq/v2rXrxIkTzU8fPny4c+fOnJwcSxUGAPkIa9DAxdXVNGBJ1qgdMa1mTZs9T+iKAEq8lwZhWlpaWFhY7n0HfX19g4ODf/nll4yMDEvVBgD50Gh1JFMSES9TqrU6ocsBKPFeGoReXl7Ozs6TJk0iIo7jdu3adfXqVS8vr/Pnz1uwPADI64sJYzw3jVRs+J94xUAnhfThw4dCVwRQsv0nCLOyslJTU3Of9uvXb/HixTzPf/HFF1OnTp06derBgwebNm1q8SIB4Jng4OCb/xzyzbppajl2JYU37tAd4zQAb+M/0yf++eefdu3aBQUFNWvWrEmTJk2aNBkwYMCff/6Zmpp6+fJloUoEgDzs7e1T1Xo+uBsdX5EmcTt48CAu4QYotP/0CD09PWfMmDFv3jxHR8cff/yxQoUKjenwPf0AACAASURBVBs3XrJkCcuyaWlpQpUIAHlIJBI5Y6QFncizak6d7pO+xIxegML7T48wICBg6NCh3t7ebdq0ISKVShUdHX348OHDhw+XK1cuMDCwWbNm06ZNw/14AQT3/ecT+8zdaLSzp3NbEp/orly5UrNmTaGLAiiRXneJtezs7OPHjx8/fnzIkCH+/v7FXdZrwhJr1gMre1lYRkZG5ZDwdIkbDVpGsRfD7q45uX8n4UBYDSyxZg1ec4m11w1C64QgtB54/7W8CRMnzYv14Hmerh1QpN97fOeSvb09DoSVQBBag9cMQixUCFBSDR0y2P7kCkq4QUNXaiNGjJsyQ+iKAEokBCFASVWtWrUm9YOpXnc6sIi/fnDnnr9L9AAPgFAQhAAl2KgP+tptm0oyJQ1allmu3vKVq4WuCKDkebMgjIuLu3HjhslkKqZqAOCNdOrQ3stJTsFd6O+5+rT4VRs2C10RQMnzBkG4YMECPz+/5OTkNWvWrF69GoMwANagb48uoiXvU6Uw6vfzhUcZWPsC4E29QRAGBAQMHDiwUaNGAwcObN68+ePHj4uvLAB4TbOmTnBgDFS5ER3+Vc0qdu/dK3RFACWMMNMnTp8+vXDhQrVa3atXr969e7/YIDU1dfbs2bdv3w4JCZk0aZJCkf89SDF9wnrgqn0Btena+8DpS3y32SRXum/96P6FEzgWgsP0CWtQjNMn4uPjo6OjC1UVEdHDhw9bt24dFhb2wQcffPrppxs3bnyxzTvvvJOWlvbRRx+dPn162LBhhX4tAFuwdP43ckdXqhhM9/7Ryj3OnTsndEUAJckb9AgPHz68cOHCwMDA7t27e3h4JCQk1K9fvxAvOWXKlEePHq1Zs4aIli9f/ttvv508efL5BidPnuzYsWNycrJEIklOTvbx8bl//763t/eLu0KP0HqgRyggnud9goLjnhiow0QmMz445fDZw5FCF2Xr0CO0BkXfI0xOTl66dGlERMTixYsHDx5sNBoLV9nZs2fDw8PNj8PDw8+dO8dxXJ4GDRo0MJfu6enp7+9/8eLFwr0WgC1gGOabqZ9IAsIpqA3vUflBQopGoxG6KIASQ1xwk395eXm5uLi0bdu2bdu2b/OSSUlJrq6u5sdubm4GgyEtLc3DwyPfBuY2SUlJ+e7KYDBkZWW1aNEid0vnzp2HDh36NuVB4ahUKqFLsGnBwcHKrxemz2tLId0yWYc5P/706ZhRQhdl09RqNS6tFxzHcRKJpMAe4RsEoZ2dXY8ePQYMGNCkSZO3uQGFvb29TqczP9ZqtUSUZwDB3t4+Jycn96lGo7G3t893VxKJRKFQTJ06NXdLUFAQBuiEgt+8gAIDAzs2rLPGWIdv0JsLbLFi7fhZUycKXZRNYxgGQ6OCMw+NFtjsDYZG7927N2HChDt37gwYMKBx48ZRUVGFq6xChQoPHjwwP37w4IG7u7tcLn9ZA47jYmNjfXx8XrY3sVjc6jleXl6FqwqgpOvYtpXs0Tma05xOrE5KST195ozQFQGUDG8QhOHh4YmJiePGjdu9e/eRI0caNWpUuJfs1avX+vXrzecwVqxY0atXL/P2TZs23bx5k4g6dep069at8+fPE9H27dvt7e0bNGhQuNcCsB09enR3S7lC782lnt8Zey/67Jsfha4IoGR4g6FRHx+f3J6ZWCwWi9/ge5/XtWvXDRs2BAYGurm56XS6/fv3m7fPnDlz/PjxAQEBLi4uP/zwQ+vWrQMDA2/evLly5UqRSFS41wKwHSKRqG2LJisenuE3TiT3iicTLqWnpz9/uh0A8vXS6RMGg2HatGlGozEwMDAsLCwoKKhoX/ju3bsajSYoKCg35LKysmQymVQqNT9NTU2NiYmpWrXqK848YfqE9cD0CWsQFxfnHxxumHiUZA7s8d9n1uWnT/pU6KJsFKZPWIO3vTGvTqfr0KFDnz59OnTocPny5f3799esWXPgwIEMwxRDtYWEILQeCEIrEdSw+XX/rhT9B9k5+Nlp71w4iQEVQSAIrUHR36F+1qxZ9+7dW7Vq1VvXVmQQhNYDQWgl/tqxs/sHo7mvrpMxR7l9ypZJvVq3bi10UbYIQWgNXjMI3+A83/Tp0z08PM6dOxcSEvJ2tQFAcWnZvJm7m+vjvXPo+kFVjnrx73oEIcCrFXzVaGpqqlqtJqL4+HiVSlXoBWUAwDKG9+/NXtxBk4/yH245duX+y9ajAACzgoPw6NGjPXr0GDx4cPv27b/88kvMZACwciMH93dx96Sfu9PGTzPSUpauXCN0RQBWreAg7Nat2969e3/66afKlSu7uLhYoCYAeBve3t5B7lJG6U6DlnHdZi9aulzoigCs2quC8PlRUIVCsWrVqo8//jgmJqb4qwKAtzJtwkcK1kRz21DM+XSjdOWadUJXBGC9XhqEJpNp8ODBO3bsyN0il8tFItHt27ctUhgAFF7z5s1dki9R+4kUMZjrMPnLuQuFrgjAehUwfSIyMvLSpUsODg4Gg2H79u1SqXTr1q0vu1+85WH6hPXA9AkrkXsgFvy0+NNtl40PLlCN1uIz66N3/Fk/NFTo6mwIpk9YgyKeR2heGtR6ItAMQWg9EIRWIvdAaLVar8B6TwavJ68AenCmW9wfW/5YInR1NgRBaA0KeWPegwcP3r1798V2ycnJmzdv1uv1RVYgABQnuVzerkVTJv0REbEZcR6uzkJXBGCl8gbhvn375s+fT0Q8z0+YMKFt27aLFi1KS0vz8/Nr0aLFokWLzHcQBADrN3fW1IDob5XfhDLbpp0+f+nvffuFrgjAGr10aDQxMXHLli1OTk4bNmw4dOhQu3btBgwY4OLiYjQan78jvLAwNGo9MDRqJV48EP61Gzzo/CM5enqs6nshcrO3t7dQtdkUDI1ag7ddYq1s2bJ6vb5Tp079+/dPSkpat27d/PnzpVLp5s2bi7pUACguOp1OzUvItx4RGSo3uXHjBoIQII9XXSyj1+vXrFkTFBRktavJoEdoPdAjtBIvHog64S2vVu9vcnB33THl3MGdvr6+ApVmW9AjtAZvfLEMx3F5viaVSgcPHvyKFHzxWwDA2kTt2DTO475iy0S1s1+DDj1X/7lR6IoArMuzIFSpVDNnzjSvr10gjUYzY8aMzMzMYisMAIqGq6trWHAtCumeM2rL4/8d+Gz2XKErArAuz4LQ0dFx5MiRY8eO3bBhg8lketk3cBy3ZcuWUaNGDR061NXV1SJFAsBb4XmeWIaIiHmDW5AC2Ii8fxUcxy1dunTNmjW+vr5hYWHlypVzdXVlGCY9PT0pKenkyZO3bt3q3bv3mDFjrOG21zhHaD1wjtBK5HsgdDpdWMuOjxjX9DsXyjg5DO7T8+vPp7BswWvuQ6HhHKE1eNuVZc6fP//333/HxcXFx8cTkbe3t7e3d8OGDSMiIgrcqcUgCK0HgtBKvOxAcBzXvf+QHeJQrkFfh53T5vUMGT54kMWrsyEIQmvwttMngoODg4ODi7oqABAGy7L3Yh5xvb4jqUJVq9uJczuHDxa6JgDrgLERAFvRvVM7h13T6fZRWeSc995pJ3Q5ANbiPz1CvV5/4MCBmJiYhISEpk2btmrVyrx97dq1d+7c8fLy8vHxadGihZ2dnRClAsBbmT7x4wpeayNP7o5sOzakSRuhywGwFv85R/jw4cNKlSr9/PPPAwcOlMvlz7fLyMj4/PPPf/rpp7t371aqVMnideYP5witB84RWonXORBjT5rsWPq+gfDXu5ViOEdoDQp594m2bduOHDkyTwoSkYuLy4IFCypXrlyUNQKAEKbUFq24zSVoMI8CgOjFIGzatOlLm7JsWFhYMdcDAMXOS0EDqrBzL2NlKACiF4PQ0dHxFa1dXFyKsxgAsJCJtdglv/7S5N0+3/6w4BULaADYgrxBKJVKcx9v27atQYMGZ86cyd0iFr90ugUAlCCbV/1mjL1yLPSzrw48+OaHhUKXAyCk/KdPPHnyZNCgQd26dZNIJNu3b9+zZ495O8MwFqwNAIpL5JET+qajqGxVddMxkYejhS4HQEj5BOGJEydCQkLWrVs3Y8aMI0eOfPXVVykpKb/99pvliwOAYtKmSUOH6CWUcl986Oc2TRsJXQ6AkPIG4c8//xwRESGXy0+fPj1z5kzzgqIDBw50d3f/9ttvsVwvQOnwvxFDP+8YFH7ic7l7ucZd+uNOMmDL8gbhpUuXPv7447Nnz9apU+f57V26dGnSpMmOHTssWBsAFBeWZSeM+9+xnRtCFE/ad363anj7T6bNErooAGH8Z0K9Tqe7evVqvXr1Xtb6wYMHZcuWfXGWoVAwod56YEK9lXijA5GUlFSzXe/U/+0jIo+FLa5HbXN3dy/O6mwIJtRbg8Isui2TyV6RgkTk5+dXBKUBgHUwGAyM9N/PtVK5Xq8XtBwAYWA6BIDtqlChQkS1csdWvp+pM1X09StXrpzQFQEIAEEIYNO2rP7txo0bt57Q6PuVn+jJSVrwtwCUMghCAFtXvXr16kQ79abxm86/6/CoSZMmWEMKbAruRwgARESeJ39a+f3nfVadq9moRVJSktDlAFgOeoQAQES0Zs1afswRDSvWKctt3LJt7IejhK4IwELQIwQAIiIHhYIyk4hImvbA091N6HIALAc9QgAgIlr50/c9P+iVpTfllK3Zvkt3ocsBsBwEIQAQEdUPDY25eoaI+h4yzb5M34YKXRCApWBoFAD+44cw0Ypb3MU0LCwMtgJBCAD/4Smnr+qJRhw3cYhCsA0IQgDIa2g1VsLSsluc0IUAWAKCEADyYhlaEi6adtaUoEGvEEo/BCEA5CPIhRlajf30FDqFUPohCAEgf5/XFZ16zO9+hE4hlHKYPgEA+ZOLaXFj0ZAtt0ca9oYEVWvfvr3QFQEUC/QIAeClfNW3Uxf0nn7D6f056yfN+ErocgCKBYIQAF5q8/Zd+lYfU8N+WX2Xbti2U+hyAIoFghAAXqqav69DzAnieXp41qtsWaHLASgWOEcIAC/VvXu3gyfO7PyxkdbeM3jaQqHLASgWCEIAeCmGYX754dtfiFJ1VHur8b0kvklZRuiiAIoYhkYBoGDuMvqlMfvBEVO2QehSAIoaghAAXss7FdlmXszY/Sl79uy5f/++0OUAFBkEIQC8ro88H64e1brHshMNugzcsHmr0OUAFA3BglCtVl+5ciU7O/tlDdLS0q5fv56VlWXJqgDgFVatWWvqNF3bcUbqkM1f/vCT0OUAFA1hgnDnzp0VK1YcPHiwr6/vhg0b8nyV5/n69etXrlz5vffeq1ChwsyZM4WoEQDycnN2lGTFExFlJjo5OgldDkDRECAIDQbDyJEjV6xYcebMmS1btowaNUqr1eZpM3HixNTU1CtXrpw6deq77747c+aM5esEgDzGjBxWN+Fvjx+byNYMbzDiS6HLASgaAgTh0aNHiahz585E1KxZMw8Pj8jIyOcbMAzTo0cPkUhERAEBAV5eXvHx8ZavEwDysLe3P3VwT9z5I3cvnlqvrXYsCetxQ2kgwDzCmJgYf39/hnk6G6lSpUoPHz58WeNDhw5lZmaGh4e/rIHJZDp37lzu0woVKpQpU6boigWAvKRSqbeUVjcT9ztsOttF7CETuiCAt1NcQTh+/PgXr3N555133n33XbVabWdnl7tRLperVKp8d3Lv3r3+/fv/8ssv7u7u+TYwGAwqlWro0KG5W/r16zd8+PC3Lh/e2MsOIliYxQ5EfSX1rCDuc8C4uYmBxST7F6jVap5Hj1lgHMdJJBKJRPLqZsUVhA0bNnzxzJ+vry8ReXp6ZmRk5G5MT08vm98ahrGxsa1atZo2bVqvXr1e9ioSicTJyenChQtFUzS8HaVSKXQJQGTBAzGnEbXeqX3vi+Wy5Gt9u3bo2a2rZV63RGAYxsHBQegqbB3HcSaTqcBmxRWEr0ivOnXqXLt2LSsry9HRMScn59y5c/PmzcvTJi4urkWLFqNHjx45cmQxVQgAb0nEkN/hr3+/Q3zwiCNzvnBQKNq3ayt0UQBvTICLZapWrdqyZcuhQ4ceP358xIgRderUCQ4OJqKVK1d26dKFiHJyclq2bOnm5ubk5LR06dKlS5deu3bN8nUCQIGijx3lO0ykCrUzw0fv2H9Y6HIACkOYeYRr16718fGZNWuWs7Pz1q1P16eoWLFiWFgYEZlMpmbNmtWpU+fcvx4/fixInQDwao3q15Me/40yE9mTaxrUDxW6HIDCYEr06dzU1NTq1aunpKQIXQhQdnY2zhFaAwsfCLVaPXnm7FMXLsmCO4qaDotsL5Zi3UYiIlKpVDhHKDjzOcICL5ZBEELRQBBaCaEOBMdTj4MmuYjWNBfhGlJCEFqH1wxCfHgDgCLAMrS2meh+Nv/lBU7oWgDeDIIQAIqGXEw724jX3OX+uIMshJIEQQgARcZdRjtaiyadNu2PL8HnXMDWIAgBoCgFODObWor7HjJeTkcWQsmAIASAIhZRllnUSNQp0hSnRhZCCYAgBICi954/O7I6++5+k8ogdCkABUEQAkCxmFqHDSvD9IoyGnHpDFg3BCEAFJeFDUVihhkVXfCqxwACQhACQHERMbS2uehsKv/dZfQKwXohCAGgGCkltLutaPF1bs1dZCFYKQQhABSvcgpmd1vRJ6dMhxJxESlYo+K6HyEAQK4gF2ZDC/H7u1LaXF+YnZo8YfTgRg0bCl0UwFPoEQKAJTTzYhxWD1mTVfmvioO7DB3/8OFDoSsCeApBCAAWokpN5Bv2Jd96WTW7njp1SuhyAJ5CEAKAhZTzcGOu7KGU+8aLu8sH1BW6HICnEIQAYCG71698X72v8bGp/cdOGXLbPx4LsIF1wMUyAGAhXl5e63772fz4+8tcxC7TwQ4iPyXu4wsCQxACgAAm1GKVEmqxx7S/vaiyI7IQhIQgBABhjKzOihhqusu0r70oyAVZCIJBEAKAYIYFsPYSarPXtLedqJYrshCEgYtlAEBIfSqx88PY5r9dqdqwVYWgelNmzRa6IrA5CEIAEFgvf9Z+w9g7HebHjTv+a9S1w4cPC10R2BYEIQAIz6RVkVcAsaLsCg3u338gdDlgWxCEACC8Ns0jHDZ9xJz4g46tuFaupdDlgG1BEAKA8JYvmvfHyFbf19Kc3Lf9iK7ckKMmPe7aBJaCq0YBQHgsy3br2tX8+Egl6nvI1HKPcVsrsbtM2LrAJqBHCADWxV5M21qLWpZjGu4w3nqCZdig2CEIAcDqMEQzg0WTa7NNdxkP43a+UMwQhABgpYZUY9c2F78fZVx5GycMoRghCAHAerUsxxzvLP72Etd57s6A0CZ1IlqdOPmP0EVBaYOLZQDAqlV2ZLY3zKg55UvDmJ2Uo+nxQbdH186KRCKh64LSAz1CALB22UkxTpXrkr0ruZbPUZZNT08XuiIoVRCEAGDtgoKC7OPPs6fWiY8uyVZrp91yVRuFrglKEQQhAFg7uVz+z/6dn1VM/KqO6e6RHURUawuuJoUiw/B8Cf6fKTU1tXr16ikpKUIXApSdna1UKoWuAmzlQOx+xI88bupQgZnXQOQgEbqa/KhUKgcHB6GrsHUcx5lMJomkgP9F0CMEgJKnYwXmSncxEdXeajyEriG8HQQhAJRIzlJaEi5a1Eg08LBpxHGTyiB0QVBiIQgBoATrUIG53F1MRLW2GqMS0DWEwkAQAkDJZu4a/txINOgIuoZQGAhCACgN2v971rDWVuNBdA3hTSAIAaCUcJLSknDR4saiwUdNI46bstE1hNeDIASAUqVdeeZyt6ddwwPx6BpCwRCEAFDamLuGvzYWDT2GriEUDEEIAKVT2/JPLygN2GTcGYsbOcFLIQgBoNRylNCScNHvTUX/O8H1OmjKyBG6ILBKCEIAKOXaeDNXuou9FFRzq3FHDLqGkBeCEABKP0cJLWgo+rO56JNTXK+DpnR0DeE5CEIAsBURZZlL3cT+jlRzi/GH3WfrNW8fFNZs+67dQtcFAsMd6gHAhijE9G2oqHU5rm3joaYRG0juOGxS1wb1QsqWLSt0aSAY9AgBwObUc8h2c3UjD39ycFf7NLh0867QFYGQEIQAYHOcnJy85CQ5tow5u5m5e/z9O0GTz5gy9UKXBQJBEAKALTqyZ+v0QM2EMndvHtp+/n3njBzy32BAHNom3KEeioaN3Bjd+uFAFNqDbP7bS9ymB9zwAHZybZGz9K32hjvUWwPcoR4A4A34KZkl4aJzXcToHdoaBCEAwDPmODyPOLQlCEIAgLx8X4hD8/Js2dnZ169fNxqNQhcIRUmYIDxx4sSwYcOGDRsWHR39imbHjx+fPHlyQkKCxQoDAMj1fBxW2mjo/evBSqFNm334dbWQxrg0oTQRIAjPnz/frl27unXrhoSEdOjQ4dy5c/k2e/LkyejRo+fOnZucnGzhCgEAcpnj8PS74sgl36QM+yul78qHoaN+WrpC6LqgyAgQhAsWLBgxYsTo0aNHjhw5evToH3/8Md9mn3766bhx48RirH0DAMKr7Mj4ODDE80TE8dz6+/xfMZwRK3iXCgIE4T///NOsWTPz46ZNm/7zzz8vtomKirp3797gwYNff7eRkZHzvptz5cqVIikSACCPH7+a7rm8a5nV/SudXzbpwyE/X+d81hvGnTRdyyjBk9BKn8zMzLnfzp49a2ZSUtJrfkux9Le0Wq1KpcqzkWEYd3d3IkpKSnJzczNv9PDwSExMzNNSrVaPGTNmy5YtDMO8+oUMBoNJp+nZvdvDe3ddNamdK3n0Xjhv8brNwcHBRfSjwOtSq9UFHi+wAByI4lMvJPjykb8TExN9fX1FIlGvqrq72cyq+2yL3SJvBX1QyfReRZPi3/dUtVotaLG2q32ziG5lWAcR0/7PNZEHDymdXQucR1gsQbh27drp06fnfSWx+NGjR0Qkk8n0+qfXI+t0OoVCkafllClT+vfvHxAQUOALSSSScfX9B/voepyK/7lLsEwscpFJ/pozveYHvViF8ul/9kpG7pD7lFhcKFssTCbTi4cSLA8HolgpFArzB3qzWgqa60lzwuhwIi29yc66LO7my4wIYOq6MRzH4UAUE16fw2myzf/xWhWnyebU2Zwmi1NnJyUm2qvTegfWIqIz8dfPTe3X7KvlZF/51TssliAcOnTo0KFDX/bVChUqxMbGmh/HxsZ6e3vnabBz506O45YtW0ZEer2+c+fO33zzTf/+/fPd20+Xk2YcvFL5YedTCUlNfdz+SVYFtmgqKVOe02Sb0pMMcXee/r7UT39rrJ2CtVc+i0mFklUomRe2sPZKRoTTk2+AZVkWHzKsAA6E5bFErctT6/JsvJpfc5fvGcW52vED/SSDAlllAV0ReIbP0ebG2/Nv2s8e/PsfEeW+Uee+aYsUSomHd4VqwckL1iZk62Ri9pqGqTd3s9jZucCXFmCJtdmzZx89enTv3r1E1KlTp4YNG06bNo2ItmzZEh4e7unp+eTJE457eg7ay8srMjIyLCzMzs7uxV3lLrGWlJQ0YkDfR7Ex4U2bzV+8RCQSvezVeYOe02RzWhWnUfHabE6j4jQqTqviNNmmrHTTk/R/N2YTw7ByB1bhwMiVrMKBVTiw5p6l3MH89Ol2uYPIya14flUlCVb2shI4EILjeIpK4Bdf1R9MYtuWZ4cHsK28S95gdU5OTr7vum/qFW+55o1cVrrpSVoRvuWePHFi2icfGQz6KV983bZd+9dZYk2AIMzMzGzRooVYLGYYRq/XR0VFubi4EJGDg8O2bdtat279fGOZTHby5Mm6devmu6tiXWu04OOnUfHabNOTdN6of93j5+hKLz9/c/PmzY9HDktLSx324dihI0cVxw9VfPD+ayVwIKyESqUySBw2PeAWXeOMPA2qwg6uxq5bsfTn5X/4VKiwYuF3Pj4+QteYP71e/36XznF3b+UQu3TVugZhYS+2ed0ehTqLWLbAt0eRk6vI0Y2RvN3qrvl5zbVGhVl022QynT59mohCQ0NzJ0jExsaWKVNGJpM93/LBgwflypV72QcTK1l0+7kBa9V/u/ZZebr5ZDKyCiXz3ADsc+OxDq0GDP+ijoevs2LUoXvfrd5UNzhE2J/rjeD910rgQFiJ3EW3eaKjifxvt7i/Dkbr98zVD1lNMedDTn939tBeoWv8D06j4nRqXqte8+f669tXjwupmKTSfRQdu+vrCS8OVz47x2TvyCoc8pxXev5dTthzTK8ZhMKUKBKJGjZsmGdjvp+P/Pz8LFLRW2GkdiKpncjZvcCWvNHwYl7y6mxj4kNOk63JSA9w9yeiMDfpyelDywT5sDIFI7N/+q/CgZUpWJk9I1OwcntGpmBl9qzcnpHbP90osy+Oj1QA8DYYoqZeTFMv0cJbd8YHtiWpgqqEX9s8Yc4lLqIsE+rBSIr8lC7Pc1o1p1XxOg2nU/M6NafVcDoNr1P/u928Rc3rNJxWzevUnE7Dyh0YuYKV2T8+c7uMnZiIXOVSXY5e7ObFVqia59KKUnbVIa4HsShGLBE5uoocXfP9atCOUwsv3KmklO5M0O09dqisuxtv0POGnBfHH0xZGXxa0vMjtJxGxRv1nEbFSKR5BiIYiR0jlrw4VMtIpIxY+urR2jxWLFuyac3q6jVqzpj9rZOTU9H9YgBKv46tmn3xVWvmyAKS2PlXr5mh5yee5q5m8DVcmPCyTKtybGNPRp7fW3KecUjeYOANOS+eqSGj3rzdpMpkROLc9wFGImUk0twxSYmjT+47A4mfbhcpnYl9emnFoOTkthGN4vTxF1PUoyZNc2jaxaK/JiHgfoRWRK/Xr1r5e1J8XO8BgypVqlSIPfBGg/kzIKdV81o1p1NzOg2vVT/9MKjTcFqV+emzD4OGHFZuz8qVjEzByhWMnT0rVzztg8rtWbnDv11PxYF/zv0ye+bMsIpHHqVfcan2+58bn3/pAkfkkpKSpn06Pi42bIdisAAAGNxJREFUdtDI0e/36fv6P1RRnbS3ERgafX2vOW6WKysr6/NJE65fudy5e48x4z95deMX70d4+vTpKQN6Lmxe+fLj7FWJ/F+rlnNatVateZCiepSuepypyc5We7Pq8qzajdE6cWrG/Per1z3/l2geImLl//0Lldvn/s2aGxTyN/Kv7Ozs6Ohof3//qlWrvuWuhGXV5wiLSikLQmFwHKdTc1oVr9U8C86cp/H5fHAu2HPcV2JoX9mTJ+q84dSeIW1YOwUjUzBSO0Zix0llUoUDI7Vj7BSseaNUxsodGKmMkdoxdvIuffq/72YM8lBOiI75fvWmkJCCz4AmJyf36NhOn50pkjts2v33izNtCuHevXvOzs65SzqUPqU4CBMTE6VSaVEdu1Urfvvuy1kiht7t+d4Xc+a+zreMGNivStKV1r5uX5+O7fW/T99p3ZzXaTidljfo+Byd+WMlr9dxOg2v0xi0atZk5LVqTq/jDTpep1l79rbJZOpXswIRddp4Zv/Uwc9OasjtWZnCaOfwwKC4prE/mSU7maVwdLCvVc6+ia99cy/WXVZgdZAPBCEUsejo6M+G9ptU1+toQpa6SuiP8+dzOVpeq+ENOk6fo0lPtRMxvF7H52j/fUfI4bQqPkdnfptoM2/trl71iWj1lUcSkbhPvSqMVMaIJKzCgURiRipj7eSMWMLI7RmxhJHYsTLFZ7+srG183Na/zNHYtCOM58JvvmQkUhJLGYkdI5EwEjtGLDUP+7xO/SaTqVuHtsbkRymanEEfjhs9bvzr/+w8zy+cN/fvnduDQ+tP+/JruVye+6Xr16//sexXz3LeI0b/z97e3rwxKirq4vlzrdu2q1mzZr57O3LkiEajadWqlVQqvXr1amJiYkREhEwmi4mJefjwYUhISG5n4saNGxqNpkqVKteuXVMqlXq9nmEYmUymVqvt7e3v3r2r1+ufnzuo0+nkcrnJZGJZ1mAwiMViFxcXDw8Pg8FgMpnkcrlarQ4MDFQoFJmZmdeuXatevXpcXBzDMDVr1rx//35cXFy9evWenwmekZFx9OjRKlWqBAYGElFsbOzhw4dr1KiRu35TVFTU5YsX27RrZ27wCklJSV9MnZQYHz9y/Cdt27V//kvXr18fM2RQenpa9/f7TJv15YvfO/KDATdPR2sNpnf7DJg684tXv9DzOK2KeOK0KiKe12p4nuN1aoPBENau8853a0lE7ODIG99OGl/Z0403GjitikxGLkfL63W80chrVbzRwOtzuBwNGY3vrIra3KO+mGUOx6SezzRN7hzByOxZO7n5ox4rUzBSGSOVsTIFYyfP4XiFkwsjt2ckMkZqx8rs78Q86vdO+4m1y1xK0zzyCPhj4+ZXlG3k6FI6fzyJj07mDyRwLlKmlTfT2JNp5sX4OJS8yRhCseqLZaAkaty48cS5i7asWx3QrvaMjz5mZTLW4dlMVUN2tv0rOyJBJ2J/v/ogyEW+LVa9avuesr4+vF7LG42cVkVGA2/I4XRa3mTgtWre/FSrzsrM9CxjR0Se9nZpN+5lHdjIG3LMjc2nSXijnjfoeYPeHIfMcxlJYnNSShip+V/Z8ZsPXNNivmxWycjxned+27+OHyuREhFjJ2NEkn8fiImIsZObz5ewMoX5ooD123acWb90dnD5bef3z5hs+G7BIvMPlZCQ0PedDh/XcH9wWT/gUNSWPZFE9MuiBXuWLWxdzn7Ikp8Wrlof9sJ1YR/0fs94/4qTVPTtDGmb9h2ObVtf2Vn+WaZp1NjxS7//ulYZ5Ucp2r+PRnt6en784cgbx6NYnrsc97iSm+PjJ9k6E1/W3u5RlqasUpGUpTZwHM+TmGEkItbI8UaO44mkIpZhyGjiTBzZiVlXudTE81l6k7ejIj5LE+Zb9paan7No8cQPRzQoq/z7Rmxtbw+RWJwpVRozUwPdHMZn6PcdO2HueMXFxXVsHtG6nMOPKeo+Yz5p1KRZn87tu/s5/pGk6Td+8gdDh//84/y/V/zc0kvRf/GPyzZsC35lR39Aj6593A2VvRQTxo6quOmvqpX8uByN+Usj+/ScVcvV36Xyx3+t21/FJyI0hNeqiYg3Gfkc7Z3YuNjTR9e0DuCJ3l2x5INATxmZiOf/G3JqIuI0KiLitNlExGnVxPOszJ5YhpXZE8MycgXDsIxMoTGYFAwnEbFE5CkXZ6alUVl3VuEgdvUksZi1Uzz930buYP5fiLWTk0jcRj/th1ORrSo4/X4rbfrPK8o0b/6KH5ZTqez+OzRarXr1xes2rV2+1Des6mcf/u8V30tEYpZC3JkQd2ZcDTLxopuZfHQyvyuW/+SUyVHCNPZkwssyrbwZfyVCsQigRwhFo8ARObVaveiHeXEP7/cdMrxho0avs88zp08P692jZQXnQ4+eLFy5Jjwi4mUtzXHIG/W5GUn/ZiRv0PNGA6/XHTlzfvOaVV809jdyfOetF44umEkmIxHxOTreZPj3gXmLljgT/ftOSkSf7zzerqJrfW8Xld44dNeF9d1CzdG7+1b8/eT0D+v5EVHnjaf3jnqXkSl6/LptaZtApZ34RFz60Qx+RvdWuXUyIrGOEXX7+petXesS0bSjd04lZES+F8oyzE/nYnbeS93atY69RLTpRoLWs9Kg1o3f+ez7v7rVnXX0Zks/j6+P3Z7cuMqRmNQnOcZ+Nct/sv+ag0TkKBPrjFxWjsHEk4tMYidis/VGld5kJ2ZdZBKOp8mNq0yNuj4yxC8mU3MhOfOLptXLO8r3P0hZej1lYt1y5ZSyyVHXVr8bQkSNfj92eEBjqYhddSVO7h80qHVjIvp1+z5ZamzP6uX0Jq7X/9u797Coyn0P4L+15j7DICODOAOCWolF6VMclYi8IKjb8pBZeMFjFzuQ6dn0aD2Vx63Pbteuvb1VVm7tqGXEIe0i6dZjcLi0QUq30fEGJYpymwGC4TIXhnV5zx+LRmIUKM01M+v3+WOeNcvF4ofzyNf3Xe/l4KkZd8XE0vbkMWH2HvbJo+fy1jyWuvmDPXNidQpZyeUfT3SQl+bGE1dvthGeIz/lHGFZ0tP9QG753xfdCwA7vr1kHj4s9c7RtKq36Tl7+xd/T5sEAB+cqtNHjF6SMJHS6IS/MUqlqaqzvr77o3dSbgeA+Qf+74ttfwkONQrzrwGA1gYBANXnmNYEAVC0RjfAQLCMx5bazv4zRCm/ROuPFP9jKBvdMAyz8913zlSc/Ne0Rb+b+8DAF3s/I7xRLnaRggZSaiUlVsLwJDGcFnLxHiMuMtsftgiRb9HpdC/+Yf0v+pJJkycfPXbi5MmTz919t8lkGuDKoXSQzp4ye8fh4qdLaloc7hXPv2RY8MzQK3ngrs/f+8NzWqXsYE37rMeejvzTq0LE3ltV9eGSR5e6mRqbUzfCbHzmz6TbeUvZpeOWppmjjceb7ONmPBicnOa5D+FYtcvp5P/WzXIquazBxSrUGnsPF6ySW7tZtUbT6WZ0Clmbmw03RaijbushwBGiktF2NwsAFAXdLK+gqR6WEEJ4ClieUAAAFE94ngABoCiKAPCEMBwvpymG4zme0BS4WE5B0y6WAwAXw8mVyq4e1nNnAkAAunrYUI2yzc3dYY5SRt4KAGHjGisvXwAAq92tDxkeGXv3D8fyk8fA9612U/TY4OS0MQePn7S0T40K/WeLI2buouGLllKa3myjaBn1U85RcjmlVI86Pf2z75tuDdEctXZ//MERc5/JUdMbMtaXf3XnMOXHl+1HPtxp/Pnz4HhC9OWVGcXnXAybsmCRad6yoX9217Lzg+xvvvnG6XROnTp1gLWo+lIoFCuznr3+b32dxuqpjPFUxngAgItdvd2nO6v4VjeZFEYlhtPJEdTdoRSNqThk2CJEN4a/jNGorq4OCQnpu27yEO3LzT2a99k98QmZK1f1bT3k5nz07paN4SNHvv7mO8JY36ampsxl6bWXaqYkJL6187+8/ze6779zXv7Pl4Dwi5c9/i/3JqxZ+bSapiZOjn9q5e8zly0JVtBB4RGfH/lSo9G8++bWnW9upYG0O13DNcqW9k6ZTKaW0zanO0yvbety8oT08LyMomQ0RRFgCSEEhLc84RmeKGV0sErh4ngeqGEqZavDNd48wqHU7fwwZ/mShQYZf6mlPUitpGSyuMRp5V+VBCtlYdFj9x88rFQqAYBhmCUPp16oPEsrVe9l594RG/v4okerTp8KCTXuzt0/ZswYi8WS8W+L62tr4++7f9t7uwZuV7W1tW157dXGutqn/uPZhPvu6/tHPM/nHThQe/ny/AULrrXkyvnz51Uqlc8uyNLPb9civJZGJylr6s3F8x1k8ojeKRmTR1DKgJr19wvgYBl0U/lLEPomzxQRQojD4ej7C9TlcjEMo9frm5qagoKChF2WVCpVV1eXTqerq6uTyWQURQnL83IcJ3y5TCbjeV54ValURqORZVmWZXU6XVdXl9lsFrZqam5uDgsLE+6p0+l4nnc6nd6/vnmex4W8f6mbH4R9NbngeAtf1kQKGsgPHWTKCOq+cCoxnE4cSamH1PoNEBiE6KbCIPQR+EH4CHGDsK8uBr5pJgWNfKmVfNdKxof0DkCdZqKDA31zDHxGiBBCCPQKSI6gkiNkAOBgobyJlDbxb53lFxdy40N+GoBqpg19Vq1gGKaqqspsNgfwjNu+MAgRQkgqdPIrocjwcKqNFDSQvef5zFLOpOmdkhGn7ZrzuznthnFgOff2n9amLZgvdtW/OewaRTcG9sj5CPwgfITvdI0OBcPDiRbyDyv5ysoX5e7sdjMk6Rno7grf8cAn+ccidGDWUio/fLiIXaMIIYSGREFDQjiVEE69MJHeflmZVehgAIBnHRz94gmuwQkWJwlWwEgtNUoHJi0VqQOzljJrqUgdjNRSIzWDfgefhkGIEELoimXpi7e//2DT7jJovbznrb/Mnd0bEzY3NDqJxdn7etZGChqI54xBBWYtZdL2f70lmArx+d3hMAgRQghdodPpviv93/r6eqPR2HfVWYMKDCoq1gAA/efq9/DwY/eVjGx0kpM/9mbkhU7SzfVPx7HBYNJQZh2M0v0G2zH+chiECCGEfoam6V+0cIGSFnpKIc4rIwHAxYLFRRodV173X4RGJ29xQq2D6BW9oSi8jtVfiUyT9iYtGodBiBBC6DekkcNYPTVWD95NSfDqcfU0JRudxOYGoce1bzoKr1FBlP4aI2AsFsuft77tdve88Punx/RZxm8AGIQIIYREM0CPq5uDVjexOOFiV29T8qwNLC6+0QG1DsLyP+tlFV5HqmHpnIfq71tDFOrD89KqjpcMZVtvDEKEEEK+SCX7qcfVeJWmZLMLrC5S5wCLkzQ44Fw7KWgkNQ3WeiqMv2c+AFhPHbrj7bP/8+/33DHYqgAYhAghhPzPCA2M0FAThkPfpiTPm6P+2thQdwpUOoPl5KeLXhulG/xWGIQIIYQCBE3Thz9+P2vdK253z193bYuLMnAcN+hXYRAihBAKHBMmTCj6Yp9wLOzKMigfmMGBEEIIiQeDECGEkKT5dxASQnp6esSuAgHDMAUFBWJXgcDlchUVFYldBYLOzs7S0lKxq0DQ2tr69ddfD3qZfwehzWaz2+1iV4GgtrZ2zZo1YleBoLKycv369WJXgeDkyZOvvfaa2FUgKCsr27Rp06CX+XcQIoQQQtcyxH0GMQgRQghJGgYhQgghSfPveYSEEEJISkqK2IVIncvlslqt+EGIrqur6+LFi/hBiM5ms+EH4QtaWlqGcpl/B2FMTMzLL78cHx8vdiFSRwi5fPny6NGjxS5E6jiOq6+vj46OFrsQqWNZ1mKxjBo1SuxCpM7tdvfdUvFaqCE+S0QIIYQCEj4jRAghJGkYhAghhCQNgxAhhJCkYRAihBCSNP8eNdrXDz/88O233yoUisTExPDwcLHLkS6O486dO9fU1DRz5kyKusq+0ug3wvN8YWFhfX39/ffff8stt4hdjnQxDHP69Gm73T516lSxa5G0oYdCgLQIt2zZkpSU9Pnnn2dnZ48bN+7LL78UuyKJ+u6774KDg5OSklJSUoayHya6gdLS0p577rljx45NmTLl4MGDYpcjUUePHtXr9bNnz05LSxO7FknbvHlzUlLSgQMHsrOzY2JiBt4VIECmT9TW1ppMJoVCAQCvvvrqoUOHysvLxS5Kiuw/ue222xiGkcsDp8vBx5WXlz/00EPV1dV6vf6jjz56/fXXT58+LXZRUmSz2XieP3PmzMKFC61Wq9jlSFffUHjllVeOHDlSVlZ2rYsDpEUYFRUl/MAAYDKZcG8msQQFBY0cOVLsKqTo0KFDs2bN0uv1ADB//vzKysqamhqxi5Iig8EQGhoqdhWofyi43e4BLg6QIPRwOp2bN29+6qmnxC4EoZuqoaEhMjJSONZqtSEhIQ0NDeKWhJAvcDgcg4aC3/Rctba2zpw50/v8pk2bkpOThWOGYRYvXhwbG5uZmXlzq5OQqqqqRYsWeZ/fu3fvhAkTbn49SMBxXN+hSXK5nGVZEetByBcwDLNkyZK77rorIyNjgMv8JgiHDRuWnZ3tfT4qKko44Dhu6dKlPM9nZ2fTdKC1dH3H6NGjr/pB4DBFcZlMpubmZuGYYZi2tjaz2SxuSQiJi2XZ9PR0QsigoeA3QSiXy++8885r/SnHcY8//nh7e3teXp5SqbyZhUmNWq0e4INAYpk+ffqqVatYlpXL5YWFhSaTCf9rgqSM47gnnniio6MjLy/P87DwWvwmCAe2adOmnJyc9PT0rKwsANBqtVu3bhW7KCnq7u7Oysrq7OwEgBUrVuj1+i1btohdlCTMmTMnNDT0kUcemTp16htvvLF27VqZTCZ2UVJktVo3bNhgsVg6OzszMzMjIiLWr18vdlFStHHjxpycnKVLlw4lFAIkCKdNm7Z9+3bPW2wUikUmk8XFxQHAjBkzAECtVotdkVTQNF1UVLRnzx6r1bpr1y7cCU8sarVa+Cfw4IMPAsDw4cPFrkii+oWCSqUa4OIAmUeIEEII/To4qAQhhJCkYRAihBCSNAxChBBCkoZBiFDg4Hl+gAUVEUJXFSCjRhFCAFBcXNxvHv3GjRtbW1sBQKFQBAUFAUBnZ6ewMciKFSuio6NFqRMhn4JBiFDg+P7775OSkvqeef7551evXp2bm1tTUyOMIGcY5uLFi+np6QsXLsQgRAiwaxShgGG324cNG+Z9/sSJEzNnzvTMo1IoFDExMQsWLBg1atTNLRAhH4VBiFCAOHjw4Ny5c/uddDqdx48fnz59uvC2trZWOAgKCjIajTezPIR8FgYhQj4hLy/P++S2bdvOnj07xDt0dHSEhIT0O3ns2LGenp5p06YJb3NycoSDJ5988tdWilCgwWeECImvsrKyoqIiNTW13/ndu3d3dXXFxsYOeoeamppbb73V+3xJSYlGoyksLDx8+HBVVZVnJSmdTnf9ZSMUGDAIERLfp59++uKLLwrHTqdTq9UKx3/84x/tdvtQ7pCfn798+XLv88XFxampqY8++mh3d/cnn3yCS18i5A27RhESGSEkODjYs0D5vHnzhOkNABAUFDSU7Y4JITRNe283ITwgTE5ONhgMJpMpKipK6CN1u90VFRU39IdAyI9hixAhkTU1NUVERHjeGo1GT6SVlZWtW7du0DuUlJQkJiZ6ny8vL+/p6fGMlPF0ve7bt2/WrFmeyxwOx/vvv19fXz9mzBgAYBimrq5u7dq1wcHBv/JHQsivYIsQIZFZLBZhzjsA7N+/3xM/FRUVMpmMoigA4Hm+uroaAOx2u/fwmaqqqvHjx3vfubi4ODIyst8OvR0dHZWVleHh4Z4zOp1u0qRJxcXFGRkZGRkZK1eu1Gg0wvdFSAqwRYiQyCIjIx9++GEh6qqrqy9cuLBhwwaDwZCdnZ2fny9ck5eX197evnfv3tTU1IaGhoKCAmG7UQBwOBxXbboRQvLz8z3jRQVnzpzJzMz03qG0pKREmInf3d0t7Kin1+tv/I+KkE/CIERIZGFhYVlZWbt27UpISMjJyamoqHjhhRfCw8Nzc3MNBoNwTXR09IULF2JiYuLi4lpbW1taWjxfftXpg/v37y8tLY2KilKr1Z5hOHa7vbm5OT4+fvLkyf2uLy4ufvbZZwFg3759y5Yt874hQgEMN+ZFyD/Mmzdvz549RqNx9erVqampnqbejh07MjMzr+fOLMuaTKYNGzbYbLa2tjbv9iJCgQ2fESLkB1iWtdlswlow5eXlcXFxwi4Tly5dGjt27HXevKKiIiYmZtWqVevWrcO2IJIgDEKE/MD58+cTEhKE49tvv/3IkSOTJk0CgPz8/H6rbP8KRUVFM2bMAACKolJSUq7zbgj5HewaRciPeTLsV7Pb7ampqcuXL1+4cKH3TESEpACDECHp4nn+s88+0+l0drt93LhxEydOFLsihESAQYgQQkjS8BkhQgghScMgRAghJGkYhAghhCQNgxAhhJCkYRAihBCStP8HCgQvx5G918gAAAAASUVORK5CYII=",
"image/svg+xml": [
"\n",
"\n"
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n"
+ ],
+ "text/html": [
+ "\n",
+ "\n"
]
},
- "execution_count": 8,
"metadata": {},
- "output_type": "execute_result"
+ "output_type": "display_data"
}
],
"source": [
- "plot(dlr.ωn/para.EF, imag(sigma_wn.dynamic[1, 1, kFidx, :]),marker=2, label=\"imag\", xlabel=L\"$\\omega_n/E_F$\", ylabel=L\"$\\Sigma(i\\omega_n, k_F)$\", xlims=(-2, 2))\n",
- "plot!(dlr.ωn/para.EF, real(sigma_wn.dynamic[1, 1, kFidx, :]),marker=2, label=\"real\")"
+ "plot(ωn_grid/para.EF, imag(sigma_ωn[:, kFidx]),marker=2, label=\"imag\", xlabel=L\"$\\omega_n/E_F$\", ylabel=L\"$\\Sigma(i\\omega_n, k_F)$\", xlims=(-2, 2))\n",
+ "plot!(ωn_grid/para.EF, real(sigma_ωn[:, kFidx]),marker=2, label=\"real\")"
]
},
{
@@ -333,7 +397,7 @@
},
{
"cell_type": "code",
- "execution_count": 15,
+ "execution_count": 10,
"metadata": {
"tags": []
},
@@ -343,521 +407,779 @@
"output_type": "stream",
"text": [
"47\n",
- "-1.8006326323142126, -1.8006232505919049 - 3.1140812186430085e-7im\n",
- "-1.8006326323142126, -1.8006232505919049 - 3.1140812186430085e-7im\n"
+ "-1.8006326323142126, -1.800623250593426 - 3.1140629700555206e-7im\n",
+ "-1.8006326323142126, -1.800623250593426 - 3.1140629700555206e-7im\n"
]
},
{
"data": {
+ "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlgAAAGQCAIAAAD9V4nPAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nOzdd3gUVdsG8Gdmtu8mm94TktATakIvoQgoSKQjKEVQQRHRgEoRqSIgisDLJwoIigioCEgPRUBQeu8kIaQS0stm+858fyzEEEJPsrB7/y4vr9nZszPPCZPce6YygiAQAACAo2JtXQAAAIAtIQgBAMChIQgBAMChIQgBAMChIQgBAMChIQgBAMChiWyy1kuXLi1ZsqS4uLhfv35du3a9t0FhYeHChQuvXbsWGRk5atQoiURinX/o0KFVq1axLDt8+PBmzZpVbdUAAGCHbDAiTEtLa926tZeXV7t27YYMGbJ169Z72/Tq1evUqVOvvPLKhg0bRo8ebZ159OjRl19+uXHjxmFhYZ07d75w4ULVFg4AAHaIqfoL6qdOnXr58uXffvuNiJYsWbJu3boDBw6UbnDy5MmOHTveunVLJpOlpqbWqFEjOTnZy8trwIABNWvWnDlzJhF9+OGHxcXFy5Ytq+LiAQDAzthgRHj48OH27dtbp9u3b3/kyJEyYXz48OEWLVrIZDIiCggICAoKOnHiBBEdOXKkXbt2JR88fPhwldYNAAD2yAbHCDMyMtzd3a3THh4eRqMxJyfHw8Oj3AZE5OnpefPmTev8kmYeHh7WmeX67LPP9u/f7+npaX3Jsuz8+fPd3NwqvC82YbFYOI6zdRWVxb57x/M8y9rnGWrWr7MMw9i6kEph35ulffdOJpOJRA9JOhsEoUwmMxqN1mmDwUBEcrm8TAOTyVTy0mAwKBQKIpJKpaU/aJ1Zrtzc3ICAgF69epUs0MfH56E/i+dFUVHRA/r+vLPj3gmCoNVq7bV3er2e4zixWGzrQiqFHW+WZO+9e5SvnjbIhoCAgOTkZOt0cnKyq6urUqks02DLli3WaUEQUlJSAgICSj5oPVk0OTnZ39//fqtwcnLy9/fv379/ZfXBpliWtddRBdl17wRBsOPesXfYupBKYcddI3vv3aOwQed79+7922+/Wcd2q1ev7tOnj3X+zp07ExMTiah79+7nz5+/fPkyEe3atUskErVo0YKI+vTps3r1aiLieX7NmjUlHwQAAHhiNhgR9u/f/6effoqIiPD29k5ISNi/f791/rhx42JiYt566y0PD48ZM2a0b9++WbNmR44cWbx4sXV/y5gxY9q3b9+qVSuTySQIwsiRI6u+eAAAsDM2CEKJRBIbG3vy5EmtVtusWbOSA4SxsbFqtdo6/dFHH/Xp0ychIWHZsmU+Pj7WmR4eHqdPnz569CjHcc2aNbPjo7sAAFBlbHP+CMuyTZs2LTPTeiCwREhISEhISJk2YrG4TZs2lVscAAA4Eoc+QAoAAHbGaDTGxsYeO3bs0T9iJ1cUAACAfbBYLL+sXXftetKg/r3r1KnzWJ81GAwRUZ3TvJpyBenR4d4/LlnwKJ+y2yBcs2bNkCFD3N3dGYbhOE4kEmm1WqVSqdPprLfwNhqNLMuaTCaZTGYwGMpcwlEuvV5vvd8NAACUyMnJSU1NrVOnjlQqffqlvTl63PoUrjiw2dKerx/e9lv16tWJSG8hnZmKTIKJp3wjGXkqNpHGLJh4yjOQiSeNibRmIe74PwmujQ095hDRujmtU7Zoo0NkH9Z7yL5P+wzCxTnBxam5gQ1aEMNJZTIZx/M8yVy9BG0BSZXGohxiODNPBoNeIJYsJoW7jxNn/n3FktatWpa7wKtXr77Ub3CxIApwVe7bsr7kpJ7HNWnm7OWr1oo4buEX0/r17vEUXSzLZDINGz32wKHDdWvVXLt8celb81SgffsPfPXdymr+vjMnfVThq7h48eLx48fr1avXpEmTp1lOdnY2wzCV9BMAgDK2bN/x1kdTBP96qlsXj+7dVnJLr9KsQVVgFLQWKjZRgZE0ZkFrJo2J8o2kNQtaMxUYqchEWjPt3HXIMPEwEWXpCht9voPrMLLASFKOFCJSiRgxS65SkrCkFJNSxEhYcpGSmCUnMSlEjKerWlyQZiAio1Yu6Cc0ltV2efjdjmxw0+0qwHynpc+b04jVtH4ijfmT5rSj936n9Mv0z080fAVNb0LvrKWVb5FcTU37Un46dZtA+Td91wyZtuavchf4zZiBV1qMo6DG3L8/9nTLfvW98U9QVVrClSnjPyp6ZyMZter57ZfuOukme7wTX5VikrBU7t1Jfl66+LsTWYZun7GnNnQq3D97wbcuD/xmJmFJKXrQ9mFdV2nx8fEtew7J7rOIy7jcNPH3w7s2P1bxD7Z9Z+yQj2fmN+ynPr3u10WzOr3Q8cmWMzJm/Ma9/5LAD3ql8/xZ08ptYzab+w4dcfjkGWeF/LXe0ct/+Y1j2XkzJr/ap+f9Frvg/7778dcN4bVrLf5ypqur65PVJgjC3r17i4uLO3fubH838rD7O8s4OTnZuoqnxfP8H39sSEm/+Wrf3qVvSPJYvcszUIFRKDBR/p2JT/p1zBi6jlTuzL8/N5DmhPYZW2QSCo1UbCatmfKNQrGJeCInMTmLGYWIlCJykZJSxChE5CQmtYQUIlKIGBcJqcSkENEnr3WL6zCVAhs6r35r5aQ3X3jhBbXkMbr57tgJf2zfxfLmb2ZNG9iv96N8xD5HhMQQiSRkNpKTBxERbya5mgwaUnuTIBDLEm8hmYqcPEivIScvIiKli0arP5ld/teC7LxCUvsQkcXZ90xqPHv9Sb49ZF/MMbiFEsOSVKUXq9Zd1RUx8od/rJRiExl5slhEHGcp89aNwwmGsB5ExNduf3DFD28fsuQbHrQoI0/F5gf1QmMiE3/3rONHKawnBTa0BDY8umuh28+mez+lljD37oPgWHIu78+jjCP5nQ3w9Nyfc/p9S351c2tGDZmzsA21K93SRUplQttZTNzda1KKGF1W6i9/nSoe8xcRLZ3fwS96tLPLfzeYlYvI+sXj4J+/7izyNnx0NPP60VnfjbBMPkYW0zsTOgS1eFEivb3r27XU14h/9u6cumZfYb+fLlzek/tOzI5ffyynM4/gzdFjN13MMrkF+0ydfeafvY+yN/7RHTx0aNbC7z3dXedOmeDn51eBSwa78dpb723LkBZ7hX21pPvpAzu9vb2t47D0QsasE/KNVGAU8o1knSgoNVF62kVCagmjlvw3YSKOzEYiYiz6um7ifjUYZzHrJCalmBQicpUwChFJH+c7f8PVSwa982Ha5oxhr/fv3fWFx+3mkvlzFs2d+VjfyewzCGVpp/RGLa3/hHLTuWWD5FKx6auO8prNio9vUpLOwAn02wemwkxL9g0yaOjGSVHaWXXWxVnjR41oU/4/V7epo9+aMMBYq4Pi4rZtm3+tXftJLmE0tGnWePWklG3TxJrMdpG1N3R7wi+YRUW6e7++7Xft3Wf0p7mt33M5/8e0kX0/6Fnx/7LXopq27rMku8HLbMaVxiHeu/uXs50VGAX+npkWngrLCU3SmUl/J9C/qekTm37B4leXSz3XINS3X+hdwZd3T6gXmshy95qKzUKe3iSIbieYwEkuZpskpcK+ZHVXrt4yetYiIhLLeM8aJJKSSFqsCnhnV7bY1ffeNWZtPF3UsA+pPCxNXt01exGz/HZnRCw53fkZiFlS3RlhSzhS3vnxW/fnEBEJwsGdfxsmHiGiG5uFvt/uC2rZlUolrovk9v2q1RJiGaJSYe8kZqzLLhmml/xlkXOMjKP0tNQ+I8Zl9V/C5CWf6fv6+X/3lfPjfmoJCQmfzJirKdZ+MTEmMjKyMlYBZaSlpSUnJ0dERDzusbc8A2XqhWw9ZeuFLD3d0lG2Xth84Jhu4hEiytAV1JqxWxf5qogltZicRBJXmcUabK5SUktILWF8FSVR91/suZQ3Mvt70fT+b0ULbtU8+PzvFm1Rq5/2YoSQkJB/Yv98miU87p4J+9w12qdPn40bN3bo0OHTTz+VSCQqlcrJySklJSUiIuL06dOenp6CIGRnZxNRQkJC7dq1DQZDcHBwaGjoA5aZnJx88eLF5s2bP81TLPR6/c6dO9Vqdfv27Z/4Pv33249x4sSJrbF7WzWL6NK58xNX+GA7Y3d9+e2KID+fL6dP8vLyqsAl5+bmdh8w9HpKeoi/7471q11cXJ5sOf2GjjhwLYMs5q6RtX66zwljycnJzbv0yGs0QHn9oCgnsTBiAGc2hBecPPrXjnLbHz16tNvICbldp8uv7B7gX7xi8XzrfDNPRXcC3sST5k7oGi1UbL4932Ah7Z3pVztE5I7eRQoX1crB73/wYXCDZlQqcfONt38VC4zEC0Slwr7IJFiXXTJMLzaT0UJEpDULBp5MZ7Zr4k7SK58RETOnrcvkfUqpWMLejmGOIWcJ0Z2stabpf2+JGboz5laJScySjCM5x5TEvDWqVSJq0rRZSre5JHf2+W3k5cN7y/wbPeWu0by8vD179gQEBLRsWf5xetuyya7RFat+GT/vWyGggTr91In9O6375HmBsvWUbRCy9ZSlE27pKNtA2XohW0+3dEKW/va0nCNvBeMpIw8Z4yElLzl5yph5r7fNGPADeVRT/zh4zRfjOkW1tH6vevreGQyGzMzMgICA5/TxI/YZhBMmTHB2dp40aZKtC6kU9nG44n6evndxcXEcxz34a01mZubevXtr1qwZHh7+xx8bpFJJz549H/BHfO/ev35YtzGyfp3333nbetbxE9iyfefImPFmXujbo/u3X33xZAspV1paWkTnXpl9/sfkpYad+r+D+/4qNgtG/nYMWwQqNBLdyVprmlrPwTMLVGQS6E4eF5nIzJPOIugt/8W89a2C3Kzcb98QPthKRMzaD9VtB7rVbaESk4wjZwkpRIyYLK5SRilhZRy5Shjrfm+1hGQcKUWMs4RkHKnE5CRmZNx/I2mrrKysiHYv5tTvI087+2bHel/OmFKBP5wKUZW/dEaeUouFFA31faFF9vu7SSxn9i+p7iKRtH/Lmn/uUvKQMR4y8pAx3nK6nXYy8pQx3nKyzpeUNyo7c/bsoHc/LCgoGj5owPSJH9mkd88mBOHzx763WjvunfUxTBV7aLDEv4cPf7FoqbeH26xPPy65K2HFqt6oRWL7yYLMyfvPmH/272cUTkUmMlio0EjFZqFIbyoyswaB01so3yhYgzbfSHoLac1CgZH0Fio2UZFJ0FuoyERKEck4UksYhYiK//75Rmae8ML7JPDKL1uPXXfUVcq4SshVSqUnlA/b35+UlPTG6I9S09JGDR8SM/qdCuz7k22WgiB8/NmM3zdtDa4W9Mv3C8vcPCtTRynFQkqxkFREycVCioZSioVkDWXrBT8lE6ikM5+0KxqxnlQekp1zYtoGDR4yxBp4XEUPuuz4l+4R2ecxQgBH06ply62VvFNx3+bfJs/+Wpunn/b7qurezne/yej1AseRWPyox4c0JtJbqNAkFJvoL5P3pyuOFRNRQYZCJhWxTLJGOGugPCPlGfi82xOCiSdXKblKGFfpXRMu1gkJjX99RFzUp9St0cxVQ5s1rt+6detyV33z5s3/LV1BvKVujRCVStW+ffsnPhP4wTZt+nPpkdSiDw6lxP3TcfD7r329PllDyRohpZiSNYJKTEEqJlDJVFNRoIpp5kmBSjbYiXzkjPUg8W75zMGjXiIn72C1aPo7n0ilz+Vex+cCghAAHklQUNCqJQsramkqManE5CFjiKjeq12P79ux68umSpl0zbJFLRuXn6ZGnvIMlGcQ8ox3TaRrhYt5lGeg5Iwcqt6SiPJqdO628nxoZnN/BfkomNL/9+QM7bu8ktr4DT52AREx3jWY9A/Uzs6h/j51w8JCgvzHjHxTr9f7+fk94BF9WVlZX877atPuv7NycsVSeWiQ/9Z1P3p6egpEN4qEi3l0MU+4kCfs3XS1qHoHYlihZuv0TeOJqI0PE6hkA5VUTcXIH/bXt/MLHVMunCgoKPDw8HjSHzM8EgQhANgYwzCrv//fQ5tJWPKWk7e8ZGBUdoTUr03Etu0zdf4RnidXbln/q8SbS9fSTa2QrqXzucLOVMrQ8jcuX8x0bUACQyFNqUZLIfuG0HJQXq22J+d3Pdl0nOjwX3MWNnGr3kCck9ggIrJ+reoD+7xivcuX2Wzes2fPkpU/7z/4b5HOILAc1e5AL0RR66E5Z7a0HDnDbeg3l/MFdylT14XquzEv+DHRb7886u3hOZzIKX7f8D5dp0U89tnmYrEYKVgFEIQAYCd+WbZ46YofE5JOvfnbynphoUTU2J3K5GVeVGidn85mBjchbT5ZTJR/k9o0pMQT1KQvhXUy715EY3ff0ubTmjGpHq/s+P6TL79ZIBiKSKQgkZiRqgSlGwW3pfw0kihILCX/ekQk+Idzl9cuasmFuTJ3XTJbq17d3374deOW+h1e6N8XDxJ/diEIAcBOSCSS0e+MeHAbV1fXHxfO/njqrOTcVO2us4LKW/j+NSaih3Bmi9B+JOkKSCyl+H8p6m26tEfgxKQrJJmavGtQ3Y6CxUxX/iJnLxJ4un6UGnSj1e9R8wFul7bMmzOhhVc5x/Dq169fv379yukuVBgEIQA4lq5dOnft0pmITCZTampqdnb2hQsXCl54e+VvQ8R+yhvf97S4Bhbm5fCe1akwkwSB5E7kU5tc/OjIL1S3I6VepKwExqSlbV+4Oaveq6N5dep3YWFhtu4WPDkEIQA4KLFYbH0AuPU54R+OeY+IcnNzExMTz168vGjpyoskmAWB9Bq6dZWSTpLMmTm6tmaQ39vj3howYECZyyHg+YUgBAD4j5ubm5ubW2Rk5PAhg/R6fZ06dZKyCijzOonETE5y/q1UZ2fnhy8FnisIQgCA8slkshs3bhAuOX/e5OTkyOXyR3/Gy9PeHRUAAOAZIQjC8NcG9GjTtFW9OsuXfPuIn8KIEAAA7IFgMh7eu1Nz5eS6F2uZeSH68ymvUKq8ZkOnTq8++IMIQgAAeIbk5+fPmjo5OSHhtbdG9OjZ694GlsJcS16mJT/LnJdlyb1lycsy52Va8rMEgy6rgFeQhYg4luE4sbzNK4rg2g9dI4IQAACeISOHvN7alN7dSzVlwoderLlRkI8556Y5J8OcfZMvzDHnZTEcxzm7sc7uIg9fkbuPNLQe6+zGqd1Ebj7eFst3nTt+/G/yLa2x95BhTvVbPMoaEYQAAGBLglFvzko3Z6ebs9LM2elXTx39um8TIuoR5Hzgl6V1enblXL0k1eooGkVxrp6ciycjvu+j0EQi0fa/Dpw+fdrV1TUkJOQRC0AQAgBAFRHMJktBtiUnw3QzyXQr2Zx905KTYSnM4ZzdOXcfkYev2Dekdlj4uss3G3s5bUrWLJz7pWvjxo+1CoZhIiIiHu8jdvk8ws6dO+/du1cul0+bNs3d3d3V1VWtVqekpLRs2fLYsWPWJ9TfunXLbDZfvXo1LCyM47hatWo1b978AcuMi4u7ePFiixYtnvJhb0lJSUql8mlupGvfZ3Lbce8q9XmENveUT6h/xtnxZklP3bvMzMwPRrwZH3ft5Z69p37+RclD6gWL2ZKf9aDM8w4S+1bj3H1Ert5U6lkfBQUF8z6fcT3+2pAR777UtdvTdu8R2GcQ/j2s46htJxt6Ox9Ny4uu5Xsqs0jKMU18XVaeSXozImTLtYxabqqzGXkWnm8X5LH3RtagBoGn8kwvDhnx4cfjy13gxj/++HLSuChf5x3J+eu27nyy2ykJgjCob6/suEsFelOfN978eNLkJ+tduVstz/PTP524e8e2xpFN5y78n0qlerKFP5ggCNevX3dzc6uk57eRXf/FQRA+v+x4s6Sn7t1rvV55UUhvX839s0MJnV6Ojg4Ptu7h5IsLRe4+Ik9/kae/yMNP5OEn8vTjXDyJeeYerGifu0YDncQSjnm9XiBDNK5F9ZfXHd7yaos9iVkDwvxGNgrcciX13caBs//ReSql1VzkMc2r963rZ+aFPou/GexpLneBi+cuX9Iu1EMhaegm/X7sW9OGvXq/f0tGLGVE5f8tOJuYYow/90PHWrwgRC/53/DIUIlIxMpVd39cwoilpeewirINLEaz8c6FogzDMnIFEf3w46pb+zb93Lba+qunJo/9YP6ChYxIxEhkj/YDeyQmk+mVLi+I8jPTi3TvT5j8xltvV+DCiejq1avHjx+vV69eo0aNKnbJAFCxLAU5pps3TGnXTemJ104c/qJXBMcwHfydr1yN6/NiJ0XjdiIPP871Wcy8ctlnEAok6M28mGMS8rRmXlCKRUkFOk+F9HxmERFZeIEYytQaik2Wlv5ul3OKiOhGgdbF3UPRKKrcBbr5bk8p1HkoJKlFBs/aDSVBteg+I2nBZBDMpvLfYkXWwT/DMAwJvFbDi0TmnIy7P24UTIbSc3itpszyLQa9kbv9YDNB4AWdlohObDvSLchFyrHdq3u8u2NLxucFgtksGPXWZoxUzlg/wnKs9E6ISkoym2Hkt4cpjFjK3PlSz8pV1u2Y4USMRBZ74lyoLmNiu1Cjhe81Y0q/MH8qldyMTMGwLLEiRiojIlYiI5GIiLEmPSOWPOD4NhHt3hX72eiR3QJVC5MK5y39sX3Hjg9o/ACJiYn/++pLTsSN+Wh8YGDg/ZoJgpCWlubh4SGTVeR3BQC7JJiMpowkU3rinf+uE8OI/UIlfiHS2o279eoz7ci+dj6q7y7nfLtumTIy0tb1Pjb73DU6Oir815PxLMOwLKNQKCUyhUTMMYLg7OKan5dr4QWGYfQ6rUFvsPAWjuNUCoW7p9f3q9eGh4eXu8CEhIRBvXsYiou8A4J+37rjyXY8CoIwoGe0JikuX2+Kfm3IhCnTnqx35e7H2Llj+7xx740M89x8oyCi3xtjx0+8a9UGnWCxEBHxFt6gvT3TeCezBUHQF9+eaTIIpttBzus01rwXzCbBZNh55NSBrRsnNA8xWfgem84e+PpTKpXcgl4r8DzxZsGgJyLeoCOLhQSe1xffaWYkIlamJJYhVsRK5UTESGUMJyKGGb1m17thHjXdVJezi1Zc13z74ZuMSMyIpbfHtZyIlciI41ipgliWlSmIYRiZkmFuRzgrdyIiEytq0TB8XJibief/L05z9PwljivnUajFxcXdOrRTGIvSivQdury4Z+cOhpiPJ3829M37jnHj4+O3bdlSJyzsxRdffJx/q7tg1+jzy9F2jVpyM003E03picb066b0G5bcDJFXgNgvROwbIvYLEfuFcs7/HRzhef63X9ddOnvmlb79mzRpUuXlVwD7DMIJEyY4OztPmjSpYhfL8zzLPu1N6eLj41Uq1dOccXO/38ndu3bFbt4U2ar1gIGvMZWwR8JoNEZ37iDX5KQV6t/9eMLwEe88wUJ4nYYEEiwm62hVMOgFi5kE4eNJk8Py41+u4fVn3K0EdcjM0W9Z09c6rhUsZsGoJ4uZN+iIt/B6HQm8oNfeGRALvE5DRFdTM749kfhNl/pE9O7O8591bx3s600iCSOWsFI5cRwrVzEs99OB47kXjr/dKOhajmZ07PkdA1vwvNBz84UDu3condSMRMrKlIxYUrJj+cqVKwNe7jK0hss/mbrIXq9PfNJvMIIgXLhwQavVNmnSpNyEfq4hCJ99Wq32+//7X1py8pC3RzZo0KBkflFutrQo27qf03Tzhik9kZFIxb7BYv9Qa+yJvAIYzj53H1rZc98q3NOnIBHVqFHj6RdSrs5dunTu0qWSFk5EEolk5/5DcXFx7u7u7u7uT7aQMsdES3y24NvBfXt+t/WKl5//utU/KJ/oZBx1cfGHDcOv5WjMgpBk5OrHfMnxJsFkJLOR1+uIt/DaIoHnLcwpGccSkc5s8VVJOYbhOMZdLKSuWeQpEwkmI6/TCGajYDSwMgWJJL8cixtZxy26pk+v2kKfH757p4YTI5ExUjmrUDFSOSORsVIZK1cxUvntlzIlK1PQPVvLF9Om7Pn9lwBneRqj3HngkETyoH3FABVu2IB+dYtTItXSN3ptXvvtQl/SmdITTWnXzQXZEp9qYr8QsV+IvGFrsV8oq3SsJ2wgCOExMAxTq1atyliyh4fHjv2H6Om+eiuVyh/W/j5n6mSO41ZvWiL3Dy632bDw1i9Gtb5wOOVKjoaRu80+kaq3CK61G9ab9kOZlry+WDAZ665bt3/5gu416Xxmoa+fv6RGA8GgEww6XqsR8jJ5o0Ew6HhdsaDXCkYdb9QLei2vK2ZEYkYqZ6UKRq5gJXKSyNctX7KtbyTLMF8cSdyxeHbnqLasQsXKVYxcdXvigYdRH8xisRw/ftzV1bV27YffUAocjWAyGpOvXjtz8uvejYjoaq5m3++rB/Z6RdEoStRtiF6mdlKrbV2jLSEIwa40adp0/fbYB7dxdXU9dPLMpUuX/P39XV1dd+3aJZFIOpZ3eg4rU5JMOeDNkcdPn31l+y4fX78la/9QBgc/SiWC0SAYdbxBJ+iKeYNOMOgFbrmJF6QcU2g0izX5+ssneF0Rr9XwumJBV8TriomIlSsZuROrULFyJSt3YuVKVlGSlHe/lP13rNFoNL7Urq2/oLmpMbR5pe+0L+Y8xo8M7JRgNBhvXDYknDfEnzOmxot9qzmrlMfS8mp7qP7ONHyzcKq6YUNrS31RkW1LtTkEITgiiURScpFGt24PuWKXZdlvvv3ucVfBSKSMRMqqXErmTPp8do8pk13kkuoNIrt+Nv/e47iCychri0pykddqeJ2G12osBTnmjKTbL3UaXqcRtBreoGPvDCUPJWfXMOdMaV1DIIpe/eO43l0kLp6sSs2p1KzKob/pOxrBoDNcv2hIOG9IOG9KT5T4h0qrN3DqPEAaEsZI5T/1/HDih++nHU99Z+K0hndSEAhBCFBlBg4e2jW6B8MwLi4u5TZgxBJO7c6pH+34K8+X5KLbP/8UHjxBRCYLb7ZYjBeP6otyeU2BRVMg6IpZlZpVqTln99vRqHbnlGpWpebUt+cQ+6hn7hw9elSn07Vt29b+TvZ5fvFajTHxgiH+vCHhvOlWikgmyXIAACAASURBVCSolrR6ffXLb0iq1Smzsz04OHjtpi22qvNZhiAEqDoSiaTCLp9gWVbpbD2pod2rNX/eEttvx1GN0Tzxi3nuQ98o3ZDXaiyFObxWwxfmWgpyeJ3GkJlqKbgzJz+LWJZzdmed3ViFilO7c85urFzFqe/McfFkZQoienPQQM3VM0oRO1viuv2vA8hCG+I1Bcaky4brlwzXTpsyU8V+odLQcHX3YZKQ8Kc50uywEIQAzz2GYX74ZV1hYaFMJrv3ZFRWoSpzf6IyeE2+RVPAawosBTm8poDXFJjSr+uvnb49p7iA4UTFUqdLB/f+3iuCiD7+O/7EhtURrdpy7j7WjITKcPbMmZGDX9NpNU1btlm6arWgyTcknDcmnDfEn7cU5EhCw6XV67v0fU8SWPPRx/RQLgQhgJ1wdn7CU95ZlUvpY5n34vVaXXaG5sfdZl7gWCajSCeKP5V7/bA59xYjEovcfDh3b1J7iNx9pJ7+Incfzs37fjcahEf34cg35zfxDFIHzfr3+IrXO3cLdZeG1pPWqK9s0VXsF3Lv9TnwxBCEAPAQrEyhDAiNmTSl+5xZHEO9Bw6OnHj7xFReU2DOzTDnZBhupZrTE00Xj5pzb1nyMlmls8jdh3P3Ebn7ity8OTcfkbv3A264XFRUtG/fvsDAwMaP+cwduyQYdLoLh/PTU/xbNCai6i4ybXhbv0+nPi+37nzuIAgB4JEMH/nO0Lfe5nm+9O1jWJVaolJLgmqzd99ZhtdqzDk3LTkZ5pybhhuXzSf+Kv0IHk7txjm7i9x9RB6+nLtPgSDu1LpFB2/55TxdVL9BT3zvnuedYDYZrp7UnjmoP39EHFSrf4/oUfsONHaXbUkp3r7ybaRg5UEQAsCj4jjuEc+RYRUqiaImBdYsPVMwmyy5t8y5tyw5GebcDP2105bDO825GevPXn/FR/x2Y39eEF5Zufzj0e+K3LwrpwfPJJ7XXzutO7Vfd+GI2D9UEdHepdc7rMJpItGxY8cSEhI+7Nz5aZ5gCg+FIASAKsKIxCKvAJFXQJn5IRv+ODf3UyLK1ZnEgjlr4VjBZBQH1JAE1pQE1BAH1hS5+9jheEgQjDcuaU8d0J05yLl5KyLaOb/8RpmLZ5o1a9asWTNbFeg4EIQAYGMv9+y1cf3v3bf8y3Di/636zTcqitdqTBlJppQ43aWjhbGrLfnZIt8QSWANSWBNcWBNsXfQc52Lpowk3ZmD2hN7GZFY3ijKc8xXIk9/Wxfl0BCEAGBjLMuuWLOu9NNdWIVKGhouDb39WDRepzHdTDKlxOmvnS7a+1vZXPQKtJ5CmZOT8/lnk1ISEwePfLdHz142608pc2ZO2/jbr6HVa8xfstRTyupOHyg+tlswGuT1W7q+9lFJB8G2EIQA8Ex4wNNdWPnduagtMqXGG1PidJeOFcausRTlif1CJYE1hn39Q1cXSx9vp08njq0WHFJyFz1b2bJ589lNa37tEPJvasaILq2X9mohbxzlNnCcpFrt53pEa38QhADwnGEVTtJajaW1bl9oweuLTSnxxtT45OTk6KYRRNTNX/nPbz81qBFi21utXjx+OMpHIWKZtkHuX5/P8p22Ghf/PZsQhADwfGNlSmnNhtKaDWs1+mHd5bQGnqo/E/PmdWYzvnhT5O4rC28uC28uCahRlYMwS35W0V/rm2ScnHD+plzE/XtL27l7D6TgMwtBCAB2YtnqtV/PnrU6/trM76e37tSJeN5w47L+4tHc1fP44kJZ3Uh5eAtZ3SY//7Jm/pwvWIbtO2TogIGvhYaGVmANlvzsor9+157cp4js0GH+r99fu77p17VRfRsMfO31ClwLVCxGEARb11DxJkyY4OzsPGnSJFsXUime5tG1zz477p0gCFqttsJuuv2M0d99Qf2zxnwrRXfxqP7SsZvXLg7feOz7l+oN23xKxDJZOpOnj+833//QqlWrBxR/v81Sq9XGx8eHhoaqVCpLQU7R7nXa0weULV5y6tj3OXrIux3/0j0ijAgBwP6JvAOdvAOdOvZNP3G0xj+D9yRmNfByFnPsGw2DRm47M2XYq4Vi1eoNf9auXVsketS/ihs3bvw05v1G3uqL2UWLRw2unhunbN7FZ+IyPAPyuYMgBAAHEt4oMpmXG9Jz83X6PnX8lp66MbtjWBM/17e2nO7zQltOppw+d17fVwc8eCHp6ekdmkUyhuIZ7eq0CnC7mFW0eMP2tbH7Hnzvcnhm4eAtADgQkUi055/Dr3w0vUDuNu9IQkqhXsyyJ9LznKSiPQOa/fJizY9Gj+rcuvnCr+bdbwnbtm1rXq+uWjC0DXQzWXgiMll4qW8wUvD5hREhADgWhUIxbPjwYcOHX7hw4eeffvxo3RpfpbiaTExEC44mjKjv272mcsbq73MKCgYNHlKrVq3Sn505ZfLq7/+vhrMkR2cc1jDone1nvzt5I1+s/G3rChv1BioARoQA4KDq1as3d95XJ6/Effnz+niRy2fHUvcl5QwI9yei86m3bu5cN6LnS9//3/+sjQVBGNAzetmiBVH+6ghfFzPPD918qsBo9mzQ/OTluLp169q0K/BUEIQA4NCUSmXTpk33Hznx5pffde/b/+uTKUtPJw2qHxjTJFik13wxeWLHVs2Tk5PbtWhWeOmEl0LcP8x/b2KWhOM0JJr3/coNW7dJJBJbdwKeCoIQAIBEIlGbNm0WLlla77VR1xT+mQZ+ycnE5n5uxQbTqVOn2zSJqGnMcpZw41vV/CD2fLGZD2wQeSP9Vv/+/W1dOFQABCEAwG0ikejd0e9v2733plvIjhs5P59P1vMCxzKNvJ2jgty1Rsvy00nEssPej9mya69UKrV1vVAxbBOE165d69GjR926dQcNGpSZmXlvA51OFxMTExYW9sILLxw8eNA6s6CgYOHChUOGDOncubNWq63akgHAUYjF4vVbd/zfqnV5epPOaFaIuI4hHvMOx4V7OicW6D6cPnvS1Om2rhEqkg2CkOf56OjoRo0abd++XS6XDx48+N42kydPPnv27J9//jl06NDo6Ojs7Gwiys7OPnXqVGho6J49e8xmc5UXDgAO5MWXXnr7vfcZhvRmy7oLaR2CPdddTG3cvsuId961dWlQwWxwi7Xdu3e/8cYbqampDMMUFBR4e3tfvHixevXqJQ2MRqO3t/fOnTubN29ORF26dOnatWtMTIz13bS0tICAgIKCAmfn+97BCLdYe37Zce9wi7XnkUQicREzRGTiBY6hbK3R1hVVPDv+pXtENhgRXrhwISIigmEYIlKr1TVq1Lhw4ULpBqmpqQUFBZGRkdaXTZo0KdMAAKBqGI1Gk0SRozcXmXi7TEGgyrugfufOnffODA8PDwwMzMrKUqv/uxefq6trmcOEWVlZSqWy5I5/Li4u586de6y1nz179ujRo8uWLbO+ZFl2165dXl5ej9eHZ5VGo7F1CZXIjnsnCIJOp+N53taFVAp7HRESUXJyskajUalURUVFtq6lUtjxLx0RyWSyh26WlRKEgiAsXrz43vmjRo0KDAx0cXG5fPlyycyioiJXV9fSzVxcXKx/L6xPrNZoNG5ubo9VQL169cLDw0eNGmV9KZPJ/Pz8HrsbzzD73o9hr70TBIHjOHvdNSoWi+01CK3sdbO0su/ePVSlBCHDMFu3br3fuyEhIT/++KN12mg0JiYmlnkemL+/P8dxCQkJNWvWJKKrV6/WqVPnsQrgOM7Z2bliHzMGAAB2yQbHCLt3737r1q0dO3YQ0fLlywMCAiIiIoho586d1oBUqVQ9e/acP38+EV25cmXHjh2vv377mZZ5eXkFBQVElJ+fn5+fX/XFAwCAnbHBTbflcvnPP/88fPhwIlKpVGvXrrXOP3HixOXLl9944w0i+vrrr/v27evl5WU2m2fPnm29763ZbLaeXOrq6tqoUSOZTJaenl719QMAgD2xzdMnunXrlpaWVlhY6OLy34NLJk+eXDIdEBBw5MiRoqKi0sc5RSJRbm5uVdcKAAB2zWaPYWJZtnQKlsvBj98CAEAVwL1GAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAoSEIAQDAodkyCDUazYMb6HQ6i8VSZqbJZDKbzZVWFAAAOBbbBOHRo0dr1KgRHBxcrVq1AwcO3NsgPz+/W7dufn5+bm5uX3zxhXXmnj17IiIinJ2d1Wp1ly5dUlNTq7ZqAACwQzYIQp7nBw0a9Mknn2RnZ8+ZM+e1114zmUxl2syYMUMsFmdlZZ07d27BggVHjhwhIkEQvv7666KiopycHBcXl1GjRlV98QAAYGdsEIT//vtvXl7em2++SUQDBw4UiUR79uwp02bVqlVjx44ViUTVqlV79dVXf/rpJyLq3Llzhw4dRCKRTCYbOHDgxYsXq754AACwMzYIwoSEhNq1a3McZ31Zp06dhISE0g3y8/NzcnLCwsKsL+vWrXv9+vUyC9m4cWO7du3utwqe53U6Xd4dBQUFFdoDAACwH6JKWu4333xTXFxcZmZUVFRUVFRBQYFCoSiZ6eTklJeXV7pZfn4+EZW0UalUubm5pRusWLFi3759J0+evN/az58/v3///sWLF1tfymSyQ4cOeXl5PUWHniEPPcnouWbHvRMEQafT8Txv60IqhV6v5zhOLBbbupBKYcebJdl772Qy2UM3y8oKQrPZfO+5ndY/AZ6enqWHaHl5eWUiyvqyoKBAqVRaG3h7e5e8++uvv3722Wd79+59QLA1bNiwbdu2kyZNqoiuPIucnJxsXUIlstfeCYLAcZx1q7Y/YrHYjoOQ7HeztLLv3j1UZQXhxx9/fL+3wsLCLl26pNfrZTKZxWI5d+7c9OnTSzdQKBTVqlU7deqUn58fEZ0+fbpkN+nGjRtjYmJiY2Pr1KlTSZUDAIBDscExwoYNG9avX3/ixIlpaWnTpk3z8fFp3bo1EW3YsKHkRNCRI0dOnz49Pj5+586dGzdutJ5Zs2vXroEDB44dO/bWrVt79uz566+/qr54AACwM5U1Inyw3377LSYmJioqql69eps2bWIYhohK71f5+OOPCwsLo6OjnZ2dV61aVbt2bSJKTU1t27ZtbGxsbGwsEUml0o4dO9qkfgAAsBuMIAj3zuV5/vjx48eOHUtPT09PT2dZNiAgICAgoGvXrkFBQVVf5eOaMGGCs7OzvR4jLCoqsuMd+nbcO0EQtFqtvR4jtO+TZex4syR7792jKDsiFARh+fLly5cv9/f3b9GiRa1atZo3b87zfG5ubmpq6ogRI0wm05gxY3r06GGTcgEAACrWXUGYkZExderU9u3b//PPPyJR+XtNNRrN8uXL33vvvblz56pUqiopEgAAoLL8l3Y5OTnff//9N998U/oiv3upVKoPP/wwMzPz888/nz59ulQqrfwiAQAAKst/Qejk5DR16tRH/JiXl9esWbPKPb4IAADwHPnv8gmJRFLmPYPB8ICo4zjufrtPAQAAnhf3vY7w5s2bDRs2HDJkiPXllStXli5dmpKSUlWFAQAAVIX7Dul0Ot3kyZN9fHysL+vUqVOjRo2ffvqpV69ebm5uVVUeAABA5bpvEIaGhu7bt6958+ZElJGRsWDBAicnp1u3brVt2xZBCAAAduOuXaPXrl3bu3evVqu1vhw6dOgPP/xgNBrff//9goKC1NTUVq1a1apVyxZ1AgAAVIq7RoRnz5598803DQZD06ZNo6Ki2rVrN3jw4JUrV9arV+/RTygFAAB4jtw1IvTx8VmxYkVqauqHH35YWFg4btw4X1/fWbNmxcXF7d27V6fT2apKAACASnLXiLB169YWi0UsFvft27dv375ElJmZuX///v3797///vvXr19v2bLl2rVrS86gAQAAeN7dFYQsy7LsXWNELy+v/v379+/fn4gyMjKOHDni4PdmBQAAO/MYV8T7+Pj07Nmz8koBAACoejZ4MC8AAMCzA0EIAAAODUEIAAAO7TGOERYXF2/bts1isbRq1apatWqVVxMAAECVeYwRYUxMzPDhw3v37n3t2rUFCxaYzebKKwsAAKBqPMaI8N13323Tpo1UKu3cuXPr1q0NBgMewwQAAM+7x0iyxo0bN27c2Dr94KfYAwAAPC8eIwi3b9++YcOG1NRUi8Xi5eXVoUOHIUOG3Ps4XwAAgOfIowbh6tWr9Xr9gAEDnJ2dBUEoKipKS0ubOXPmzJkzK7U+AACASvWoQejs7Dxo0KAyM/fv3280GjEoBACA59ejnjWalJSUlZVVeo5Gozlz5gxSEAAAnmuPOiIcMGDAhAkT9u/fb7FYLBaL2WyOjIz88ssvK7U4AACAyvaoQejp6fnDDz8QkUaj4Xne2dm5MqsCAACoIg8KwsuXL4vF4tDQ0NLPZlKpVJVfFQAAQBW5bxAmJyeHhYV17NjxxRdfZFnWyclp0KBBSqWyKosDAACobPcNwqCgoEuXLtWpU4dhGCK6dOlSdHT0jBkz2rRpU4XlAQAAVK4HnTVat25dawoSUVhY2NKlS3v27FlcXFwlhQEAAFSFx7jptr+/f0FBwdWrVyuvGgAAgCr28LNGv/rqq6SkpJdeeumnn36KiIgIDw+vgrIAAACqxsODMCYm5sqVK8eOHbt+/frmzZulUmkVlAUAAFA1Hh6EHMeFh4eHh4f7+PhER0cfO3aM47gqqAwAAKAK3PcYYXFx8Zw5c4xGY8mcF1988dy5c+fOnauSwgAAAKrCfUeESqWyS5cukyZNql+/fpcuXZydnVesWKFUKgMDA6uyPgAAgEr1oF2jERERERERmZmZ//zzz40bN0Qi0enTpz08PKqsOAAAgMpWNgjNZrNIdNdMLy+vXr16ZWdnq1QqmUxWhbUBAABUurLHCLt16zZ//nzr9LFjxzZv3qzRaIiI5/lvv/12y5YtVV0gAABAZSobhAsXLuzSpQsRJSUlffnll7Nnz/bz8xs2bNjly5djYmJkMllKSoot6gQAAKgUZXeN1q1b1zohl8sHDRrUs2fPhISEVatWDRs2zGKxBAUF/frrr1VeJAAAQGW57+UTXl5eLi4u8+fP53l++vTp169f3759++rVq/38/KqyPgAAgEr1oLNG27dv365dO71eb32Jm6sBAID9+W9EmJ6efvz48TJvMwwjl8vL/eTu3bvxJAoAAHje/ReEfn5+p0+fXrx4sclkevBntFrtzJkzCwsL8ZxeAAB43t21a3TEiBF///13dHR027ZtO3To0KRJE4lEUvKuVqs9fvz4rl27zp079/nnnzds2LDKqwUAAKhgZY8RRkVFtWzZMjY2dtmyZa+//rper1er1TzPazQauVzeqlWrgQMHzpw5k2Uf40GGAAAAz6xyTpYRi8Xdu3fv3r07EVkslszMTJZlvb29q7w2AACASveQxzBxHOfr61s1pQAAAFS9cvZwGgyGuLi4pKSkkjlnzpxJS0uzWCxVWBgAAEBVuCsIL1265O/v//LLLx84cKD003fNZvOqVauCg4NDQkKstx4FAACwD3ftGtVqtb6+vrGxsWWeQd+kSZMmTZr06dOnXr16GBcCAIA9KbtrdODAgWVSsEStWrUiIyMrvyQAAICqUzYIfXx8HtA6ICCgMosBAACoamWDUKFQlEzn5uauXLmy9L5QPJgXAADsTNkgZBjGOrFz58769eu/9957Y8aM0el0t1vjOnoAALAv5QRbcXHxqFGjunXr5uvre+LEiUmTJsXExOTk5FR9cQAAAJWtbBAeO3ascePGy5Yt+/TTTw8fPhwWFubv7z937tzJkycnJiZW4Irj4uL+/vvvoqKi+zXIycnZv39/SkpK6ZnJyckHDx48deqU0WiswGIAAMBhlQ3C2bNnMwxz6NChmTNnisVi60y1Wr1gwYKFCxdevny5QtY6ZsyYDh06zJgxo2bNmkePHr23wZYtW2rXrj179uzIyMh58+ZZZ86fP79t27ZTp0594403atSoceHChQopBgAAHJpQSnx8fExMTHFxsVAei8UyZcoUrVZb7ruP7vTp066urhkZGYIgWLPt3hWFhIT8/vvvgiBcuXJFLpffvHlTEITSq3777bf79+9/v1WMHz9+1qxZT1nnM6uwsNDWJVQiO+6d9eb1tq6isuh0OqPRaOsqKosdb5aCvffuUdw1Iqxevfr8+fNLnzhaGsuy06dPv99zeh/d+vXrX3rpJetdvIcOHXro0KGMjIzSDU6ePJmTk9OzZ08iql27dmRk5ObNm4mo9Ko9PDxEoofcKBUAAOChbJAlKSkpwcHB1mk3NzcnJ6eUlJTS1y+mpKQEBASU5FxwcHDJkcL4+Ph58+ZlZmbm5eX99NNP91uFRqPJzMzcvXu39SXLslFRUfe7UcBzh+d5nudtXUVlsePeCYJgx73jeZ5hGDvunb12jey9d49ysUOlBGF6evqoUaPunT9jxowGDRro9frSz/uVyWQll2dYlWkglUq1Wq11WqVSRUZGpqenr1y58tixY9WqVbtfAWfPnr1x44b1pUQiqVmzppub29N161mh0+nsJtTvZce9EwRBp9OVXKFkZ/R6PcdxJScW2Bk73izJ3nsnk8keuvuwUoJQrVa/9dZb98738/MjIh8fn5KLMSwWS25ubpknPZVuQETZ2dl16tQpeWvEiBFEFB4ePn78+H79+pVbQK1atZo0aTJp0qSK6M0zRxAElUpl6yoqix33ThAElmWVSqWtC6kUIpHIjoPQjjdLsvfePYpKCUKlUml9rm+5mjdvPnfuXOv04cOH3d3dS/aUWjVs2DArKysxMTEkJMRsNv/777/jxo0rsxA7/v4CAABVyQbHCPv06fPpp59+8MEH7dq1mzJlypgxY6zfIocMGRIYGDhr1ix3d/ehQ4cOGTJk3Lhx69evDw0Nbdu2LRFNmjTJ29s7ICAgOTl53rx5Y8eOrfriAQDAztjglmlSqfTgwYMikWj9+vUxMTETJ060zu/WrVubNm2s04sWLerVq9e6deuCg4N37NhhndmpU6cbN278+uuviYmJP/7440cffVT1xQMAgJ1hBEGwdQ0Vb8KECc7OzvZ6jLCoqMjJycnWVVQWO+6dIAhardZejxHa98kydrxZkr337lHgJtoAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQEIQAAODQRDZZa0ZGxsaNGwVB6N27t4+Pz70NTCbThg0b0tLSoqKimjRpUubd2NhYuVweFRVVJcUCAIA9s8GIMCkpqUGDBidPnjxz5kyDBg1u3Lhxb5uXX3558eLFWVlZL7/88urVq0u/FRsb27Nnz88//7yKygUAALtmgxHhokWLunbtunz5ciKyWCwLFy785ptvSjc4cODAhQsXEhIS5HJ5y5Ytx44d+9prr7EsS0SFhYUfffTR6NGjz549W/WVAwCA/bHBiDA2NvaVV16xTkdHR8fGxt7boFOnTnK5nIheeuml5OTk+Ph461vjx48fNWpUQEBAVRYMAAB2zAYjwvT0dF9fX+u0n59fenr6vQ38/Pys0xKJxN3dPT09vVatWvv37z979uz//d///e9//3vwKm7evHnq1ClBEKwvWZYdOXKkUqms0H7YjMFgkEgktq6isthx7wRBMBgMIpFtDsxXNoPBwHEcz/O2LqRS2PFmSfbeO7FYbN2h+ACV8ju5YsWKKVOmlF2TSFRyOLAkogRBYBimTEuGYUoalLTRarXvvffe77///tAuWen1+vz8fOu0VCp9/E4AAIBDqJQgHDBgQLdu3crMLAk8Pz+/W7duWaczMjJKRoclfH1909LSrNNGozE3N9fX13fHjh15eXnTpk0jori4uJs3bw4aNKjMeTSll1C7du1JkyZVUIeeLUaj0Y6j3Y57JwiCxWKx495xHCcWi21dSKWw482S7L13j6JSglChUCgUivu926VLly1btvTu3ZuItm7d2qVLF+v8uLg4f39/hULRpUuX119/Xa/Xy2SyXbt2BQYG1qhRQ61Wr1q1ytpy48aNR44cGT58eGUUDwAADsUGhyvGjBnTrFmzESNGMAyzadOmY8eOWec3btx448aNnTt3bteuXVhY2IsvvtimTZsffvhh3rx5LMt6e3t7e3tbW168eDEuLq5jx45VXzwAANgZOgFMLgAAGZ1JREFUGwRhcHDwuXPnNmzYQETTp08vuaB+7dq1jRo1IiKGYbZv3/7HH3+kpaVt3br13gvqu3Xr1rhx4youGwAA7JJtTmDz8fEZNWpUmZnR0dEl02KxeMCAAff7eM2aNWvWrFlZxQEAgCPBvUYBAMChIQgBAMChIQgBAMChIQgBAMChIQgBAMChIQgBAMChIQgBAMChIQgBAMChIQgBAMChIQgBAMChIQgBAMCh2efDsgEA7IBer9fpdJW9Fo1GYzabK3stT49lWbVaXRlLRhACADyjmjdvnpSUxLLYdUdEpNPp5syZ88EHH1T4khGEAADPqOLi4pMnT1avXt3WhTwTJk+erNFoKmPJ+KIBAAAODUEIAAAODUEIAAAODUEIAAAODUEIAAAODUEIAGC3NBrN+fPn8/LybF3IMw1BCABgh9asWVO/fn1nZ+cGDRq4ubkFBwcvWLDAYrE82dISExMvX75csRU+OxCEAAD2Zty4ca+//vqFCxcEVkSeISSSJiUlxcTE9OnT58luIrN27dpFixZVeJ3PCFxQDwBgV7Zs2TJ//nziRNRvLkW9TWIZ8RY6+Qf9/N6ff/65YMGCjz766HGXGRMTw/N8ZVT7LMCIEADArixcuJCIKHoyvfA+iWVERCxHTfvTkO+IaNGiRYIgPO4yN2zYsHTpUiK6dOnSsGHDFi1aVKtWrSZNmhw+fHjr1q0RERHh4eG///67tfHmzZvbt28fHBzcsmXLTZs2WWcKgjB79uywsLAmTZqsXbs2Ojq6YnpbETAiBAB4bhSayPKwgdm/hw8TEbUaWvaNiJ4kd05JSbmQkBIQGPTghbhK73qZlJSUkpJCRAUFBWvWrAkJCTl69OiPP/746quvduzYMTY29uzZs3369HnppZecnJzc3NyWLl0aHBx86tSpXr16NWjQIDQ09Oeff169evW2bducnJwGDx588ODBx+p4pUIQAgA8HzQmCltv1pofOJ4zG3RaLRGRs1fZt1iOVB6kK2zzaxYX5PuAZYhZ+qubKNyVKfddlUo1efJklmWHDRs2duzYqVOnenp6durUyd3dPS4uLiIiok2bNleuXNm+fbvJZAoODv73339DQ0PXrVsXExMTEhJCRBMmTEAQAgDAY1OJKXXgQ/9oi13HuuTn51NWIvnUuusdk44KbhLRtbcCvL3FT1yGp6en9YEYcrmciLy8bieuQqEoLi4mogkTJmzfvr1bt26urq56vT4nJ4eIMjMz/fz8rC39/f2feO2VAccIAQDsygsvvEBEFPt12Tf++paMunr16nl7e1fe2i0Wy6JFi2JjY+fMmTN+/HiFQmGdHxQUdPXqVev0lStXKq+AJ4AgBACwKxMnThSJRHRoJa0YTmkXyKSnzAT6YxJtmExEU6ZMqdS1cxynVqv37dtXXFz8/fffnzx50jp/1KhRX3755aZNm/bs2fPFF19Uag2PC0EIAGBXIiMjly9fLpFI6PBqmhZBo5zp07q08ytG4KdMmdKvX78nWGZISEjdunWJSK1Wt2zZ0jqTZdlOnTpxHGd92aJFCxcXFyL65ZdfFi5c2Lhx4/j4+M8++ywoKIiIOnXqtGzZsrVr165du/aTTz5xd3evmN5WBBwjBACwN0OHDm3cuPHcuXN3796dlZXl7Ozcrl27sWPHtm/f/skWOHDgQOtEWFjYypUrrdNisXj37t0lbZYvX26d6Nix49GjR8sswWKxdOrUqXv37kQ0ZcqUZs2aPVkllQFBCABghxo0aPDLL78QEc/z1nNbbEun09WpUycsLCwrK8tsNpdcX/gsQBACANit7Ozs+Pj4gICAgIAA21aiUqni4+Nv3LihVCoDAwNtW0wZtv+aAAAAFW758uWhoaGenp4tW7YMDAz09vaeNWuWyWR6sqVdu3bt7NmzT1mSTCarU6fOs5aChCAEALA/77777ttvv52YmCgXcbXdVU5SUWZm5uTJk6Ojo41G4xMscP369d99912F11lGTExMo0aN3NzcqnjHKXaNAgDYlT/++OO7774Tc+zn7ev2D/MXsYxAtOd65tjdF2JjY7/++uuJEyc+7jLHjx//BHcofVxhYWE9evR49913nyytnxhGhAAAdmXx4sVE9HHLGq/VCxCxDBExRJ1Dvf6/vTuNauIK+wB+Q8CAQZBEkBCFoigBioJBQTwU6CuLqICKWKtVPBRwr6221uW49u3RilVPe1qlLhREaLWtULCAGBYVEC1UEaiCC0KQQFjEEEJgMu+HeZvmAEUBk2HI8/s0SW7If0gmT+bOzL1Hfd4mHh1ESYuPjz969ChCqLS0dMmSJYcPH544caKdnV1ubu6lS5d4PJ61tXVcXBzR+NKlSzNnzuRwOE5OTklJScSdGIbt2bPH2tra0dHx7NmzHh4evV8lIiLCy8tr1KhRg131QYI9QgAAoAystRHHXjGhYGFBAUJoMc+ix/3/Y206hqFbV1f36K/bb1n2f6COpssaj2j/jjVaV1dHDLotkUhSUlLc3NyqqqrOnTu3fPny+fPn//nnn3fv3vX39w8ODjYyMrK0tLx8+TKXy713756/v7+Li4uNjU1sbGxycnJhYSGTyVy+fHlJScmg/gFqAYUQAACoAe/saPx+J+p3Zl15Nybr7EQIsQ167lfRaTSW/qiXnd3Vp//X0JzV3yvREPvDfXrmVn0+aGxsvG3bNhqNtnz58nXr1u3evZvJZLq7u5uamlZWVvL5/FmzZpWUlAgEgo6ODg6Hk5+fb2Njc/HixS1bthCju23fvj07O3tA665WUAgBAIAaaAwD8x0/vLIZ69u05ubmJ63tU1iGqvd3dGP17TKE0LR9p01NTQcdY9y4cTQaDSGkr69P3CTuNzAwkEqlCKGPPvooPz8/ODjYxMREV1e3paUFISQWi5XDc5ubmw/61dUBjhECAMCI4uPjgxD69vaTHvfHFD/t7FY4OTkNpQq+EoZhMTExqampu3btWr9+Pe2f/lVra+uKigpiuaysTH0BBgH2CAEAYETZuXPnb7/99uvfdbJuLHLGW5NNmHUS2YX7tXF3nyGE9u/fr9ZXp9Pp48aNS01NXbhwYUJCwr1794jh2TZu3BgaGmphYWFkZHTo0KE+n1tcXNzc3CyRSEpLS1ksFp/PNzExUWtaAuwRAgDAiDJt2rT4+HgDA4MrVaLgn285nhL4JeT/ePcZTUfn0KFDgYGBg/ibtra2Tk5OCCE2m+3v70/cqaOjs3TpUl3d/9+h8vPzI7pJk5KSEhISfH19W1tbDx06NGXKFISQp6dnYmKiQCBIT0/fvXu3skNV1ZUrV2JiYvh8/oMHD2JiYp4/fz64/8BA0TRwaYjmff7550ZGRjt37iQ7iFq8fPlyzJgxZKdQlxG8djiOS6VSJpNJdhC1kMlkdDpdT2/w070OZ2R9LG1sbDIyMiZPnjyI5z548ODrr7/OyMh4/vw5m8329PTcsmWLq6vrGw/5mjAMk0qlxL/xs88+q62tvXDhwoD+wu7duw0MDHbt2vXGs0HXKAAAjEC2tranTp0iO8W/ZDIZj8ebNGmSSCQyMTH5+eefyU70LyiEAAAA1I7JZFZXV9fX1zOZTM0c+Xt9UAgBAABogq6uLumTYPQJTpYBAACg1aAQAgAA0GpQCAEAAGg1OEYIAADD1/3791tbW8lOMSw8f/580qRJ6vjLUAgBAGCYmj179sGDB9X9KgqFQkeHGr2Dymv53ywohAAAMEzFx8dr4FVG8CgWr4kavwIAAAAANYFCCAAAQKtBIQQAAKDVoBACAADQaiOzEIrF4ubmZrJTqEtubi4xDfSIlJ6eTnYEdWlpaSksLCQ7hbqUl5c/ffqU7BTqkp2d3dHRQXYKtcBxfARvdK9pZBbCsrKy+/fvk51CXfbt2zdS1669vT0sLIzsFOpy8+bNo0ePkp1CXeLi4i5fvkx2CnXZs2ePcnb1Eaatre3DDz8kOwXJRmYhBABo3oic3BRoAyiEAAAAtBoUQgAAAFqNNiJ7M3g8nlgsdnZ2JjuIWhQVFfF4PCMjI7KDvHkKhSI3N9fb25vsIGohFotra2udnJzIDqIWDx8+ZDAYVlZWZAdRi1u3btnb24/I4VcwDLt+/bqXlxfZQdRl0aJF69ev77/NyCyEWVlZtbW1w3MGyKF79uwZl8ul0+lkB1GLJ0+eWFtbk51CLTo7O5uamiwsLMgOohZNTU16enoj8vcZQqi6unrChAmw0VGRtbX15MmT+28zMgshAAAA8JrgGCEAAACtBoUQAACAVoNCCAAAQKtBIQQAAKDV6Pv27SM7w5AoFIpr167l5eUZGhqyWKw+25SVlaWlpb18+ZJy53ZjGJaZmXnz5k0TE5OxY8f2biASia5evVpSUsJgMNhstuYTDsXTp0+Tk5NFIpG1tXU/E2SXl5eXl5e/9dZbGoz2BhQXF6enp3d1dXG53D4bENeKCASCpqYmLperq0uZWbJfZ6MrLCwUCAQikcjKyooqs58TWlpaioqKurq6/muDkkqlv//+e0lJCYfDGT16tIbjDZFYLL59+7ZCoTAxMen9KIZhBQUFOTk59fX1VlZWI/Us2T7gFBcSEjJ9+vSIiAg2m52cnNy7wblz50xNTaOiomxtbdeuXav5hIOGYZifn5+Li0t4eDibzc7KyurRIDc3d+zYsYGBgStXrjQ2No6OjiYl5+BkZmayWKzw8HAXF5eAgACFQtFnM6FQaGpqymazNRxviI4cOcLlcqOioqysrPbu3du7gUQi8fLycnBwWLNmjbe3d1pamsYzDt6SJUuUG11KSkrvBmvWrJk6derHH3/s4eHh7Ozc3t6u+ZCDEx4ePmrUKGNj461bt/bZ4MWLF/b29r6+vu+//76ZmdnDhw81nHAoli5dqq+vP2bMmH379vXZgHi/wsLC+Hy+nZ2dWCzWcEKyULsQFhQUmJmZvXjxAsfxhIQEBweHHg2I3+Pp6ek4jjc0NDCZzMrKShKCDkpGRoaVlZVUKsVx/OTJk7Nnz+7RoKGhobW1lVi+evWqvr5+Z2enplMOlqur68mTJ3Ecl0qllpaWvcs8YdGiRdu2baNWIWxraxszZkxJSQmO45WVlQYGBo2NjT3afPrpp35+fl1dXWQEHJL8/HwzM7O2tjYcxxMSEt5+++0eDcRiMULoyZMnOI53dXVZWlr2WSyHp+rq6s7OzoiIiP8qhMeOHfP29iZ+t23atCk8PFyzAYfk6dOncrl82bJl/1UIq6qqiAUMw2bNmnX48GENpiMTlbosektNTfX19SWu4Q0ODv77778fP36s2qC4uLi9vd3HxwchZGpq6uHhkZaWRk7WgUtNTQ0ICDAwMEAIhYSEFBQUEF8xSqampsbGxsQyh8Pp7u7GMIyEoAPX0NBw69atkJAQhJCBgUFAQEBqamrvZhcuXGAwGIGBgRoPOCS5ublmZmbECDI2NjY8Hi8zM7NHm8TExE8++aS8vDwvL49ak2oRGx0xxkpwcHBFRcWTJ09UGzAYDCaTSaxUd3e3XC6nUKe9paXlqFGj+mmQmpq6ZMkSGo2GEAoJCenzcztsWVlZ6enp9dNAeeG5jo7O+PHj5XK5RnKRjzKHJfokFAqVw8eMHj3axMREKBROmjRJtQGHw1EeouByuUKhkISggyIUCmfMmEEss9lsfX19oVA4bty4PhsfOHBgxYoVRNUc/oRCob6+vvL7kcvl3rt3r0cbsVh88ODBnJychw8fajzgkKh+LFFfnzqZTCYUCg8dOsRgMLq6uqqqqrKysmxsbDSedDD63OhUxyUxNDS8ePFiSEiIg4NDRUXF1q1b3d3dSQr75gmFQuVBXy6X29DQIJfL+6+dVFRUVJSXl3fs2DGyg2gItfcIMQwjfpoRdHV1u7u7+2lAp9N7NBjOMAxTPcugn/D79++vqKg4fvy4pqIN1eus2oYNG7Zv3z5+/HjNRnsDXvmxJHqwnZyc/vjjj6ysLD8/v71792o85iC9zkZ3/Pjx6dOnh4aGBgcHx8TE1NXVaTymuqh+dOl0Oo7jCoWC3Ehv3OPHj0NCQr755ptXjkw2YlC7EHI4nIaGBmK5u7u790COHA6nsbFReVMkEnE4HI1GHALVtZNIJO3t7X0OU/nVV18lJSVdvXq1z9NKhycOhyOVSiUSCXGz9/tSXV2dkpKSn58fFRUVHR3d3t4eFRVVX19PRtgBU33jEEIikajHG2dsbDx69GjlMMdeXl4Ummn5lRtddnb2X3/9df78+aVLl37xxRf29vanT58mI6laqK6+SCRisVj6+vrkRnqzampq5s6du2PHjg8++IDsLJpD7ULo5eV17do14gepQCAwNzcn+pckEglxiMLZ2bmrq6u4uBgh1NHRcf36dQrNbODl5ZWVlUX83szIyODxeObm5gihly9fdnR0EG1OnDhx+vTpa9euUWvPycLCYurUqcSRM4VCkZWVRbwvxBcrQojNZsfFxfn4+MydO5fP5+vp6c2dO5fJZJKc+/W4u7s/evSopqYGIdTU1FRSUuLp6YkQkslkbW1tRJt33323srKSWK6srKTQAPFeXl4CgUC50XE4HGK/QbnR0el0DMOUh6vlcjm1Lp/oDcMw4mOJEPLy8srIyCCWMzMzR8CkDd3d3c3NzcSyUCj09vZet27dunXryE2laWSfrTMkGIbNnDkzMDDw6NGjEydO/O6774j733vvvU2bNhHL+/fvnzp16rFjx+bOnevv709e2AHr7Ox0cHAIDQ2Njo42NzePj48n7p8/f/6uXbtwHBcIBDQaLSAgIPIf9fX1pEYegB9//NHc3Dw6Onrp0qWOjo5yuRzH8YKCAoRQj0sp8vLyqHXWKI7jGzZscHZ2Pn78uJub26pVq4g7T5w4MXPmTGI5Pz/f1NT0yy+/3L9/P5vNzsvLIy/swBAbXVBQ0H9tdB0dHY6OjvPmzfv+++/XrVtnbGz86NEjUiMPwOXLlyMjI3k83rRp0yIjI4kTzu/evYsQkkgk+D/X82zevPnAgQPGxsaFhYVkRx6ApKSkyMjIyZMn8/n8yMjI3NxcHMdv3Liho6NDNHBzc7OwsFB+n8TGxpKaV3MoP/uERCKJjY19/vy5p6enr68vcWd2dra+vv7s2bOJm8nJyYWFhdbW1qtXr2YwGOSFHbDW1tbY2NjGxkZfX19irwIhlJGRwWazXVxcHj9+nJWVpdp+2bJlyvNIh7/s7OysrCwzM7OwsDAidkNDQ3JyckREhGqz+vr6zMzMVatWkRRzMBQKxU8//XT37l17e/sVK1YQFyaXlZVVVVUFBQURbcrKyi5dusRgMBYtWmRra0tq3oF55UYnlUoTExMfP35sZmYWGhpKoeMRt2/fLikpUd50dXWdPn16c3PzL7/8smbNGmLQg5qamvPnz3d2di5ZssTR0ZG8sAOWn5+v2gnv4eFhZ2cnEolSU1PDw8MRQklJScpOC4SQra2t8mtnZKN8IQQAAACGgtp99wAAAMAQQSEEAACg1aAQAgAA0GpQCAEYUQQCAVVG2gNgmIBCCAA13Lx5c/HixU5OTv2f4Pbw4UPl7DnffvvtO++8s3XrVo0EBICqoBACQA1z5sxZvHgxi8VSHeGsB5FIpDr94caNG9lsturouwCA3qAQAkAZ169f7/+6LmJqCOVNHMdv3LihJZeCATBoUAgBoIzs7Oz+x/TCMEx1yIjS0lIcx+3t7dWeDAAqo/Y0TABoD6FQ+OzZs1mzZiGEUlJS6urq7OzsVPf2ioqK+Hy+6lNycnLeeecdHR0diURy6dKlhoaG8PBwCs0OCIBmwB4hANSQk5Pj6urKYDBOnTrl4uLy66+/xsfHqza4c+dO70Lo6elZX19//vx5YrKn27dvazY1ABQAQ6wBQA0REREsFmvKlCmLFi1is9mlpaVWVlZGRkbEo3K5PDExcfXq1cr2CoXCzMzszJkzbW1tK1euxHG8qKjI1dW1n3NtANBOsEcIADXk5ORkZmbKZDJi6GdHR0dlFUQIpaWl+fn5qbYvLS2VSqU//PADi8XCcVxHR8fNzQ2qIAC9wTFCACigpqZGKBS2tLTk5uY6OjpWVFT0mJ1RJBIR01UqZWdnL1iwIDExcfPmzcnJyTExMaqPYhh29uzZ5uZmQ0NDPT09mUwmk8nCwsLMzMw0sT4ADCdQCAGggJycHHd3dwaD4evr29TUVFdXx2Awampq5syZgxASiUQ9poknnuLr60un0xcsWLBjxw6EUFJSUmhoKDFNLp1OX7hw4YwZM2pqaogL8I8cOSIWi6EQAi0EXaMAUEBOTg5x4QSO43Q6ncvlpqamurq6Eo9euXKlR7+oQqHIy8vz9vZGCOno6FhYWOA43tjYqDpZ/PXr1729vel0OnGiAJ/Pt7S01NgaATB8QCEEgAImTJgQEhKCEKLRaGfPnj137hyfzycOFiKEuru7e8w4LZPJgoKCeDweQsjHx8fNze3MmTMrVqxQbaMsrhcvXlQoFO7u7oaGhppZHQCGFThrFABqu3PnDo1G63HhxOtwcHBYu3atnp5efn5+XFycOrIBQAlwjBAAartz587atWsH+qyGhoa2trZNmzYhhDgcjhpyAUAZ0DUKAIXJ5XIDA4NBPFHZL4oQCgoKepOZAKAaKIQAUFhtbe28efMG8cRr1665u7u/8TwAUBEcIwRA6yQnJ4tEIj09PTs7Ozc3N7LjAEAyKIQAAAC0GnSNAgAA0GpQCAEAAGg1KIQAAAC0GhRCAAAAWg0KIQAAAK32f2eQKMm8ESNVAAAAAElFTkSuQmCC",
"image/svg+xml": [
"\n",
"