From 0b9a6d3f5c28ea571d352e388957aa6e5e62149c Mon Sep 17 00:00:00 2001 From: L-Nafaryus Date: Mon, 11 Oct 2021 16:30:43 +0500 Subject: [PATCH] Fix: faceCentered fillets radius --- anisotropy/samples/faceCentered.py | 30 ++++----------- playground/geometry.ipynb | 61 ++++++++++++++++++++++++------ 2 files changed, 58 insertions(+), 33 deletions(-) diff --git a/anisotropy/samples/faceCentered.py b/anisotropy/samples/faceCentered.py index 6ed70d4..fda75d5 100644 --- a/anisotropy/samples/faceCentered.py +++ b/anisotropy/samples/faceCentered.py @@ -4,6 +4,7 @@ from anisotropy.samples.structure import StructureGeometry from math import pi, sqrt +import logging class FaceCentered(StructureGeometry): @property @@ -22,34 +23,15 @@ class FaceCentered(StructureGeometry): @property def thetaMax(self): - self.Cscale = 0.8 return 0.13 @property def fillets(self): if self.direction == [1.0, 1.0, 1.0]: - #C1, C2 = 0.3, 0.2 - #Cscale = 0.9 #C1 + (C2 - C1) / (self.thetaMax - self.thetaMin) * (self.theta - self.thetaMin) - #delta = 0.thetaMin = 0.01 - thetaMax = 0.13 - - L = 1.0 - r0 = L * sqrt(2) / 4 - radius = r0 / (1 - structure["theta"]) - - C1, C2 = 0.3, 0.2 - Cf = C1 + (C2 - C1) / (thetaMax - thetaMin) * (structure["theta"] - thetaMin) - delta = 0.012 - fillets = delta - Cf * (radius - r0)012 - - return self.Cscale * (2 - 1 / (1 - self.theta)) * self.r0 #delta - Cf * (self.radius - self.r0) + return self.filletsScale * (2 * sqrt(3) / 3 - 1 / (1 - self.theta)) * self.r0 else: - C1, C2 = 0.3, 0.2 - Cf = C1 + (C2 - C1) / (self.thetaMax - self.thetaMin) * (self.theta - self.thetaMin) - delta = 0.012 - - return delta - Cf * (self.radius - self.r0) + return self.filletsScale * (2 * sqrt(3) / 3 - 1 / (1 - self.theta)) * self.r0 def build(self): ### @@ -213,12 +195,16 @@ class FaceCentered(StructureGeometry): ### # Shape ## + #poreCell = self.geo.LimitTolerance(poreCell, 1e-12) + #grains = self.geo.LimitTolerance(grains, 1e-12) + self.shape = self.geo.MakeCutList(poreCell, [grains]) self.shape = self.geo.MakeScaleTransform(self.shape, oo, 1 / scale, theName = self.name) - isValid, _ = self.isValid() + isValid, msg = self.isValid() if not isValid: + logging.warning(msg) self.heal() ### diff --git a/playground/geometry.ipynb b/playground/geometry.ipynb index 68eb32b..5a95305 100644 --- a/playground/geometry.ipynb +++ b/playground/geometry.ipynb @@ -36,7 +36,7 @@ }, { "cell_type": "code", - "execution_count": 55, + "execution_count": 10, "id": "60e986af-be50-46ac-b8e6-73e12b0175b9", "metadata": {}, "outputs": [], @@ -57,12 +57,13 @@ "C1, C2 = 0.3, 0.2\n", "Cf = C1 + (C2 - C1) / (thetaMax - thetaMin) * (theta - thetaMin)\n", "delta = 0.012\n", - "fillets = 0.9 * (2 - 1 / (1 - theta)) * r0#Cf * (r0 - radius) #(radius - r0)" + "fillets = (2 * sqrt(3) / 3 - 1 / (1 - theta)) * r0 #Cf * (r0 - radius) #(radius - r0)\n", + "fillets2 = (sqrt(2) - 1 / (1 - theta)) * r0" ] }, { "cell_type": "code", - "execution_count": 56, + "execution_count": 11, "id": "dfa331c2-a2b0-4c81-8113-91ba38c35c42", "metadata": {}, "outputs": [ @@ -72,7 +73,7 @@ "" ] }, - "execution_count": 56, + "execution_count": 11, "metadata": {}, "output_type": "execute_result" }, @@ -99,7 +100,7 @@ }, { "cell_type": "code", - "execution_count": 64, + "execution_count": 12, "id": "b1ef64d1-ba45-4079-bac6-cbc1468174db", "metadata": {}, "outputs": [ @@ -109,13 +110,13 @@ "" ] }, - "execution_count": 64, + "execution_count": 12, "metadata": {}, "output_type": "execute_result" }, { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAEGCAYAAAB/+QKOAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAAx00lEQVR4nO3deVxU9f7H8ddhV1wAlwEVSZMUBRUtl2RJXBAQFZcs0yzXuuFaLpXZ1duiZSmauZFWmmWRWjIuJSaokSuJC6aYKC7ghgIqg4zn94fJTxIVhGGG4fN8PHpchzkz8/5c7u3tOd855yiqqqoIIYQQ/2Jh7ABCCCFMkxSEEEKIQklBCCGEKJQUhBBCiEJJQQghhCiUlbEDlJY///wTW1tbY8d4IJ1OZ/IZi8pcZjGXOUBmMVWmPotOp6Nly5aFPmc2BWFra4uHh4exYzxQUlKSyWcsKnOZxVzmAJnFVJn6LElJSfd9Tg4xCSGEKJQUhBBCiEJJQQghhCiUFIQQQohCSUEIIYQolBSEEEKIQklBCCGEKJQUBBCdeJZj6VnGjiGEECbFbE6Ue1R5+lt8oE0iPUvHoHZujO3sjkNlG2PHEkIIo6vwexBWlhasG+XDc0+58nV8Cs/M2sry+BTy9LeMHU0IIYyqwhcEQI0qtrwf5oV2tC8eztV456dDhMzdzo7ki8aOJoQQRiMFcRcPl2qsHN6WhQNbcf1mHi9E7mTE13s4eemasaMJIUSZk4L4F0VR6Obpwq/j/JkQ2JjtyRfp8mkcMzYcIVuXZ+x4QghRZqQg7sPO2pLXOjbitzeeIbRFHRbGHueZj7fy/Z5Ubt1SjR1PCCEMTgriITTV7Pjk2Rasfa0Drk6VmBiVSM/5O9iTctnY0YQQwqCkIIqopasDq199mjn9W3IhS0ffhfGM/jaBs1duGDuaEEIYhBREMSiKQi/vumx5w5/RAY3YdCiNgE+2MmfzUW7k6o0dTwghSpUUxCOobGPF+K6NiXndn04eGuZsPkanT7by8/6zqKqsTwghzIMURAnUc6zM/AGtWDWiHY72Noz+NoF+C+M5cPqqsaMJIUSJSUGUgrYNa/BzuA8zentx4uI1eszfzsSo/ZzPyjF2NCGEeGRSEKXE0kLhuTb1+W3CMwz3bciahDMEzIplYexxdHmyPiGEKH+kIEpZNTtr3gr24Jdx/rRr6MSMDUfoOjuOXw6lyfqEEKJcMWhBxMXFERgYSJcuXVi8ePE9z589e5ZBgwbRq1cvQkNDiY2NBeDmzZtMmjSJ0NBQgoKCWLRokSFjGkSDmvZEDn6Kr4e0wcbSghHL9/L2r2kkn5fLigshygeDFYRer2f69OlERkai1WqJjo4mOTm5wDYLFiwgKCiItWvXMnv2bKZNmwbAxo0byc3NZd26daxevZpVq1Zx+vRpQ0U1KL8narF+jC//DW3KsUs6us3ZxnvRh8nMuWnsaEII8UAGK4jExETc3NxwdXXFxsaGkJAQYmJiCmyjKArZ2dkAZGVlUbt27fyf37hxg7y8PHJycrC2tqZKlSqGimpw1pYWvNShAUvCXOn3ZD2+2HGCgFmx/CCX7RBCmDCD3TAoPT0dZ2fn/McajYbExMQC24SHhzN06FBWrFjBjRs3WLZsGQCBgYHExMTg4+NDTk4Ob775Jg4ODg/8PJ1OR1JSUqnPUZrsuMmLHna0q1WXBbsuMiEqkcitf/Fq2xo0rmln7HjFkpOTY/L/fReFucwBMoupKs+zGPWOclqtlrCwMIYMGUJCQgITJ04kOjqaxMRELCws2LZtG5mZmQwYMICnn34aV1fX+76Xra0tHh4eZZi++JKSkvDw8MADCO2gsibhDB9uOMK49Wd5trUrE7o1pmYVW2PHLJI7s5R35jIHyCymytRneVB5GewQk0ajIS0tLf9xeno6Go2mwDZRUVEEBQUB4O3tjU6nIyMjg+joaHx9fbG2tqZGjRq0atWKAwcOGCqqUVhYKPRpXY/f3vBnuG9Dftx3mo6ztrJ0+wluyt3shBAmwGAF4eXlRUpKCqmpqeTm5qLVagkICCiwjYuLC/Hx8QAcP34cnU6Hk5MTLi4u7Ny5E4Dr16+zf/9+GjZsaKioRlX1n6/FbhzrR0tXB6ZHHyZk7jZ+l7vZCSGMzGAFYWVlxdSpUxk2bBjBwcEEBQXh7u5ORERE/mL15MmT+f777+nRowfjx49nxowZKIrCCy+8wLVr1wgJCaFv37707t2bJk2aGCqqSWhUuwpfD2nD4kGtuXFTz4DInfznm72ckavFCiGMxKBrEP7+/vj7+xf42ZgxY/L/3KhRI7777rt7Xmdvb8/cuXMNGc0kKYpC12bO+D1Ri8Vxf/P51mS2HDnPq/6NGOnfEDtrS2NHFEJUIHImtQmys7ZkdCd3Yl5/hk5NNMzefJTOn8aySc7GFkKUISkIE1bXoRLzX2jFymFtqWxjycjle3lx6S6Sz2cbO5oQogKQgigHnm5Uk/WjfXk3tCl/pl6h25w43tceJkvOxhZCGJAURDlhZWnByx0a8Nsbz9C3dT0it5+g46xYovaelrOxhRAGIQVRztSsYsuMPs1Z+58O1HOsxBs/7KfPwt9JPH3F2NGEEGZGCqKcauHqwOpXn2ZWvxakXr5Bz/k7mPxjIpeydcaOJoQwE1IQ5ZiFhULf1vXY8oY/w3waELX39tnYX8enoJfDTkKIEpKCMAPV7Kx5O6QpG8f64lWvOlN/OkSPz7az92SGsaMJIcoxKQgz0qh2VVYMbcv8Aa24lJ1LnwW/M+GH/VyUw05CiEcgBWFmFEUhpLkLMa/7M9L/zr2xt7JcDjsJIYpJCsJM2dta8WaQR/5hp3fksJMQopikIMzcncNOnw3wzj/sNDFqv3zbSQjxUFIQFYCiKHRvXif/sNPqfWfoKIedhBAPIQVRgdx92Mmz7u3DTj3nb2ffKTnsJIS4lxREBdSodlW+GdaWec97cyFLR+/Pf2dSlJxkJ4QoSAqiglIUhdAWdYh5/RlG+t2+5WnAJ7Es/+OkHHYSQgBSEBVeFVsr3gz2YMMYX5q6VOOdtQfpOX87CXLYSYgKTwpCAOCuqcrK4W2Z+89hpzA57CREhScFIfIpikKPfw47jbjrsNMKOewkRIUkBSHuUcXWireCPVg/xhcPl6pMWXuQXvN38GfqFWNHE0KUISkIcV9PaKry7fB2RDzXkvTMHMI+v31J8cvXco0dTQhRBqyMHUCYNkVR6NmyLgFNajM35hhLd6Sw4WAaL7asTuPGKhYWirEjCiEMRPYgRJFU/eeS4utH+9LEuSrz4i/SZ+HvHDp71djRhBAGYtCCiIuLIzAwkC5durB48eJ7nj979iyDBg2iV69ehIaGEhsbm//ckSNH6N+/PyEhIYSGhqLTybdpTEFj56p8N6Idr/vU4tSl64TO2870dYfJ1uUZO5oQopQZ7BCTXq9n+vTpLFu2DI1GQ9++fQkICKBRo0b52yxYsICgoCAGDBhAcnIyI0aMYMuWLeTl5TFhwgQ+/vhjmjRpQkZGBlZWcjTMVCiKQufHq/JigDczNx1h2e8n0B44y9TuzQj2ckZR5LCTEObAYHsQiYmJuLm54erqio2NDSEhIcTExBTYRlEUsrOzAcjKyqJ27doA7Nixg8aNG9OkSRMAHB0dsbS0NFRU8YiqV7bmgzAvVr/6NDXsbXlt5T4GL9tNysVrxo4mhCgFBvtreXp6Os7OzvmPNRoNiYmJBbYJDw9n6NChrFixghs3brBs2TIATpw4gaIoDB06lMuXLxMcHMzw4cMf+Hk6nY6kpKTSH6QU5eTkmHzGorp7Fjvgo841iP7Lmq8SLtLl01j6N3egn6cDNpamvTdhrr+T8k5mMQ1GPW6j1WoJCwtjyJAhJCQkMHHiRKKjo9Hr9ezdu5eoqCgqVarESy+9hKenJ+3bt7/ve9na2uLh4VGG6YsvKSnJ5DMWVWGzeDaDlzvnMD36MCv+PMeO07n8r6cnPu41jZTy4cz9d1JeySxl50HlZbBDTBqNhrS0tPzH6enpaDSaAttERUURFBQEgLe3NzqdjoyMDJydnXnqqadwcnKiUqVK+Pn5cejQIUNFFaVIU82O+QNa8fWQNqiqysAvdjLq2wTOZ+YYO5oQopgMVhBeXl6kpKSQmppKbm4uWq2WgICAAtu4uLgQHx8PwPHjx9HpdDg5OeHj48PRo0e5ceMGeXl57N69u8DitjB9fk/UYuNYP8Z0cmfTwTQ6fRLLlztOyCU7hChHDHaIycrKiqlTpzJs2DD0ej19+vTB3d2diIgIPD096dSpE5MnT2bKlCl8+eWXKIrCjBkzUBSF6tWr89JLL9G3b18URcHPz49nnnnGUFGFgdhZWzKuyxP08q7L1J8O8t91h4nad5r3e3nRwtXB2PGEEA+hqKpqFn+lM/XjfFA+MhZVcWdRVRXtgXNMX3eYC9k6XmhbnwmBTaheydqAKR+uIv9OTJnMUnYelE/OpBZl4u77Yr/09GOs3HmKTp9sZU3Caczk7yhCmB0pCFGmqtpZ825oM34O96GuY2XGrdrPgCU7ST6fbexoQoh/kYIQRuFZtzqrX32a93p5cujsVYIi4vh40xFu5OqNHU0I8Q8pCGE0lhYKA9u5EfP6M4Q2r8P8347TZXYsW46kGzuaEAIpCGECalW15dP+Lfl2eDvsrC0Z8uUeRi7fw9krN4wdTYgKTQpCmIz2j9dg/WhfJgQ2JvboBTp/GssX20+Qp79l7GhCVEhSEMKk2FhZ8FrHRvw6zp+2DZz4X/Rhen2+gwOn5b4TQpQ1KQhhklydKrP0paeYP6AV6Zk6es6X+04IUdakIITJUhSFkOYubB7vz4C29Vn2+wm6fhrLr4dlEVuIsiAFIUxe9UrWvNfLi6hXnqaqnTXDv769iH3uqixiC2FIUhCi3Gjt5kj0aB8mdWtC7NELdPk0Ti4AKIQBSUGIcsXa0oJXn3mcX8b6413fgf+uO0zvz3dw6KwsYgtR2qQgRLlUv0Zlvh7ShrnPe3Pmyg16fLaD97WHuSaL2EKUGikIUW4pikKPFnWIGf8Mzz7pypJtJ+g6O07OxBailEhBiHKvemVrPuztxQ+vtKeyze0zsf/zzV7S5S52QpSIFIQwG0895oT2nzOxNyedp/MnsSyPT5FFbCEekRSEMCt3zsT+ZawfLVwdeOenQ/RZ8DtJ5zKNHU2IckcKQpilx2ras3xoG+b0b0nq5et0n7edDzckcT1XFrGFKCopCGG2FEWhl3ddNo/3p2+reiyK/Zuus+P47a/zxo4mRLkgBSHMnqO9DTP7NmfViHbYWlnw8rLdhK/cx/ksWcQW4kGkIESF0bZhDdaP8WV8lyf45VA6nT6JZf1fmdySRWwhCiUFISoUWytLRndyZ+NYXzzrVGfeHxd5dlE8yeezjB1NCJNj0IKIi4sjMDCQLl26sHjx4nueP3v2LIMGDaJXr16EhoYSGxt7z/Pe3t588cUXhowpKqCGtaqwcnhbxneoxbHz2QRHbGduzDFy8+TmRELcYbCC0Ov1TJ8+ncjISLRaLdHR0SQnJxfYZsGCBQQFBbF27Vpmz57NtGnTCjw/Y8YMfH19DRVRVHCKotClUVU2j/cn0NOZT389Svd529h7MsPY0YQwCQYriMTERNzc3HB1dcXGxoaQkBBiYmIKbKMoCtnZ2QBkZWVRu3bt/Oc2b95M3bp1cXd3N1REIYDb98Se97w3Xwx+kqycPPou/J3//nxIbk4kKjwrQ71xeno6zs7O+Y81Gg2JiYkFtgkPD2fo0KGsWLGCGzdusGzZMgCuXbvGkiVLWLp0KUuXLi3S5+l0OpKSkkpvAAPIyckx+YxFZS6z3D1HHeCzEGe+3HeZr35PQbs/lfB2tWhTr7JxQxaRufxOQGYxFQYriKLQarWEhYUxZMgQEhISmDhxItHR0Xz22WcMHjwYe3v7Ir+Xra0tHh4eBkxbcklJSSafsajMZZbC5mjdHF46mcHkHxN5NyaNHi3qMDW0KTWr2BopZdGYy+8EZJay9KDyMlhBaDQa0tLS8h+np6ej0WgKbBMVFUVkZCQA3t7e6HQ6MjIy2L9/P5s2bWLWrFlkZmZiYWGBra0tAwcONFRcIQq4c3OihVv/5rPfjhF37AJTQprSp1VdFEUxdjwhyoTB1iC8vLxISUkhNTWV3NxctFotAQEBBbZxcXEhPj4egOPHj6PT6XBycmLlypVs2bKFLVu2MHjwYEaOHCnlIMqcrZUlYzq7s360L4/XqsIbP+xn0Be7OHXpurGjCVEmDFYQVlZWTJ06lWHDhhEcHExQUBDu7u5ERETkL1ZPnjyZ77//nh49ejB+/HhmzJghfzsTJsddU5UfRrbnf708+TP1Cl3nxLIk7m/y9PKVWGHeDLoG4e/vj7+/f4GfjRkzJv/PjRo14rvvvnvge4waNcog2YQoDgsLhUHt3OjsUZt31h7k/fVJ/Lz/LDP6eNGsTnVjxxPCIORMaiGKwaV6JZa8+CTzB7Ti3NUceny2gxkbjpBzU2/saEKUOikIIYpJURRCmrsQ889VYhfGHqfbnDh+P37R2NGEKFVSEEI8ouqVrZnZtzkrh7VFBQYs2cnEqP1cvX7T2NGEKBUPLYiPPvqI7Oxsbt68yeDBg2nXrh0//fRTWWQTolx4ulFNNo314xX/x/lx3xk6fRqLNvEcqipXiRXl20MLYseOHVSpUoWtW7dSt25dfv31V7l4nhD/YmdtyeSgJvz0Wgecq9vy2sp9DP96L+eu3jB2NCEe2UMLIi/v9vVotm7dSrdu3ahatarBQwlRXnnWrc7a/3Tg7WAPtidfoMuncSyPT5F7Tohy6aEF0bFjR7p168ahQ4do3749ly9fxtbWtC85IIQxWVlaMNyvIb+M9ce7vgPv/HSI/ovj+ftCtrGjCVEsDy2I0aNH89133/Hjjz9ibW2NnZ0dCxYsKItsQpRr9WtU5ushbZjVrwV/pWXRLWIbC2OPywl2otx4aEH0798fBwcHLC0tAahcuTLDhw83eDAhzIGiKPRtXY/N4/3p2LgWMzYcIezz3zmSlmnsaEI81H3PpL5w4QLp6enk5ORw+PDh/G9kZGdnc+OGLLwJURy1q9mxcGBr1h9IY+pPBwmdt53/PNOI1zo2wsZKvm0uTNN9C2L79u2sXr2atLQ0Pvzww/yfV6lShfHjx5dJOCHMyZ0T7No/XoPp6w4REXOMjQfT+Khvc1q4Ohg7nhD3uG9BhIWFERYWxqZNmwgMDCzLTEKYNSd7G+Y8501oizq8veYgYZ/vYLhvQ8Z1eQI7a0tjxxMi30P3bVu1asVbb73FsGHDAEhOTuaHH34weDAhzF0nDw2/jPej/1OuLIr7m6CIbew6cdnYsYTI99CCePPNN/Hx8eH8+fMAPPbYY3z99dcGDyZERVDNzpoPe9++XEferVs8uyieqT8dlPthC5Pw0ILIyMggODgYC4vbm1pZWeX/WQhROu5cruPlDo+x/I+TBM6OI+7oBWPHEhXcQ/9NX7lyZTIyMvJv5PPnn3/K2dRCGEBlGyveDW1G1CvtsbW24MWlu5jwg1z8TxjPQ28YNHnyZF599VVOnTrFc889R0ZGBhEREWWRTYgKqbWbE+tH+zI35hiL4v4m9ugF3uvlSddmzsaOJiqYhxZEs2bNWLFiBSdOnEBVVRo0aIC1tXVZZBOiwrKztmRityYEe7kwISqREcv30r25C9N6NKNGFbnUjSgb9y2IX375pdCfp6SkANC1a1eDBBJC/D/PutX5ObwDC7ceZ96WZHYkX+S/PZrRo0UduX+7MLj7FsRvv/32wBdKQQhRNqwtLRjVyZ1AT2cmRCUy5rs/Wbf/LO/18sK5up2x4wkzdt+CuPvsaSGE8T2hqcrqV59m2Y4TzPrlL7rMjmVKiAfPPukqexPCIO5bEMuWLXvgC19++eVSDyOEeDBLC4Vhvg3p7KFh0o+JTPrxAOv2n+PD3l7GjibM0H0L4tq1a2WZQwhRDI/VtOfb4e1YuesUH65PInBOHC+2dKBxYxULC9mbEKXjvgURHh5e4jePi4vj/fff59atW/Tr148RI0YUeP7s2bNMmjSJrKws9Ho9b7zxBv7+/uzYsYNPPvmEmzdvYm1tzYQJE2jfvn2J8whhTiwsFAa2c6Njk9q8tfoAC3ddYN+FP/i4b3PcatgbO54wA/ctiCVLljB8+HD+97//FXp8c8qUKQ98Y71ez/Tp01m2bBkajYa+ffsSEBBAo0aN8rdZsGABQUFBDBgwgOTkZEaMGMGWLVtwdHRkwYIFaDQajh49ytChQ9m2bVsJxhTCfNV1qMSXLz/F3HW7idybQbc525jUrTEvtn9M9iZEidy3ICIjIxk+fDiurq5Ur1692G+cmJiIm5sbrq6uAISEhBATE1OgIBRFITv79m0Ys7KyqF27NgBNmzbN38bd3R2dTkdubi42NjbFziFERaAoCl3dq/KsvxeTfzzAf9cdZsPBND7u24L6NSobO54op+5bEDVq1CA9PZ3Vq1ezfPny/BsGFVV6ejrOzv9/5qdGoyExMbHANuHh4QwdOpQVK1Zw48aNQhfGN23aRNOmTR9aDjqdjqSkpGJlLGs5OTkmn7GozGUWc5kDbs9y5WwKk9pVwbumyuLdl+g6eytDWjsR0rgaFuXom07m9nspr7PctyCef/55XnrpJVJTU+ndu3f+z1VVRVEUYmJiSvzhWq2WsLAwhgwZQkJCAhMnTiQ6Ojr/YoDHjh1j1qxZLF269KHvZWtri4eHR4kzGVJSUpLJZywqc5nFXOaAgrM0bQrP+t1g8uoDfL7zAvsuqHzctwWuTuVjb8Jcfy+m6EHldd+L9Q0aNIgNGzbQp08fYmJi8v/ZsmVLkcpBo9GQlpaW/zg9PR2NRlNgm6ioKIKCggDw9vZGp9ORkZEBQFpaGuHh4cycOZP69es/9POEEAXVcajEVy8/xYzeXhw8k0ngnDiWx6dw61bxjgaIiuuhV3OdNm3aI72xl5cXKSkppKamkpubi1arJSAgoMA2Li4uxMfHA3D8+HF0Oh1OTk5kZmYyYsQIXn/9dVq3bv1Iny+EuL028Vyb+mwa50drN0fe+ekQL0TuJPXydWNHE+WAwW7sYGVlxdSpUxk2bBjBwcEEBQXh7u5ORERE/h7I5MmT+f777+nRowfjx49nxowZKIrCihUrOHXqFPPnz6dnz5707NmTS5cuGSqqEGavrkMlvh7Shg97e3HgzNXbexN/nJS9CfFAilrc1WcTZerH+aB8ZCwqc5nFXOaAos9y5soNJkUlsj35Ih0a1WBmn+bUczSttYmK+Hsxlgflk1vDCVHB1HWoxPKhbfggzIs/T10hcHYc3+w8WexvKgrzJwUhRAWkKAoD2t5em2hZ34G31xxk0Be7OJ0haxPi/0lBCFGB1XOszIqhbXk/zJOEU7fPwl6585TsTQhACkKICk9RFF5o68bGsX60cK3OW2sO8OLSXZy5csPY0YSRSUEIIQBwdbq9N/FeL0/2nswgcHYc3+6SvYmKTApCCJFPUW5fIXbTWD+a16vOm6tlb6Iik4IQQtzjzt7E/3o2y9+b+E72JiocKQghRKEsLBQGtX+MTWP98KxbjcmrDzB42W7OXZW9iYpCCkII8UCuTpVZOawd03s2Y/eJy3SdHceahNOyN1EBSEEIIR7KwkLhxfaPsXGsL401VRm3aj+vrtjHpWydsaMJA5KCEEIUmVsNe1aNbM+bQU3YcuQ8XWfHselQ2sNfKMolKQghRLFYWiiM9H+cdaN8cK5ux8jlexn//Z9cvXHT2NFEKZOCEEI8ksbOVVnznw6MDmjET3+epducOLYdu2DsWKIUSUEIIR6ZjZUF47s2ZvWrT1PZxpJBX+zinbUHuZ6bZ+xoohRIQQghSqyFqwPa0b4M9WnAip0nCY7Yxt6Tl40dS5SQFIQQolTYWVvyTvemfDu8HXm3VPotjOfDDUno8vTGjiYekRSEEKJUtWtYg41j/ej/lCuLYv+mx7wdHDxz1dixxCOQghBClLoqtlZ82Ls5y156iozrufSav4N5McfI098ydjRRDFIQQgiD6dikNr+M8yPYy4VPfj1KnwW/k3w+29ixRBFJQQghDMqhsg1zn/dm/oBWnLp8nZC521i6/QS3bsmlOkydFIQQokyENHdh0zg/fBrVZHr0YQZE/kHqZbnFqSmTghBClJnaVe2IHPwkH/VpzsEzmQRFbGPVbrmMuKkyaEHExcURGBhIly5dWLx48T3Pnz17lkGDBtGrVy9CQ0OJjY3Nf27RokV06dKFwMBAtm3bZsiYQogypCgKzz7lyoYxvnjWrcakHw8w9Ks9nM/MMXY08S8GKwi9Xs/06dOJjIxEq9USHR1NcnJygW0WLFhAUFAQa9euZfbs2UybNg2A5ORktFotWq2WyMhIpk2bhl4v36UWwpzcuYz4u6FN2ZF8ka5z4li3/6yxY4m7GKwgEhMTcXNzw9XVFRsbG0JCQoiJiSmwjaIoZGff/kZDVlYWtWvXBiAmJoaQkBBsbGxwdXXFzc2NxMREQ0UVQhiJhYXCyx0asH6ML2417Bn1bQLhK/eRmSN/ITQFBiuI9PR0nJ2d8x9rNBrS09MLbBMeHs66devw8/NjxIgRTJkypcivFUKYj8drVeHHV9rzRtcn2HQojVd/Pk3sUbnwn7FZGfPDtVotYWFhDBkyhISEBCZOnEh0dPQjvZdOpyMpKamUE5aunJwck89YVOYyi7nMAeYxSycXeCy4DjNj0xm8dBehTaoxpLUTdlbl9/s05fn3YrCC0Gg0pKX9/41E0tPT0Wg0BbaJiooiMjISAG9vb3Q6HRkZGUV67b/Z2tri4eFRihOUvqSkJJPPWFTmMou5zAHmM4sH4Fr9ED+dUFi64wSHL+mZ098br3rVjR3tkZj67+VB5WWwWvby8iIlJYXU1FRyc3PRarUEBAQU2MbFxYX4+HgAjh8/jk6nw8nJiYCAALRaLbm5uaSmppKSkkLz5s0NFVUIYWJsLC2YGtqUFUPbck2nJ+zzHXy2RS7VUdYMtgdhZWXF1KlTGTZsGHq9nj59+uDu7k5ERASenp506tSJyZMnM2XKFL788ksURWHGjBkoioK7uztBQUEEBwdjaWnJ1KlTsbS0NFRUIYSJ8nGvycaxvkxZe5BZvxzlt78uMPvZltSvUdnY0SoERTWTM1RMfTcOykfGojKXWcxlDjDvWVRV5ef9Z5my9iC3bqm8G9qMfk/WQ1EUI6YsGlP/vTwoX/ld+RFCVBiKotCzZV02jvXDq151Jv6YyMjle7mUrTN2NLMmBSGEKDfqOlRi5bB2vB3swda/LhA4ZxtbjshX4A1FCkIIUa5YWCgM92vIT+EdqFnFhiFf7uHtNQfkPtgGIAUhhCiXPFyqsfa1Dozwa8jKXacImbudP1OvGDuWWZGCEEKUW3bWlrwV7MHKYe3Q3dTTZ8HvzNl8VL4OW0qkIIQQ5V77x2uwYawfoc1dmLP5GH0WxnPi4jVjxyr3pCCEEGaheiVr5jznzbznvTlxIZvgiG18s/Ok3GuiBKQghBBmJbRFHTaN86OVmwNvrznIsK/2cCFLvg77KKQghBBmx6V6JZYPacvU7k3ZlnyRbnPi+PWwfB22uKQghBBmycJCYYhPA6JH+aCpZsfwr/cwKSqRazr5OmxRSUEIIczaE5qqrH2tA6/4P873e1MJitjG3pMZxo5VLkhBCCHMno2VBZODmrBqRHv0t1T6LfydT375i5vyddgHkoIQQlQYbRo4sXGsL2He9Zi3JZm+C+M5eUm+Dns/UhBCiAqlqp01nzzbgvkDWuV/HTZq72n5OmwhpCCEEBVSSHMXNoz1o1nd6rzxw35GfZvA1es3jR3LpEhBCCEqrLoOlfh2eDsmBDZm48E0giLi2Pn3JWPHMhlSEEKICs3SQuG1jo2IevVprK0seH7JH8zaJAvYIAUhhBAAtHR1QDvalz6t6vHZb7KADVIQQgiRr4qtFR/3kwXsO6QghBDiX2QB+zYpCCGEKIQsYEtBCCHEfd29gG1TARewpSCEEOIh/r2A3a+CLGAbtCDi4uIIDAykS5cuLF68+J7nP/jgA3r27EnPnj0JDAzkySefzH/uo48+IiQkhKCgIN57770Ku0gkhDAN9nctYP9dQRawrQz1xnq9nunTp7Ns2TI0Gg19+/YlICCARo0a5W/z1ltv5f95+fLlHD58GIB9+/axb98+fv75ZwAGDBjArl27aNu2raHiCiFEkYQ0d8G7vgPjVv3JGz/sZ+tf53k/zIvqlayNHa3UGWwPIjExETc3N1xdXbGxsSEkJISYmJj7bq/VaunevTsAiqKQm5vLzZs38/+zZs2ahooqhBDFUsehEivvWsAOjthmlgvYBtuDSE9Px9nZOf+xRqMhMTGx0G3PnDnD6dOnadeuHQDe3t60bdsWHx8fVFVl4MCBPP744w/8PJ1OR1JSUukNYAA5OTkmn7GozGUWc5kDZBZjCHCGukEufBR3nueX/MGzng680NIRKwslf5vyMkthDFYQxaHVagkMDMTS0hKAkydPcvz4cWJjYwEYMmQIe/bsKbBG8W+2trZ4eHiUSd5HlZSUZPIZi8pcZjGXOUBmMRYPoEubPKatO8R3e05z5ApEPNcStxr2gOnP8qDyMtghJo1GQ1paWv7j9PR0NBpNoduuX7+ekJCQ/Me//vorLVq0wN7eHnt7e3x9fUlISDBUVCGEKBF7Wys+6mt+C9gGKwgvLy9SUlJITU0lNzcXrVZLQEDAPdsdP36czMxMvL29839Wp04ddu/eTV5eHjdv3mT37t0PPcQkhBDGFtLchY1j/fC86wzs7Fy9sWM9MoMdYrKysmLq1KkMGzYMvV5Pnz59cHd3JyIiAk9PTzp16gTc3nsIDg5GUf7/mF1gYCB//PEHoaGhKIqCr69voeUihBCm5s4C9qK443z6y1F2HbdkgUNdWrs5GjtasSlqed8H+oepH+eD8pGxqMxlFnOZA2QWU5RwKoNXv97Fhet6xnd5glf8H8fyrgVsU/Cg/67lTGohhDAQ7/qOfBZaj2AvFz7e9BcvLt1JemaOsWMVmRSEEEIYkL2NBXOfa8lHfZqz7+QVgiK28duR88aOVSRSEEIIYWCKovDsU66sG9WB2lVtefnL3bwXfZjcPNO+6J8UhBBClJFGtauy9rUODG7vRuT2E/RZ8DsnLpruRf+kIIQQogzZWVsyracniwe1JjXjOt3nbmNNwmljxyqUFIQQQhhB12bOrB/tS7O61Rm3aj/jv/+Ta7o8Y8cqQApCCCGMpM4/d60b29mdtQln6D5vOwfPXDV2rHxSEEIIYUSWFgpjOz/Bt8PbcSNXT9jnO/hi+wmTuEyHFIQQQpiAtg1rsGGML880rs3/og8z9Ks9XMrWGTWTFIQQQpgIR3sbFg9qzfSezdiefJGgiG38fvyi0fJIQQghhAlRFIUX2z/G2v90oIqdFS9E7mTWpr/I05f9ORNSEEIIYYKa1qlG9Cgf+rWux2e/JdN/8R+czrhephmkIIQQwkRVtrl9n4mI51ryV1oWwRHb2HDgXJl9vhSEEEKYuJ4t67J+tC8NalXh1W/28daaA+TcNPx9JqQghBCiHKhfozI/jGzPSP+GrNx5ih6fbedoepZBP1MKQgghygkbKwveDPLg6yFtuHwtl9B52/lm50mDnTMhBSGEEOWM3xO12DDGjzYNnHh7zUE+/fWoQT7HYLccFUIIYTi1qtry1cttWLnrFE72Ngb5DCkIIYQopywsFAa2czPc+xvsnYUQQpRrUhBCCCEKJQUhhBCiUAZdg4iLi+P999/n1q1b9OvXjxEjRhR4/oMPPmDnzp0A5OTkcOnSJfbs2QPA2bNnmTJlCufOnUNRFBYvXky9evUMGVcIIcRdDFYQer2e6dOns2zZMjQaDX379iUgIIBGjRrlb/PWW2/l/3n58uUcPnw4//GkSZN45ZVX6NChA9euXcPCQnZ2hBCiLBns37qJiYm4ubnh6uqKjY0NISEhxMTE3Hd7rVZL9+7dAUhOTiYvL48OHToAYG9vT6VKlQwVVQghRCEMVhDp6ek4OzvnP9ZoNKSnpxe67ZkzZzh9+jTt2rUDICUlhWrVqhEeHk6vXr2YOXMmer3hrzsihBDi/5nEeRBarZbAwEAsLS0ByMvLY8+ePaxduxYXFxfGjRvH6tWr6dev333fQ6fTkZSUVFaRH0lOTo7JZywqc5nFXOYAmcVUledZDLYHodFoSEtLy3+cnp6ORqMpdNv169cTEhKS/9jZ2RkPDw9cXV2xsrKiU6dOBdYnyis7OztjRyg15jKLucwBMoupMvVZdLr739bUYHsQXl5epKSkkJqaikajQavV8sknn9yz3fHjx8nMzMTb27vAazMzM7l8+TJOTk7s3LkTT0/PB35ey5YtS3sEIYSo0AxWEFZWVkydOpVhw4ah1+vp06cP7u7uRERE4OnpSadOnYDbew/BwcEoipL/WktLSyZNmsTgwYMBaNas2QMPLwkhhCh9imqo68QKIYQo1+TkAiGEEIWSghBCCFEoKQghhBCFkoIQQghRKCmIUhIXF0dgYCBdunRh8eLF9zyfm5vL2LFj6dKlC/369eP06dMA7Nixg969exMaGkrv3r2Jj48v6+gFPOocd5w9exZvb2+++OKLsop8XyWZ5ciRI/Tv35+QkBBCQ0Mf+F3xsvCos9y8eZNJkyYRGhpKUFAQixYtKuvo93jYLLt37yYsLIymTZuycePGAs+tWbOGrl270rVrV9asWVNWkQv1qHMkJSUV+N/W+vXryzJ28aiixPLy8tROnTqpp06dUnU6nRoaGqoeO3aswDYrVqxQ33nnHVVVVTU6OlodM2aMqqqqeujQITUtLU1VVVX966+/VB8fnzLNfreSzHHHqFGj1FGjRqmRkZFlFbtQJZnl5s2bavfu3dWkpCRVVVX18uXLal5eXpnmv1tJZvn555/VsWPHqqqqqtevX1c7duyopqamlmn+uxVlltTUVDUpKUmdMGGCumHDhvyfZ2RkqAEBAWpGRoZ65coVNSAgQL1y5UpZj6Cqasnm+Pvvv9UTJ06oqqqqaWlpaocOHdSrV6+WZfwikz2IUlCUCxNu2bKFsLAwAAIDA4mPj0dVVZo2bZp/hrm7uzs6nY7c3NwynwFKNgfA5s2bqVu3Lu7u7mWe/d9KMsuOHTto3LgxTZo0AcDR0TH/MjDGUJJZFEXhxo0b5OXlkZOTg7W1NVWqVDHGGEDRZqlXrx5NmjS55wrO27dvp0OHDjg4OFC9enU6dOjAtm3byjJ+vpLM0aBBAx577DHg9hUnnJycuHz5cllFLxYpiFJQlAsTpqen4+LiAtw+ibBq1apkZGQU2GbTpk00bdoUGxvD3ID8YUoyx7Vr11iyZAnh4eFlmvl+SjLLiRMnUBSFoUOHEhYWxpIlS8o0+7+VZJbAwEAqVaqEj48PHTt2ZMiQITg4OJRl/HtyFvUinqX52tJWWlkSExO5efMm9evXL814pcYkLtYn4NixY8yaNYulS5caO8oj+eyzzxg8eDD29vbGjlJier2evXv3EhUVRaVKlXjppZfw9PSkffv2xo5WbImJiVhYWLBt2zYyMzMZMGAATz/9NK6ursaOVuGdP3+eCRMmMHPmTJO9341ppipninJhQo1Gw7lz54DbV6vNysrC0dERgLS0NMLDw5k5c6ZR/yZRkjn279/PrFmzCAgI4KuvvmLRokWsWLGiTPP/O+ejzuLs7MxTTz2Fk5MTlSpVws/Pj0OHDpVp/n/nfNRZoqOj8fX1xdramho1atCqVSsOHDhQpvn/nbOoF/EszdeWtpJmyc7OZuTIkYwbN86kryMnBVEK7r4wYW5uLlqtloCAgALbBAQE5H/rYtOmTbRr1w5FUcjMzGTEiBG8/vrrtG7d2hjx85VkjpUrV7Jlyxa2bNnC4MGDGTlyJAMHDjTGGEDJZvHx8eHo0aP5x+53795d4E6IZa0ks7i4uOTf1vf69evs37+fhg0blvkMdxRllvvx8fFh+/btXL16latXr7J9+3Z8fHwMnLhwJZkjNzeX1157jZ49e9KtWzcDJy0h466Rm4+tW7eqXbt2VTt16qR+/vnnqqqq6pw5c9TNmzerqqqqOTk56qhRo9TOnTurffr0UU+dOqWqqqrOnz9fbdGihdqjR4/8fy5evFju5rjb3Llzjf4tJlUt2Sxr165Vg4OD1ZCQEHXmzJlGyX+3R50lOztbHTVqlBocHKwGBQWpS5YsMdoMdzxslv3796u+vr5qixYt1DZt2qjBwcH5r/3hhx/Uzp07q507d1ajoqKMkv+OR51j7dq1atOmTQv8f/7w4cNGm+NB5GJ9QgghCiWHmIQQQhRKCkIIIUShpCCEEEIUSgpCCCFEoaQghBBCFEoKQohiyszM5JtvvgFg586djBw5slivX716tdEuESFEcUhBCFFMmZmZfPvtt4/8+jVr1nD+/PlSTCSEYch5EEIU07hx44iJiaFBgwZYWVlRuXJlHB0dOXr0KM2aNWPWrFkoisLBgweZMWMG169fx9HRkQ8//JB9+/bx5ptvUrt2bezs7Fi1ahWRkZH89ttv6HQ6vL29mT59OoqiGHtMIeRMaiGKKzU1VQ0JCVFVVVX/+OMPtVWrVuq5c+dUvV6vPvvss+ru3bvV3NxctX///uqlS5dUVVVVrVarTp48WVVVVR04cKCamJiY/34ZGRn5f37jjTfUmJiYshtGiAeQq7kKUULNmzfPv/RzkyZNOHPmDNWqVePo0aO8/PLLANy6dYtatWoV+vqdO3cSGRlJTk4OV65cwd3dvcjX9RHCkKQghCihu+/fYWlpiV6vR1VV3N3dWbVq1QNfq9PpmDZtGj/++CMuLi7MmzfP6Lc3FeIOWaQWopjs7e25du3aA7dp0KABly9fJiEhAbh9b+hjx47d8/o7ZeDo6Mi1a9fYtGmTAZMLUTyyByFEMTk6OtKqVSu6d++Ora0tNWvWvGcbGxsb5s6dy3vvvUdWVhZ6vZ7Bgwfj7u5OWFgY7777bv4idb9+/ejevTs1a9bEy8vLCBMJUTj5FpMQQohCySEmIYQQhZKCEEIIUSgpCCGEEIWSghBCCFEoKQghhBCFkoIQQghRKCkIIYQQhfo/TcHgo6JLt0wAAAAASUVORK5CYII=\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAEGCAYAAAB/+QKOAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAA3EElEQVR4nO3deXxU9b3/8dcsmck+kwSYhGQyCSZAMGEJiyBoahBQIlUg1Ou1ltILeKl6sfXa296fN9Z4LbUVFbXFBU1tUW8LSimmLdUIDTuRbViCrEMSJAOShTXb5Pz+OJMJgRACyWQm4fN8PM4D5izJ52uAt+d8z/f71SiKoiCEEEJcRuvrAoQQQvgnCQghhBCtkoAQQgjRKgkIIYQQrZKAEEII0Sq9rwvoLDt37sRoNPq6jDbV1tb6fY3t1VPa0lPaAdIWf+XvbamtrWXo0KGtHusxAWE0GklJSfF1GW0qLi72+xrbq6e0pae0A6Qt/srf21JcXHzVY/KISQghRKskIIQQQrRKAkIIIUSrJCCEEEK0SgJCCCFEqyQghBBCtEoCQgghRKu8GhCFhYVMmjSJCRMm8Pbbb1/1vNWrVzNgwAB2797t2ffWW28xYcIEJk2axLp167xZJuz6I9j/BKcOQGOjd7+XEEJ0E14bKOdyucjNzSUvLw+LxUJ2djaZmZkkJSW1OO/cuXP8/ve/Z8iQIZ59hw4dIj8/n/z8fJxOJ7NmzWL16tXodLrOL7SxEdb+Aiod6mdDGMQMhr7DIGao+mtkP9DKzZYQ4ubitYCw2+3YbDasVisAWVlZFBQUXBEQixYtYs6cObz77ruefQUFBWRlZWEwGLBardhsNux2O8OGDev8QrVaeHwbfPMVfL0Dvt6p/rr1HXDVqucYwyFmCPQd2jI0NJrOr0cIIfyE1wLC6XQSHR3t+WyxWLDb7S3O2bt3L+Xl5XzrW99qERBOp7PFHYXFYsHpdLb5/Wpra9scMn5tWggcDv2GQz+gsQFj9VECK4sJqthPYOV+jCVvom2sB8AVEEpNxEAuRg6kJmIgNZEp1If0bTM0ampqOlij/+gpbekp7QBpi7/qzm3x2VxMjY2N/PKXv2TBggWd8vW8MxdTGvDt5o8NdXCqGL7eie7rHYSc2EnIwT+Bq049HmhW7zIufTxljveEhr/PyXI9ekpbeko7QNrir/y9LW2Fl9cCwmKxUF5e7vnsdDqxWCyez+fPn+fAgQN873vfA+DUqVPMmzePxYsXX/Nan9Eb1EdNMUNg+Ex1X0MdnNynPpY6sVP9deMb4L7TICjCExhhSm/oGwbhsfJ4Sgjh97wWEGlpaTgcDkpLS7FYLOTn57Nw4ULP8bCwMLZs2eL5/Mgjj/CTn/yEtLQ0AgMDeeqpp5g1axZOpxOHw8HgwYO9VWrH6A3uu4ahzfsaasG597LQeI24xgbY8DMI6aOGRmy6+mvfdAjt7aMGCCFE67wWEHq9npycHGbPno3L5WL69OkkJyezaNEiUlNTGT9+/FWvTU5O5t5772Xy5MnodDpycnK88waTt+iN6j/+senN++prOLrlUxINlXB8uxoaB/8BKOrx8DiIdYdF32HqFmT2RfVCCAF4uQ8iIyODjIyMFvvmz5/f6rl/+MMfWnyeN28e8+bN81ptXS4gkJqoW+HSZ5G15+DELvfbU9vV4Che1Xw88pZL7jTS1ddvDSFdX7sQ4qbUYxYM6paMoZAwVt2aXKhQH0s13WWUbII9y9VjGi30Hth8hxGbDpZU9Y5FCCE6mQSEvwmOhFsy1a3JWWfzXcbXO+DAatj5gXpMGwCWW1v2Z/QeCDr50QohOkb+FekOwiww4B51A1AUqC5Vw+L4djU4di+HL99TjwcEq29axQ5394UMB7NN3pwSQlwXCYjuSKNRx1eY42HQ/eq+xkaoONLcl3F8W8vR4EGR7sC4JDRCevmuDUIIvycB0VNotdArSd0Gf0fd56pXx2gc3+betsPhAlDcExKa49Wg6OsOjJghar+IEEIgAdGz6QKaB/aN+IG6r+nNqePb3Hcb22DvCvWYRgu9U9TXbZvuNvoMUr+OEOKmIwFxs2ntzalzp1o+mtr/V9ixVD2mD4TowS0fT0X2803tQoguJQEh1FHc/SepG6id4FXHmh9LHd8O29+HLYvV44FmrOYBcCID4kZIf4YQPZQEhLiSRgMRCeqWOl3d52qAU/vVO42yL9Ef2QjrXmruz4hIcN9hjFBDI3owBAT6qAFCiM4gASHaR6eH6FR1S/8eR4uLSbklXl0/4/iXUPYllGyGPR+r52sD1HObAiN2uDoyXBZeEqLbkIAQN84QcmV/xpkT7kdT7tDY9REUvaMeCzS1vMuQR1NC+DUJCNG5wmMg/D5IuU/93OiCU19dEhrbWj6aMtvcYSGPpoTwNxIQwru0OrAMUrf0R9R9decvezS15ZJHU3p1fqm4ERA3Ut1keVchfEICQnS91h5NnS1Xw8LzaOr/oGiJeiwoojksmh5NBZp8U7sQNxEJCOEfwqLVx1KXP5oqK3JvX8LBz1DXz9BA7wEt7zJ6D1TvVoQQnUYCQvinSx9NNS3vWlOtjsko+1INjUsH9BlC1UF8TYERO0JW6ROig7waEIWFhbzwwgs0NjYyY8YM5s6d2+L4Rx99xIcffohWqyU4OJjnn3+epKQkysrKmDx5MomJiQAMGTKE3Nxcb5YquoNAE9xyl7qBOqCv4khzYJQVwYZF0NigHo9IaPloypKmLhErhGgXrwWEy+UiNzeXvLw8LBYL2dnZZGZmkpSU5DlnypQpPPTQQwAUFBSwYMEC3n33XQDi4+NZuXKlt8oTPYFGA1G3qNuQB9V9dRfUuaaaAsOxHnYvU4/pjOra4U2BETdSDRkhRKu8FhB2ux2bzYbVagUgKyuLgoKCFgERGto8c+jFixfRyJsqoqMMwWAbo25Nqo+37MvY+g5segOApKDesHsMWEeB9TZ1YkNZoU8IwIsB4XQ6iY6O9ny2WCzY7fYrzvvggw/Iy8ujvr6e999/37O/rKyMBx54gNDQUJ588klGjBjR5verra2luLi48xrgBTU1NX5fY3t1u7ZoB0D8AIj/LrjqCaw+RNDpPRhO7iK05EsMxX8BoFEbQE3EAC5GpXKxVxoXo9JoCO7j4+Lbp9v9TNogbfEPPu+kfvjhh3n44YdZtWoVixcv5sUXX6RPnz6sWbOGiIgI9uzZw2OPPUZ+fn6LO47LGY1GUlJSurDy61dcXOz3NbZX92/LYGAaxcXFRKakqMu6lm1FW7qV4LIigo/8GQ78n3pqeBxYR0Kc+y4j2j/7Mrr/z6SZtKXrtBVeXgsIi8VCeXm557PT6cRisVz1/KysLH7+858DYDAYMBjUv4CpqanEx8dz9OhR0tLSvFWuuNmFWSBliroBNNRB+W4o2wqlW9XHU03rZugDIWaoGhrW29TgCLv6n20huiuvBURaWhoOh4PS0lIsFgv5+fksXLiwxTkOh4OEhAQA1q5di81mA6CiogKTyYROp6O0tBSHw+HpyxCiS+gNEDdc3UbPU/ed+bo5LEq3wpa3YOPr6jFzfPMdhnWkOhpcFloS3ZzXAkKv15OTk8Ps2bNxuVxMnz6d5ORkFi1aRGpqKuPHj2fp0qVs2rQJvV5PeHg4L774IgBFRUW89tpr6PV6tFotzz33HGaz2VulCtE+4X3h1gfUDaChFk7YoXSLeqdxbCPsWa4e0wc1j8uw3qZ2gsvEhKKb8WofREZGBhkZGS32zZ8/3/P7Z555ptXrJk2axKRJk7xZmhAdpze6HzONbN5XXaYGRmmRGhqbfgMbXlWPRd4C8aOb35jqNUCmPxd+zeed1EL0KKY4dWtaaKn+ojoxYekW9bHUgdWw8wP1WKBJfSwVf5saGLHD1XmqhPATEhBCeFNAUMtxGU2jv0s2u0NjC3zxmXpMo1PfkPLcZYwGU6zvahc3PQkIIbrSpaO/hz2s7rtY2bwiX+kW2P572PKmeiw8rvkOwzpKnS5EJ39tRdeQP2lC+FpQBCRPUDcAVz0496jrZJRuabmUa0Cw+igqfrT7FduREGT2WemiZ5OAEMLf6AKg7zB1G/3v6r7qMvcdxlYo3QzrXgbFhTr1+UCIvw2TzgoWI0QkygJLolNIQAjRHZjiIC1b3QBqz6nLuDYFxp4V9K2thq3PQ6jFfYcxWv01erA8lhI3RP7UCNEdGUOhX4a6ATQ2cmRLPv30J92PpTbBPvdsyAEh6uy18e7AiBsJxjDf1S66DQkIIXoCrZZac5I6VcjIf1P3VR9X7y5K3IFR+GtQGkGjVUd6x49pDo3wvr6tX/glCQgheipTLJimN4/JqDnjnibEHRg7/gBb31KPmeObA8M62r2Eqwziu9lJQAhxswgMh6Tx6gbq21Ll9uY7jMNrwP5H97mm5j6M+NHQNx0CAn1Xu/AJCQghbla6APWV2djhMOaHzYP4mu4wSjbDwdXucw3qW1Xxo9U7DettEBzp2/qF10lACCFUlw7iG/qv6r7zp1sGxqbfqut+A/QZpIaF7Xb1Vxn13eNIQAghri4kCgZOVjdQ55Y6vg2ObYKSjeojqS/VdeQx25rDwnY7RCXJeIxuTgJCCNF+AUGQME7dAFwN4NzdHBgHP4NdH6nHQnq7H0ndrs5FJdOEdDvy0xJC3DidvnnUd1M/xjcH1bBoCo3iVeq5hjB1PinbGDU0YodLx7efk4AQQnQejQZ691e34d9X91UfV/swjm1Uf/3if9X9OoP6dlRTYMTfpr49JfyGVwOisLCQF154gcbGRmbMmMHcuXNbHP/oo4/48MMP0Wq1BAcH8/zzz5OUlATAW2+9xfLly9FqtTzzzDPccccd3ixVCOEtptiW04RcqFA7vpsCY+PrsP4VQKMO4LONIUxvA2sUhPbxaek3O68FhMvlIjc3l7y8PCwWC9nZ2WRmZnoCAGDKlCk89NBDABQUFLBgwQLeffddDh06RH5+Pvn5+TidTmbNmsXq1avR6XTeKlcI0VWCI2HAveoGUHfePd25+y5jx1Li6i/Axv8HUcmQMBZsY9WOb1Ocb2u/yXgtIOx2OzabDavVCkBWVhYFBQUtAiI0NNTz+4sXL6Jxv/FQUFBAVlYWBoMBq9WKzWbDbrczbNgwb5UrhPAVQ0jLeaVc9Rzd9GcSNV+DYwPsWQHbfqceM9vUsEhwB4bMXOtVXgsIp9NJdHS057PFYsFut19x3gcffEBeXh719fW8//77nmuHDBnS4lqn09nm96utraW4uLiTqveOmpoav6+xvXpKW3pKO6CHtSUkieLAVIicCENdGKsPEXxqJ8GndhC8/6/od30IQH1Qby70HqZufYZSF5bgd4HRnX8uPu+kfvjhh3n44YdZtWoVixcv5sUXX7yhr2M0GklJSenk6jpXcXGx39fYXj2lLT2lHdDT25IKPKD+trERvjkAx9YTcGwjJscGTCX/UI8F91LvLJruMvrc6vM5pfz959JWeHktICwWC+Xl5Z7PTqcTi8Vy1fOzsrL4+c9/fkPXCiFuIlot9BmobiNnN08RcmyD2ofh2ADFf1HPDTS5x2HcrgZG9BAZi3EdvPZfKi0tDYfDQWlpKRaLhfz8fBYuXNjiHIfDQUJCAgBr167FZrMBkJmZyVNPPcWsWbNwOp04HA4GDx7srVKFEN3ZpVOEpH9P3VdVooZFU2gc+Ju63xDqHosxVh3s1zcd9Abf1e7nvBYQer2enJwcZs+ejcvlYvr06SQnJ7No0SJSU1MZP348S5cuZdOmTej1esLDwz2Pl5KTk7n33nuZPHkyOp2OnJwceYNJCNF+5nh1G/Iv6uez5S0D44vn1f36IDUwEu5QAyM2HfRG39XtZzSKoii+LqIz+PtzPugeNbZXT2lLT2kHSFuuy/nT6ihvxwZwrFenCwGvBIa//1zaqk8exgkhbj4hUerqeylT1M8XKtz9F+vVbY17tLc+8LLAGH5T3WFIQAghRHAkpNynbtAyMI6thzW/AJTmwLC5JyyMG9GjA0MCQgghLtdaYJRsct9hrIO1C/AERtzI5juMHhYYEhBCCHEtwZEwMEvd4LLAWN9KYLjvMGJH+LTsjpKAEEKI63V5YFysVKc399xh/BJQQGckPioVTt6j3mXEDu9Wr9VKQAghREcFRbRcee+SwNDt/wzWvKDuDwhWF1FKuAMSMyDGvwfu+W9lQgjRXV0SGEdtj5BiszTfXRwthILn1POM4e5R3ndA4h3qqns+nhrkUhIQQgjhbcGRMOjb6gZw7mRzWBxdBwf+ru4PilBHeSfeqW69B/p08kEJCCGE6GqhfSB1urqBuuqeY50aFo5C2P+puj+kt9rZnXgnJNypTifShYEhASGEEL5milWnBWmaGqTS4Q4L913G3hXq/rC+6qOohDvU0IiwebUsCQghhPA3EQnqlv6IOlvt6cPqncXRQjhUAPY/queZ49U7i1FzoO/QTi9DAkIIIfyZRgO9ktRtxA/UwDi1391/UQhf5YMuAPq+2unfWgJCCCG6E40G+qSo222PqoHhpX4J/3mfSgghxPXzYqe1BIQQQohWefURU2FhIS+88AKNjY3MmDGDuXPntjiel5fHsmXL0Ol0REZG8otf/ILY2FgAUlJS6N+/PwAxMTG8+eab3ixVCCHEZbwWEC6Xi9zcXPLy8rBYLGRnZ5OZmUlSUpLnnJSUFD7++GOCgoL48MMP+fWvf82rr74KQGBgICtXrvRWeUIIIa7Ba4+Y7HY7NpsNq9WKwWAgKyuLgoKCFueMHj2aoKAgAIYOHUp5ebm3yhFCCHGdvBYQTqeT6Ohoz2eLxYLT6bzq+cuXL+fOO+/0fK6trWXatGl85zvf4fPPP/dWmUIIIa7CL15zXblyJXv27GHp0qWefWvWrMFisVBaWsrMmTPp378/8fHxV/0atbW1FBcXd0W5N6ympsbva2yvntKWntIOkLb4q+7cFq8FhMViafHIyOl0YrFYrjhv48aNvPnmmyxduhSDwdDiegCr1cqoUaPYt29fmwFhNBr9emFw8P/Fy69HT2lLT2kHSFv8lb+3pa3w8tojprS0NBwOB6WlpdTV1ZGfn09mZmaLc/bt20dOTg6LFy8mKirKs7+6upq6ujoAKioq2L59e4vObSGEEN7ntTsIvV5PTk4Os2fPxuVyMX36dJKTk1m0aBGpqamMHz+eX/3qV1y4cIH58+cDza+zHj58mGeffRaNRoOiKMyZM0cCQgghuphX+yAyMjLIyMhosa8pDAB+97vftXpdeno6q1at8mZpQgghrkFGUgshhGiVBIQQQohWSUAIIYRolQSEEEKIVklACCGEaJUEhBBCiFZJQAghhGjVNQPiV7/6FefOnaO+vp6ZM2cyevRomYZbCCFuAtcMiA0bNhAaGsratWuJjY3ls88+49133+2K2oQQQvjQNQOioaEBgLVr13LPPfcQFhbm9aKEEEL43jUD4q677uKee+5h7969jBkzhoqKCoxGY1fUJoQQwoeuORfTf/zHfzB79mzCwsLQ6XQEBgayePHirqhNCCGED13zDuLBBx/EbDaj0+kACA4OZs6cOV4vTAghhG9d9Q7i1KlTOJ1Oampq2LdvH4qiAHDu3DkuXrzYZQUKIYTwjasGxPr16/nkk08oLy9nwYIFnv2hoaH8+Mc/7pLihBBC+M5VA2Lq1KlMnTqV1atXM2nSpK6sSQghhB+4Zh9Eeno6//3f/83s2bMBOHToEMuWLWvXFy8sLGTSpElMmDCBt99++4rjeXl5TJ48mSlTpjBz5kyOHz/uObZixQomTpzIxIkTWbFiRXvbI4QQopNcMyB+9rOfMW7cOE6ePAlAQkICv//976/5hV0uF7m5uSxZsoT8/Hw+/fRTDh061OKclJQUPv74Y1atWsWkSZP49a9/DUBVVRVvvPEGf/rTn1i2bBlvvPEG1dXVN9I+IYQQN+iaAVFZWcnkyZPRatVT9Xq95/dtsdvt2Gw2rFYrBoOBrKwsCgoKWpwzevRogoKCABg6dCjl5eWA2v8xduxYzGYzJpOJsWPHsm7duutunBBCiBt3zXEQwcHBVFZWotFoANi5c2e7RlM7nU6io6M9ny0WC3a7/arnL1++nDvvvPOq1zqdzja/X21tLcXFxdesy5dqamr8vsb26ilt6SntAGmLv+rObblmQPz0pz9l3rx5lJSU8C//8i9UVlayaNGiTi1i5cqV7Nmzh6VLl97w1zAajaSkpHRiVZ2vuLjY72tsr57Slp7SDpC2+Ct/b0tb4XXNgLj11ltZunQpR48eRVEUEhMTCQgIuOY3tVgsnkdGoN4VWCyWK87buHEjb775JkuXLsVgMHiu3bp1a4trR40adc3vKYQQovNcNSD+8Y9/tLrf4XAAMHHixDa/cFpaGg6Hg9LSUiwWC/n5+SxcuLDFOfv27SMnJ4clS5YQFRXl2T9u3DhefvllT8f0+vXrZeyFEEJ0sasGxJo1a9q88FoBodfrycnJYfbs2bhcLqZPn05ycjKLFi0iNTWV8ePH86tf/YoLFy4wf/58AGJiYnjzzTcxm8388Ic/JDs7G4DHHnsMs9l8nU0TQgjREVcNiEtHT9+ojIwMMjIyWuxrCgOA3/3ud1e9Njs72xMQQgghut5VAyIvL6/NC2fNmtXpxQghhPAfVw2I8+fPd2UdQggh/MxVA+Lxxx/vyjqEEEL4masGxDvvvMOcOXN4/vnnPYPkLvXMM894tTAhhBC+ddWAWLJkCXPmzMFqtWIymbqyJiGEEH7gqgERFRWF0+nkk08+4Q9/+INnwSAhhBA3h6sGxEMPPcT3v/99SktLmTZtmme/oihoNJorJt4TQgjRs1w1IB555BEeeeQRnn32WZ577rmurEkIIYQfuOa83RIOQghxc7r2wg5CCCFuShIQQgghWiUBIYQQolUSEEIIIVolASGEEKJVEhBCCCFaJQEhhBCiVV4NiMLCQiZNmsSECRN4++23rzheVFTE1KlTGTRoEH//+99bHEtJSeH+++/n/vvv59///d+9WaYQQohWXHUkdUe5XC5yc3PJy8vDYrGQnZ1NZmYmSUlJnnNiYmJYsGAB77333hXXBwYGsnLlSm+VJ4QQ4hq8FhB2ux2bzYbVagUgKyuLgoKCFgERFxcHgFYrT7qEEMLfeC0gnE4n0dHRns8WiwW73d7u62tra5k2bRp6vZ65c+dy9913X/P84uLiG663K9TU1Ph9je3VU9rSU9oB0hZ/1Z3b4rWA6Kg1a9ZgsVgoLS1l5syZ9O/fn/j4+KuebzQaSUlJ6cIKr19xcbHf19hePaUtPaUdIG3xV/7elrbCy2sBYbFYKC8v93x2Op1YLJbruh7AarUyatQo9u3b12ZA3ChFUfjuu1sASI01MTjWTFqsCWtkUKsr6QkhxM3CawGRlpaGw+GgtLQUi8VCfn4+CxcubNe11dXVBAUFYTAYqKioYPv27cyePdsrdSoKDLdFsvark7y3/ij1LnVhJFNQAIPjTO7QUH+Ni5DQEELcPLwWEHq9npycHGbPno3L5WL69OkkJyezaNEiUlNTGT9+PHa7nccff5wzZ86wZs0aXn/9dfLz8zl8+DDPPvssGo0GRVGYM2dOi87tzqTVavjxhP78eEJ/6hoaOeA8i72smt3Hq9h9vJol6454QiMiOEANjDgTabFm0uJM9DUFSmgIIXokr/ZBZGRkkJGR0WLf/PnzPb8fPHgwhYWFV1yXnp7OqlWrvFlaqwx6LanuuwVQH2fVNrj4qtwdGmXV7D5ezVv/PEJDoxoaUSEGT2g0/RodLqEhhOj+/LaT2l8Y9ToGx5kZHGf27Kupd7G//Cy7y9S7DHtZNb9d+w0ud2j0CjWQFmsiLU7tzxgcZ8ISHuijFgghxI2RgLgBgQE6hlrNDLWaPftq6l3sO3HGc5exu6yafx44iDsz6B1mJNGkY+xxPYPj1NCICjX6pgFCCNEOEhCdJDBAR3p8BOnxEZ59F+oaKHaHhr2smi+PnOTVggMo7tCINQe5w8LMkDgTqXEmwgMDfNQCIYRoSQLCi4INeobbIhluiwTU942t/ZLZ477D2FVWhb2smr/taX4duF+vENIuCY1BfcMJNsiPSQjR9eRfni4WatQzul8Uo/tFefZVXajDXlaN3R0YW45UsHLn1wBoNdDfEqb2ZVjV0BgQHYZRr/NVE4QQNwkJCD9gDjZwZ//e3Nm/t2ffyTM1zaFxvJqC/SdZtq0MAINOy8CYMPXxVKyZwVYTSb1D0etkTishROeRgPBTfcIDuXtQIHcPUkeUK4pCWeVFdh93P5oqrWbljq9ZurkEgKAAHbf2DVcfTVlNDIkzY4sKltdthRA3TAKim9BoNFgjg7FGBjM5LQaAxkaFo6fPex5N2cuq+XDrMd7b0AiAOTjA05cxJE690+gTJq/bCiHaRwKiG9NqNdzSO5RbeocydZg6dXqDq5EDznPsKqtiV2kVu8qq+c2aQ57XbfuaAhliNatbnDoaPNQofwyEEFeSfxl6GL1Oy6C+4QzqG85Do9TR4BfqGtj79RlPYOwqrfK8OaXRQFLvUHdgmBhiNTMwOhyDXvozhLjZSUDcBIINekYmRDIyIdKzr+J8HfayKnaVqn0aa/afZPklneApfcMZ6g6MIVYziVEhaLXSnyHEzUQC4iYVGWLgWwP68K0BfQC1E/x41UVPYOwqrWLZtjLe33QMgDCjnsHWpvEZZkJqGnxZvhCiC0hACEDtBI+LCCYuIpiswWonuKtR4dDJS/szqninsHmiwuh/nFSnHIlXpx1JizURIv0ZQvQY8rdZXJVOq2FAdBgDosP4zgh1bfGaehd7vz7DP778ihN1RnaWVvH3vWp/RtOgvqZ5qobGm0nuE4ZOHk0J0S1JQIjrEhigY7gtguALJs8yiqfP1bKrrIqdJVXsdE8d8n9FpQCEGHSkxZkYao1gqNXMsHizzGwrRDchASE6LCrUSOZAC5kDmwf1Hf3mPDtLqzzbu+ubF16KDg+UR1NCdANe/VtZWFjICy+8QGNjIzNmzGDu3LktjhcVFfGLX/yCr776ipdffpl77rnHc2zFihUsXrwYgHnz5jF16lRvlio6kUajoV/vUPr1DmVaujo+o2k69J0lzaFx+aOpYe7AGGKVR1NC+AOvBYTL5SI3N5e8vDwsFgvZ2dlkZma2WDo0JiaGBQsW8N5777W4tqqqijfeeIOPP/4YjUbDtGnTyMzMxGQyeatc4WWtTYd++lwt9rJqdrgD46+7y/loa/OjqcFx6iOpYfHq46neYbJ+hhBdyWsBYbfbsdlsWK1q52ZWVhYFBQUtAiIuTv2/S6225aCs9evXM3bsWMxmMwBjx45l3bp13Hfffd4qV/hAVKiRuwb24a6Bza/aXv5o6u1L3pqyRgYxzBrhCY1BMTKgTwhv8lpAOJ1OoqOjPZ8tFgt2u/2Gr3U6nW1eU1tbS3Fx8Y0V20Vqamr8vsb28mZbUoIgpb+Oh/pHUdsQwaGKOvafqmH/qVo2HnTyl13qVOgBWg1JUQYG9g5kYG8jA3sF0jtEd10TFMrPxD9JW/xDj+kZNBqNnrdq/FVxcbHf19heXdmWoZd9PlF9kZ0lVeworWJHSSV/PVDNin3VAPQJM3ruMIZZ1bmm2lpwSX4m/kna0nXaCi+vBYTFYqG8vHmlNKfTicViafe1W7dubXHtqFGjOr1G0T3FmIKISQviXvestvWuRvafOMuO0kp2lKihsXqvesep02oYGK12gDc9nkrsFSLToAvRDl4LiLS0NBwOB6WlpVgsFvLz81m4cGG7rh03bhwvv/wy1dXq/xWuX7+eH//4x94qVXRzATotaXEm0uJMfG+Muq9pbIYaGFX8+ZK1M8zBAeqYDGsEUVwgNrFe1gIXohVeCwi9Xk9OTg6zZ8/G5XIxffp0kpOTWbRoEampqYwfPx673c7jjz/OmTNnWLNmDa+//jr5+fmYzWZ++MMfkp2dDcBjjz3m6bAWoj0uH5vhalQ4fOocO0oqPaHxzwMHUBT4n4JykvuEqm9Z2dQ3rfr1kskJhdAoiqL4uojO4O/P+aB71NhePaEtZ2vqWbXBzjeEs72kku3HKjnjnoTQFBTAsHgz6fERDLdFMMRq9vt1M3rCz6SJtKXrtFWff/+JF8KLwgIDGNY3mJSUZEBdoe/IN+fYfqyK7SWVbDtWydqvTgHNg/ma7jDSpS9D3AQkIIRw02o1JPUJI6lPGN8ZqY7fqb5Yz87SKrYfq2R7SSWrdn7Nh1vUvozIEAPDrGbSbWrn95A4s0wZInoU+dMsRBtMQQFk9O9NRv/eQPMU6E2PpLaXVFKw/yTQ/MaU2pehPp6KjwyWuwzRbUlACHEdLp0CvWlJ16oLdewoUR9LbS+p5JPtZfxhs7rQUq9QA8Pc/RjDbRGkxZoIDND5sglCtJsEhBAdZA42tJgyxNWo8FX5WU9gbD9WyWf71HEZAToNqbEmRrgDI90WQZ8wmf5c+CcJCCE6mU6rYVDfcAb1Dee7o20AfHOulu3H1I7vbccqeX/TMd5ZdxSA+MhgRrjDYkRChMxkK/yGBIQQXaBXqJGJt0Yz8VZ1jrHaBhd7jp9h27EKth2rpPDgKT7ZcRxQ1/8eZotgeLwaGN3hFVvRM8mfOiF8wKjXefolQJ3JtqTiAtuOVfLlMfWx1KsF6kA+rQZSYsI95w+3RRBrDpLOb+F1EhBC+AGNRoMtKgRbVIhnkaWmV2y3OSrYVlLJ8m1l/H6T2vkdHR7YIjAG9Q33Zfmih5KAEMJPXf6KbYOrkf3lZz39GNuOVZK/+wQAgQFakiMNZJRoGZGg9mfI/FKioyQghOgm9DotqbEmUmNNzLw9AVCnPm8Ki/X7v2bxPw/jWqOg0cAASxgjEyIZkRDBiIRIYs1Bvm2A6HYkIIToxmJMQdw3OIj7BvelOElLfL9kdpZW8aWjki+PVbQYk9HXFMiIhEhGJkQw3BbJgGh5W0q0TQJCiB4kxKhnbFIvxib1ApofS33pqKDoWCVbjp72rMgXZtSTbovwBMZQq5kggwziE80kIITowS59LPX9sYkoikJZ5UW+PFZBkaOSLx0VvPQPdUJCvbZ5EN8I96OpXqFGH7dA+JIEhBA3EY1GgzUyGGtkMFOHqW9LVV2oY3tJJUWOSrY5Kvn95mMsWa8O4uvXK4ThtghPX4bMYHtzkYAQ4iZnDja0WFxJHcRX7b7DqOSzYifLtpUBEBViYESCGhijEiMZFBOOXqf1ZfnCi7waEIWFhbzwwgs0NjYyY8YM5s6d2+J4XV0dP/nJT9i7dy9ms5lXXnmFuLg4ysrKmDx5MomJiQAMGTKE3Nxcb5YqhHBTB/FFMtwWCRnN62QUOSopclRQ5KjwrPkdYtC5+zHUwBhqNctkhD2I1wLC5XKRm5tLXl4eFouF7OxsMjMzSUpK8pyzbNkywsPD+eyzz8jPz+ell17i1VdfBSA+Pp6VK1d6qzwhRDtduk5G0wy25dU1bHVUUHRUDYxXPldHfQfoNAyOM7sDQ+38NgXJeIzuymsBYbfbsdlsWK3qwitZWVkUFBS0CIgvvviCxx9/HIBJkyaRm5tLD1kBVYgeLdoUyLeH9OXbQ/oCUH2hni+PVXhCY8m6I7z5z+bxGKMSIz13GZZwmb22u/BaQDidTqKjoz2fLRYLdrv9inNiYmLUQvR6wsLCqKysBKCsrIwHHniA0NBQnnzySUaMGNHm96utraW4uLiTW9G5ampq/L7G9uopbekp7QDft6Uv8EACPJAQSU2Dma9O1bLnZA17nTX8qajEM01ITJie1D6B3GoJJNUSRN8w/RUd375uS2fqzm3xy07qPn36sGbNGiIiItizZw+PPfYY+fn5hIaGXvUao9Ho1wuDg/8vXn49ekpbeko7wP/aMuyS39e7Gtn39RmKHBVsdT+W+uzwOQB6hxkZ6e74HpkQSUpMOAe+2u9XbekIf/u5XK6t8PJaQFgsFsrLyz2fnU4nFovlinNOnDhBdHQ0DQ0NnD17loiICDQaDQaDAYDU1FTi4+M5evQoaWlp3ipXCOFFATotQ6xmhljNzL6jH4qicPjUObYcberHqOSvu9V/L8KMegb0CmC808CoxEgGx5kIkDelfMJrAZGWlobD4aC0tBSLxUJ+fj4LFy5scU5mZiYrVqxg2LBhrF69mtGjR6PRaKioqMBkMqHT6SgtLcXhcHj6MoQQ3Z9G09zx/fBt6qJKx6suUnRU7cdYt/8EL/59PwBBATrSbWZuS4ySN6W6mNcCQq/Xk5OTw+zZs3G5XEyfPp3k5GQWLVpEamoq48ePJzs7m6effpoJEyZgMpl45ZVXACgqKuK1115Dr9ej1Wp57rnnMJvN3ipVCOEHYs1BxA6L5YFhsRQX6+lt7UfR0Qq2uLemN6UMOi1DrWZGJUZyW79I0uMjCJEFlbxCo/SQ14b8/TkfdI8a26untKWntAN6fluqL9SrfRgONTD2HK/G1ah4pgi5LVF9S2pEgn+9WuvvP5e26pPYFUJ0C6bgAO4eZOHuQWpf5rnaBra7JyDcerSCvA0O3io8gkYDKdHh3NYv0h0aUUSGGHxcffckASGE6JZCjXru7N+bO90LKtXUu9hRUuUJjI+2lpC3wQFAcp9QbuunhsVtMhaj3SQghBA9QmCAjjG3RDHmligA6hoa2X28is1H1FdrV2w/ztLNJQAkRAUzKjGS0f2iGN0vir6ymFKrJCCEED2SQa/1zCn12F3q2hj7Tpxh69EKNh9R55P605fqJITxkcGM7ieBcTkJCCHETUGv0zI4zszgOHUsRmOjwv7ys2w+cprNR05fERi3Nd1h3BJ10y7XKgEhhLgpabUaBvUNZ1DfcH4wLvGKwLh0mnNrZBCjE6NuusCQgBBCCFoPjK+cVw+M25oCo18kcRHBPq7eOyQghBCiFVqthpSYcFJiwpk1tmVgbDlSQUGxk+XuwIiLCPL0X9yWGIk1smcEhgSEEEK0Q2uBceDkWTYfPs3mywIj1hzkubvordTjv8Pk2iYBIYQQN0Cr1TAwOpyB0eF83x0YB0+e8zySWvPVST7e7n4kteYbxvSL4vZbejHmlqhuMw5DAkIIITqBVqthQHQYA6LDmHl7gicwVmzcy5Fz+hZvSfXrFeIZszG6XxS9Qo0+rr51EhBCCOEFTYFxf4qJlJQUXI0KxSfOsOnwaTYdOc3KnV/zwRZ14N4AS1hzYCRGYQr2j7mkJCCEEKIL6NyTCqbGmphzZz8aXI3sPl7NpiOn2XT4NP9XVMLvNjrQaGBQTDi3uwNjZEIkYYG+CQwJCCGE8AG9Tsuw+AiGxUfww28lUdfQyK6yKjYeOs2mI9/w/sZjvLPuKDqthrRYE2NuieL2W6IYYYskyNA162FIQAghhB8w6LWeZVfnk0xNvYvtxyrZdOQ0Gw+f5p3CIyxee5gAnYahVjNjbunFmH5RDIv33gJKXg2IwsJCXnjhBRobG5kxYwZz585tcbyuro6f/OQn7N27F7PZzCuvvEJcXBwAb731FsuXL0er1fLMM89wxx13eLNUIYTwK4EBOm5P6sXtSb14Cjhf28CXxyrZePgbNh8+zRtfHOS1goMY9Vrm353MD7+V1Ok1eC0gXC4Xubm55OXlYbFYyM7OJjMzk6Sk5kYsW7aM8PBwPvvsM/Lz83nppZd49dVXOXToEPn5+eTn5+N0Opk1axarV69Gp5NlBoUQN6cQo56M/r3JcE9vfqamnq1HKth85DSJUSFe+Z5eWwncbrdjs9mwWq0YDAaysrIoKChocc4XX3zB1KlTAZg0aRKbNm1CURQKCgrIysrCYDBgtVqx2WzY7XZvlSqEEN1OeKC6gNIz9w3i3rQYr3wPrwWE0+kkOjra89liseB0Oq84JyZGbZherycsLIzKysp2XSuEEMK7ekwndW1tLcXFxb4uo001NTV+X2N79ZS29JR2gLTFX3XntngtICwWC+Xl5Z7PTqcTi8VyxTknTpwgOjqahoYGzp49S0RERLuuvZzRaPTrhcHB/xcvvx49pS09pR0gbfFX/t6WtsLLa4+Y0tLScDgclJaWUldXR35+PpmZmS3OyczMZMWKFQCsXr2a0aNHo9FoyMzMJD8/n7q6OkpLS3E4HAwePNhbpQohhGiF1+4g9Ho9OTk5zJ49G5fLxfTp00lOTmbRokWkpqYyfvx4srOzefrpp5kwYQImk4lXXnkFgOTkZO69914mT56MTqcjJydH3mASQogu5tU+iIyMDDIyMlrsmz9/vuf3RqOR1157rdVr582bx7x587xZnhBCiDZ47RGTEEKI7k0CQgghRKs0iqIovi6iM+zcuROj0T/nVBdCCH9VW1vL0KFDWz3WYwJCCCFE55JHTEIIIVolASGEEKJVEhBCCCFaJQEhhBCiVRIQQgghWiUBIYQQolUSEJ2ksLCQSZMmMWHCBN5+++0rjtfV1fHkk08yYcIEZsyYQVlZGQAbNmxg2rRpTJkyhWnTprFp06auLr2FG21Hk6+//pphw4bx7rvvdlXJV9WRtuzfv58HH3yQrKwspkyZQm1tbVeWfoUbbUt9fT3/9V//xZQpU7j33nt56623urr0K1yrLUVFRUydOpVBgwbx97//vcWxFStWMHHiRCZOnOiZ6NNXbrQdxcXFLf5s/fWvf+3Ksq+PIjqsoaFBGT9+vFJSUqLU1tYqU6ZMUQ4ePNjinKVLlyr/8z//oyiKonz66afK/PnzFUVRlL179yrl5eWKoijKV199pYwbN65La79UR9rR5IknnlCeeOIJZcmSJV1Vdqs60pb6+nrlvvvuU4qLixVFUZSKigqloaGhS+u/VEfa8pe//EV58sknFUVRlAsXLih33XWXUlpa2qX1X6o9bSktLVWKi4uVp59+Wvnb3/7m2V9ZWalkZmYqlZWVSlVVlZKZmalUVVV1dRMURelYO44cOaIcPXpUURRFKS8vV8aOHatUV1d3ZfntJncQnaAjy6sOGjTIs9ZFcnIytbW11NXVdXkboGPtAPj888+JjY0lOTm5y2u/XEfasmHDBgYMGMDAgQMBiIiI8Olswh1pi0aj4eLFizQ0NFBTU0NAQAChoaG+aAbQvrbExcUxcOBAtNqW/zytX7+esWPHYjabMZlMjB07lnXr1nVl+R4daUdiYiIJCQmAuiZOZGQkFRUVXVX6dZGA6AQdWV71UqtXr2bQoEEYDAbvF92KjrTj/PnzvPPOOzz++ONdWvPVdKQtR48eRaPR8G//9m9MnTqVd955p0trv1xH2jJp0iSCgoIYN24cd911Fz/4wQ8wm81dWf4Vdd7ocsL+tBRxZ9Vit9upr68nPj6+M8vrND1mydHu7uDBg7z00ku89957vi7lhrzxxhvMnDmTkJAQX5fSYS6Xi23btrF8+XKCgoL4/ve/T2pqKmPGjPF1adfNbrej1WpZt24dZ86c4V//9V+5/fbbsVqtvi7tpnfy5EmefvppXnzxxSvuMvyFf1bVzVzP8qpAi+VVAcrLy3n88cd58cUXffp/Eh1px65du3jppZfIzMzk/fff56233mLp0qVdWv/ldd5oW6Kjoxk5ciSRkZEEBQVx5513snfv3i6t//I6b7Qtn376KXfccQcBAQFERUWRnp7O7t27u7T+y+u83uWEO+PaztbRWs6dO8ejjz7Kj370o6tOlOcPJCA6QUeWVz1z5gxz587lqaeeYvjw4b4o36Mj7fjwww/54osv+OKLL5g5cyaPPvoo3/3ud33RDKBjbRk3bhwHDhzwPLsvKioiKSnJF80AOtaWmJgYtmzZAsCFCxfYtWsX/fr16/I2NGlPW65m3LhxrF+/nurqaqqrq1m/fj3jxo3zcsWt60g76urqeOyxx7j//vu55557vFxpB/m2j7znWLt2rTJx4kRl/Pjxym9/+1tFURTl1VdfVT7//HNFURSlpqZGeeKJJ5S7775bmT59ulJSUqIoiqL85je/UYYMGaJ8+9vf9mzffPNNt2vHpV577TWfv8WkKB1ry5///Gdl8uTJSlZWlvLiiy/6pP5L3Whbzp07pzzxxBPK5MmTlXvvvVd55513fNaGJtdqy65du5Q77rhDGTJkiDJq1Chl8uTJnmuXLVum3H333crdd9+tLF++3Cf1N7nRdvz5z39WBg0a1OLv/L59+3zWjrbIdN9CCCFaJY+YhBBCtEoCQgghRKskIIQQQrRKAkIIIUSrJCCEEEK0SgJCiOt05swZPvjgAwC2bNnCo48+el3Xf/LJJz6bIkKI6yEBIcR1OnPmDB999NENX79ixQpOnjzZiRUJ4R0yDkKI6/SjH/2IgoICEhMT0ev1BAcHExERwYEDB7j11lt56aWX0Gg07Nmzh1/+8pdcuHCBiIgIFixYwPbt2/nZz35Gnz59CAwM5I9//CNLlixhzZo11NbWMmzYMHJzc9FoNL5uphAyklqI61VaWqpkZWUpiqIomzdvVtLT05UTJ04oLpdL+c53vqMUFRUpdXV1yoMPPqicPn1aURRFyc/PV376058qiqIo3/3udxW73e75epWVlZ7f/+d//qdSUFDQdY0Rog0ym6sQHTR48GDP1M8DBw7k+PHjhIeHc+DAAWbNmgVAY2MjvXv3bvX6LVu2sGTJEmpqaqiqqiI5Obnd8/oI4U0SEEJ00KXrd+h0OlwuF4qikJyczB//+Mc2r62treW5557j448/JiYmhtdff93ny5sK0UQ6qYW4TiEhIZw/f77NcxITE6moqGDHjh2Aujb0wYMHr7i+KQwiIiI4f/48q1ev9mLlQlwfuYMQ4jpFRESQnp7Offfdh9FopFevXlecYzAYeO211/jf//1fzp49i8vlYubMmSQnJzN16lSeffZZTyf1jBkzuO++++jVqxdpaWk+aJEQrZO3mIQQQrRKHjEJIYRolQSEEEKIVklACCGEaJUEhBBCiFZJQAghhGiVBIQQQohWSUAIIYRo1f8HyIqpzNZ/4B0AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] @@ -127,22 +128,60 @@ "source": [ "seaborn.lineplot(\n", " data = {\"theta\": theta, \"fillets\": fillets}, \n", - " x = \"theta\", y = \"fillets\"\n", + " x = \"theta\", y = \"fillets\", legend = \"full\"\n", + ")\n", + "seaborn.lineplot(\n", + " data = {\"theta\": theta, \"fillets2\": fillets2}, \n", + " x = \"theta\", y = \"fillets2\"\n", ")" ] }, { "cell_type": "code", - "execution_count": null, + "execution_count": 13, "id": "f587e6da-c921-4533-8077-260f122636da", "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "(array([0.14459953, 0.13429238, 0.1237727 , 0.11303387, 0.10206896,\n", + " 0.09087075, 0.07943172, 0.06774402, 0.05579944, 0.04358943,\n", + " 0.03110503, 0.0183369 , 0.00527525]),\n", + " array([0.40411255, 0.3938054 , 0.38328573, 0.3725469 , 0.36158198,\n", + " 0.35038378, 0.33894475, 0.32725704, 0.31531246, 0.30310245,\n", + " 0.29061806, 0.27784993, 0.26478828]))" + ] + }, + "execution_count": 13, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "fillets, fillets2" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "52d9bd9a-b2ed-48aa-8f5c-dd980ea1ff3a", + "metadata": {}, + "outputs": [], + "source": [] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "8f04f541-29a8-42cd-99d6-d0d1cd766f7b", + "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { - "display_name": "Python 3 (ipykernel)", + "display_name": "Python 3", "language": "python", "name": "python3" }, @@ -156,7 +195,7 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.9.6" + "version": "3.9.7" } }, "nbformat": 4,