{ "cells": [ { "cell_type": "markdown", "metadata": { "editable": true, "slideshow": { "slide_type": "" }, "tags": [] }, "source": [ "# Running simulations in parallel\n", "\n", "In this tutorial we show how the number of threads can be specified when simulations are executed in parallel. We start by defining a simple system object." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import discretisedfield as df\n", "import micromagneticmodel as mm\n", "\n", "import oommfc as mc\n", "\n", "region = df.Region(p1=(-50e-9, -50e-9, 0), p2=(50e-9, 50e-9, 10e-9))\n", "mesh = df.Mesh(region=region, cell=(5e-9, 5e-9, 5e-9))\n", "\n", "system = mm.System(name=\"mindriver_steps\")\n", "\n", "system.energy = mm.Exchange(A=1e-12) + mm.Zeeman(H=(0, 0, 1e5)) + mm.Demag()\n", "system.dynamics = mm.Precession(gamma0=mm.consts.gamma0) + mm.Damping(alpha=0.5)\n", "system.m = df.Field(mesh, nvdim=3, value=(1, 0, 0), norm=1.1e6)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We are now going to drive the system using `TimeDriver` for $2\\,\\text{ns}$ and save the magnetisation in 200 steps. When we pass those arguments to the `drive` method, we can also specify `n_threads`. `n_threads` is a positive `int` which specifies how many threads are going to be used in simulation." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Running OOMMF (ExeOOMMFRunner)[2023/10/18 12:39]... (1.2 s)\n" ] } ], "source": [ "td = mc.TimeDriver()\n", "td.drive(system, t=2e-9, n=200, n_threads=4)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We can finally plot the magnetisation." ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "editable": true, "slideshow": { "slide_type": "" }, "tags": [] }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA5EAAADaCAYAAAA/pcaNAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy81sbWrAAAACXBIWXMAAA9hAAAPYQGoP6dpAAB8g0lEQVR4nO3dd3gU5d7G8e/MbN9NsumBhBo60puiYkfUY0FBOCpWBGxYj92jKCrYULFgOwqIYkVfz1Es2BuIKKJSpAuEkt6zbd4/NgnEtM1mW/D34doryczszsPu7O7c8zRF13UdIYQQQgghhBAiAGq0CyCEEEIIIYQQou2QECmEEEIIIYQQImASIoUQQgghhBBCBExCpBBCCCGEEEKIgEmIFEIIIYQQQggRMAmRQgghhBBCCCECJiFSCCGEEEIIIUTAJEQKIYQQQgghhAiYhEghhBBCCCGEEAGTECmEEEIIIYSIOevWrWPMmDH07NmTTZs2Bf045eXlzJo1i1GjRtG/f39OOeUUXnvttRCW9O/HEO0CCCGEEEIIIcSBFi1axKxZs0hISGjV4/h8PqZNm8Zvv/3GDTfcQHZ2NkuWLOHf//43JpOJsWPHhqjEfy8SIoUQQgghhBAxY8WKFcyePZs777yTnJwcnnjiiaAf6/3332f58uU8+uijnHTSSQAMHz6cXbt28dNPP0mIDJKESCGEEEIIIUTMcDqdLF68mD59+jB37txGt/vqq6+YN28ev/32G4qi0KtXLy677DJGjRpVu80777xDRkYGY8aMqXPf+fPnh638fwfSJ1IIIYQQQggRM3r06EGfPn2a3Obzzz/n0ksvxW63M3fuXB599FESEhKYOnUqX3zxRe12q1evZvDgwSiKEu5i/61ITaQQQgghhBCiTXnggQfo0aMHTz75JEajEYDDDz+cU089lTlz5nDUUUdRXFxMcXExGRkZLFq0iAULFrBz507S0tI477zzuOCCC9A0Lcr/k7ZJaiKFEEIIIYQQbUZOTg6bNm1i9OjRtQESwGAwcPTRR7N27VoqKyspLy8H4MMPP+Sjjz7i1ltv5bnnnuOwww5j9uzZPPzww9H6L7R5UhMphBBCCCGEaDP27NkDwNy5cxvtM7l3716sVisAbrebZ555BovFAsBhhx3G3r17mT9/PpMnTyYpKSkyBT+ISIgUQgghhBBCtDkXXXQRp59+eoPr0tLS0DQNTdPo27dvbYCscfjhh/Pll1/yxx9/MGLEiEgU96AiIVIIIYQQQgjRZrRr1w4Ar9dL7969m9y2W7du5Ofn11vu9XoB6jSHFYGTPpFCCCGEEEKINiM9PZ3s7Gw+/PBDXC5XnXXPP/88r7zySu3fJ598MmvWrOGPP/6os90XX3yB1WqlV69eESnzwUZCpBBCCCGEECJm7NixgzVr1rBmzRr27t0LwMaNG2uXuVwurr/+evbt28dFF13EV199xffff8/MmTN58MEHqaioqH2s8847j44dOzJ58mQ++OADvvvuO2655RaWL1/O5MmTsdls0fpvtmmKrut6tAshhBBCCCGEEAA333wzS5YsaXT9smXLyMrK4quvvmLevHn8/vvveDwesrOzufDCCznjjDPqbL9v3z4eeughvvjiC0pLS+nSpQvnn38+48ePD/P/5OAlIVIIIYQQQgghRMCkOasQQgghhBBCiIBJiBRCCCGEEEIIETCZ4kMIIYQQQggRUb7dPaJdhDrUjA3RLkKbIiFSCCGEEEIIEVE+fNEuQh3SPLNl5PkSQgghhBBCRJRb98bULVDr1q1jzJgx9OzZk02bNgV0nxUrVnDeeecxcOBABg8ezJQpU1i3bl2wT11MkBAphBBCCCGEiChfjP0LxKJFixg/fjylpaUB/z9//PFHLr74YqxWK08++SSPPvooxcXFnHfeeezYsSPYpy/qpDmrEEIIIYQQIqLcemw1Z23OihUrmD17NnfeeSc5OTk88cQTAd3v0UcfJSUlhSeffBKTyQTAIYccwrHHHsvTTz/NvffeG85ih43URAohhBBCCCEiyoseU7fmOJ1OFi9ezLhx4wL+PxYWFvLDDz9wwgkn1AZIgKSkJA4//HA++eSToJ67WCAhUgghhBBCCBFRPvSYujWnR48e9OnTp0X/xw0bNqDrOj161B+Jtlu3bhQWFpKTk9Oix4wV0pxVCCGEEEIIEVFuvfng1tbl5eUBkJiYWG9dzbK8vDzatWsX0XKFgoRIIYQQQgghREQF0oS0rauqqgLAaDTWW1ezrLKyMqJlChUJkUIIIYQQQoiIch/8GRKLxQKA2+2ut65mmdVqjWiZQkVCpBBCCCGEECKivCjRLkLYpaSkAFBQUFBvXW5uLgCpqakRLVOoSIgUQgghhBBCRJTvb1AT2bNnTzRNY/369fXWrV+/ntTUVNLS0qJQstaT0VmFEEIIIYQQEeVCjalbOMTFxTFy5EiWLl1ap+/jnj17+O677zjppJPCst9IkBAphBBCCCGEiCifrsTUrTk7duxgzZo1rFmzhr179wKwcePG2mUul4t33nmHPn368PXXX9fe7/rrr6ekpIQrrriCb7/9ls8++4xp06aRkJDA1KlTw/b8hps0ZxVCCCGEEEJElAst2kVokSeeeIIlS5bUWTZ9+vTa35ctW4bP58Pr9eLz+WqX9+7dm5deeolHHnmEyy+/HE3TOOyww5gzZ05tn8m2SNH1v8EkLUIIIYQQQoiYsXxbl2gXoY4RnbZEuwhtitRECiGEEEIIISLq7zA668FMQqQQQgghhBAioty6xJC2TF49IYQQQgghRERJTWTbJiFSCCGEEEIIEVFuvW0NrCPqkhAphBBCCCGEiCivzDTYpkmIFEIIIYQQQkSU9Ils2+TVE0IIIYQQQkSUV5c+kW2ZhEghhBBCCCFERElz1rZNQqQQQgghhBAioqQ5a9smr54QQgghhBAioqQ5a9smIVIIIYQQQggRUVIT2bbJqyeEEEIIIYSIKB9SE9mWSYgUQgghhBBCRJRXl4F12jIJkUIIIYQQQoiIcutatIsgWkFCpBBCCCGEECKiZIqPtk1CpBBCCCGEECKipCaybZNLAEIIIYQQQoiI8ulqTN0OJu+88w45OTlNbvPVV1/xwgsvBL2Pg+sZE0IIIYQQQsQ8L0pM3Q4mt9xyC2vWrGlym02bNjFv3ryg9yHNWYUQQgghhBAR5fZJDAmldevWsW7dOgB0XWfFihWUl5c3uG1lZSVvvPEGHo8n6P3JqyeEEEIIIYSIKJknMrS++uorHn74YRRFQVEUFi1a1OT2uq4zZsyYoPcnIVIIIYQQQggRUW6fDKwTSpdeeiljx47l559/5sorr2TMmDF069atwW01TaNDhw6ceOKJQe9PQqQQQgghhBAiomSKj9BLSUnh+OOPZ9iwYZx33nkMHTo0bPuSECmEEEIIIYSIKJ8uzVnDZeHChWHfh4RIIYQQQgghRETJPJHhlZ+fz5dffsmePXtwu90NbqMoCldccUVQjy8hUgghhBBCCBFRUhMZPj/88APTpk2jvLwcXdcb3U5CpBBCCCGEEKLNkJrI8Hn44Yepqqri9NNPZ8CAAVgslpDvQ0KkEEIIIYQQIqKkJjJ81q5dy3nnncfNN98ctn1IiBRCCCGEEEJElE+X0VnDxWg00rt377DuQ0KkEEIIIYQQIqLcEiLDpnfv3vz5559h3Ye8ekIIIYQQQoiI8ulqTN0OJldffTWvvvoqGzZsCNs+pCZSCCGEEEIIEVGeNhjc1q1bxyOPPMKPP/6Ix+OhX79+TJ8+neHDhzd6n0mTJrFixYoG102bNo1rr7025OX8/fffGTFiBGPHjuXQQw+lS5cumEymetspisK//vWvoPah6E2N+yqEEEIIIYQQIfbP76dEuwh1vHros02u3759O2PHjqVr165cfvnlWCwW5s+fz9dff82iRYsYMGBAg/ebNGkSZWVlzJgxo966tLQ00tPTQ1L+A/Xq1QtFUZqc3gP8IXLt2rVB7UNqIoUQQgghhBAR5fG1rSk+nnrqKbxeL8888wxJSUkADBkyhNGjRzNnzhxeeumlRu9rt9vp169fhEoK999/f9j3ISFSCCGEEEIIEVE+2s4UH7qu88knnzBy5MjaAAlgMpkYPXo0CxcupLi4mPj4+CiWcr+xY8eGfR8x2xh53bp1jBkzhp49e7Jp06YG10+ZMoUhQ4YwYMAAzjvvvEbbGwshhBBCCCFih09XYurWlF27dlFSUkL37t3rrevevTs+ny+sg9i0xrZt2/jmm28oKioK6ePGZIhctGgR48ePp7S0tMH127dv59xzz6WgoICHHnqIefPm4XA4uPjii1m9enWESyuEEEIIIYRoCY9PjalbU/Ly8gBITEyst65mWc02DSkoKODmm2/m6KOP5pBDDuGUU05h0aJFrXj2mvftt98yZswYxowZw+TJk+uE3PPPP58vvviiVY8fcyFyxYoVzJ49mzvvvJMJEyY0uM2BbZKPOeYYDjvsMB5//HFSUlKYM2dOhEsshBBCCCGEaIlo1zy2pCbS5XIBNDjCqdFoBKCysrLR++/YsYP27dvz8MMPM3fuXLKysrj77rt54YUXWvEMNm7NmjVMmTKFPXv2MHLkyDrr8vPz2bhxI1dccQU///xz0PuIuT6RTqeTxYsX06dPH+bOnVtvfVtrkyyEEEIIIYSoqy1N8WE2mwFwu9311tUETKvV2uB9586di8FgwOFw1C47+uijmTBhAo8//jgTJkyosy4UnnvuOZxOJ6+++ioJCQl1piBJSkri3XffZfz48Tz//PM88cQTQe0j5l69Hj160KdPn0bXt+U2yUIIIYQQQoi2VROZmpoK+Gvx/io3N7fONn/ldDrrhURFUTjuuOOorKxk48aNwTx9TVq1ahX//Oc/6dChA4pS//+WmprKxIkT+fHHH4PeR8zVRDantW2ShRBCCCGEENHVXHCLJRkZGSQmJrJ+/fp669avX4/RaKRHjx4N3tfn8+Hz+TAY6saumuavNbWcoVRYWEjHjh2b3CYzM5OSkpKg99HmQmRr2yTX8O2uX5MpwKv7ol2EOnw0PUlqJLl1b7SLUOvMziOb3yhCVIc92kWoy9DmPtZElCmxdMyY63+3RYseQ2XxxTfcTCxaHn/rmWgXoVaqGjsn4mYldt5LRqVtzUEYSeZ2m6NdBAC8zQxmE2tOPPFElixZwr59+2prHcvLy/noo48YNWoUdnv986Ht27dz8sknM2nSJG666aba5V6vl08++QSn00m3bt1CXtaEhAT27dvX5DabNm3C6XQGvY+29erRujbJQgghhBBCiOjzocTUrTmXX345VquVadOm8fnnn/PNN99w+eWXU1FRwXXXXQf4Bwjt06cPr732GgAdO3bkhBNOYP78+Tz88MN89913LFu2jClTprBhwwZuuOGG2kqwUBoyZAivvfZaozWNy5cvZ8GCBQwZMiTofcTOJaMAtaZNshBCCCGEECL62lpNZHp6Oq+88goPPvgg1113HbquM3DgQBYsWFBbm6jrOl6vF59vf8u+2bNn06dPH9544w1efPFFjEYjffr04emnn+bYY48NS1mnTp3KxIkTGTduHMceeyyKovD+++/zxRdfsHLlSlavXo3RaGTatGlB76PNhcjWtEkWQgghhBBCRF9b6hNZIzs7m3nz5jW6fsSIEfUyislk4tJLL+XSSy8Nd/Fq9e3bl6effprbb7+dF198EYBXX321dn1GRgb33nsvvXv3DnofbS5EQnBtkoUQQgghhBCxQW+DIbItOeKII/j4449ZuXIlf/zxB2VlZTgcDnr06MHQoUPRtNb1G465ELljxw4KCgoA2Lt3LwAbN26kvLwcgJ49e3L55ZezdOlSpk2bxlVXXYXRaOS5556r0yZZCCGEEEIIEZu8PgmR4WY0GjnssMM47LDDQv7YMRcin3jiCZYsWVJn2fTp02t/X7ZsGVlZWc22SRZCCCGEEELEpkAGsxGts2vXLvbs2YPH42l0m2HDhgX12DEXImfNmsWsWbOa3a65NslCCCGEEE3RdZ2i0j9BUVAVDVXRUFQNVTEc8Hv1ckVFUdrWQCBCxLK2NrBOW7J7926uuuoqfv3112a3Xbt2bVD7iLkQKYQQQggRCYqiUFq5j7Vb3wtwexVF0Ti0nwejScFoBKNR8d9MNb+DyaTUXW+CxCSVqVc6SE6R+QuFANBjZyrwg84999zDmjVryMzMpH///lgslpDvI+gQuWbNGlauXMmOHTvIz8+nuLiY+Ph4kpKSyMrKYujQofTr1y+UZRVCxBBd1/HgwuP14NIrcahOjKq5RY/h9lVhUEwoSmw0afHpPva5tmFRHVg0BybFErWy6bpOnnsnJsWMRYvDqJij+jwVuHdjUExYtTgMSujntGoJj+4mz7UDu+bEpsWjxsCk4nsrt+LV3cQZkrEZEmKiTAWVuyis2kWCOYN4UxoG1RTtIqHrOpvzv8NuSibRmoXZEP2B8LJSh1DlKmHzrs+b3VbXfei6j6IigMDPgMf8w8LkaRIgReh5PDput47V2vZq9WRgnfBZsWIFxx9/PHPnzg3buUOLQqTH4+HVV19lwYIF7NixA/B/ITREURQ6dOjABRdcwIQJEzAYpNJTiFil6zpePLipwqVX+X9ShVt3HfD7/p9uqnDjQkdHKVboaR2B05zW4v3ucW/j9/KvMSs2zKoNs2rHotpq/7ao1csVW0TCpqqolHjy+Ln8I//faFg0B1Y1DovmwKI6sNb+jMOi2sMWFhRFocpXxo8l/wNAU4xY1bjq/cdh1fy/+8sWF/bA6/ZVsaL4XQBMihWbFl9dhnj/72ocNi0es2oP++tkUIzsc23n58qPAAWrGoddS8BmSMCuJVSHywQsqiNiwdtucPLNvsX48KKg4jAk4jAkEWdMxmFIJs6QhEWLi+iFgARzGr/mfsT6gq8ABYcxiQRzBgnmDJymDBym5IiHXUVRMGk2Vu96B6A6THYg0daBJGsHLMa4iJanRtf2R1HlLmHnvh9D+rjtszT+PTOeo44NfS2AED/9UsWsxwpZ9EzLv39jgTRnDR+Xy8Vxxx0X1u+cgJPdL7/8wo033sjWrVuxWCwcfvjhHHHEEXTp0oWkpCTi4+MpLi4mPz+fLVu28NVXX/Hjjz9yzz338PLLLzN79mz69+8ftv9INO3a7aF9hoRkEVt0XaeUwv2h8IDw99dlLqrQ8TX/oH9hwsLAuONJNKQHVcYscw8KPbvZ6fqDSm8ZePc1uq2GoTpo+kNlbcBUbVgUe+3vWitb6WfbBpPvzqHAvQsfXsq9RZR7i8Dd8Pb+sFsTKveHTIsWh1V1tCr8tjf3IM+1g5yqjXh1N6XefEq9+Q1uq2I4IFT6y1MTOq1qHCbV1qovkzRzJzpY+vJn5W+49ApcngoKPXsaKIe2P1yqB4TM6mWhqsXsZT+MfNdOKnwlVPiKqfAVg/vPemWxVYdKuyGh+vcEbJoTkxrak3q7wUn3uENZX/INOj5KPHmUePLIqfyjdhtNMRJnSD4gXPp/mlRrSMtSQ1UM9Ek+lh/2vAXolLrzKHXnsbP0t+r1GvGmNH+wNGXgNGdgNSSEPehmJfRne8EqSl37KHPlUebKY0fRzwBYjU4SrVkk2TqQaO2I1Rj+8oA/3PbudAoudyn7CuvPQ91SmgbHnWgmo53G8m9dfPuVC92n49PB5/PfdB/4fHrt3z6f/3P7wPXe6m30mm30+vezGxX+dV0cw4ZEv6b572Znjof2GVrEW4mUl/uY+XAhjz1bxH23J2EwtM0aPWnOGj4dO3akyN9kImwUvbGqxAO8+eab3HXXXTidTi6//HLGjh2L1dr8l15FRQVvv/02Tz/9NIWFhcyYMYOzzjorJAVvLd/u7iF5nLx8L1Nu2Mtb/2kXkseLNq/e8iARTr4WNBcKN7fujXYRap3ZeWRA2+30bWat9wf0MDyPiUoq/bSRWONSWvU4Xt3D9yXvNRqOWsqgmKsDpj9YWjT/zzRTZyyaI6DHqPSW8V3Bm7j0ylaXp6YG0V+j6TigRjOOBENqszVBHp+Lbwve8oekVlBQq4NlTcCtDphaHBY1Dotqa3bQEK/u4fuCtyn1FgRdDpNixarFVQfLeGxaHFY1PqhazAJ3DisK/y+ochgVS3WgTMBuSsSuJWI3OLEZEtCU4C5E6LqP5XlLKHTvbtH9TKrVHy6NycRZ0nCYUnAYkzGooQncq/e9T05ZYMHIqFqqQ2U6CY4OxFsywtLkNK9sKyt3vNbsdmaDw19TmZSNM64TdktKyE/Ydd1HcVkOecWbyC3cQFHZjlY9Xv+BRmbMSqBTFwNTL8jnh+9dISppw/uaN8dJj+7BHys+n47LBW6PjtsNLpeOx1O9zK3j9viXud3+ZR6PjssNbpdO924GevbYv29zkO+dQHi9Ort2e9m6zcPWPz3s3OVl0gQ7me0b3qcxDLXsXq/ODz9V8f4n5bz/cTmXTopn6oXxId9PU5Z9WcGVN+WyZZuHOIfCxpUdSYhvWY2eud3mMJWuZXq9fXe0i1DHujP/He0ihMzrr7/OggULeOuttzCbW9bVKFABhchevXoxbtw4brnlFuz2ln+ZlJeXc++99/L2228HPQJQqIUiRFZU+Dhxwi6cTpX/W9A+BKUKDV3X+WZFJcu+KufGKxJb1E4+3CHS69UpK9eJjwusTG0tROq6TmWlTkVl3Z+VVfWX1flZUX/bhHiFm65JIDmp/hdhoCESYJ9vJ794v8VH6EJwJ7UX3dT+qIqK6mj9CWaZt4jvit/F21h1XyskGNLobh9GsimrRffb59rOqqIPQl4egGRjFtn2ISQaMwLavsi9j+WF7wRVW9wcs2qni20gWZZeAYWnEk8e3xcsCenxVENFI8vSi272YQH3r91QupwtFT+HtBwWzVEbKpPNWaSZuwQcXEo9BXy777WQPD9WQwJxxhQcpuTan3ZjYouboFZ6Svlq53y8enBhxmpMIMHSrvYWZ0kPSf/KVTveYl/Zxhbdx2iwkRjXqfbmsKYHNWJqZVURecWbyCveRH7xZtyeihY/xl/ZHQrX3hTHxPNsaJr/eCkr9XHJefmsXhXazzajCa66Lo6LpthpZ6r7/3/9rXIWLCqvDX77Q2B1KHSDy63jdvnXeYM4VFUVrr3KwfVXx2Ey7X9vtCZE6rpOXoGPrds8bNvuD4rbtnvZut3Dtj89/LnTg7v6aYxzKCx4JoXjj268QiNUIbKk1MeyLyv430flLF1Wzr48/+fwtAvjePS+1l1EbYm8fC833Z3Py6+X1i6bPiWeB+5KbvFjxUqI7BljIXL9QRQi161bx3/+8x9+/PFHJkyYQOfOnRsNk0cddVRQ+wjo3X7PPfcwfvz4oHYAYLPZuPfeexk4cGDQjxFrfD6dC6bv4ZsfKrn60oRoFweAzdvcvPxmCQvfKGbLdg8fvd4+Kh2ty8p9bN7mYcs2N1u2u9m81cOW7W42bfVfNXv9+YyAQ2RDdu/1sHaDG49Hx+v1XxX1eDngdx2PBzxeHZ+X2t89nv3rvV7wVG9vMChcdlE86amBfeF8s7yS2+4ubDQUhsKJx1mYeXtygwGypVLVTIZwDD95v8RD666Iaxjoq40gXe3Q6nLV8OlePLqbVGMWu91bQva4Di2J7vZhpJo6BRwAdN1HubeYEm8+pZ58TIoVl976k8saLQmPuq5T5SujxJtPmacAq+ag3Nu62sgDWVQ7XWyDyLL0CiiUuHyVlHkLKPUUYtMSQlZzXCPV1JGutsE4jc03jfb4XJR7iyjzFoVlyoVKbykGxUR7Y09SzU0fP17dTYW3lApvCZXeEio8JZg1OxUheK0qPEVUeIrYW7EJAINqpmPcALomDK9XS6nrOl7dhctbQZW3HLfP/9PlrcDlK8es2Sj3BPf+r3AXUeEuYnfJuuolCg5zygHBsj0OcwqqovqbY+pevD4XHl8VHp8Lj89V5++a34N57dyecvYWrGVvgf+CtEGz4HR0rA2VcbZ2qGr949nrdVFQso284o3kFW2irDI3qOeiMSeebOHWu+JJy6i7b7tD5dn5SVw4MY+1vzU+T1tL9O1n5P5HEujes+Hax3FjrXz1TRWvvh66z64Ddems8fRjiQwNovlsWbnPHxC3e9j2p5dt2z1sqQ6J27Z7KC1r/js0s53GmwtSOaRP+Jrvbtvh5v2PKvjfx+V8+V0Frr+8dY490sJDd7c8vAVD13Vef6eMG/6dVxtgwR/kL78ksrWgoab72mYz3LbgjDPOQFEUdF1nzpw5TW4b1ik+WhMgw/E4seDGu/N4679lAGR3jt5IhUXFXt54r5SX3yjhq+X7m97dfm0ixx5hC8s+dV1n914vW7Z52LTNzZZt7trQuHmbhz37Gr6secwRFl57Lh1nQuuCUUKcyuzHC/n869Y3NRzc38Tzj6UGHCABDh9hYeJZdv51R/BN+hpjsSjcd4eTyReEdkAQp5rCMOVYfvAswxNkbZ+deAYYjsCuBP+l5dFdlHjyKfbmUuzNp9ibR5m3IKTNbW1aAt1sQ8kwZzf6HPoDWjml3nxKPP7AWOrNp9RTEJYatmRjZnV4rN/sXdd1XHoFpZ6C6nLs/+kJsuaoKRbVQVfbIDItPeuFR13XqfSVUeYtoMxTSKm3gDJvIWWegpA07W1IurkrXW2DiDfUvaLv0d21/VHLPUWU+4op8/j/DmWw/6tEUzu62AeTau4EgFuvpMJTQmV1UKzwllDpKfH3xfSW4PaF53k5kEWLI9XaBYcpBa/uYmPht/5w6C3H5av+6a0Iy7HbMJ3Sqn2UVu1jd8m66lpjpTocusJSY94Yj7eS3KIN5BZtAEBTjSQ4OuB0dMRhS6e8Mo/8ok0UlG5HD0OXhHaZGnfcE88xxzfexzY+QWXuc4mcOSaX4uLgP+uMRrj8GgeXTHNgNDb+/aCqCo8+6MTjgTfeDu175aLzbdx1ezx2W9MXAHLzvfx3aUVtYKypTdyX27pjY8AhRt6Yn0q7EI9DUdNM9X8f+5up/rau8e/JHtlGFj2TFpF+iNt3eJh+Sy5Ll9V/HceeYqdzh+iOlN1aMjpr+NSEyHAK+l344Ycf8tlnn7Fv3z7c7obfbIqiMH/+/KALF6vmPl/InGcKa/+OdIj0eHQ++bKcBW+U8O7Ssnq1X0eOsHDHdUkh3++u3R7OmbqHn9a4qGhhjdt54x3MezC1TrOXYFmtKm+9lM6p5+zm2xVVQT2GwQC3XOPkxqucTX4ZN2baxXGoKlx/W+iCZL8+Rv7zZAq9eoT2eKrSK9jl28Iu35agA2S60pE+2rAWDYpS5Sun2OsfXKTYm0eJN5/yVvbra4pFsZNtHUR7W2/UA2o43L6q2mBWUhsW83HrwR07LZFkzKTbAeHR5aukzFNAqbfggLIU4A5TQDuQVY2ji20QmZYegEKFt7hOSCz1FlLmLcSrh75J8V8pKLQzd6OjtR+qolLuLSLXtWN/aPQWUeUrD3s5DhRvSCXemAaKzvbyNawv+ZZKbwlePTS1R61R6S3hz9Jfol2MBnmraxZjha7rVFQVUFy2C483fO8rTYMLJtu54loHtkYCVe4+L599UsWnH1Xy3ddVVLXiI6f3IQbuf9hJz96BfQZrmsLcR5x4vLDk3dYHyfR0lccfcnLcMYENSJWcqLIv18tDc0P3mT/meAsvPpWCwx6a1gclpT4++cJf2/jhAc1Um5LoVHl7fjqJzvCOauz16jz9YjF3ziqgrLzh862rp7btWkiQgXXCadasWWHfR1Ah8umnn+bxxx9vdHqPGrEy91soLXm/lGv/XbcZTKRC5Jq1VSx4vYRX3i5h996Gr6gmJaosfDI9LFfI2mcYePaRVM68cA8bNgV+onn7dU7uuD4xpMeD3abyzsIM/vHP3axY1bJv5j49jTz/WCqD+7euo/GUC+PQVLjmltYHyaumxnHnTU7M5tA8Rz7dR66+i52+zeTpOUHX9Cko9FAH0kHt0WStXoWvpDoo5lFcHRrDWVt0IJNioatlIJmm7lT4Stldtam2ZrHEk0+lr7T5BwmxREM72pm7oSgqe6q2sKlsFaXe/IgHIwCzYiPF1BGTaiHP9SfbKtZQ7i2KaG3Rfgp2LQGDYiLfvYtdVX80f5cIKfbso9jT+OjAIrZpqhmvrwqf7qGiKrjPZFUxkBjXiaSEbHbs/aHRx4m3t+elN1307lv/u3/zRg/LPqrk048rWb3K3eqTZIMBpk13MOWKpmsfG/PP8VY++bSSkpLgC3L6Pyw8eL+TpMTAw5uiKPxregK9exq59Kq8gJqpNmXqRQ5mz0is7WsarK1/unn/o/LqZqqVNFIH0iCDAV59Lo1uXcN7zvfrWheX3ZDLDz81fm5z6FAzwwe3/WljdJniI2JKS0spLy/HbrcHNb5NQ4IKke+88w6pqancdNNN9O7dO2yj/sSa71ZWcN4Ve+p8KagqdMoKf4h8flERV926r167/L96YU46HTLDV54e2Sa++V8mJ03MYeXPTYc3gwHmPZjK+RNCP+/Xrt0ebrorn7UbAr8Crihw3WUJ3HGDE4ul9R9cZeU+khI1OnXQ2PZncM2kMtI1nn00mWNGhebLoEwvrq11dNG6q/AmLPTXDidRTa1d5tN9lFFMiV6w/1ZUFJZml81RMZBoSMOgmNhRtY71FcvDMgptSxgVMwoqBZ4cCjw5US2LgoKOTpVezs6qdc3fISJ0yryF0S6EOGgoUP2e9/qCq+ZzWNNIjs8mOcE/+qumGvHpXjbt/LTetppqolvWcXRIG0bvvs/59+vV+eUnd21w3LIpdE1me/QycP8jTvoc0oIWIFU6X35dxX/fr+SDjyrJLwj+YlFCgsLsmQmcdYY16IvA/zjRxrL/MzDhon1s3d7y50ZR4P47nVw+OTRzrO7Z6+XVt8uaDGiNefS+ZI4+PDzT8QBUVvq4/7FCHn6yCE8zDSCunhobY3G0ltREhldubi5PPPEEn3zyCXl5ebXL09PTOemkk7jiiitwOAIbtb4hQYXI3bt3869//YtTTjkl6B23NX9sdnH6+Tn1mo52yjKEpIlmcyafm8CYY+zc9VAeL75a0uA20ycncNqJoR+K/UC6rnPx1XubDZDxcQqvPZfOcaPC0y/TGa/yv4/LKa8I7BOoa2cDzz+aysjhoQlruq5zxIm72bg5+KZup55kZe4DSa0ePMere9ij/8lO32YK9dDUpNRM32FWrOT4tlLg20sJhZTqhfjCVIOloGJWrFTqZQFt78NDnmdXWMoC/uafDkMShe49ATc3DVcTWZNixWFIpMxbRJUvsOcnXIFaw4Dd4KTSV44rCjWrB1JQsKhxmFQLRZ69US0L+C8iWDX/dC57q7ZGuzgYVSuaolHpjXxtfEM0xYRBM2FQTWiqGYNqIr98OwR9rLb8fkaDrTY0JsV3xWKq3ySwrCIXn6/uZ3taYm96djwJiyker8/NZx9XsuzjSj7/pIq8Vvbz+ytNgylXOph2lSPg84vPvqjktTcr+PCT1tU61jjqSBNzH06kffvWN9t0Jqgcc6SVFxe17Di0WhReeDKZU8eE7jxixBALX7zXjlfeKuX2ewsabdX1V1ddGs/k88LXfPSr7yu4/IY8/tjcfNVo544GTgvhcxJNMrBO+OzZs4fx48ezd+9ejEYjnTp1wmq1Ul5ezs6dO3nxxRf54osvWLx4MfHxwR3bQYXIdu3aYTCEby6gWLMv18sp5+aQ18AVva4R7A/pTFAbDW+D+5mZdXv4h5pWFIUuHZt+7Tu013hnYTv69Q7fyGk2m8roY6y8837zJ7HTLozj3tuTmh0IoCUURWH0MVY2bm440DfFZlV44J5Ezp/YsjnxDqTrOsV6Pjv1zez2bcNL4GFWRQWURgfhOHD6DoAc31by9JbNe9ccDSNxhiTitWTitWTitGQcmpMtlb+wsXJVSPfVHKNiIc6Q5J/4XfP/dGiJGFQTHt3Np7kvRawsBsWEQ0sizpCE3ZBYXZ5ETKoVXdf5PG9hxMpiVCw4DE7/VBeaE4fB/9Oi+q9afpG/KEIlUbCqDmw1czpqCdi0eGxaAlYtDlXR+LPid4pKwx0iFSyq3T/HplY912b1zVI932bNtBf5rl1hCZFG1YJJs2FSrZg06wG/2/x/q/6fZs2GUbWgKCq5FdtYueftoPfZUPCruR34t6aaMNT7e/8yTTXV+7wrKN/BivLwHkeKopHo6EhSQjbJ8dnE2ZqfEqSkfH8rAospgV6dTibBnkVu0Qb2Fqwnr3gTn14SXN/h+HilycF1uvUwMOsRJ337t+zc4sOPq3hzSeu7EFjMcOdt8VxyoR1Vbf0J/m9rXRwxZnezNWt/lZaq8vpLqQwZGPqWbqqqcN74ONJTNcaev6fZsp14rJX77wj9OBMAhUVebru3gBdeDvxc4srJ8a1u1hszpCYybJ566iny8/P597//zVlnnVWn1Wh5eTmvvPIKjzzyCM8++yw33HBDUPsIKgmec845vPHGG5x55pmYTOELCrHA59OZ8XA+6akazniVH3+pG+KyO0UuRDrsKocPt7BmresvyxVemZcesv50zblycgJPvFCMr4GLrwP6mnh3YQbtQzxyWkNOO8nOO++XY7cpDXY8z2qv8cwjqRw3KjzNT8aeauOpF1oWIgcPMPH83GS6Zwd33OTle3nt7TK+9yyllKIW3deBk0y1K+3UTqzwfEI5dcve2PQdcUpiq0KkSbEQpyUTb9gfGG1qfIMBOte9s94yBYVMUw8qfKXkeeqvD5SKAYch0R8Yq4Oaw5CESWm8qVaea0eDfQftWiJl3uD7wmqKEYeWiMOQiKM6KDq0JMyqrdGyFHn2NtjP1K45qfKVB92k2KI6sGtO7AYnjgMCo0lt/H1T5N7bYI2oAXP1NDItPzOwVAdFe3VYrAmKNi2+2SlIdldtavH+/kpFazAY+n+Pw6LZA56fMadiQ0DbGRUzRtWKWbXuD4ZGR51AWBMQjaq1zoBRgTJrdpIsHSiq2tPIPJEKTlMGabZsEi1ZGA8Igqol+AtdgdhW8ENYHtduSamuacwmKa4zmtay85Tisl2AQqf0Q+maeTSVrmK++Pkhgj3jze5u4LjRZo4dbaGqSueCs+tPj6OqMPkyO1dcE4cpiO/yf5xs4fmXAmul0JgB/Y3Me9xJ926hO6/p08tIt65G1m0IPHT36mHkzQWpdOoQvvOIdz8oY+Lkvc02p+zdw8iCp8I3Euu6P9zYbQoD+pr45XdXs+VJiFe5YGLouwhFi4zOGj5ff/01//znPznnnHPqrbPZbEyePJnt27fz8ccfRzZEnn/++RQVFTFmzBhOO+00MjMzMRob/tA544wzgipYrFBVhSfu9/cJm/FQXm2INBj88w926xLZkVkfviuFTVvdHDbUwt0P+09in5qdRveukQvznTsYeeiuZEYMMXPEP3bVfuideIyVV55JJ84RmY7Sp5xg48M3M3C74R//rBtyJp3t4KG7k0mID19Zhg8xkdVeY8eu5pvDKApcf2U8t16fENTgCCt/quKJZ0t4b2l5s/1iD2TASIbaiUy1K3HsH9yoirphpKnpO+KUxID3Z1Xj/IFRSyLe4A+MZqXxYHQgt6+KQm/d2qR2pmyyLYOxa/H8UfFjgCHSP3iLo7pWsSYs2tS4Fs9Lt8+1rc7ficZ2dLMNxaRa+abg9Wbvr6JhNyRWB8Yk4qp/WtSWT+Hy15otu5ZItm0wGeaurCx6n/wGAngNBaU6oDkPqFVMxG5wtmjE3Rp7qurO52lWbXS2DiDL2pvvC95utN+jRbX/JSD6Q6NVi0cLcpLyKl8F+e66/U+TjVl4dDdFnj21y4yK2R8M1bi//HRgNTkxqcH3+zqQT/eyu2JjnWUm1UaauTNpli5YNIc/MKqWhkOpObSf5XGmFHokHs73OYvrLLcbEsmK60d7R2/MWsNN4/QwBUiXt4Kc4t/ZUxpY2G6OUbOSFN+V5OraRou5df3FFEVhRJ8pxNv9oyrbLcmYjQ6q3IFdNFRVGDTEyLGjLRw72kLnLvuP7cUL6we9rt007n/ESf+Bwb/2h40wkZKsktvACKMOh8JhI0x8vKzhlkyaBtde5eD6q+OC+n5qiqIoXDkljitvCGxe2aOPMLPw2VScCeE9jzj+KCudOhjYur3xasjkRJW35qeH9Tzi0KEWDh3q72aTX+Dl4un7GpzOo8Yl58ZF7BwrIiREhs3evXvp27dvk9sMGDCAd955J+h9BPWt/dNPP/HKK69QUFDAM8880+A2uq6jKEqbD5EHmj7ZicWs8sb/lfLAncmcfM4uukawJhLAYlFZujiT4hIfdz9cwIUT4zj3rMhflbpycgKVlb7aADn5vDgevy8lIvMm1UiIVzlqpJVnXto/hHh6qsaTDyTzjzD3DQX/BYYz/mHj0y8rOfMfNh59urjBEeiy2ms8NzeZIw4Nvj/m+j/cvPV/gfc/S1RSyVSzSVOy6p2ce3VPneavzU3f0VCIVFCwE0+ckki8Jb26hjGper644Pj7N/qfvzRjJ7pZBxOn7W9CFKfVL4dZsfmbfGr+oBanJWI3pwQdSA6k6zr7qrYD+8Njkqk94A8KCmptLaWC6g+u1U1ha37atPigJlRvSE2gdWhJZNsGk27uWht6HJqTfPdOVAzYDQnYtUQcmrM2wAZSm9cSe11bgZppQwbS3tKj9jl3GtMxqdYDQuL+GkUtiMDabFmqtuA/bhQyzNl0sQ4g3phCgTsHj8/tb3qqOmqbmjZE0UJX47Gvaltt39g0c2c62QeQZMqM6mjlf5asAUBVNDJsPciKO4REc+TK5NO9FFbsIq9sC7nlWymubN2gU4qikmDPIjmhG8nx2cTb24XsfQbQo8OYOs+NoqikOnuyY9/KRu9jscDho/y1jUcfZyYpueH326Y/9n/2qipceKmd6dfHYba07rXQNIWTTrSw8BX/94Qzwf/3P062ctQRZn79zd1giMzuqvH0Y4kMHhS+C9ETxtqZMauQfbm+2gvwDTlvgp3HZiVFZJwJu03l6YdSOOls/wXov5bLaITFL6RF9Bxv2ZcVTQZITYPLLm7703rUIc1Zw8ZkMlFYWNjkNhUVFY1WAgYiqG/O+++/n8LCQg4//HD69OnztxmdNdGpcdNVidx4pRNFUXhsZmrE54isER+ncvLxNh6/N7X5jcOkZlCbmbck8a8rE6J2krR+o7+ZzJn/sPH4/SmkNPLlHQ4zbnFy/50KJaU+Zj5Uv3npWafZeHRWUquvqp7xDxv/uqOAktLGP3FNWGivdiFT7YpNafzCQlX1qK2BTN8BYMNBopKGXYknTnESRyIOJaE2NKiW0AT2XPcOkg3t6WYditNQ/7iO11LpYOpVHRaTsGuJmBoKrSEIkOBvPmo3OBlgO742PNZQFY2ejsP8TXUNSdi0hJCGtL8q9xYDCgPiTyDd1KXe69XZNoBOtv5Y1dCMYNiUMk8hOjqHxB1NO3O3ev/vQ+KODuv+/2qfazsdLH3pbOuPTdt/glUzN2ek7a3cSidbfzrZ+2MzRH8ERbevijJ3Pr2TjqG9vRdGLfLTAhRV5PDDn6+06jFs5qT9TVTjO2MI4/+jofdQamL9EJmYpHLM8f7gOPJIM1Zr8++9jdUhsnNXjfsedjJoSOjC27ixVlQVTj3ZwuGHmevUKjY0CNwlF9i48/Z4bNbw1mxZLArXXxlPbp6PiePsDD2q/kWEf9+YwA3TG+7mEC7HHGHlikviSU5SMZsUbrt3fxeFJ2ancOSh4RuJtSEDDzGT3cXApi0Np+yzTrXTIfMgG49EBtYJm+7du7N06VLOP/98VLX+e9zr9fLf//6X7t27B70PRW9usscGDBw4kAkTJnDLLbcEveNo8+0O/kmr8zg+PSSdz4NRWuYL2aS7Nbx64KPM7dnn4bOvK5k4NvjhgZvjC+Ay1aRpezl1jI3xp4e3/45bb7zZ6prfXIwcvb9JbZxD4eF7k5h4VmBNOQMx/cb8eqPbKSikKO3JVLuSrLQLqM9UoS+X1d6v603fESzVEZoQWeotwNFAbWOLhWjQL4/uRsMQE/PdxlJZ3L4qDIoxpDU/reH2VbWqBhxACeFAcV7d07qa8BA3Z/XpPhSUoI4dPURl8elePt34OF5f4G3xDaqZJFsnUuxdSLZ1xpISnYsCNXw+D5//9AAmo4O0xF7c+8R6Bg4xtniAk6OH7+GkU61c/a84LK2sfayRGsA5yMxZxTz6hP/7IyNdZe4jTo45KvRB3NzMsZ+z20OPIftH1TaZYN6cZMafEfrWQ8YWXNh7531/H0nwTwV2X5gG0mlObp6X8Rfv4bsfqnjo7iRuu7eAqir/edA3H7RnyIDQVNqY220OyeO0Vqf/PBDtItSx7eIbo12EkHnrrbe47bbb6NevHxMmTKB79+61o7Nu2LCBV199lXXr1jF79mxOO+20oPYR1DedzWajd+/eQe3wYBOtAAmEPEC2VHqqIawBMlBPP5wS9edi89b9Vw5HDDXx/NwUOjczim1LXXCOvTZE2oirHiSnM2alZVdLjYqJQw0ntvh+4RaSABlCwfQVDJdYKktrA1uoxVp5QtGUOpSCGYwn9GXQSLZ1Ym/pH01speC0tCPZ3oUUexfiLXUvioVnYqHAqaqBkf2uxGz01/YPGb6l+Tv9hdut8/ATiQwZHvkBCTdV10SOPd3KAzMTSEyMznFx4CB4iU6VV/+TwuEjIl87/lc1NXynnGDjnluj912UkqzxwWsZfPhZBaefZOenNS4WvVHK4SPMIQuQMUWas4bNWWedxdq1a3n55Zf59ddfG9xmypQpQQdICDJEHn300Sxfvvyg6u8o2q5oB0iAzds8aBrcfG0CN1wVH5a+oYMHmLjhqnhOONbCg+OPCbpWqqHBc4QQIpySbZ3rhUirMYFkWxdS7J1JsnWKSlPblmhoTsmWMBqVqARIgNw8H88+mciZp0f34mF5dYjs2tnAmwtSgx6pPNQ6Zhk4pLeRl55Mjfr0GRaLyukn+Wtmp10Yz6I3SrlmavSbxoeDIs1Zw+r2229n7NixfPjhh2zatIny8nLsdjvdu3fn5JNPblVTVggyRN52221ce+213H333YwbN47MzMxGp/qwWmOrtkOIcNB1+PDtdEYMDd+VQkVRuPNmZ+3vQgjRViTbu6ApJpLsHUmxdSHZ3hmbMVE+yyLkpecSSU6K3HgBjSmv0Bkx1MTiF1NJiYHy1EhJUnl7fkbMjXw6bJCZs8+wc/IJDY+g3OZJTWTY9e3bt9lRWoMVVJ/I4cOHA1BS0vRw14qi8PvvvwdXsjALVZ/Ig01L+kRGQiB9IiOlqT6RNaMRR8qZnUdGbF/NCVWfyJAJYf828fcQyj6RrRbiPpGtEao+kTV8ujfoAah88bF1QfrxtxoemT4aAukTGSnN9Yncu89LfJwasv6gTWlJn8hYVlnpw2IJbbiNlT6Rnec9FO0i1LF1WnDzJcay/Px8Vq1aRW5uLsXFxTidTtLS0hg0aBAJCa2r4Q7qm9PhcNT5GQ2TJk1ixYoVDa6bNm0a1157bYRLJP7O5Gq6EEI0LZwjGIu2IS1VjoGWCnWAjCmxU09w0PF6vdx///0sXrwYr9dfCXJghYfRaOTCCy/kuuuuC3ofQYXITz/9NOgdhlLfvn2ZMWNGveVpaWlRKI0QQgghhBAiIBIiw+bZZ5/l5ZdfJiUlhWOPPZb09HQcDgclJSVs376dL774gueee46EhAQuueSSoPYRQ214Ws5ut9OvX79oF0MIIYQQQgjRAjKwTvi8/fbb9OvXj/nz52Oz1e9TW1payvnnn89rr70W+RDp8Xj46aef2LNnDx5PwxOjAjKCqxBCCCGEEKIuqYkMm927dzN16tQGAyT4uyROnDiRmTNnBr2PoELkhg0bmDJlCnv27Gl0m5p2txIihRBCCCGEEAdS2mCIXLduHY888gg//vgjHo+Hfv36MX369NpBR0N9v2AlJibi8zU9WKbP5yM5OTnofQQVIu+//352797NsGHDGDBgABZLdOZ2Kigo4Oabb+b7778nNzeXTp06cc4553Duuec2e9+mRtqMNJXYqc6PpdFQIbZepyq98Rr3SNObqP2PNF9FZbSLUIdqjaG55jQZRKJBaow9L7E0MFYMHTO6LXYmN/c4YmfUWoC4GDr7NSuxMdcixNaIqLFUFtEIPYY+ewOwfft2zj33XLp27cpDDz2ExWJh/vz5XHzxxSxatIgBAwaE9H6tcdxxx/H5559z9tlnN7rNl19+yejRo4PeR1AhcvXq1Zx11lnce++9Qe84FHbs2MHo0aN5+OGHKS4uZvHixdx9991UVlYG3b5XCCGEEEIIEWaxcy0mIE899RRer5dnnnmGpKQkAIYMGcLo0aOZM2cOL730Ukjv1xrXX38911xzDVdccQUTJ06kR48eOBwOKisr2bRpE4sXL8btdjN16lQqKirq3NdqDWxKpaBCpK7rDBs2LJi7hszcuXMxGAx1phk5+uijmTBhAo8//jgTJkyI6hQkQgghhBBCiIYpsTU1eZN0XeeTTz5h5MiRtUEQwGQyMXr0aBYuXEhxcTHx8fEhuV9rjRgxorY5a1Ozahx++OF1/lYUhd9//z2gfQQVIrt168a+ffuCuWvIOJ3OessUReG4445j9erVbNy4kYEDB0a8XEIIIYQQQohmtKGayF27dlFSUkL37t3rrevevTs+n48NGzYwdOjQkNyvtdLT00P6eA0JKkROmTKF++67jzPPPLNVHTJbw+fz4fP5MBjq/hcqK/39s8zm2OnLIYQQQoi/H7dbp6DAR0GBj/x8nYJCH4cdaiIx8SCeQF6IALWlmsi8vDzAP2DNX9Usq9kmFPdrraZqH0MlqBCpqipDhw7lxBNP5JRTTqFz586YTA13eg9kkJuW2r59OyeffDKTJk3ipptuql3u9Xr55JNPcDqddOvWLeT7FUIIIcTfk9frxuMpx+0u55tvqigo0MnPrw6ItUHRR0GBTmH1spKS/VUtVqvC7NnxEiCFqNGGBtZxuVwADeYdo9E/uFVNRVYo7tcWBBUir7jiChRFQdd1XnvtNZQGRrarmeIjHCGyY8eOnHDCCcyfPx+DwcDIkSMpLy/nlVdeYcOGDcycObP2hRFCCCGEaIyu6+TmrqWqqgi3uxy3u6z65/6bx1OO1+uqvc+EiS3bR9euGs8+46RXLzk3EaJGW6qJrGnh6Ha7662rCYoNDUgT7P1CweVy8dNPP7Fv3z48TYzqH+x0jK0KkdE0e/Zs+vTpwxtvvMGLL76I0WikT58+PP300xx77LFRLZsQQggh2gZFUTAaraxb9yYeT+hrBE4+2czDDyUQFyc1kELU0Yb6RKampgKQn59fb11ubm6dbUJxv9ZatWoVV1xxBYWFhY1uU1PhF9EQedVVVwW1s1AymUxceumlXHrppdEuihBCCCFaoMJVhNnoQI2Rufyczi4MGnQpq1e/hMtVEpLH1DS47dY4Lr3UFvUL70LEohiabrVZGRkZJCYmsn79+nrr1q9fj9FopEePHiG7X2s98MADFBYWMnLkSPr06ROWsWICCpGPPfYYV199dat3FqrHiXVVVTpms3xhCNHW+HQfqhIbtQU+3YeOD00J6lpfyHl8LgxqbEz47tO9uH2VmNTYOTkvdu/Dp/uwGxIwqpZoFwfwH0N/5H+NSbVgNyVhMyZhMybERHDbW7iODTs/wmJyYjMnHXBLxGZOxmpyoqqRLafD0Y7Bg6eyevWLVFS0bqCLtDSVp59yMmJEbLxnhIhJbag5K8CJJ57IkiVL2LdvX23tYXl5OR999BGjRo3CbreH9H6t8ccff3Deeedx2223hfyxawR0dvLcc8+xZs0aZs6cSUZGRot3kpOTw+23387y5csP6hDpcunc/VAB4061M7CfjA4rRGO2etai48OhJOBQnFgUe0yEgfUVyynw7CZOS8KhJRGnJRKnJWFSrBEvn4LCquIPqfCW4DAk4tASsWtJOAyJ2LWEiIfLHVXr2Fz+E3bNiV1z+stj8P9uVR0oEQzfCiprSj6j0L0bu+bEpiX4fxoSsGv+5yfSgdegmPgmdzFe3Y1RtVSXx+l/jgz+321aAgY1cn3iVEUlwZzBz3v+r3aZgorVmIDdmITdmIjdlFT7u0mLXCjvmDqc3OI/yCvZTIWrgLySTX/ZQsFqclaHyiTMJelYrcnYbMlYLImoaniOf6s1iQEDLmTlyieDbto6YriRp55ykp4e/bAuRCxrSzWRAJdffjlLly5l2rRpXHXVVRiNRp577jkqKiq47rrrAFixYgUXXnghd955JxMmTAj4fqFmNpvp27dvWB67RkCfwk899RTXXXcdJ5xwAuPGjeOSSy4hKyur2fvt2LGD559/nrfeeguz2czTTz/d6gLHqrUbXFx05T50He65pf4wvtGk6zpeLxgM0T9JF03LL/CxebOHzVs8bKq+xTlU7r4jPmz9aXy6Dw/uOjcvngP+9lQvq/u3BzdGTPTyjcChJrRon0lqOitcS9GrO0RoGLArCThUf6j0/+7ETGTDW3frEL4peptd3o11lpsUS51QGaclYdecYQ1yiqLQ1zGKbwrfpNxVzF62HbgWmxaPQ0vEUR0sHdXhKVy1TJ0s/dhbtY0CTw6Fnj111qlotUGuNmTWhCYl9KFJURQOiTuGb/Nfo9izj2JP/XmLTaoNu5aATXNi1xL8QU5zYtPiw/Ic2QwJ9Io/gt+KPsPtq6TQt5tC9+5625lVe51guT9ghqdcGY4etC/rw65S/+TROj7K3QWUuwv467NmUM11wqXN1w6bORmbOQktxKFNURT6djqd79bOw+2taGALnQpXQXXA3Ay5de6NxeKsDZVW64G3pKACpq7rFBVtJSfnR/buXYPPV38QjECcdZaF/v2M/PKLm8xMH5mZGvHxSkxcJBMi5rSxEJmens4rr7zCgw8+yHXXXYeu6wwcOJAFCxbUzgrhP+f24vP5WnS/UDviiCNYtWpV0P0dA6Houh7QS/jnn39yyy23sHLlShRFoWvXrhx++OF06dKF5ORk4uLiKCkpIS8vjy1btvDNN9+wefNmdF1n2LBh3HfffXTo0CFs/5GWqsrpGpLH8fl05r1YzK33FlBZqfPYfclMvTC+RY+hEp4vly3b3Sx6s5Qduzw8OTsFTWt+P74IvKP37PPy85oqjjzUgs3WdDBy696wlwf8r2Nevg9VheSkhk/kqvTGR7ZqiZISH5u2eNi8xesPi9WhcfMWDwWFdZ//Y48y88IzicQ56j5PEzuMbHIfefpuSihqIBDWD4W+INuTpNKePgzFZIkL6v4b3avZ4v21yW0MGP2hsjpcOqrDpUlpvLmgam1dU8J97j9ZVfpRAFsq2NX4A2ot/SHTojr2nzBqrQ8Ff1b8zu9lXwe0rYKCTUuoDpeJOAxJOLREbFpCSJrplnuL+bbwLbx64CfYFtVxQO1ldcg0JmNSW3+BYHflJlYXB/JaHUjBqsXVlsmmJeAwJ2EzOOu+dkHQdZ2fCt5nb+WWIO7tvzBgM1UHOGMidmMiNmMiVkNcq2p63d5Kvtkxn0pP8H39/LWCydgtKdjNydgsydjNyZiNca16zvYUrmX1ljeCvn99+wOm1ZqM3Z5KRsYQDIaGWwdVVRWze/dP5OSsbFUTVodD4ZFHEjj+ODNjx+ax+hdPnXWZ7TXaZ6pkZmpkttfIzKr+mamRnq5iNAb/HCZU17pXVel8t6KK9DSVtFSNpEQ14uHVGANNpWvEUllCbcn7pZw+xo6qBvf6qhl/hLhEwek1Y060i1DHujuvjXYRQiY/P5+LLrqII488kvHjx5OZmYnBEOKLgYGGyBoff/wxCxYs4KeffsLj8TQ6vYfBYGDQoEFccMEFHH/88SErcKiEIkTu2u1hyrW5fPKF/yqq1aKw9eeOJMS37Ms+lCGyqNjH2/8t5eU3S/nq+0qSE1VWfpJFZrvADpxQh8i8fC+rfqnix9UuVq2u4sfVVeTs8fL8Y6mcc5aj2fuHIkS6XDo5e7zs2u0lJ8fDzhz/77uqf+bs9rBrt5cTj7Py3OPJOOwNv34tCZHlFT62bvUeEBC91cHRw959gYW2CeOsPPqgs8GTi+ZCZJVewSq+oozigMscKAWF7vSnA91QFAUlyM7aPt3LctdSSvXCFt/XiBmH6g+VNbWWDiUBo2JudYgE+KXsc3Jcf21eFxgDRhyG6lBpSK79Pdjmlbqus7L4f+S7dwV1f/A3YbTXhMvqYOnQErFq8S0Ol39WruX30q+CLksNg2LyNz01OLFriTiqf7a0TL8ULSOnakOrywOgKQZ/zaXhwKaoiSQY0wI+Ga/ylvPNvldx+RqqXQuOgobd6PQHTGMSnRIGYTG07OJNXvk2fsgJZVjz01QTNnMyHVOH0T5pQItCi8tTTnllHr//+T9KK/eGtFyqaqBdu6F07HgkFkvd1kE+n4e8vPXk5PxIXt56Wlsd0rOngeeeddK1q/97dutWDyednFdnjsimywrp6dUBs+bWXiMzU6V99e9N1WYmHPDZcvOdBTz5vP9igdEI6aka6Wkaaaka6WkqGWka6akaadU/M9I10lJUrNbQtHSJpeAWrrL4fDq/rnOR3dmIvZkL4aGm6zozHspnzVoXb/2nXdCPEzMh8q4YC5F3HTwhsqKigvvuu48333yzye0UReH3338Pah8tjqQnnHACJ5xwAqWlpXz//fds3bqV/Px8SkpKiIuLIykpic6dO3PooYficDQfEtqqJf8r44obc8kv2B8IzjrN3uIAGQoej86yryp4+Y0S3l1aTmXl/i+uFx5LCzhAtlZRsY9Vv1SxanUVq35x8ePqKrZurxu8VJWAA2QgSkp97MzxkpPjZWeOh5zd/mC4M8cfDnfmeNmX23xou/HqeG67IaFFV/VcLp1t2z1s2uKtU5u4abOHXTmt6y1+7VUObr0x+Cv8ZsXKEP0ofuIrSihsVVkOZMVOP0YQryS1+rFURaOv8bA6zVoD5aaKAt8eCqjbrNKElThP4v6aOC2ReC2lxUGpu3Uoe13b8dLyJm0e3BR69vibfFbtX25V44gzJOHQkok3pJBm6hhQ7ZKOj67WgeS7cwj2ZFfHR6m3gFJvAbg21y5X0bBrTrraBpFu6tLs8ebTfSQbM7FrTsq8hUGVpYZHd1Hk2UPRAU1jNcVIJ2s/utgGN9t3UNd1qnxlpFu6srtqI3oIRmjw6h5KPLmUeHJRUMm09cJpSm/2edF1HbevkjJvIWWeQuIMKeS5/mx1eWofHy9l7kKc1kyy4vo1GyA9Phfl7sLqW0Ht7wpqSJ6n/RRS47vTOX0k8baGT2h9Pg/lVfmUVeVRXplHWVVe7e8NN2NtHYPBQmbmoWRljcRkqvs9U1a2h5ycH9m9+yfc7rKQ7G/sWAuzZ8XXaVXTubOBB2bHc9nlRQE9hs8HOTk+cnJ8rFzZ8GeOw6EcEDD9gfOoo8z061f3fXLP7U5W/eLiuxVVuN2wY5eXHbuavxibEK+QXh0sa3+ma6SnqqSnaWSk+YNncqIadO2Xruu4XFBZpVfffFRV6gf8rVNVpVPZzLKqKp1Dh1o469TQD0bSELdbZ9WaKr76voKvllfw9fJKrrokgbv+lRyR/ddwuXSm3LCXhW+U8OSs0E8JEQ1taZ7ItmbGjBm8++67ACQnJ2MyhX6cgKDThcPhiMkaxnArLvFx/R15LHy9tN66i88JrllfsNasdfHy6yW8uqSU3Xvrf0lcMzWBk4+3hWXfpWU+fv51f+3iqtUu/tjc9Am3qsILj6fyzzNbHiBXra7iv0srqsOhl127PezK8VJS2roryBaLwrxHkjjr9JZ/Gb3xdgX/urUQl6v5bQOlKDB7ZgIXX9D6L0eTYq4Okl9TROtGGgRII4s+DAmqj5uu67ipokIvo1Ivo0Ivrf2pouElNE2FXVSQ56kg35NDurEzDs1ZL0D6dB9VPn85Kn1lVPhKqfSVVd/8v7v1qkb2ELwKXwketwuHlkiCIaU2QPp0r78c3hIqfCV/+VlKlS80J7t/pWEg09KTjtZDsGv7+7R6dDcV3hLKvcVUeIsp9xVT7vXfKn0lLQ78gTCrdjrZ+pNl6Y1R3V+z7Q+K5ZR7iyj3FlLmLfL/7imi3FuML0THzYEMiokO9kPoZB+ARav7PvT4XJR5/EGxvDowlnuKKPMU4gnDMQOgKSY6JgygU8IQLIb9n50eXxVlfwmJNaGxyhueY6aGqhjITB5Ep7RDsZkT0XWdSlcxZVW5lFf6A2NZZS7lVflUuAqJRMcnk8lBVtbhZGaOwGDY3xrB46lk79415OSspLi4ZcHebE7A661qcIAdoxFm3BXPpEkNN8s+9VQr33zr4uWXQxOUS0t11q/3sH69h+7dNa69xkHfvvVP4YxGhQXzUjj8xJyAW70AFBXrFBV72LCx6feUwQDnT3Qw41YnzoTALs69+0EZF125j/KK1h8HJhPcfXMSY08Jz7kNQEWFj+U/VfLV95V8tbyC736orFP2Ky9J4M4bWn8htSWKir2Mu2Q3n37tP55OPCZ8//+IamN9ItuSzz//nEGDBjFnzhzS09PDso/YGDu+jfhmeSUXT9/Htj/rf8j27mHksGHhH5F1zz4Pi5eU8vIbpaz+rfH0MmSAmZm3hOZDrqLCxy+/u1i12sWP1TWN6/5w42vBFSRVhf88nsrEIAIkwCG9TTzzYimvvBG6k6P2GRqLX0xlUP/grs6cO9FGrx4GLpyST87u1l9Os5jhmScSOeUka6sfq4ZBMTJAH8lyPqGK4E5mVFR6MJBMGq+p0nUdF5V/CYhl1aGxlAq9DB/h79+qYaC9qRspxg7o+Cjy7GOPa2udsFilVxDpby676iTV1BGrFk+Vr4wN5Suo8JZS6SuhMkwhsTFmxUampReJhnRceiU5VRv9YbH65tJDXzvUGIeWRGfrAJLMmVR6S9lbtYUyb2F1aPTfvCHqi9wcs2qjk30g7W09cfsqKHTtptxT6C+Pp5AyTxEuX3lEygJg0mx0iOtPkrUDVd5ydhSv8QdGTyHlroKQNpcNlEG1kJ7Yh3hbO6rcJWzc9am/VrEqD2+QA9G0lsWSSMeOR5KRMQRN81/gas0gOYqikZrah3bthpKYmM333z9UL0SazQm88YbG4EF1vzsqK3X++MPDunUe1q1z88cfoT12u3Xzh8d//MPS5BgHGeka8+el8I+z9+IN4cfuocPMzLg1gZHDW9Zl4PST7Dx2v48p1+bSsg5UdfXuYWT+k6n07xvac62iYi/f/lDJV9/7axlX/FyJu5FDZtL4OObcnRLRvqZ/7nTzj/Ny+HWd/5yve1cjXTpGbpTncJKayPCpqqpi3LhxYQuQICEyIC6XzsxHCnjoiaJGg9PF57ZucIGmVFb6eO+jcha9WcqHn5U3+6UQ51B4+ek0TKbWl0fXdZ5/uYRbZ+Y3+qHanNYGSACTSWHenCSyOxu458HAmgg1ZeggE6++kEpGK4dgHzLYxLIPUrl4WgHfLw++StKZoLDopSRGhOhChK7rFLCPHLazlx1B1/TZcNCPQ4lTnHh0D6V6EZWUUUE5lZT7f6+qoDJCIbExKhoaBry4+dO1jj9d66JWlgOpaPjwUuYrpKyyMKplMShmDIoRt6+KzRWroloWqxqHRYvDo7tYW/Y13tLoBBAAo2LBYUxCQWF7+S9sKPk2amUBf5Nek2bF43OxqfB7NhV+H9XygL/puaoY8Pgq2Zm3ip2tb9jQanZ7Oh07HkVaWr/a+SRbM0iOw9GOdu2GkJ4+EKPRX8vjb35Zt9VRYmI3evc+m6TEZ1m6tJK11YFx3ToPW7Z4W3RxNVDZ2f7weOqpTYdH8LeWmv9KKT+vcWGzKq1urQPQt7eRu252cuJxlqDPcyadHYeiwKXXBBckr7gknpm3Joas76bPp/PAEwW89d8yfv6tKqDX7bQT7Tz/SFrQzXmD8dOaKk6dtIucPfu/X0cffZDUQoLURIZRr169yM3NbX7DVpAQ2Yz1f7i46Kp9rPql8YBgMhGyPn4HWvVLFc+/XMwb/1dGUXHg30xPP5hKdufQXKVSFIWrLk3gyEMtXHjlPtb90bKTvVAEyAPLcuM1CXTpbGDatXlBNyOdeKaNuQ8mY7GE5osgLVVjyeJkLpySz4cft7xJW1amxusvJ9Gje+tfs1K9iBy2s5vtQdc81sigI70YjKF6Goti8lhFAwOqxMCXgA9vVENsY2KpTB69KmxNLluqwudvshsL3HolBa7gBy0KNa/upsITvVDdEJ/uxRfGkbJVxYCOjh7APuLjO9Kp01EkJ/dEUVR8Pi/79v1GTs6P5OdvQNcD/640GCykpw+kXbshxMVl1lvv8VTi8+2/ANep0zF06XIcy5c/whFH5ge8n2B17apxzTUOTj+t+fBYQ1XhtnsKW1XjV6NzR43b/+Vk/Bm2Vgcnl0snNVmjV3cjazcEfnxnpGk892gKJ4Q4OKmqwuRzE/htvYtVa5r/XDz2CCuvzkuP6FRpHywrY+LU3ZSW1X0xD6YQ2dbmiWxLbrzxRq6//noOPfRQ+vfvH5Z9SIhshK7rPDu/hJvvzqeisumj/IyT7Y1OC9EanbIMJCdpeL2Bv8suPieOs08PfaAd2M/Mtx+056a783luQWAnf6oKL85NZcLY0JVH13Xat9M4pLeJVatbliIVBWbc6uSay0JXa6zrOt985+LxJ0v59IuWn6Af0sfAqwuSaZcR/PFTpVeym+3ksJ3SEAyio6LRi0G0o1Od58lCZAYxEKKl/LXQkWn22hwVDZvBSaW3NCZCu1mzY9SsaIoBj89NmTt61YgWUwL26mlCbOak6qlCUjAbHHz+68N4mmhmk5TUnY4dj8Lp9DerLyvbS07OyqAGyUlMzKZdu6GkpPSpbQLbEJfLP7q1wWCld+/xpKT08v8/LElUVIQvRHbpsj88tjS0OOwq3bMNzfZtbEpaqspN1yRw4TmOkLRoAvj48wrOunBP8xse4LQxNp56MIWU5PCMspqSrLHwyQx6Zudz54ONv57DB5l5+8V2WCyRGzjxuZeLuOLmffVanhmNcPTI0HV5iTppzho2S5cupWfPnkyYMIHu3buTmZnZ4OA6iqLw6KOPBrUPCZEN2L3Xw9Trcvnw08Bqci4+NzwD6iQnadxzcxLXTEngkXlFPPhEYZPb9+5h5JG7wzda2GdfV/LN8voDDDREVeGlJ1I5+4zQBtq7Zxfx0NyWT1sR51B44YkUTjohNB++Pp/OBx9W8tiTpaz6Obhag1FHmJj/XBJxcS3/YvLqHvayi91sI4+WfTE3xU48/TgUh1J/rlMLobn6qaBgxoZVsWNRHGiKxg5v88ONGzFjU+Io0kPXPENFw6zasKoOLKq9+ub/vdxbzLqKyDYlNCgmrGocVi0Oq+qo/ulv+pnv2sn68siWx6Lasarx/vkLtXisWjw2NR6TauXrgtciWtNaM/2GTUvAriVgq7kZEih272NV0fsRK4uCilWLr54KJAG7IRG7IaHOXJOf7f4PnghcZVdQsBoSsBmd2IyJ1T+dVHhKWJv7CVXesrAPtHMgg2apM4+kf27JZGzmJLRGRtzNK96Ex9vwd0u6sw9d0kdibd8Vj6eSnJyVQQ+S067dEDIyhmC1JjZ/B8DlKsXhaM8hh5yD1bp/jAG7PZ2Cgo0t2n8gOnf2h8czTm95eDzQoP6moEJkfJzCNZfFc9nkuEanugrWcaMsxMcpFAcw7YndpvDQ3clc+M/WzdsaiM++KeeuhxoPkH17mvjfovb15moOF13XuWNWPvc/XtDg+iOGW0P+2kST1ESGz/z582t/37BhAxs2NDwVVmveYxIi/2LZlxVccPlecvMDuzyS3cXAUSNbPy9dUzxenVWrm76ibbEoLJqXXmeY8VBb/HYpv69vPjCFK0AC/GOMtcUhsksnA6+9mELvnq0f3tjt1nn9rQrmPl3Kxk3BX+kdN9bK4w87W3yV95vvqli0uJwveQ9vC0/gTZjx4m20xqY9nenJQDSl4Y8FVVEx69Zmm8kqKFgUe3VI9AfFmt+tih0ztjojpm73NN1/0aE46aj1IkPrzE7vHxR5Ag2RCmbF2mBArPndpDTex+cX1+eNPrJRMaPrOh5aVhtuUIxY1XismgNLbViMw6o5sKpxdUYm/at1ZY331dMwBjUdiYqGVYvDplYHxJqwWF3Gxo6FXZV/hCVAqhj8ZTAkYK8OjDWh0aTaGn2t9lRtbnB5a1m1OP98kdr+eSNthoRm57IsdRdQFcJBeBRUrMYE7AeExJrAaDXEozYwH15xVWjnXaxTHkXFZkqqDopJ2CwptYHRZGj8dWrMnsK6nwGKotI+aQCd00ZiMydRWLadtWvfbPUgOYFMqXMgqzWZwYOn1qutdDgymryfwQDZ2QZ69TJgMSu89nrTn5mdOmlcc7WDsWNbFx5rDOxn4rW3Az/+zGaYelEc110RH5ZWVQAWi8opo228+lbTFzSGDTLz4txUunWNzMAxR46wckgvE2vW1v8s79rJwNLF7UlKjNzcl5WVOiccZeOQ3iZeWlzMx1/UPXYOpqasIAPrhNOCBQvCvo+gQuQPP/xAhw4dyMho/IP0119/ZdGiRdx///1BFy4ajj3SwhfvtWflz/6pKz78tIL1Gxv/0rron+EbUKdGeYXOymZC5Jx7kjmkV+jngDnQXTclsuT9MjxNZCdN8wfI8WFoUgswZKCZoYNMrPwpsJP3Iw8zs/DZlJB9MXq9cO/s4hYNnf5XV13m4I5b4oLqY/LG2xW8/lbgfR1VNNJoTwadSCKNb1laL0RqGOjFYNopHZt9PAs2XFRhxYYFGxbsWA3xB4REB2YsLTpZy/FubXB5ippJR60XSer+efosyv4mtUbMWBRbdVC1YVFsmBU7Nou/Nsis2Fo8P2QNj+5mr2tbveUOLZFO5r60M2WzqvQj8j05ddYbMNbWIlrUOP9zE2BIbEqlt5QCd919KShkmLPpZO3Htoo15FQ1XDNiUMz+YKbG7Q+K1aHRotqD+vza7dpUb1myMZPO1v78UvJpk1OkqGjVYdUfyA6sWTQHUR6f7mVP1ZY6ywyY8RBYU1KzaqsNiHaDE5vm/2k1xDcaopuzp7LltVQKqj88m5Lq1SpaDE2H1obYjM4Wl+GvzAbH/hrFA2oWrebEoN9bf6XrPvYW+UOkqhrpkDyYTmmHoSgKu/JWszP/Z8qrWtZ0tKFBcoJhsTgbXG637x/tMDNTpVcvI716GejV0x8cs7MNtRcIH3+8/nRgNTp11Lj6ajtnnmkNaV+7I0damHy+g4H9TTz6VDEbNzf8pa2qMGminVuuTSCzffjrFMaf5mDTFg8D+5l4dn7dbjGqCjdf7eSWa5wYjZHrd2gwKDw6M4UJl+5mwCFmln3p/35tl67x4WuZtM+IbF2L1apy1Egruq7z4JP+2khVpXbgnxOOOrhCZCyMqXCwGj58eNj3EdS7Y9KkScTHx/Pggw9y1FFHNbjNzp07eeedd9pciFQUhewuRrK7GBlznI0PPmn8hN1ggEkTwj83ZJeORp6bk8rZlzTcbHHcqfaIzFGZ3cXIxefG8ez8Euw2hfIKvU7n/XAHyBpTLoxj5U/N9+uZfL6DB+5ODOkXksWiMPUSO/fMarxf6KEjTA2O1KoocO9d8Uy5JPjn5+yzrLz8avNXmJNII4NOpNG+dl5HXdfr1SI6SKAfh2JXAjt+BnEEGoY6J/qKIfgRZct8xRTr+19LFY1MLZsOWk/sav0mtUlqBiNNp2JRbI3XmBpa3zJgr2vrAWFbIc3YiU7mPiQaMmr/71nmnqQaO/qbnar+4FgvJGqhuXix+4CaNoNipoOlNx0tfWvnMYw3pFLpLavX7NSmxQcdXBvj9lWR69oB7A+yna39iTekAGDTEij25NaWxabWbX5qUe0oWuhOzPJcO2v7HlrUODrbBpBp7cW3ea9T4fO3WjAq5upwmFCvSaxBNaEYQlvTsLuifsgG/wAyNi2huhz+prg2g788Fs3uv/hiDU3LFoNqwqLFUeltug+7phiwGZOwmxKxG5Oqb4nYTYmozsCafbZGYdmf+HQvXTNGkZUyhOKynfz+5//IK97YojlJDxwkx+FoH9aLuw5HBkveTqJnTwPx8U2H6U+W1b+Y0bGjxtXT/eExHIFpwCEm5tyfRHGJj2tubjiAn3GKldtvdNKzW+Smijj5BBsnn2DjnffL6oTILp0MvDg3lUOHhrdVV2OOHmljyw+dee3/Sln2ZQWJTpWli9vTtVP0ptFQFIUv3snirgfzaZeu8cGycn5b72JA3/BWFkSchMiwc7lcrFixgs2bN1NRUYHdbic7O5thw4ZhMLTuuzjoe5eXl3PZZZdx9dVXM3Xq1FYVIlZ9+W0Fm7Y2Xgt56ok20lMj08zhjJPsXDU5npJSnfYZGvc9WghAl44Gnn4wNWJzFt16rROjQeHG6U4OOfzP2uHDNQ3mP5nKuNPCGyABxv7DxoLFpYw7zcau3V4eeKxu81ZNg4fuSWTyBeEJ1hedb+fRJ0spOaBvh6rCGadauepyB6kpKocMqRv4zWZ4+rFETvtH6/pkHjrcRFamxo6d9ZsSOkggg45k0BGLUn8/LirrnJRlkU13+qM10BSuMTWBNFR2+7YCYMZGR0MPMrVuGJXGQ49BMYa8DA3Z5dqIQTGRZepJR3NvrFr9Y6mdKTvs5agtT9Uf2LQEOln60d7Svd5z0Nnaj87WfhEpy17XVlRFI8vcl07Wfli1uu/5wfFjMCqmFjcdDNaeqk04tCS62AeRYc6ubdrZN/6o2n6UJjVyJ6dlnkI8uos+CUfj8VVhUq3VYTG4mtbWsJsSqazwn6wbFBNWo5NESyYOU3J1UEzCrDXe7ywSLc1MBjtDss9ld+HvfL/uWVye8AySE0qqamDYsOZP5nNzvfz00/5ziA4dNKZPtzPurPCEx7/66rvKei2Hjh1l4c6bExg8IPzzWjdm05b9z8n5Exw8fE9yxPocNsZmU0lMULHbFP63qD2H9Ire81PDYVd56C7/xbmTj7Pz8LyCiE4vEgnSJzK8lixZwqxZsygu9p8n67pe+3mfmprK3XffzdFHHx304wcdIu+44w4++OADHn30UX7//XdmzZqF1XoQjRgFnDrGzv8WZzBpmr+P5InHWvn0q4ra+RLDNaBOYx68KxlVVVj6qb8mymCAl59OI6GZK6GhlJFm4OF7/IP3mE0KJehoGix4Ko2zTo3M6J0Wi8IHb/qbE/37vsI66xKdKi8/m8Kow8N30hgfr3LheXbmPl2KxQznTLRx+RQHnTv5307b/6z7rZ2QoLDwhSRGHtr6LyVVVRh/ppU5c/1NpMxYaoNjnOJs8r6V+I8bDQN9GEq6ktXq8rSGrutU+ErpZzycNLVjyJrHtZZX95Bh6sog0wlBN2cMJY/uprt9GCnGDhENII2xavEclXhOozWckQxsAFnWPvSNO7rec5Nsis7xbVQtHJl2Xkwczz2Tj/b3XzQ40dToH8sH8nir2F3wGzvzf6aobEeL7rt/kJzBdQa8iTWfflqFrkNWlsr06Q7Gj4tMeKzd/xf7BysaMtDEXTc7OfrI6NT2HWjjFjeJTpUnH0jhzH/EzqjfGWkGlrzUjhGDo/8c/VXvHiYeujMl2sUIOQmR4fP1119z6623YrfbOfXUU+nWrRsWi4Xy8nLWr1/P559/zpVXXskrr7wS9BQgQX+rJCYm8p///IcHHniAl156ia1bt/Lkk0+SlRXdE9NQO+YIK98uzeTO2fnMnZXCDf/O46VXS+mYZeC4UZENzTVXoPr38V8BvffWJIYNit6HndmsRDxA/lVR0f5r5b16GHntxRS6hmiOzKZMnWzHaIRLL7aTmlK3Jq/qgJas7dupvP5yMr16hq5MZ59lZfHcZDLoSBJpAQeLSsqJJ5FDGIFNCX+NcSAOMY2MdhHq0RQDWeae0S5GLYNiJNXUfH/VSEkytot2EepwGtOb3yiCIh2imxJvTot2EerQdZ3Csj/ZmfcTuwt/j+ggOdGwbp2H2bPjGT/OGrKpMlpi2ReVdM82cOfNTk47yRoTF6EAOnUwsnJZJpntYuvCxqFDYue92xBnQuQG+IkYGVgnbF588UU6duzIyy+/TGpqar31O3fu5LzzzuO5555j7ty5Qe2jVe9gVVW5+eab6du3L3fccQfjxo1jzpw5HHbYYa152JjTMcvAi3P9X8bXXZbA/MWlXPhPR9SaFbRL15h0toOrpyREZf817DaVBU8lRy1AAhRWh8gxx1t44YkU4oOYLiMYGekat95Yv88egKvKf2mtd08Dry1Mpn370H7wd+9mpK8yrMX3c+BkKMfERA0JtG5YaSFE21HlLmFX/i/szPuZ8qqWzVMZZ00nM3kQKZ2HtmqQnGi4/fbgBlALhaJiH9ddGc854+whHbQnFG6+2hntIogYITWR4fPrr78yefLkBgMkQGZmJmeffTYLFy4Meh8huQxUU016xRVXcOmll3L99dfTvn37UDx0zOnRzcS4U+2cH4EBdRqjKArPPJQa9bbx/5mbytCB0e03UFTs45rL47jrZieaFhtflFUuncMPM7Hg+SQSEmIjsAEBD54jhBChsiv/F37b9m7LBsnRzLRL7Edm8iDirP7BrNzG2K4lakg0v6MT4lXOnxgbLU6EaJSEyLApKysjLa3p1ijt27ev7S8ZjJC1Jejduzdvv/0211xzDQ888ACdOnUK1UPHnLmzk6PerCAWrixGO0AC3HFjAkNioBwH6tzJwOsvJ2M2R/81EkKIaHLaOwQcIJMcXchMHkiasxeaGr2RMYUQkSHzRIaP0+lk69atTW6zfft2nE5n0PsIaTWJ0+nkP//5DxdccEGzBQ+FdevWMWXKFIYMGcKAAQM477zzWLFiRdj3G+0AKfaLtQAJkJSoSoAUQgjAZk4kydG50fUWYzxdM0ZxRJ+rGNp9Eu2S+kmAFOJvQtFj63YwGTZsGK+++ipr165tcP0vv/zCyy+/zIgRI4LeR1A1kcuWLSM5ObnBdTX9JIcOHcq6deuCLlhztm/fzrnnnkvXrl156KGHsFgszJ8/n4svvphFixYxYMCAsO1bCCGEEIHJTB5EfunW2r8VRSMtoSeZyYNIjuvSJgbJEUKEgdREhs1ll13GZ599xllnncWgQYPo1q0bNpuNsrIyNmzYwC+//ILVauXyyy8Peh9BhcjMzMxmtzn++OM5/vjjg3n4gDz11FN4vV6eeeYZkpL8Q3wPGTKE0aNHM2fOHF566aWw7VsIIYQQgUlz9sKww4zF6CQzeSDtkvphMrStQXKEEKF3sNX+xZIePXrwwgsvcNddd/Hjjz/y448/1ll/yCGHMGPGDLKzg5/zOrbGVw6Qrut88sknjBw5sjZAAphMJkaPHs3ChQspLi4mPr7h0TOFEEIIERmaamRkr8swG+NkVGYhRC3FJykynIYMGcJ7773H1q1b2bhxI+Xl5djtdnr06EGHDh1a/fhtMkTu2rWLkpISunfvXm9d9+7d8fl8bNiwgaFDh0ahdEIIIYQ4kMUkF3WFEH8hGTLs3G43WVlZdO7cuXZZYWEhLpcLk8nUqsdukyEyL88/z1RiYmK9dTXLarZpjLnd5tAXTIRcLA2bE0uDpX8s/QiEECKKbol2AYRo86Q5a3i9+eabPPjgg8yfP59evXrVLv/444+ZPXs2t956K2eeeWbQj98me7O7XC6ABhO00egf1a2ysjKiZRJCCCGEEEIERvHF1u1g8sUXX3D77bejqiqaVndWiY4dO+JwOLjtttv49ttvg95HmwyRZrO/fsrtdtdbVxMwrVZrRMskhBBCCCGECJAeY7eDyHPPPUefPn34+OOP63X/GzFiBB988AF9+vRh7ty5Qe+jTYbI1NRUAPLz8+uty83NrbONEEIIIYQQIrZEu+bxYK6J/O233zj77LNxOBrujGW1Whk3blyrpmNsk30iMzIySExMZP369fXWrV+/HqPRSI8ePaJQMiGEEEIIIURzpE9k+KiqisViaXKb1g6s0yZrIgFOPPFEvv32W/bt21e7rLy8nI8++ohRo0Zht9ujWDohhBBCCCFEo3Q9tm4HkZ49e/Lhhx82ut7r9fL222/TrVu3oPfRJmsiAS6//HKWLl3KtGnTuOqqqzAajTz33HNUVFRw3XXXRbt4QgghhBBCiEYcbE1IY8nZZ5/NzTffzCWXXMLYsWPp2LEjZrOZoqIi1q9fz5tvvsmGDRu45557gt5Hmw2R6enpvPLKKzz44INcd9116LrOwIEDWbBgQatStRBCCCGEECK8JESGzxlnnMGGDRt48cUXGx2B9fzzz2fcuHFB70PR9YOs/la02oUXXsh3333HggULGDFiRO1yr9fLs88+y5IlS9i1axfJycmcfPLJXH311c22uxYHl5KSEp566in++9//UlBQQPv27Rk1ahRXXHFFnflb5ZgRB1q3bh2PPPIIP/74Ix6Ph379+jF9+nSGDx8e7aKJGPDtt98yd+5cfv/9d8xmM926dWPq1KkcddRRtdvIMSQa8sMPPzBp0iSGDRvGwoULa5fL8RLbDh//cLSLUMc3b1wfkscJ9txn+fLlnH/++Y2u/+233zAYWlb/t2nTJv773/+yYcMGysrKcDgcZGdnc9JJJ9WZOzIYbbYmUoTH22+/zXfffdfgupkzZ/L6669z5ZVXMnz4cP744w8eeOABtm7dytNPPx3hkopocbvdTJ48mU2bNjF9+nR69+7Nr7/+yuOPP86KFStYsmRJ7ZxEcsyIGtu3b+fcc8+la9euPPTQQ1gsFubPn8/FF1/MokWLGDBgQLSLKKLo008/5bLLLmPUqFHMnTsXXdd56aWXmDJlCo8++ignnXSSHEOiQS6XizvuuIO/1onI8dIGHKT1WK0995kxYwZ9+/att7ylARIgOzubq6++usX3C4SESFErLy+P2bNnc/TRR/P555/XWffnn3/y6quvcvHFF3PZZZcBMGTIEHw+HzNmzGDlypUMHTo0CqUWkfbpp5/y888/M2vWLMaOHQvAsGHDcLvdPPzww6xYsYLDDjtMjhlRx1NPPYXX6+WZZ54hKSkJ8B8Po0ePZs6cObz00kvRLaCIqjlz5tC5c2eeeuopjEYjAMOHD+foo49m4cKFnHTSSXIMiQY99dRTFBcXc8ghh9RbLsdLbDsYR2cNxblPly5d6NevX0jLVVlZyfPPP8+ZZ55J+/btQ/KYbXZ0VhF69913H06nk/POO6/eumXLlqHrOqeeemqd5aeccgqKovDJJ59Eqpgiyrp168a9997LCSecUGd5TbOInJwcQI4ZsZ+u63zyySeMHDmy9mQO/MOLjx49muXLl1NcXBzFEopo0nWdyy67jBkzZtQGSPDPY9apUyd2794tx5Bo0IYNG3j++ee5/vrrsdlstcvleGkboj0vZDjmiYzVc5+KigqefPJJdu7cGbLHlBApAPjiiy/473//y4wZMxqcN2b9+vUoikJ2dnad5QkJCaSlpbF27dpIFVVEWXZ2NuPGjas3ge3mzZsB6NixIyDHjNhv165dlJSU0L1793rrunfvjs/nY8OGDVEomYgFiqJw8sknc+ihh9ZZ7na72bZtGx07dpRjSNTj8/n497//zaBBgzjrrLPqrJPjpY2I9pQeYZjiI5bPfUI9DI6ESEF5eTkzZsxg7Nix9b7Ea+Tl5eFwOBoMmImJieTn54e7mCKG7dmzh3nz5tGrVy+GDBkCyDEj9svLywOoM+hSjZplNdsIUWPu3LkUFhZyzjnnyDEk6nn11Vf57bffmDFjRr11cry0DdGueQxHTWQozn0++OADzjrrLAYNGsTw4cO5+uqr2bZtW2gKGELSJ/Ig5Xa72b59e5PbxMXFkZaWxmOPPUZZWRk33nhjo9tWVVU1+IYAMBqNFBUVtaq8IvpacswcKC8vj2nTpuFyuXjwwQdRFAWQY0bs53K5ABo8HmqaL1ZWVka0TCK2LV68mGeffZYzzzyT0aNHs3LlSkCOIeG3e/duHn74YaZMmULXrl3rrZfPnLahLfWJDPQcKRTnPuvWrWPKlCmkpqbyyy+/8OSTT/LTTz/xv//9j7i4uKDKHw4SIg9Se/bs4eSTT25ym7Fjx3LuueeycOFCZs6cWaffwF9ZLBbcbneD61wul0zXcBAI9JiZNWtW7d/bt29n8uTJFBQU8Pzzz9OjR4/adXLMiBpmsxmgweOh5mTParVGtEwidj3xxBPMnTuXU089tXYibDmGxIFmzJhBWloaU6dObXC9HC9thK/tpMhAz5Fac+4zaNAgvv76a5KTk1FVf2PRwYMH06VLF6ZMmcILL7zANddcE1B5t2/fXtu9CPzvibFjx5KSkhLQ/QMhIfIglZWVxfr165vcxuPxMG7cOAYNGsSJJ55IWVkZsP/qXGVlJWVlZdjtdlJSUigpKcHlctW7wpKXl9fquWZE9AVyzBxo3bp1XHzxxVgsFl599VW6detWZ70cM6JGamoqQIPNeHJzc+tsI/7e7rzzThYvXszkyZO54YYbals2yDEkanz44Yd89tlnPPPMM7jd7toTdq/XC0BZWRnJycmAHC+xLlRNSCMh0HOk2267LehzH5PJ1OBxeeSRR2KxWFizZk3A5R09ejSHHnooZ599NieccAI2m437778/4PsHQkLk39ju3btrO/gOHjy43vopU6YA/k7CvXv3Rtd1/vjjjzpz1+Tm5pKbm9vgfDbi4LVz504mT55McnIyzz//POnp6fW2kWNG1MjIyCAxMbHBL+D169djNBrr1GKLv6c5c+bw2muvcdttt9WbcFuOIVHjs88+Q9f12nOUvxo8eDBXXnmlHC9twUE4T2Rrz30aCp8ejwe3211bwx6IgQMH8v3337N8+XISEhI444wzGD9+fL0Bf1pDQuTfWFpaGosWLaq3fO3atcycOZPbb7+d3r17A3DCCSdw//338+6779Z5A7z77rsAnHTSSZEptIg6r9fL9OnTMZlMvPjii402jZBjRhzoxBNPZMmSJezbt6/2Smt5eTkfffQRo0aNwm63R7mEIpo++eQT5s2bxw033FAvQNaQY0gATJs2jXHjxtVbPnPmTABuv/122rdvT25urhwvMa4t1UQGqjXnPrfeeivvv/8+n376aZ0uZsuWLcPr9TJ8+PCAy7F48WJ27tzJe++9x//+9z9eeukl5s+fz8CBA5kwYQJjxoxpdbciCZF/YyaTqcEJT2uahPTo0aN2fXp6OhdffDEvvPACiYmJDB8+nN9++43HH3+cM888szZsioPfkiVL+PXXX7nuuuvIycmpnReyRmJiIllZWXLMiDouv/xyli5dyrRp07jqqqswGo0899xzVFRUcN1110W7eCKKPB4Ps2bNIisrixEjRjTYZKtnz55yDAkAOnfuTOfOnestrxlwpOa8RY6X2KcchDWRLTn3ueCCC9iyZQtffvklAOeeey7vv/8+F110EVOnTiU1NZXVq1fz9NNPk52dzfjx41tUlszMTKZNm8a0adNYv3497733Hu+//z4333wz9957L6eeeirjx48P+nxM0UM9aYho85YvX87555/PggULGDFiRO1yXdd54YUXeP3119m1axdpaWmcccYZXHbZZXUmiBYHt5tvvpklS5Y0uv7AwXfkmBEH2rRpEw8++CArVqxA13UGDhzIddddR79+/aJdNBFFO3bs4Ljjjmtym2XLlpGVlSXHkGjUpEmTAFi4cGHtMjleYttxx4a2j15rLfv0lpA8TqDnPpMmTWLz5s188803tct+//13Hn/8cX788UfKy8tJS0vjuOOO48orr8TpdIakfMuXL+fRRx/l559/BmDIkCFceumlHHXUUS16HAmRQgghhBBCiIg67pgYC5GfhSZExqo9e/bUNm+tGRNl4MCBbNu2jcLCQk499VTuvffegC/yS3NWIYQQQgghRGRJPVbY+Xw+Pv30U958802+/vprPB4P8fHxTJo0iYkTJ5KdnU1FRQWPP/44L774IklJSdx8880BPbaESCGEEEIIIUREKW1onsi2Ztu2bbz55pssWbKEvLw8dF2nf//+TJw4kVNOOaXOSK9Wq5WbbrqJvLw83n33XQmRQgghhBBCiNh0MI7OGitOPPFEFEXBarUyfvx4/vnPfzY7gM6RRx7Je++9F/A+JEQKIYQQQgghIktqIsOmZ8+eTJw4kdNOOy3g6WwGDRrEQw89FPA+JEQKIYQQQgghIupgnOIjVtTMS9kSWVlZZGVlBby9hEghhBBCCCFEZEmIbNMkRAohhBBCCCEiSvFKiGzLJEQKIYQQQgghIktqIts0CZFCCCGEEEKIyJKayDZNjXYBhBBCHBx+//13Bg0axL/+9a+o7P/LL7+kT58+zJkzJyr7F0IIEThF12PqJlpGQqQQQohWKyoq4vLLLycrK4uZM2dGpQyjRo3iyiuvZN68eSxdujQqZRBCCBEgXY+tm2gRCZFCCCFa7d577yUnJ4cZM2ZgNpujVo6pU6fSs2dP7rzzTvLz86NWDiGEEM3w+WLrJlpEQqQQQohW+e2333j33Xc54ogjGDx4cFTLomka119/PYWFhcybNy+qZRFCCNEEX4zdRItIiBRCiL+plStX0rt3b/75z3+i/6Upz5tvvknPnj256aabmn2chQsXAjB58uQ6y99++2169uzJSy+9xA8//MCkSZMYMmQI/fv3Z+LEiaxYsaLO9pMmTaJnz57k5+cze/ZsjjzySPr378+pp57KsmXLAHjvvfc4/fTTGTBgAMceeywzZ87E7XbXeZxRo0aRmZnJG2+8QXl5eYufFyGEEOGn+HwxdRMtIyFSCCH+poYOHcoll1zCqlWreOWVV2qX5+bm8sADD5CZmckdd9zR5GPous6XX36J3W5n2LBhDW6zevVqpk+fzvDhw7njjjuYOHEia9asYerUqeTm5tbbfubMmWzdupWrr76aSy65pPb3//znP8ydO5czzzyTW265BbPZzMKFC5k/f36d+yuKwhFHHEF5eTmff/55y58YIYQQ4efTY+smWkSm+BBCiL+x6dOn8/XXX/PII49w3HHHkZGRwcyZMykpKeGpp57C4XA0ef9169aRl5fHsccei8HQ8FfK0qVLWbJkCb169QLgjDPOoKioiHfeeYcvv/ySM888s872xcXFPP/887V/5+Xl8dprrzF37lw+/vhjUlJSAOjZsycTJ07ks88+q1cLesQRR/Daa6+xcuVKTj755BY/L0IIIcJMBrNp06QmUggh/sZMJhMPPvggbrebu+66i08//ZQPPviAKVOmMHTo0Gbvv2XLFgC6dOnS6DaHHXZYbYCsMWDAAAD27NlTb/vx48fX+bt3794AHHvssbUBEqBPnz4A7Nu3r95jdO3aFYAdO3Y0+38QQggRBdEeSEcG1mkVqYkUQoi/ue7du3P99ddz3333sXz5cvr27cuVV14Z0H0LCgoASE5ObnSbzp0711tmsVgA8Hg89dZlZmbW+btmtNfGljf0GElJSQAyQqsQQsQqaULapklNpBBCCM444wzMZjPl5eWMHTsWo9EY0P2Ki4sBiIuLa3Sblk75YTKZGlweaJkOLE9JSUmL9i2EECJCfN7YuokWkRAphBCCu+++G4/HQ3Z2No8//jg5OTkB3S8+Ph6IvbBWU56mwq0QQogoivZAOjKwTqtIiBRCiL+5999/n//+979ccsklzJ07l8rKSm688UZ8AfQRSUxMBPyD38SSmvLUNGsVQggRY3Q9tm6iRSRECiHE39iePXuYMWMG2dnZXHXVVWRnZ3PllVeyYsUKXnjhhWbvXzOgTs0AO7GipjxZWVlRLokQQogGRXsgHRlYp1UkRAohxN+UruvceuutFBcXc99999X2Rbzkkkvo27cvjz32GGvXrm3yMXr16kVycjLLly9vcICbaPn6668BAhphVgghRBREOzRKiGwVCZFCCPE39corr/D1119z4YUXMnDgwNrlBoOB++67D4AbbriBqqqqRh9DURRGjRpFWVkZP/zwQ7iLHBBd1/n666+x2WwcffTR0S6OEEKIhkQ7NEqIbBVF16URsBBCiOD99ttvnHnmmRxxxBEBNYENty+++IIpU6Zw4YUXcsstt0S7OEIIIRpwUuq0aBehjg/2zYt2EdoUmSdSCCFEq/Tt25fTTz+dd999l1WrVjF48OColcXr9fLII4/gdDqZNi22TlCEEELsp3tlWo22TJqzCiGEaLVbb72Vdu3aceeddzbZ/DXcnnnmGdatW8eMGTNqR44VQggRg6I9GquMztoqEiKFEEK0mtPp5KmnnmLHjh3cfvvtUSnDV199xRNPPMG0adMYM2ZMVMoghBAiQNHuAyl9IltF+kQKIYQQQgghIupExwXRLkIdH5bOj3YR2hTpEymEEEIIIYSILKnHatMkRAohhBBCCCEiSwbWadMkRAohhBBCCCEiSvdJTWRbJiFSCCGEEEIIEVm6DGbTlkmIFEIIIYQQQkSUzBPZtsnorEIIIYQQQgghAibzRAohhBBCCCGECJiESCGEEEIIIYQQAZMQKYQQQgghhBAiYBIihRBCCCGEEEIETEKkEEIIIYQQQoiASYgUQgghhBBCCBEwCZFCCCGEEEIIIQImIVIIIYQQQgghRMAkRAohhBBCCCGECNj/A9C4Q8z5vl3xAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "system.m.sel(\"y\").mpl(figsize=(10, 3))" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.10.12" }, "widgets": { "application/vnd.jupyter.widget-state+json": { "state": {}, "version_major": 2, "version_minor": 0 } } }, "nbformat": 4, "nbformat_minor": 4 }