diff --git a/notebooks/DCinverse.ipynb b/notebooks/DCinverse.ipynb index 2cf8b20b..f0be58c0 100644 --- a/notebooks/DCinverse.ipynb +++ b/notebooks/DCinverse.ipynb @@ -1,6 +1,7 @@ { "metadata": { - "name": "" + "name": "", + "signature": "sha256:789aaa72085626eec282d4a2b06df533d7a48163114ea5f5bee10382f0f44bf2" }, "nbformat": 3, "nbformat_minor": 0, @@ -13,29 +14,22 @@ "input": [ "from SimPEG import *\n", "import simpegDC as DC\n", - "%pylab inline" + "from pymatsolver import MumpsSolver" ], "language": "python", "metadata": {}, - "outputs": [ - { - "output_type": "stream", - "stream": "stdout", - "text": [ - "Populating the interactive namespace from numpy and matplotlib\n" - ] - } - ], + "outputs": [], "prompt_number": 1 }, { "cell_type": "code", "collapsed": false, "input": [ - "cs = 25.\n", - "hx = [(cs,7, -1.3),(cs,21),(cs,7, 1.3)]\n", - "hy = [(cs,7, -1.3),(cs,21),(cs,7, 1.3)]\n", - "hz = [(cs,7, -1.3),(cs,20)]" + "cs = 12.5\n", + "nc = 500/cs+1\n", + "hx = [(cs,7, -1.3),(cs,nc),(cs,7, 1.3)]\n", + "hy = [(cs,7, -1.3),(cs,int(nc/2+1)),(cs,7, 1.3)]\n", + "hz = [(cs,7, -1.3),(cs,int(nc/2+1))]" ], "language": "python", "metadata": {}, @@ -46,11 +40,29 @@ "cell_type": "code", "collapsed": false, "input": [ - "mesh = Mesh.TensorMesh([hx, hy, hz], 'CCN')" + "mesh = Mesh.TensorMesh([hx, hy, hz], 'CCN')\n", + "print mesh" ], "language": "python", "metadata": {}, - "outputs": [], + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + " ---- 3-D TensorMesh ---- \n", + " x0: -541.97\n", + " y0: -416.97\n", + " z0: -548.22\n", + " nCx: 55\n", + " nCy: 35\n", + " nCz: 28\n", + " hx: 78.44, 60.34, 46.41, 35.70, 27.46, 21.13, 16.25, 41*12.50, 16.25, 21.13, 27.46, 35.70, 46.41, 60.34, 78.44\n", + " hy: 78.44, 60.34, 46.41, 35.70, 27.46, 21.13, 16.25, 21*12.50, 16.25, 21.13, 27.46, 35.70, 46.41, 60.34, 78.44\n", + " hz: 78.44, 60.34, 46.41, 35.70, 27.46, 21.13, 16.25, 21*12.50\n" + ] + } + ], "prompt_number": 3 }, { @@ -73,7 +85,7 @@ "cell_type": "code", "collapsed": false, "input": [ - "nElecs = 31 \n", + "nElecs = 21\n", "x_temp = np.linspace(-250, 250, nElecs)\n", "aSpacing = x_temp[1]-x_temp[0]\n", "y_temp = 0.\n", @@ -82,13 +94,136 @@ "language": "python", "metadata": {}, "outputs": [], + "prompt_number": 5 + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "txList = DC.Examples.WennerArray.getTxList(nElecs,aSpacing)\n", + "survey = DC.SurveyDC(txList)\n", + "print len(survey.txList)" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "63\n" + ] + } + ], + "prompt_number": 6 + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "fig, ax = plt.subplots(1,1, figsize = (6.5,5))\n", + "indz = 22\n", + "dat = mesh.plotSlice(sigma, grid=True, ax = ax, ind=indz)\n", + "ax.set_xlim(-300, 300)\n", + "ax.set_ylim(-300, 300)\n", + "cb = plt.colorbar(dat[0])\n", + "ax.set_title('Depth at '+str(mesh.vectorCCz[indz])+' m')\n", + "\n", + "txLocs = np.array([[tx.loc[0][0], tx.loc[1][0]] for tx in survey.txList]).flatten()\n", + "ax.plot(txLocs, txLocs*0, 'w.', ms = 3)" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "metadata": {}, + "output_type": "pyout", + "prompt_number": 8, + "text": [ + "[]" + ] + }, + { + "metadata": {}, + "output_type": "display_data", + "png": "iVBORw0KGgoAAAANSUhEUgAAAZwAAAFRCAYAAABANxtTAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xt4FFWeN/Bvh7QLhlsGpRnSkWi6IQmEJiMkOF4WhBgT\nl34ZUCagTICwZsNgZJyRm8vaOpLLODu+QgY3+HJ3N0THwbDatKADiitJFALu2DzS0USSzkW5BMFB\nEjrn/QPpSYfuTnXodFXK7+d56un07/xO1TkJTx2q6lSVRgghQERE1MvC5G4AERH9MHDAISKikOCA\nQ0REIcEBh4iIQoIDDhERhQQHHCIiCgkOONSnxMTE4N1335W7GUTUAxxwqMdiYmJw4403YvDgwYiM\njMSdd96JkpISBOvWrgULFmDNmjUeMY1GA41GE5T1x8TE4C9/+UvA9Xbu3In4+HgMHDgQBoMBH3zw\ngbvsjTfewNixYzF48GCMHTsW5eXlPtczcOBADBo0yL2Eh4cjLy8PAFBXV4ewsDCP8rVr1wbeSSIF\nCZe7AdR3aTQavPnmm7j33ntx/vx5HDhwAI8//jgqKyuxefNmuZvXLY1GE/DguG/fPqxcuRKvvvoq\nkpOT0dTU5F7HV199hYcffhh//vOfkZaWBqvVioceeghffvklbrrppmvWdeHCBffP3377LUaMGIE5\nc+Z45HzzzTdBG2CJZCeIeigmJka8++67HrGqqioRFhYm/vrXvwohhPjuu+/Er3/9a3HLLbcInU4n\n/uVf/kVcvHhRCCHE/v37RVRUlMjPzxc33XSTiImJEf/5n/8phBCipKREaLVaccMNN4iBAwcKs9ns\n3ubvf/97MX78eDFkyBDx85//XHz33Xde21dTUyOmTp0qhg0bJm666Sbx8MMPi9bWViGEEI888ogI\nCwsTAwYMEAMHDhTPP/+8pD7fcccdYvPmzV7L/ud//kcMHz7cI3bzzTeLioqKbte7detWERsb6/5e\nW1srNBqNuHz5sqR2/eM//qP413/9V/HTn/5UDBw4UMyYMUN8/fXXYt68eWLw4MFi0qRJoq6uTtK6\niHoLT6lRUE2aNAl6vd59mmnlypWoqanBsWPHUFNTA6fTiWeffdad39LSgtOnT6OxsRHbtm3Do48+\nCofDgUcffRQPP/wwVqxYgfPnz7tPTQkh8Nprr+Htt99GbW0tPvnkE2zdutVne5566ik0NTXh+PHj\nqK+vh8ViAQDs2LEDt9xyC958802cP38ev/nNb7rtm8vlwuHDh/HVV1/BaDQiOjoajz32GL777jsA\ngMlkQnh4ON588024XC688cYb6N+/P8aPH9/turdt24Zf/OIX18RHjRqF6OhoLFq0CKdPn/a7jrKy\nMrzyyitwOp34/PPPcccddyA7OxtnzpxBfHw8nnnmmW7bQdSbOOBQ0I0cORJnzpyBEAIvv/wy/vCH\nP2Do0KEYOHAgVq1ahZ07d3rk//a3v4VWq8U999yDBx54AGVlZQCuDC6iyykvjUaDvLw8jBgxApGR\nkZgxYwaOHj3qtR2xsbGYNm0atFotbrrpJvzqV7/Ce++91+N+tbS0oL29Ha+//jo++OADHD16FNXV\n1XjuuecAABERESgpKcHPf/5z9O/fHw8//DBKSkowYMAAv+v98ssv8f777yMrK8sdu/nmm/Hxxx/j\n5MmTOHz4MM6fP4+HH37Y5zo0Gg0WLlyIW2+9FYMHD0Z6ejpGjx6Ne++9F/369cNDDz2E6urqHved\nKBg44FDQNTQ04Ec/+hFOnTqFv/3tb7j99tsRGRmJyMhIpKen49SpU+7cyMhIjx3yqFGj0NTUBAA+\nr12MGDHC/fOAAQM8roV01tLSgszMTOj1egwZMgTz58/v9iihs/T0dPcF+9LSUnc7H3vsMeh0Ogwb\nNgxPPPEErFYrAODIkSN49NFHcfDgQbS3t+O9995DdnY2jh075nc7O3bswN13341Ro0a5YxEREfjJ\nT36CsLAwDB8+HMXFxdi7dy++/fZbn+vR6XTun/v374/hw4d7fPf1eyIKFQ44FFQfffQRGhsbcddd\nd2HYsGEYMGAA7HY7zp49i7Nnz6K1tRXffPONO//s2bP429/+5v7+5ZdfYuTIkQB8Dzid+ctZvXo1\n+vXrh7/+9a84d+4cduzYgY6ODkl1AWDPnj04f/48zp8/j7lz5yIyMhJ6vd5n/rvvvovJkyfjJz/5\nCQBg4sSJSElJwTvvvON3O9u3b/c4uvGnc/v94UQDUiIOOHRdrp7y+uabb/Dmm29i7ty5mD9/PsaO\nHYuwsDD88z//M5YtW4avv/4aAOB0OrF3716PdTz99NNob2/HwYMH8dZbb+Ghhx4CcOV/7F988YWk\n7Xtz4cIFREREYPDgwXA6nXj++ec9ynU6HT7//POA+rtw4UKsX78eX3/9Nc6ePYsXXngBM2bMAHDl\nGs7BgwfdRzTV1dU4ePAgTCaTz/V9+OGHaGxsdPf5qqqqKnz22Wfo6OjA6dOnkZeXh6lTp2LQoEE+\n19X5d+Hv90IkFw44dF1mzJiBwYMH45ZbbkFBQQF+/etfY8uWLe7yoqIiGAwGTJ48GUOGDEFqaipO\nnDjhLr96LWbkyJGYP38+SkpKMHr0aABAdnY27HY7IiMjMWvWLK/b93dfztNPP40jR45gyJAhmDFj\nBmbPnu2Ru2rVKjz33HOIjIzEH/7wB0n9XbNmDSZNmoTRo0cjISEBt99+O5566ikAwH333Yfly5dj\n1qxZGDRoEB588EE89dRTmD59OgAgPz8fGRkZHuvbvn07Zs+ejYiICI/4F198gfT0dAwePBiJiYkY\nMGAASktL/batc9+8/V541EOyk2t63MWLF0VycrIwmUwiPj5erFy5UgghxOnTp8X06dOF0WgUqamp\n4uzZs+46+fn5wmAwiDFjxoi3335brqZTkOzfv1/o9Xq5m0H0g7Bnzx4xZswYYTAYRGFhodecxx57\nTBgMBjF+/Hhx5MgRd3zhwoVi+PDhYty4cR75/vbX3sh2hNO/f3/s378fR48exSeffIL9+/fjgw8+\nQGFhoft/wdOmTUNhYSEAwG63o6ysDHa7HTabDUuWLJF8PpuI6IfM5XJh6dKlsNlssNvtKC0txfHj\nxz1yrFYrampq4HA4sHHjRuTm5rrLFi5cCJvNds16fe2vfZH1lNqNN94IAGhra4PL5UJkZCR2797t\nvoCalZWFN954AwBQXl6OuXPnQqvVIiYmBgaDAVVVVbK1nYKDp3mIel9VVRUMBgNiYmKg1WqRmZl5\nzWOXOu97U1JS0NraiubmZgDA3XffjcjIyGvW62t/7YusA05HRwcmTJgAnU6HqVOnYuzYsWhpaXFP\n79TpdGhpaQEANDY2eswQ0uv1cDqdsrSbgmPKlCk4efKk3M0gUj2n04no6Gj3d2/7Tyk5XfnaX/si\n67PUwsLCcPToUZw7dw5paWnYv3+/R3l3D2rk/46JiLondV8pvNxoHcg2ustXxMM7hwwZggceeACH\nDx+GTqdDc3MzRowYgaamJvfNa1FRUaivr3fXaWhoQFRU1DXr+pFGg7MhazkR0bViY2NRU1Nz3esZ\noNHgux7UGzhwIM6fP+/+3nX/WV9ff809ZVL3sZ352l/7ItuAc+rUKYSHh2Po0KG4ePEi9u3bh6ef\nfhpmsxnbtm3DihUrsG3bNsycORMAYDabMW/ePDzxxBNwOp1wOBxITk6+Zr1nAVi+/9nS6Wd/usvz\nVd4b8f0Apsqw3euJdy3ri33oHGf7rz9+vetSSh+klHvLswR4f5cv3wF4rgf1/rXLUyUmTpwIh8OB\nuro6jBw5EmVlZddMszebzSguLkZmZiYqKiowdOhQj6dXeONrf+2LbANOU1MTsrKy0NHRgY6ODsyf\nPx/Tpk1DUlIS5syZg02bNiEmJgavvvoqACAhIQFz5sxBQkICwsPDsWHDBp5SIyLV0wZhHeHh4Sgu\nLkZaWhpcLheys7MRHx+PkpISAEBOTg4yMjJgtVphMBgQERHhcT/d3Llz8d577+H06dOIjo7Gs88+\ni4ULF2LlypVe99c+2xGEvvRIYmIijhw5ck38Rz/6kc9HgaxevRqrV6/u7aYREalOeno60tPTPWI5\nOTke34uLi73W9XXTsb/9tTeKuIZDfxcjdwOCIEbuBlynGLkbcJ1i5G5AEMTI3QAFUdNOWk19UYVb\n5W5AEPT1PrD98lNDH4IlGKfUlIIDDhGRgqlpJ62mvhARqQ6PcIiIKCTUtJNWU1+IiFRHTUc4GtH1\nWQZ9HO/NISIlCMauVaPRYGcP6mUGafvBpsojHEunT4vPLM98f3m+yhmXf9uMKzOuxDYFGpda7i1P\nSr5UajrCUeWAQ0SkFhxwiIgoJNS0k1ZTX4iIVEdNRziyvoCNiIh+OHiEQ0SkYGraSaupL0REqqOm\nU2occIiIFExNO2k19YWISHV4hENERCGhpp20mvpCRKQ6PMIhIqKQUNNOWk19ISJSHR7hEBFRSHDA\nISKikFDTTpqPtiEiUjBteOCLNzabDXFxcTAajSgqKvKak5eXB6PRCJPJhOrq6m7rHjt2DHfccQfG\njx8Ps9mM8+fP++0LBxwiIpVzuVxYunQpbDYb7HY7SktLcfz4cY8cq9WKmpoaOBwObNy4Ebm5ud3W\nXbx4MX73u9/hk08+wc9+9jM8//zzftvBAYeISMHCwwNfuqqqqoLBYEBMTAy0Wi0yMzNRXl7ukbN7\n925kZWUBAFJSUtDa2orm5ma/dR0OB+6++24AwPTp0/H666/77QsHHCIiBdP2C3zpyul0Ijo62v1d\nr9fD6XRKymlsbPRZd+zYse7B57XXXkN9fb3fvnDAISJSsGAc4Wg0GknbEkIE1LbNmzdjw4YNmDhx\nIi5cuIAbbrjBf18CWjsREYWUr0kAnR1ov7L4EhUV5XH0UV9fD71e7zenoaEBer0e7e3tPuuOGTMG\nb7/9NgDgxIkTeOutt/w3VKgMAC5cuHCRfQnW/kwMD3zpuv329nZx2223idraWnHp0iVhMpmE3W73\nyHnrrbdEenq6EEKIQ4cOiZSUlG7rfvXVV0IIIVwul5g/f77YsmWL3/6o8gjH0unT4jPLM99fnq9y\nxuXfNuPKjCuxTYHGpZZ7y5OSL1kQ9tLh4eEoLi5GWloaXC4XsrOzER8fj5KSEgBATk4OMjIyYLVa\nYTAYEBERgS1btvitCwClpaX44x//CACYPXs2FixY0NtdISKiXhOkvXR6ejrS09M9Yjk5OR7fi4uL\nJdcFrty3k5eXJ7kNHHCIiJRMRXtpFXWFiEiFvExz7qs44BARKZmK9tIq6goRkQqpaC/NGz+JiCgk\nVDR2EhGpEK/hEBFRSKhoLy3bKbX6+npMnToVY8eOxbhx47Bu3ToAwJkzZ5CamorRo0fjvvvuQ2tr\nq7tOQUEBjEYj4uLisHfvXrmaTkQUOuE9WBRKtgFHq9XihRdewKeffoqKigr88Y9/xPHjx1FYWIjU\n1FScOHEC06ZNQ2FhIQDAbrejrKwMdrsdNpsNS5YsQUdHh1zNJyIKjX49WBRKtgFnxIgRmDBhAgBg\n4MCBiI+Ph9Pp9HgnQ1ZWFt544w0AQHl5OebOnQutVouYmBgYDAZUVVXJ1XwiotDgEU5w1dXVobq6\nGikpKWhpaYFOpwMA6HQ6tLS0AAAaGxs9nm7q7X0ORESqo6IBR/amXbhwAbNnz8aLL76IQYMGeZRp\nNBq/73HwVba/08+1AG4NQjuJiHyp/f5zv9+sHpJ9Lx08sh7htLe3Y/bs2Zg/fz5mzpwJ4MpRTXNz\nMwCgqakJw4cPB+D9XQ1RUVFe1zv1+wXgYENEve/qfqbzvidoVHQNR/P9OxdCTgiBrKwsDBs2DC+8\n8II7vnz5cgwbNgwrVqxAYWEhWltbUVhYCLvdjnnz5qGqqgpOpxPTp09HTU3NNUc5Ut9sR0TUm4Kx\na9VoNBDmHtTbHZztB5tsA84HH3yAe+65B+PHj3cPEgUFBUhOTsacOXNw8uRJxMTE4NVXX8XQoUMB\nAPn5+di8eTPCw8Px4osvIi0t7Zr1ajQaj3dSWK7JuFZ3eb7KGZd/24wrM67ENgUal1ruLc+CIA44\ns3pQ78/KHHBkOzt41113+ZzW/M4773iNr169GqtXr+7NZhERKYuCT5EFShGz1IiISP1UNP+BiEiF\nVLSXVlFXiIhUSEV7aRV1hYhIhVR0DYcDDhGRkqloL62irhARqZCK9tKcpUZEpGRBepaazWZDXFwc\njEYjioqKvObk5eXBaDTCZDKhurq627pVVVVITk5GUlISJk2ahI8++shvVzjgEBEpWRAebeNyubB0\n6VLYbDbY7XaUlpbi+PHjHjlWqxU1NTVwOBzYuHEjcnNzu627fPly/Pa3v0V1dTWeffZZLF++3G9X\nOOAQESlZEI5wqqqqYDAYEBMTA61Wi8zMTJSXl3vkdH41TEpKClpbW9Hc3Oy37o9//GOcO3cOANDa\n2urz+Zadu0JEREoVhL200+lEdHS0+7ter0dlZWW3OU6nE42NjT7rFhYW4q677sJvfvMbdHR04NCh\nQ37bIduz1HoLH95JREoQtGep9eBpXpp8z+2//vrrsNlsePnllwEAr7zyCiorK7F+/Xp3zowZM7By\n5UrceeedAIDp06ejqKgIdXV1PutOnz4dv/zlL/Gzn/0Mr732GjZu3Ih9+/b5bJcqj3AsnT4tPrM8\n8/3l+SpnXP5tM67MuBLbFGhcarm3PCn5kknYSx+ou7L40vX1LvX19R4vtPSW09DQAL1ej/b2dp91\nq6qq3M++fPDBB7F48WK/7eQ1HCKiPm5KDGCZ8velq4kTJ8LhcKCurg5tbW0oKyuD2ez53gOz2Yzt\n27cDACoqKjB06FDodDq/dQ0GA9577z0AwF/+8heMHj3abztVeYRDRKQaQdhLh4eHo7i4GGlpaXC5\nXMjOzkZ8fDxKSkoAADk5OcjIyIDVaoXBYEBERAS2bNnity4AbNy4Eb/85S9x6dIlDBgwABs3buzt\nrhARUa8J0l46PT0d6enpHrGcnByP78XFxZLrAleOnLpOPvCHAw4RkZLxWWpERBQSKtpLq6grREQq\npKK9tIq6QkSkQjylRkREIaGivTSfNEBE1AuC9qSB/9eDeouDs/1gU9HY+XeWTp8Wn1me+f7yfJUz\nLv+2GVdmXIltCjQutdxbnpR8yXhKjYiIQkJFe2kVdYWISIVUtJdWUVeIiFRIRXtpPryTiIhCQkVj\nJxGRCqlo0gCnRRMR9YKgTYve3YN6Zk6LDhlLp0+LzyzPfH95vsoZl3/bjCszrsQ2BRqXWu4tT0q+\nZCraS6uoK0REKqSiU2occIiIlExFe2kVdYWISIVUtJdWUVeIiFRIRafUOEuNiKgXBG2WWkUP6k3m\nLLWQsXT6tPjM8sz3l+ernHH5t824MuNKbFOgcanl3vKk5Eumor20irpCRKRCKtpLq6grREQqpKJr\nOLI+S23RokXQ6XRITEx0x86cOYPU1FSMHj0a9913H1pbW91lBQUFMBqNiIuLw969e+VoMhFRn2Sz\n2RAXFwej0YiioiKvOXl5eTAajTCZTKiuru62bmZmJpKSkpCUlIRbb70VSUlJftsg64CzcOFC2Gw2\nj1hhYSFSU1Nx4sQJTJs2DYWFhQAAu92OsrIy2O122Gw2LFmyBB0dHXI0m4godMJ7sHThcrmwdOlS\n2Gw22O12lJaW4vjx4x45VqsVNTU1cDgc2LhxI3Jzc7utu3PnTlRXV6O6uhqzZ8/G7Nmz/fdFyKy2\ntlaMGzfO/X3MmDGiublZCCFEU1OTGDNmjBBCiPz8fFFYWOjOS0tLE4cOHbpmfQC4cOHCRfYlGAAI\n4Qh86br9Dz/8UKSlpbm/FxQUiIKCAo+cnJwcsXPnTvf3MWPGiKamJkl1Ozo6RHR0tKipqfHbH8Vd\nw2lpaYFOpwMA6HQ6tLS0AAAaGxsxefJkd55er4fT6fS6DkunT4vXjGvz/eX5Kmdc/m0zrsy4EtsU\naFxqubc8KfmSBeEajtPpRHR0tPu7Xq9HZWVltzlOpxONjY3d1j148CB0Oh1iY2P9tkNxA05nGo3G\n7301vOeGiFQvCHtpqftK0cN7d0pLSzFv3rxu8xQ34Oh0OjQ3N2PEiBFoamrC8OHDAQBRUVGor693\n5zU0NCAqKsrrOvZ3+rkWwK292F4iotrvP/f7zeohCXvpAx9eWXzpuv+sr6+HXq/3m9PQ0AC9Xo/2\n9na/dS9fvoxdu3bhyJEj3bZTcW/8NJvN2LZtGwBg27ZtmDlzpju+c+dOtLW1oba2Fg6HA8nJyV7X\nMfX7BeBgQ0S97+p+pvO+J2gkTBKYcg9gWfn3pauJEyfC4XCgrq4ObW1tKCsrg9ls9sgxm83Yvn07\nAKCiogJDhw6FTqfrtu4777yD+Ph4jBw5UlJXZDN37ly89957OHXqFKKjo/Hss89i5cqVmDNnDjZt\n2oSYmBi8+uqrAICEhATMmTMHCQkJCA8Px4YNG3hKjYhUTwThGk54eDiKi4uRlpYGl8uF7OxsxMfH\no6SkBACQk5ODjIwMWK1WGAwGREREYMuWLX7rXlVWVoa5c+dKagefpUZE1AuCsWvVaDRoPxd4Pe2Q\n4Gw/2BR3DScYLJ0+LT6zPPP95fkqZ1z+bTOuzLgS2xRoXGq5tzwp+VK5VLSXVlFXiIjU53K/nlxq\nV+ZN8RxwiIgUzBXek910W9DbEQyKm6VGRETqxCMcIiIFc/VTz+OiOUuNiKgXBGuW2tdiYMD1btZc\n4Cy1ULF0+rT4zPLM95fnq5xx+bfNuDLjSmxToHGp5d7ypORLdVlFL8RR5YBDRKQWLhXtptXTEyIi\nFXLxCIeIiEJBTQMOJw0QEfWCYE0acAh994ldGDUNnDQQKpZOnxafWZ75/vJ8lTMu/7YZV2ZciW0K\nNC613FuelHypOGmAiIhCgpMGiIgoJNR0DYcDDhGRgqlpwOGz1IiIKCQ4S42IqBcEa5ZalRgXcL1k\nzV85Sy1ULJ0+LT6zPPP95fkqZ1z+bTOuzLgS2xRoXGq5tzwp+VJx0gAREYWEmq7hcMAhIlIwDjhE\nRBQSHHCIiCgk1PSkAc5SIyLqBcGapbZHTAm4XrrmwDXbt9lsWLZsGVwuFxYvXowVK1ZcUy8vLw97\n9uzBjTfeiK1btyIpKanbuuvXr8eGDRvQr18/PPDAAygqKvLZLlUe4Vg6fVp8Znnm+8vzVc54aLZh\ncUd7awvyxpXVmr7xbyIUcanl3vKk5EsVjFNqLpcLS5cuxTvvvIOoqChMmjQJZrMZ8fHx7hyr1Yqa\nmho4HA5UVlYiNzcXFRUVfuvu378fu3fvxieffAKtVouvv/7abzt44ycRkYK50C/gpauqqioYDAbE\nxMRAq9UiMzMT5eXlHjm7d+9GVlYWACAlJQWtra1obm72W/ell17CqlWroNVqAQA333yz375wwCEi\nUrDL6Bfw0pXT6UR0dLT7u16vh9PplJTT2Njos67D4cD777+PyZMnY8qUKfj444/99kWVp9SIiNQi\nGDd+Sr22Heh1p8uXL+Ps2bOoqKjARx99hDlz5uCLL77wmc8Bh4ioj/v0wCl8euC0z/KoqCjU19e7\nv9fX10Ov1/vNaWhogF6vR3t7u8+6er0es2bNAgBMmjQJYWFhOH36NIYNG+a1HZylRkTUC4I1S61U\nzAy43lzNGx7bv3z5MsaMGYN3330XI0eORHJyMkpLS6+ZNFBcXAyr1YqKigosW7YMFRUVfuuWlJSg\nsbERzzzzDE6cOIHp06fj5MmTPtulyiMcS6dPi88sz3x/eb7KGQ/NNjhLre/FldimQONSy73lScmX\nKhiz1MLDw1FcXIy0tDS4XC5kZ2e7BwwAyMnJQUZGBqxWKwwGAyIiIrBlyxa/dQFg0aJFWLRoERIT\nE3HDDTdg+/bt/ttx3T0hIqJeE6wbP9PT05Genu4Ry8nJ8fheXFwsuS4AaLVa7NixQ3IbOOAQESkY\nnxZNREQhwWepERFRSKhpwOEsNSKiXhCsWWrFIjvgeks1m/jGz1CxdPq0+MzyzPeX56uc8dBsg7PU\n+l5ciW0KNC613FuelHyp1PS0aFUOOEREasFJA0REFBJquobT5x7eabPZEBcXB6PR6Pe9C0REahCM\np0UrRZ+aNOByuTBmzBiP9zJ0fTwDJw0QkRIEa9JAgVgWcL1Vmv/bNycNrFu3DvPnz0dkZGQo2uNX\n5/cyAHC/l6HzgANw0kCo4729DU4a6HtxJbYp0LjUcm95UvKlUtOkgW5PqbW0tGDSpEmYM2cObDab\nrKOmlHc69KZ/6+jAv3V0MF/GbXR0/Bs6Ov7tB5OvtL+xEv9NhKIPcnIhPOBFqbodcNauXYsTJ05g\n0aJF2Lp1K4xGI1avXo3PP/88FO3zIPV02f7vFwCo7bXWEBFdcXU/03nfEyxquoYjaSgMCwvDiBEj\noNPp0K9fP5w9exYPPvggpk+fjueff7632+gm5Z0OADD1+8/3ANwaxO0/GxbYHIsfWn4othEW9uwP\nKl9pf2Ml/psIRR+6c3U/03nfEyxKHkAC1e2A8+KLL2L79u0YNmwYFi9ejN///vfQarXo6OiA0WgM\n6YAzceJEOBwO1NXVYeTIkSgrK0NpaWnItk9EFGpqGnC6naX29NNPY9GiRRg1atQ1ZXa7HQkJCb3W\nOG/27NmDZcuWud/LsGrVKo9yzlIjIiUI1iy1X4vfBlzv3zVrFDlLrU9Ni5ZCo9FwllqI4729Dc5S\n63txJbYp0LjUcm95FgRvwFkmCgKu9381qxQ54Ch3OgMREanqlBoHHCIiBeOAQ0REIaGmGz854BAR\nKZiSb+QMlConDRARyS1YkwayRXHA9TZplnLSQKhYOn1afGZ55vvL81XOeGi2wVlqfS+uxDYFGpda\n7i1PSr5UarqG0+deT0BE9EMSrEfbSHm1S15eHoxGI0wmE6qrq7uta7FYoNfrkZSUhKSkJNhsNr99\nUeURDhGRWgTjCMflcmHp0qUer3Yxm80eT9q3Wq2oqamBw+FAZWUlcnNzUVFR4beuRqPBE088gSee\neEJSO3iBukmiAAAYn0lEQVSEQ0SkYJfRL+Clq86vdtFqte5Xu3S2e/duZGVlAQBSUlLQ2tqK5ubm\nbusGcq2IAw4RkYIF4/UEUl7t4iunsbHRb93169fDZDIhOzsbra2tfvvCWWpERL0gWLPUZorAH1D8\nhmaux/Zff/112Gw2vPzyywCAV155BZWVlVi/fr07Z8aMGVi5ciXuvPNOAMD06dNRVFSEuro6n3W/\n+uor3HzzzQCANWvWoKmpCZs2bfLZLlVew7F0+rT4zPLM95fnq5zx0GyDs9T6XlyJbQo0LrXcW56U\nfKmkXMM5deBTnD7wqc9yKa926ZrT0NAAvV6P9vZ2n3WHDx/uji9evBgzZszw206eUiMiUjAp12yG\nThmPWMtc99JV51e7tLW1oaysDGaz2SPHbDZj+/btAICKigoMHToUOp3Ob92mpiZ3/V27diExMdFv\nX1R5hENEpBbBeNJAeHg4iouLkZaW5n61S3x8PEpKSgAAOTk5yMjIgNVqhcFgQEREBLZs2eK3LgCs\nWLECR48ehUajwa233upen892XHdPiIio1wTrxs/09HSkp6d7xHJycjy+Fxd7f6qBt7oA3EdEUnHA\nISJSMDU9aYCz1IiIekGwZqlNEXsCrndAk85nqYWKpdOnxWeWZ76/PF/ljIdmG5yl1vfiSmxToHGp\n5d7ypORLxdcTEBFRSKjp9QTq6QkRkQqp6RoOBxwiIgVT04DDSQNERL0gWJMGJohDAdc7qrmDkwZC\nxdLp0+IzyzPfX56vcsbl3zbjyowrsU2BxqWWe8uTki8VJw0QEVFIcNIAERGFhJqu4XDAISJSMA44\nREQUEmq6hsNZakREvSBYs9T0whFwvQaNkbPUQsXS6dPiM8sz31+er3LG5d8248qMK7FNgcallnvL\nk5IvlZpOqfEFbEREFBKqPMIhIlILNR3hcMAhIlIwVwcHHCIiCoHLl9Uz4HCWGhFRLwjWLLWB334d\ncL0LETdzllqoWDp9Wnxmeeb7y/NVzrj822ZcmXEltinQuNRyb3lS8qVyqegIR5UDDhGRWnDAISKi\nkLjcrp4BR5b7cF577TWMHTsW/fr1w5EjRzzKCgoKYDQaERcXh71797rjhw8fRmJiIoxGIx5//PFQ\nN5mISBYdrvCAF29sNhvi4uJgNBpRVFTkNScvLw9GoxEmkwnV1dWS6/77v/87wsLCcObMGb99kWXA\nSUxMxK5du3DPPfd4xO12O8rKymC322Gz2bBkyRL3ha/c3Fxs2rQJDocDDocDNptNjqYTEYXW5X6B\nL124XC4sXboUNpsNdrsdpaWlOH78uEeO1WpFTU0NHA4HNm7ciNzcXEl16+vrsW/fPowaNar7vggZ\nTZkyRRw+fNj9PT8/XxQWFrq/p6WliUOHDonGxkYRFxfnjpeWloqcnByv6wTAhQsXLrIvwQBA4POO\nwJcu2//www9FWlqa+3tBQYEoKCjwyMnJyRE7d+50fx8zZoxoamrqtu6DDz4ojh07JmJiYsTp06f9\n9kdR13AaGxsxefJk93e9Xg+n0wmtVgu9Xu+OR0VFwel0+lyPpdOnxWeWZ76/PF/ljMu/bcaVGVdi\nmwKNSy33liclX7LL13+rh9PpRHR0tPu7Xq9HZWVltzlOpxONjY0+65aXl0Ov12P8+PGS2tFrA05q\naiqam5uviefn52PGjBm9tVkiIupC6v2JIoB7dy5evIj8/Hzs27dPcv1eG3A6N0KqqKgo1NfXu783\nNDRAr9cjKioKDQ0NHvGoqCif69nf6edaALcG3BIiIulqv//c7zerhy5LyKk6AHx0wGdx131rfX29\nx1kjbzlX97/t7e1e637++eeoq6uDyWRy599+++2oqqrC8OHDvbZD9qdFdx4RzWYzdu7ciba2NtTW\n1sLhcCA5ORkjRozA4MGDUVlZCSEEduzYgZkzZ/pc59TvF4CDDRH1vqv7mc77nqC5LGH5yRQgx/L3\npYuJEyfC4XCgrq4ObW1tKCsrg9ls9sgxm83Yvn07AKCiogJDhw6FTqfzWXfcuHFoaWlBbW0tamtr\nodfrceTIEZ+DDSDTfTi7du1CXl4eTp06hQceeABJSUnYs2cPEhISMGfOHCQkJCA8PBwbNmxwHwpu\n2LABCxYswMWLF5GRkYH7779fjqYTEYWWlCOcboSHh6O4uBhpaWlwuVzIzs5GfHw8SkpKAAA5OTnI\nyMiA1WqFwWBAREQEtmzZ4rduV1JO2/FZakREvSAYu1aNRgNU9GA9kzV8llqoWDp9Wnxmeeb7y/NV\nzrj822ZcmXEltinQuNRyb3lS8iVzBXNl8lLlgENEpBpBOKWmFBxwiIiUjAMOERGFBAccIiIKCRUN\nOJylRkTUC4I2S628B+v5P5ylFjKWTp8Wn1me+f7yfJUzLv+2GVdmXIltCjQutdxbnpR8yVR0hCP7\nkwaIiOiHQZVHOEREqtEudwOChwMOEZGSqejGT04aICLqBUGbNLCtB+vJ4qSBkLF0+rT4zPLM95fn\nq5xx+bfNuDLjSmxToHGp5d7ypORLpqJJA6occIiIVIMDDhERhQQHHCIiCgkOOEREFBIccIiIKCR4\nHw4REYUE78NRLt6HQ0RKELT7cJ7pwXqe5n04IWPp9GnxmeWZ7y/PVznj8m+bcWXGldimQONSy73l\nScn/IVLlgENEpBqcNEBERCGhogGHrycgIlKy9h4sXthsNsTFxcFoNKKoqMhrTl5eHoxGI0wmE6qr\nq7utu2bNGphMJkyYMAHTpk1DfX29365wwCEiUjJXD5auq3C5sHTpUthsNtjtdpSWluL48eMeOVar\nFTU1NXA4HNi4cSNyc3O7rbt8+XIcO3YMR48excyZM/HMM8/47QoHHCIiJbvcg6WLqqoqGAwGxMTE\nQKvVIjMzE+Xl5R45u3fvRlZWFgAgJSUFra2taG5u9lt30KBB7voXLlzATTfd5LcrvIZDRKRkQbiG\n43Q6ER0d7f6u1+tRWVnZbY7T6URjY6Pfuk899RR27NiBG2+8ERUVFX7bwSMcIiIlC8I1HKn3J/bk\n3p21a9fi5MmTWLBgAX71q1/5zeURDhGRkkl50kDTAaD5gM/iqKgojwv69fX10Ov1fnMaGhqg1+vR\n3t7ebV0AmDdvHjIyMvy3U6gMAC5cuHCRfQna/my+CHzpsv329nZx2223idraWnHp0iVhMpmE3W73\nyHnrrbdEenq6EEKIQ4cOiZSUlG7rnjhxwl1/3bp14pFHHvHbH1Ue4Vg6fVp8Znnm+8vzVc64/Ntm\nXJlxJbYp0LjUcm95UvIlC8I1nPDwcBQXFyMtLQ0ulwvZ2dmIj49HSUkJACAnJwcZGRmwWq0wGAyI\niIjAli1b/NYFgFWrVuGzzz5Dv379EBsbi5deesl/O66/K0RE1GuC9LTo9PR0pKene8RycnI8vhcX\nF0uuCwB/+tOfAmoDJw0QEVFI8AiHiEjJVPR6Ag44RERKpqJnqXHAISJSMg44REQUEnzFNBERhQSv\n4RARUUio6JSaLNOin3zyScTHx8NkMmHWrFk4d+6cu6ygoABGoxFxcXHYu3evO3748GEkJibCaDTi\n8ccfl6PZREShF4SnRSuFLAPOfffdh08//RTHjh3D6NGjUVBQAACw2+0oKyuD3W6HzWbDkiVL3A+T\ny83NxaZNm+BwOOBwOGCz2eRoOhFRaAXpBWxKoBGiB48HDaJdu3bh9ddfxyuvvIKCggKEhYVhxYoV\nAID7778fFosFo0aNwr333ut+6c/OnTtx4MAB/Md//Mc165P6VFQiot4UjF2rRqMBknqwnmpNULYf\nbLJfw9m8eTPmzp0LAGhsbMTkyZPdZVffx6DVaj2eThoVFQWn0+lznZZOnxafWZ75/vJ8lTMu/7YZ\nV2ZciW0KNC613FuelHzJFHyKLFC9NuCkpqaiubn5mnh+fj5mzJgB4Mp7FG644QbMmzcvqNve3+nn\nWgC3BnXtRESear//3O83q4c44HRv3759fsu3bt0Kq9WKd9991x3z9T6GqKgoNDQ0eMSjoqJ8rnvq\n95/vgYMNEfW+q/uZzvseupYskwZsNhuef/55lJeXo3///u642WzGzp070dbWhtraWjgcDiQnJ2PE\niBEYPHgwKisrIYTAjh07MHPmTDmaTkQUWiqaNCDLNZzHHnsMbW1tSE1NBQDccccd2LBhAxISEjBn\nzhwkJCQgPDwcGzZscE8C2LBhAxYsWICLFy8iIyMD999/vxxNJyIKLd74eX0cDofPstWrV2P16tXX\nxG+//Xb87//+b282i4hIeXgNh4iIQoIDDhERhYSCr8kEigMOEZGS8RoOERGFBE+pERFRSHDAISKi\nkFDRNRxZbvwkIiKJXD1YvLDZbIiLi4PRaERRUZHXnLy8PBiNRphMJlRXV3db19+rZrzhgENEpGSi\nB0sXLpcLS5cuhc1mg91uR2lpqfvp+1dZrVbU1NTA4XBg48aNyM3N7baur1fN+MIBh4hI5aqqqmAw\nGBATEwOtVovMzEyUl5d75OzevRtZWVkAgJSUFLS2tqK5udlv3dTUVISFhbnrdH7mpTcccIiIVM7p\ndCI6Otr9/eqrX6TkNDY2dlsXuPKqmYyMDP8NESqDnh2AcuHChUtQl+Dtz0QPFs/t/+lPfxKLFy92\nf9+xY4dYunSpR84//dM/iQ8++MD9fdq0aeLjjz+WVPe5554Ts2bN6rY/qpylZun0afGZ5ZnvL89X\nOePyb5txZcaV2KZA41LLveVJyQ+uA98v3nV99Ut9fb3HSy295Vx9PUx7e7vfut5eNeMLT6kRESma\nlPcR3AngqU6Lp4kTJ8LhcKCurg5tbW0oKyuD2Wz2yDGbzdi+fTsAoKKiAkOHDoVOp/Nb19erZnxR\n5REOEZF6XP+dn+Hh4SguLkZaWhpcLheys7MRHx+PkpISAEBOTg4yMjJgtVphMBgQERGBLVu2+K0L\n+H7VjM92XHdPiIioFwXnzs/09HSkp6d7xHJycjy+FxcXS64L+H/VjDcccIiIFE09z7bhgENEpGjq\nebYNBxwiIkXjgENERCHBU2pERBQS6jnC4X04REQUEjzCISJSNJ5SIyKikFDPKTUOOEREisYjHCIi\nCgke4RARUUjwCIeIiEKCRzhERBQSPMIhIqKQ4BEOERGFBI9wiIgoJHiEQ0REIaGeAUcjhBByNyKY\nNBqN3E0gIkIwdq1X9mev9KDmI0HZfrCp8gjH0unT4jPLM99fnq9yxuXfNuPKjCuxTYHGpZZ7y5OS\n/0OkygGHiEg91HNKjQMOEZGicZYaERGFhHqOcGR5AduaNWtgMpkwYcIETJs2DfX19e6ygoICGI1G\nxMXFYe/eve744cOHkZiYCKPRiMcff1yOZhMRyeByD5Zr2Ww2xMXFwWg0oqioyGtOXl4ejEYjTCYT\nqquru6372muvYezYsejXrx+OHDnSbU9kGXCWL1+OY8eO4ejRo5g5cyaeeeYZAIDdbkdZWRnsdjts\nNhuWLFninmmRm5uLTZs2weFwwOFwwGazydH0XlcrdwOCoK/3ge2Xnxr6EDztPVg8uVwuLF26FDab\nDXa7HaWlpTh+/LhHjtVqRU1NDRwOBzZu3Ijc3Nxu6yYmJmLXrl245557JPVElgFn0KBB7p8vXLiA\nm266CQBQXl6OuXPnQqvVIiYmBgaDAZWVlWhqasL58+eRnJwMAPjFL36BN954Q46m97o6uRsQBHVy\nN+A61cndgOtUJ3cDgqBO7gYoyvUf4VRVVcFgMCAmJgZarRaZmZkoLy/3yNm9ezeysrIAACkpKWht\nbUVzc7PfunFxcRg9erTknsh2Deepp57Cjh07MGDAAFRVVQEAGhsbMXnyZHeOXq+H0+mEVquFXq93\nx6OiouB0OkPeZiKi0Lv+azhOpxPR0dHu73q9HpWVld3mOJ1ONDY2dltXql47wklNTUViYuI1y3//\n938DANauXYuTJ09i4cKFWLZsWW81g4ioj7v+IxypN8T3+s2iQmZffvmlGDt2rBBCiIKCAlFQUOAu\nS0tLExUVFaKpqUnExcW54//1X/8lcnJyvK4vNjZWAODChQsX2ZbY2Nig7B97uv2BAwd6rOfQoUMi\nLS3N/T0/P18UFhZ65OTk5IjS0lL39zFjxojm5mZJdadMmSIOHz7cfX8C6n2QnDhxwv3zunXrxCOP\nPCKEEOLTTz8VJpNJXLp0SXzxxRfitttuEx0dHUIIIZKTk0VFRYXo6OgQ6enpYs+ePXI0nYioz2lv\nbxe33XabqK2tFZcuXRImk0nY7XaPnLfeekukp6cLIa4MUCkpKZLrTpkyRXz88cfdtkOWazirVq3C\nZ599hn79+iE2NhYvvfQSACAhIQFz5sxBQkICwsPDsWHDBveh4IYNG7BgwQJcvHgRGRkZuP/+++Vo\nOhFRnxMeHo7i4mKkpaXB5XIhOzsb8fHxKCkpAQDk5OQgIyMDVqsVBoMBERER2LJli9+6ALBr1y7k\n5eXh1KlTeOCBB5CUlIQ9e/b4bIfqHt5JRETKJMu06GDo6zePPvnkk4iPj4fJZMKsWbNw7tw5d1lf\naD/g/6avvtKHzqTcGKcEixYtgk6nQ2Jiojt25swZpKamYvTo0bjvvvvQ2trqLvP1t5BLfX09pk6d\nirFjx2LcuHFYt24dgL7Th++++w4pKSmYMGECEhISsGrVKgB9p/2yCsoJQhl888037p/XrVsnsrOz\nhRB/vw7U1tYmamtrRWxsrPs60KRJk0RlZaUQQsh+HWjv3r3C5XIJIYRYsWKFWLFihRCi77RfCCGO\nHz8uPvvss2suGPalPlx1+fJlERsbK2pra0VbW5vX89RK8f7774sjR46IcePGuWNPPvmkKCoqEkII\nUVhY6Pff09V/d3JpamoS1dXVQgghzp8/L0aPHi3sdnuf6sO3334rhLhyfSMlJUUcPHiwT7VfLn32\nCKev3zyampqKsLArv/6UlBQ0NDT0qfYDvm/66kt9uErKjXFKcffddyMyMtIj1vmmvaysLPfv1dvf\n4up9b3IZMWIEJkyYAAAYOHAg4uPj4XQ6+1QfbrzxRgBAW1sbXC4XIiMj+1T75dJnBxzgys2jt9xy\nC7Zu3eo+rG1sbPS4SbTzzUtKvXl08+bNyMjIANA3299VX+yDr5ve+oqWlhbodDoAgE6nQ0tLCwDf\nfwulqKurQ3V1NVJSUvpUHzo6OjBhwgTodDr36cG+1H65KPpp0ampqWhubr4mnp+fjxkzZmDt2rVY\nu3YtCgsLsWzZMvesCqXorv3AlRtgb7jhBsybNy/UzZNESh/UQE1vitVoNH77o5S+XrhwAbNnz8aL\nL77occYCUH4fwsLCcPToUZw7dw5paWnYv3+/R7nS2y8XRQ84+/btk5Q3b9489xFCVFSUxwSChoYG\n6PV6REVFuU9bXY1HRUUFt8FddNf+rVu3wmq14t1333XHlNR+QPrfoDOl9UGKrm2ur6/3+F+p0ul0\nOjQ3N2PEiBFoamrC8OHDAXj/Wyjhd97e3o7Zs2dj/vz5mDlzJoC+1wcAGDJkCB544AEcPny4T7Y/\n1PrsKTWHw+H+uby8HElJSQAAs9mMnTt3oq2tDbW1tXA4HEhOTsaIESMwePBgVFZWQgiBHTt2uP+h\ny8Fms+H5559HeXk5+vfv7473lfZ3JTrNru+LfZg4cSIcDgfq6urQ1taGsrIymM1muZslmdlsxrZt\n2wAA27Ztc/9eff0t5CSEQHZ2NhISEjwea9VX+nDq1Cn3DLSLFy9i3759SEpK6jPtl5WsUxauw+zZ\ns8W4ceOEyWQSs2bNEi0tLe6ytWvXitjYWDFmzBhhs9nc8Y8//liMGzdOxMbGiscee0yOZrsZDAZx\nyy23iAkTJogJEyaI3Nxcd1lfaL8QQvz5z38Wer1e9O/fX+h0OnH//fe7y/pKHzqzWq1i9OjRIjY2\nVuTn58vdHJ8yMzPFj3/8Y6HVaoVerxebN28Wp0+fFtOmTRNGo1GkpqaKs2fPuvN9/S3kcvDgQaHR\naITJZHL/+9+zZ0+f6cMnn3wikpKShMlkEomJieJ3v/udEEL0mfbLiTd+EhFRSPTZU2pERNS3cMAh\nIqKQ4IBDREQhwQGHiIhCggMOERGFBAccIiIKCQ44REQUEhxwiIgoJDjgEAH46KOPYDKZcOnSJXz7\n7bcYN24c7Ha73M0iUhU+aYDoe2vWrMF3332HixcvIjo6GitWrJC7SUSqwgGH6Hvt7e2YOHEiBgwY\ngEOHDv1gHyFP1Ft4So3oe6dOncK3336LCxcu4OLFi3I3h0h1eIRD9D2z2Yx58+bhiy++QFNTE9av\nXy93k4hURdEvYCMKle3bt+Mf/uEfkJmZiY6ODvz0pz/FgQMHMGXKFLmbRqQaPMIhIqKQ4DUcIiIK\nCQ44REQUEhxwiIgoJDjgEBFRSHDAISKikOCAQ0REIcEBh4iIQoIDDhERhcT/B4CDZAn4sGs7AAAA\nAElFTkSuQmCC\n", + "text": [ + "" + ] + } + ], + "prompt_number": 8 + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "expmap = Maps.ExpMap(mesh)\n", + "m2to3 = Maps.Map2Dto3D(mesh,normal='Y')\n", + "imap = Maps.IdentityMap(mesh)\n", + "problem = DC.ProblemDC(mesh, mapping= imap )\n", + "problem.Solver = MumpsSolver\n", + "problem.pair(survey)" + ], + "language": "python", + "metadata": {}, + "outputs": [], + "prompt_number": 9 + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "datasyn = survey.dpred(sighalf*np.ones(mesh.nC))\n", + "survey.makeSyntheticData(sigma,std=0.01,force=True)\n", + "plot(np.log(np.c_[survey.dobs,datasyn]))" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "metadata": {}, + "output_type": "pyout", + "prompt_number": 10, + "text": [ + "[,\n", + " ]" + ] + }, + { + "metadata": {}, + "output_type": "display_data", + "png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAEACAYAAABfxaZOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xt4VNW9//H3nsmEAEkmJEAykEAUEnIhhHCRYhVDcaIQ\noaBwvFChWulprRY8SkXPr21qjxBFi6C0nkOPLV6KFGsEESJJYRSkgSMBQUBBSCCQTLiEBHJPZvbv\nj20CmBuZJOzMzPf1PPuZ25pZ34n4mTVr9l5bUVVVRQghhFcw6F2AEEKI60dCXwghvIiEvhBCeBEJ\nfSGE8CIS+kII4UUk9IUQwou4HPolJSVYrVaio6NJSUmhtLS0xbYOh4OkpCSmTp3qandCCCE6gcuh\nn56ejtVq5ciRI0yaNIn09PQW2y5fvpy4uDgURXG1OyGEEJ3A5dDfsGEDc+fOBWDu3Ll88MEHzbY7\ndeoUmzZt4pFHHkGOAxNCCH25HPrFxcWEhoYCEBoaSnFxcbPtnnjiCZYuXYrBID8fCCGE3nxae9Bq\ntWK325vc//zzz191W1GUZqduNm7cSP/+/UlKSsJms3WsUiGEEB2numjYsGFqUVGRqqqqWlhYqA4b\nNqxJm2eeeUYNDw9XIyMj1bCwMLVXr17qgw8+2OzrDRkyRAVkk0022WRrxzZkyJB2ZbfLob9w4UI1\nPT1dVVVVXbJkifr000+32t5ms6l33XVXy4Xgcindwm9/+1u9S3CZO9euqlK/3qR+fbU3O12eaF+0\naBFZWVlER0ezdetWFi1aBEBhYSGpqanNPsed9t45cgSWLoXf/Ab27AH5DVoI4QlandNvTXBwMNnZ\n2U3uHzBgAB999FGT+2+77TZuu+02V7vrdHY7VFVpYe50apfFxfDhh7B+PVy6BNOmQUAA3H8/1NTA\njBnadsstYDTq/Q6EEKL9XA59d1NVBZ98Aps3Q2YmnDunBbrBAIqibWYzTJ4Mb70Fo0drjwG88AIc\nOgTvvw/z50NFBTz1FMydC35+Wpvk5OQW+66rg7ffhs8+g9tuA6sVwsK6/j1fq9ZqdwdSv76kfvei\nfDsnpDtFUTp9P36nE7Zsgddfh61bITFRC/XJk7XrruxFqqqwfbv2QZCbC7/8Jfz85xAU1LRtQ9j/\n139BZKT2zWH7dq2WiAgt/KdPh+9/X/vQEUKI9mpvdrpF6NfWwv/9H+zfr42sAwK0LTBQC9v+/aFP\nn8shXlICf/kL/OlPWptf/ALuuaf5YO6IAwfgxRdh40aIiYHBgy9vqgrLlmnXf/tbmDDh8vPq6+Hz\nz+Hjj2HNGi3wH3kE5syBfv06t0YhhGdz69BfsHkBEeYIIgIHUXoynNzdfuR+buTAAQORg4wMjzdQ\nX6dQUQHl5QoVFQplZXD+nEJlJQQHa1thocKkH8CDD8LIkQ3TN103lD5/Ho4fh9OnobBQu7x0CR54\nAL73vcvtmvtTO1WVzz+HtWshOxtuuUVl4kRIGKFyww0qBgOoqDhVJw6nQ7tUHVTXV3Pq4ilOlp3k\nZNlJ8koKuC1iEv9x66Nd9j6FEN2PW4d++icvkf35Sf7v6wJqepwiIKiGXr2d+PVygKKFXUO5Kt9e\nNt4GpwPq61V8fMBwxQ+t3eQtAs1/+Cho9zmdUF4O1dUKNTXgcCj4+Sr08AOTjxEfgwGTjxGT0YCv\nsQem6nDqzg6iJC+CM1VFDEw4SsHS9df7LQkhdOTWoR8aqpKYCE88ASkprs25e5Jz57RpoL174cwZ\nbdrqwgVtq63VfpcYOxbGjIGPDm7jD7m/o+QPNr3LFkJcR+0N/W61984//wnx8XpX0X307Qt33qlt\nbcm1B1KjlHV9UUIIt9atxtIS+K4bEGymziChL4RoXbcKfeG6AcFmHD4X9S5DCNHNSeh7iIF9A3Ga\nyrrVj9ZCiO5HQt9DhAT1ANXIpeoqvUsRQnRjEvoeQlFAqTFTcFbm9YUQLZPQ9yDGejNF52VeXwjR\nMgl9D2Jymim6ICN9IUTLJPQ9SA81kOJSCX0hRMsk9D2In2Lm7CUJfSFEy1wO/ZKSEqxWK9HR0aSk\npFBaWtpsu9LSUmbOnElsbCxxcXHk5OS4XKxoXS+jmfPlMqcvhGiZy6Gfnp6O1WrlyJEjTJo0ifT0\n9GbbzZ8/nylTpnD48GH2799PbGysy8WK1vn7mCmplJG+EKJlLof+hg0bmDt3LgBz587lgw8+aNKm\nrKyM7du38/DDDwPg4+OD2Wx2tUvRhoAegZRWSegLIVrmcugXFxcTGhoKQGhoKMXFxU3a5OXl0a9f\nPx566CFGjRrFvHnzqKysdL1a0aqgHmbKaiT0hRAta3WVTavVit1ub3L/888/f9VtRVGaXSe+vr6e\n3NxcXnvtNcaOHcuCBQtIT0/nueeea7a/tLS0xuvJycled+7KjurTy8zRskN6lyGE6EI2mw2bzeby\n811eTz8mJgabzUZYWBhFRUVMnDiRr7766qo2drud8ePHk5eXB8COHTtIT09n48aNTQvpgnPkepsF\nq9bx3uG1nPrDe3qXIoS4TtqbnS5P70ybNo3Vq1cDsHr1aqZPn96kTVhYGBERERw5cgSA7Oxs4mX9\n5C7TLyCQKlWmd4QQLXM59BctWkRWVhbR0dFs3bqVRYsWAVBYWEhqampju1dffZXZs2eTmJjI/v37\nefbZZztetWhWf7OZGiT0hRAt61anS+wmpbitD3MOM3PdDGpe/qrtxkIIj3DdpndE9zMg2Ey9UUb6\nQoiWSeh7kPB+Zpy+EvpCiJZJ6HuQfuZeYKylsrpO71KEEN2UhL4HMRgUlNpATp2T9XeEEM2T0Pcw\nxjozp8/JFI8QonkS+h7GR06kIoRohYS+h5ETqQghWiOh72H8FDNnL8qcvhCieRL6HqaXwcy5chnp\nCyGaJ6HvYfx9zJRUSOgLIZonoe9hAnwDuSAnUhFCtEBC38OY/cxcrJE5fSFE8yT0PUyfXmbK62Sk\nL4RonoS+hwnpbaaiXkJfCNE8CX0P0zcgkEqnhL4QonkS+h5GO5GKzOkLIZrncuiXlJRgtVqJjo4m\nJSWF0tLSZtstWbKE+Ph4EhISeOCBB6ipqXG5WNG2sCAztQYZ6Qshmudy6Kenp2O1Wjly5AiTJk0i\nPT29SZv8/HxWrVpFbm4uBw4cwOFw8O6773aoYNE6OZGKEKI1Lof+hg0bmDt3LgBz587lgw8+aNIm\nMDAQk8lEZWUl9fX1VFZWMnDgQNerFW0a0DcQp0lCXwjRPJdDv7i4mNDQUABCQ0MpLi5u0iY4OJgn\nn3ySQYMGMWDAAIKCgrj99ttdr1a0KSwoEHzLqa5x6l2KEKIb8mntQavVit1ub3L/888/f9VtRVFQ\nFKVJu2PHjvHKK6+Qn5+P2Wxm1qxZvPPOO8yePbvZ/tLS0hqvJycnk5ycfA1vQVzJx2iE+l4Univn\nxoGBepcjhOhkNpsNm83m8vMVtT2nUb9CTEwMNpuNsLAwioqKmDhxIl999dVVbdauXUtWVhZ//vOf\nAXjrrbfIyclh5cqVTQtp5xndRct8fhXOttn/4tbECL1LEUJ0sfZmp8vTO9OmTWP16tUArF69munT\npzdpExMTQ05ODlVVVaiqSnZ2NnFxca52Ka6RyRlIYYnM6wshmnI59BctWkRWVhbR0dFs3bqVRYsW\nAVBYWEhqaioAiYmJzJkzhzFjxjBixAgAfvrTn3ZC2aI1vk4zZ8pkX30hRFMuT+90Npne6TyhT97J\nT0fM5/dzJ+tdihCii1236R3RffWUE6kIIVogoe+BevsEcl5CXwjRDAl9DxToa6a0Wub0hRBNSeh7\nILOfmbJqGekLIZqS0PdAfXqauSQnUhFCNENC3wMF9w6UE6kIIZoloe+B+vqbqXLKnL4QoikJfQ/U\n32ymGhnpCyGaktD3QGFBZmoVCX0hRFMS+h7IEhJInZxIRQjRDAl9DzQwxIzTJHP6QoimJPQ9UKjZ\nDD3KqK6WtYyEEFeT0PdAfqYeoBooLqnWuxQhRDcjoe+hDPWBnD4n8/pCiKtJ6HsoU72ZohKZ1xdC\nXE1C30P5qmbspTLSF0JczeXQX7duHfHx8RiNRnJzc1tsl5mZSUxMDFFRUbzwwguudifaqQdmzlyU\n0BdCXM3l0E9ISCAjI4MJEya02MbhcPDYY4+RmZnJoUOHWLNmDYcPH3a1S9EOvQyBnJXQF0J8h4+r\nT4yJiWmzze7duxk6dCiRkZEA3Hfffaxfv57Y2FhXuxXXqLePmQsVMqcvhLhal87pnz59moiIiMbb\n4eHhnD59uiu7FN8K8DVzoUpG+kKIq7U60rdardjt9ib3L168mKlTp7b54oqitKuYtLS0xuvJyckk\nJye36/niMnMPM2VlEvpCeBqbzYbNZnP5+a2GflZWlssvDDBw4EAKCgoabxcUFBAeHt5i+ytDX3RM\nUE8zhefkW5UQnua7A+Lf/e537Xp+p0zvqGrzh/uPGTOGo0ePkp+fT21tLWvXrmXatGmd0aVoQ3Cv\nQMrru2ZOX1WhqAiyskC+TAjhXlwO/YyMDCIiIsjJySE1NZXJkycDUFhYSGpqKgA+Pj689tpr3HHH\nHcTFxXHvvffKj7jXSd8AM1WOzkvkr7+GX/0K7rgDwsIgIUG7nZIC5eWuvabTCZcuQWGh9vquvo4Q\n4topakvD9OtMUZQWvzGI9lux6WN+vellyl7b0uHXOnYMJkyAOXPgllsgMREGDtQemzcPTp6EjRvB\n17ft19qzB158ETZvhooK6NkTAgLA3x9KSuDBB+Gxx2Do0A6XLYRXaG92yhG5HiosyExNJ5xIpahI\nG83/+tewZAmkpkJ4OCiKtr3+OvTurX0gOBzNv4aqwpYtMGkSzJgB48dDXh7U1Wmj+6IiOHoU9u4F\nPz/t8alTITtbe64QovNI6HuosD6B1Bs6Nqd/4YI2nfPww/CznzXfxscH1qyB4mL45S+vDumTJ+EP\nf4CkJHjySfjxj7VvDQsWQEgIGL7zr2/QIEhPhxMnYNo0ePRRePbZDr0FIcR3yPSOhzp29jRDXxiL\n+lKhS8+vqNBG+OPGwcsva6P61ly8CMnJYLVqc/7r1sGRIzB9Otx3nzbKb+cevJw7B9//Pjz+uDbl\nI4Roqr3Z6fIRuaJ7u3wiFW3KpD1qa2HWLIiKgpdeurawDgzU5unvvx8iI+E3v9GC3mRyqXwA+vaF\nzEztdwSLBe65x/XXEkJoZKTvoVRVxZBmouDfqwgfcO3JW1sL//ZvWtCvW6dN3+ht715tmukf/4Bb\nb9W7GiG6F/khVwDaPwRDXQCF5y9d83MaRviKAmvXdo/AB+03gXfegZkz4eBBvasRwr1J6HswH4eZ\nwpJr24OnpkYLVaNRC/xr2f3yerJatR+FJ0/W9ukXQrimm4zlRFfwdZqxX2g79BsC39cX3n23Y/Pw\nXWn2bK3W226DjAxt104hRPvISN+D9cBMcSvrJFRVaUspTJsGPXp078Bv8PDD8MYbWs3r1+tdjRDu\nR0b6HqynIZCDRy+yc+fl++rqYNcuLexzcmDECO2Aq4ULu3/gN5gyBTZtgh/+UDuwq6VjCIQQTcne\nOx5s/NLZnMuZTL+iHzXepyjaD6NWq7ZfvdmsX30ddewY3Hkn3H67NtefkACDB18+6Ku2FnJzYedO\n7QNu/nxtv38hPEl7s1NC34M9tukxdpzcwciwkfTr1Y++vfoS0iuEnj498TH4XLUZFAMGxaDt9aMY\nUFBQFOWqS7h8joTv3r5Sw2OtUbn6v/WV/+0bHlNVtcn1hkun6sSpOiktVXnvfQfH8+s5UVBPZXU9\nAyPqwVTByfPnCLKco0/4OS5Un2d68P/jf359swt/SSG6Lwl90ai4vJhdp3dxrvLcVVuNo4Z6Z33j\nVueouypIVVW7fmXIwuVg/u7tK303zBvaXcuHw5VtrvxQ+e71735AGRRD44eXWm+i/KIPPvQkbnA/\nBgT1pW+vvizbvB7fs2PZ9fIiV/6UQnRbckSuaBTqH8q0YXL+AoBPPz9Hdv4JvcsQQney947wClFh\nFi7UNT31pxDepkOhv27dOuLj4zEajeTm5jbbpqCggIkTJxIfH8/w4cNZsWJFR7oUwiXxgy2UK0Wy\nVLPweh0K/YSEBDIyMpgwYUKLbUwmE8uWLePgwYPk5OSwcuVKDh8+3JFuhWi3oaEWVP8iOb2j8Hod\nmtOPiYlps01YWBhhYWEA+Pv7ExsbS2FhoZw2UVxXAwIt4F9EQQEEBeldjRD6ua5z+vn5+ezdu5dx\n48Zdz26FIMA3AMXg5Ku8a1+ATghP1OZI32q1Yrc3/QFs8eLFTJ069Zo7Ki8vZ+bMmSxfvhx/f//2\nVSlEBymKQm/VwsETRcwiQO9yhNBNm6GflZXV4U7q6uq45557+NGPfsT06dNbbJeWltZ4PTk5meTk\n5A73LUSDPj4WvrEXAdF6lyKEy2w2GzabzeXnd8rBWRMnTuSll15i9OjRTR5TVZW5c+cSEhLCsmXL\nWi5EDs4SXeympbPwPXYPO16/T+9ShOg01/UkKhkZGURERJCTk0NqaiqTJ08GoLCwkNTUVAA+++wz\n3n77bbZt20ZSUhJJSUlkZmZ2pFshXBLex0LRpSK9yxBCV7IMg/AaT61fwqq3LlD23ot6lyJEp5HT\nJQrRgmEDLJRjx+HQuxIh9COhL7zGoD4WfPoU0czOaEJ4DQl94TUsARaMZu0ALSG8lYS+8BoWfwuO\nnkWcPKl3JULoR0JfeI2QXiHUGy9x/ESN3qUIoRsJfeE1DIqBAEN/vi6USX3hvST0hVfp62fh+BnZ\nV194Lwl94VUGBFg4VSqhL7yXhL7wKpEhFs5USegL7yWhL7zKDf0sVPkUUVWldyVC6ENCX3iVAQEW\nevWXffWF95LQF17FEmDBN9guoS+8loS+8CoWf+1cuXKAlvBWEvrCq4T5h1HrK6EvvJeEvvAqof6h\nVCpnOVEgS20K7yShL7yKr9EXf58gjhWd1bsUIXThcuivW7eO+Ph4jEYjubm5rbZ1OBwkJSW160Tq\nQnSVsN4WTpTIvvrCO7kc+gkJCWRkZDBhwoQ22y5fvpy4uDgURXG1OyE6TXiQBXt5EXKiNuGNXA79\nmJgYoqOj22x36tQpNm3axCOPPCKnQxTdQkSQBYO5iJIS155/9iwUyRcF4aa6fE7/iSeeYOnSpRgM\n8vOB6B4s/hYCB7i2B8/x4zBuHCQkwKpVyLcF4XZ8WnvQarVib+bccosXL76m+fmNGzfSv39/kpKS\nsNlsLhcpRGeyBFjw6/c1BQWQlHTtzzt0CO64A/7zP+GWW+DHP4Z//EML/4gI12rZsQPeeAPMZggK\nurxFREBcHISGgsyKis7UauhnZWV16MV37tzJhg0b2LRpE9XV1Vy8eJE5c+bw5ptvNts+LS2t8Xpy\ncjLJyckd6l+I5lj8LRgCP2nXSD83F1JTYelS+NGPtPv+9S948UUYNQrS0+H++6Fnz2sP6exs7TnP\nPqvdLi3VvklcuAAnTsDBg+B0Qmys9gEwbRpMngw+rf5fKzydzWbr2CBa7aDk5GT1888/b7OdzWZT\n77rrrhYf74RShLgm209sVwf9bry6cOG1td+xQ1X79VPV999v/vH9+1X1e99T1R49VNXXV1X791fV\nYcNU9dZbtec4nU2fs3mz9pqfftp632fOqOonn6jqihWqevPNqmqxqOrTT6vqV19dW+3C87U3O12e\naM/IyCAiIoKcnBxSU1OZPHkyAIWFhaSmpjb7HNl7R3QHYf5hVBlbX3StpgaysuCJJ2DGDHj7be2y\nOQkJ2qi/uhrKyuCLL+CDD+Cpp+DXvwarVRu1N/jwQ5gzB9avh1tvbb3Wfv1gwgR4/HH47DPYulX7\nHSE5GSZO1PoToj2Ubz8pdKcoiuzdI66LitoKgl8IYXRmFTs/uzwQqanRwv3DD2HbNoiPhylTYNYs\nGDbMtb7q6+H11+G557SpnNGjYeFC2LgRxo51/T3U1cEvfgGXLsHf/ibz/t6svdkpoS+8UsDiQAL/\n9wSnj/UBYMsWeOwxGDpUm7NPSYG+fTuvv7NntVH/+vWwaVP7fkBuSVWV9sHx1FPaj8rCO0noC3EN\nol8dRn56Bl/viGPhQu2H2hUr4K67urZfVe3cUfmXX2rTPJ99Btdw2IzwQO3NTtl5XnilAQEWAgcW\nkZQEw4drc+5dHfjQ+dMww4fD734H992nTU8J0RYJfeGVLAEWZv97EXv2QFqatqulu/r5z2HQoMu7\nfgrRGgl94ZUs/hbCY4oYMkTvSjpOUeB//xf+/nfIzNS7GtHdSegLr2Txt1BU7jkL6ISEwH//Nyxa\npHcloruT0BdeyRJgwV7edIkRdzZpEhw5ou3VI0RLJPSFV/K0kT5Ajx7a8QT79+tdiejOJPSFVwrz\nD6PokmeFPmgHf7VxTiPh5ST0hVeyBHjeSB+0xd/27NG7CtGdSegLr9THrw819TVU1lXqXUqnGjVK\nRvqidbJIq/BKiqIQ5h/G01lPEx4YTpBfEGY/M4E9AvE1+uJr9MVkMGEymjAZTBgUA0aDUbtUtEtF\nUVBQUBRFu43S+NpXXr+SqqqoqI3XAVRUnKqz8TFV1W47VId26dQu65x11DnqqHPWUeuopdZRy8Wa\ni5RWl1JWXUZZTRmDA6L46quHqKnR5viF+C5ZhkF4rY+OfMQ++z7Kasooqy6jtKaUizUXqXXUXhWu\ndY46nKqzSRA3BHRDaEPzof5dzX0oNHxoNDz23Q8Yg2LAZDQ1fhg1fDAF9gjE3MOM2c8MwLpD6/D/\ncx5//as2vy88n6y9I4SXqnPU4b/En1lHL3LbLT2YN0/visT1IGvvCOGlTEYTg8yDGDTyuMzrixZJ\n6AvhQaKCozBHfiN78IgWuRz669atIz4+HqPRSG4rw4rS0lJmzpxJbGwscXFx5OTkuNqlEKINQ4OH\n4uxzlC+/1E60IsR3uRz6CQkJZGRkMGHChFbbzZ8/nylTpnD48GH2799PbGysq10KIdoQFRzFyfKj\nDB4Mhw/rXY3ojlwO/ZiYGKLbOGtDWVkZ27dv5+GHHwbAx8cHs9nsapdCiDZEhUTxzYVv5CAt0aIu\nndPPy8ujX79+PPTQQ4waNYp58+ZRWelZB8MI0Z0MDR7K0fNHZTkG0aJWD86yWq3Y7U1XIly8eDFT\np05t88Xr6+vJzc3ltddeY+zYsSxYsID09HSee+65ZtunpaU1Xk9OTiY5ObnNPoQQl0UGRWIvtzP8\n5mree89P73JEF7DZbNhsNpef3+H99CdOnMjLL7/MqFGjmjxmt9sZP348eXl5AOzYsYP09HQ2btzY\ntBDZT1+IThH9ajRv37WeHyTEUlYGRqPeFYmupMt++i11GBYWRkREBEeOHAEgOzub+Pj4zuhSCNGC\nocFDsdcexWKBr7/WuxrR3bgc+hkZGURERJCTk0NqaiqTJ08GoLCwkNTU1MZ2r776KrNnzyYxMZH9\n+/fzrJzIU4guFRUcxTcl8mOuaJ4swyCEh3lt92scPHOQyIN/wm6HZcv0rkh0JVmGQQgvFxUcxdES\n2YNHNE9CXwgPMzR4KN+UfENSEuzdC06n3hWJ7kRCXwgPMzhoMPZyO73N1QQHwzff6F2R6E4k9IXw\nMD4GHwYHDeb4heNXTfGoKlRWQmEh1NfrW6PQj4S+EB6oYYpn7Fj4xS8gNBT8/CAkBEaOhBtugLQ0\nKCjo2jpKS+Gjj+DFF+Hkya7tS1wbOV2iEB4oKjiKo+ePsmABTJ8OQUHa5vftQbr798P//I/2AXDz\nzTBvHowYoX049OzZ/v6cTigqgvx8OH4cdu2C7du16zfdBDfeCElJMHcuPPss9O3bqW9XtIPssimE\nB1q5eyUHzhzg9bteb7VdRQX8/e/w5ptaQNvtWuiHhkJYGFgsly8tFu2xoiJtiqjh8uRJbTObtW8Q\nkZHaqRpvvVU7UbvJpPVlt8Pvfw9r18L8+fDEE+Dv3/V/C08np0sUQvDxNx+zdOdSsudkt+t5qqpN\nydjtWqjb7ZevFxVBVZUW/gMGaJvFAhERWtD36nVtfRw7po328/O1bwSiYyT0hRAcv3CciasncmLB\nCb1LaVZtrTbddO7ctX9YiObJwVlCCAaZB1FcXkx1fbXepTTL1xdiYuDAAb0r8T4S+kJ4oCt32+yu\nGg4eE9eXhL4QHqphD57uauRI2LdP7yq8j4S+EB6qYQ2e7kpG+vqQ0BfCQzUcoNVdjRgBX34pRwdf\nbxL6QnioqJDuPdIPDNR2+/z2HEviOpHQF8JDdfc5fdDm9WWK5/pyOfTXrVtHfHw8RqOR3FYW7V6y\nZAnx8fEkJCTwwAMPUFNT42qXQoh2iDBHcKbiDFV1VXqX0qKkJPkx93pzOfQTEhLIyMhgwoQJLbbJ\nz89n1apV5ObmcuDAARwOB++++66rXQoh2sHH4ENkUKTstimu4vKCazExMW22CQwMxGQyUVlZidFo\npLKykoEDB7rapRCinRrm9eP7x+tdSrMadttUVVAUvavxDl06px8cHMyTTz7JoEGDGDBgAEFBQdx+\n++1d2aUQ4gpD+wxl16ldnLp4irLqMhxOh261qKpKZV0lxeXFfFPyDUWXirBYwMcHTp3SrSyv0+pI\n32q1Yrfbm9y/ePFipk6d2uaLHzt2jFdeeYX8/HzMZjOzZs3inXfeYfbs2c22T0tLa7yenJxMcnJy\nm30IIVr2gxt+wFNZT/Hm/je5VHOJiroKehh70NPUE1+jLyaDSbs0mjAqRowGIwbFgFHRLhVFQUFB\n+XYYrqCgoqKq6lWXDqcDp+rEoWqXdY46ah211Dm1y1pHLRW1FZiMJgJ8A+hl6oWiKOTNz2sc7UdE\n6PzHchM2mw2bzeby8zu84NrEiRN5+eWXGTVqVJPH1q5dS1ZWFn/+858BeOutt8jJyWHlypVNC5EF\n14Tock7VSWVdJTX1NY1h3BDM3w1uh9PRGOxA4/UrPwgUFAyKQfug+PYDw6AYGj9MGj5QfI2+9Db1\nxmQ0NdbR54U+HPvlMV7+fV969oTf/EbPv4z7am92dspJVFrqMCYmht///vdUVVXh5+dHdnY2N910\nU2d0KYRTHmcZAAANgUlEQVRwgUEx4O/rj7+vvgvZGxQDSWFJ7Cncw8iRd7B2ra7leBWX5/QzMjKI\niIggJyeH1NRUJk+eDEBhYSGpqakAJCYmMmfOHMaMGcOIESMA+OlPf9oJZQsh3N1oy2hyi3JlD57r\nTNbTF0Lo4m8H/sb7h9/n7zPfw2zWztcbFKR3Ve5H1tMXQriF0ZbR7Cnag8GgrcMjB2ldHxL6Qghd\nRIVEcb7yPOcrz8syy9eRhL4QQhcGxUCSJUnm9a8zCX0hhG4apnhk4bXrR0JfCKGbUZZR5BblMnw4\nHD0K1d3zlL4eRUJfCKGbhpG+nx8MHQoHD7b9nJoacDpbflxVtWUdDh8Gh36rTnRbnXJwlhBCuCI6\nJJozFWe4UHWBpKQ+7NsHo0c33/bcOVi6FP70J6irg8GDITISbrgB+veHvDwt6L/+Gnr31razZ+Gm\nm+Dmm2H8eEhOBj+/6/kOux8Z6QshdGM0GBkZNpLcolxGj4YVK+APf4Ddu7VgBygpgf/8Txg2DC5d\n0r4NnDsH778Pjz8O8fHa6H7iRHj1VThxAoqK4JtvtO2Xv9Re6+mn4dln9X2/3YEcnCWE0NWCzAUM\nCBjA/DG/4r334LPPtO3YMRg1Cg4dgrvv1oJ/8GDX+8nJgZ/9zPN2DW1vdkroCyF09dYXb7Hx6EbW\nzrx6AZ7SUm3EP3Qo3Hhjx/uprYXgYDh9Gszmjr9edyFH5Aoh3MroAaPZU7inyf1BQZCS0jmBD+Dr\nC2PHws6dnfN67kpCXwihq2EhwyiuKKa0urTL+7rlFtixo8u76dYk9IUQujIajCSGJrK3qOuPzrr1\nVti+vcu76dYk9IUQuhtlGcWeoqZTPJ3te9+D3FxtX39vJaEvhNBdw0FaXS0wUNv18/PPu7yrbktC\nXwihu5Z+zO0K3j6v73LoL1y4kNjYWBITE7n77rspKytrtl1mZiYxMTFERUXxwgsvuFyoEMJzxfSN\nofBSIWXVzedIZ/L2eX2XQz8lJYWDBw/yxRdfEB0dzZIlS5q0cTgcPPbYY2RmZnLo0CHWrFnD4cOH\nO1SwEMLz+Bh8GBE6gn32rj9y6pZbtN02W1u/x5O5HPpWqxWDQXv6uHHjOHXqVJM2u3fvZujQoURG\nRmIymbjvvvtYv36969UKITzWaMtoPj3xKWXVZdQ76zvtdR1OB5dqLlFcXkyto5awMAgJubbF3TxR\npyy49sYbb3D//fc3uf/06dNEREQ03g4PD2fXrl2d0aUQwsOkDElh3ofzWLpzKRV1FRgVI719e9PD\n2AMfg0/jZjQYMShXj1dVVaXeWX/VVuuopbKuklpHLb1MvVAUhZ8k/YRX7nylcV4/IUGnN6ujVkPf\narVit9ub3L948WKmTp0KwPPPP4+vry8PPPBAk3aKorSrmLS0tMbrycnJJCcnt+v5Qgj3NXXYVOzD\ntLxRVZVaRy0VdRXUOmqbBHpzyw6YjCbtQ0Ex4mPwwWQ00dvUGz8fPxRFYZ99H/e+dy+v8Aq33grZ\n2fDzn1/vd9lxNpsNm83m8vM7tPbOX//6V1atWsU///lP/JpZrzQnJ4e0tDQyMzMBWLJkCQaDgaef\nfrppIbL2jhCiCzlVJ5aXLex6ZBe1ZyK5/XY4eVLvqjruuq29k5mZydKlS1m/fn2zgQ8wZswYjh49\nSn5+PrW1taxdu5Zp06a52qUQQrjMoBiw3mgl61gWUVHaAVqeEPrt5XLoP/7445SXl2O1WklKSuLR\nRx8FoLCwkNTUVAB8fHx47bXXuOOOO4iLi+Pee+8lNja2cyoXQoh2ShmSwpbjW1AUbS+ea9l1s7oa\nDhzo+tquF1laWQjhNYouFRH/x3jOLjzLqyuMfP21diauKzkc2knas7Phn//U1uH//vfh21nqbkeW\nVhZCiBZYAiyEB4bzeeHnTUb6J0/CM8+AxQJz52pn33r8ce18u9018F0h58gVQniVlCEpbDm2hWe+\nP46TJ2HDBnjzTdi6FebM0c7aFRWld5VdR0b6Qgiv0jCv7+OjLcnwq19p59c9cQJeecWzAx9kTl8I\n4WWq6qro/1J/Tv/HafyUQHx8wODGw1+Z0xdCiFb0NPVkfPh4tuVtw9fXvQPfFV72doUQAqw3Wtly\nbIveZehCQl8I4XUa5vW9kYS+EMLrJIQmcKnmEscvHNe7lOtOQl8I4XUMigHrEG1JBm8joS+E8Eop\nN3rnFI+EvhDCK91+4+1szdvaqSdscQeyn74Qwmslvp5IbN9YgnsGYzKYGk/U4lAd1DvrqXPUUees\nI8gviBetL+pdbrPam52yDIMQwmu9PeNt/nXqX9Q56rSQd2qXRsXYeFIWk8GE2c+sd6mdRkb6Qgjh\nxuSIXCGEEC3qUOgvXLiQ2NhYEhMTufvuuykrK2vSpqCggIkTJxIfH8/w4cNZsWJFR7oUQgjRAR0K\n/ZSUFA4ePMgXX3xBdHQ0S5YsadLGZDKxbNkyDh48SE5ODitXruTw4cMd6bZb6siJivXmzrWD1K83\nqd+9dCj0rVYrhm9XKxo3bhynTp1q0iYsLIyRI0cC4O/vT2xsLIWFhR3ptlty53847lw7SP16k/rd\nS6fN6b/xxhtMmTKl1Tb5+fns3buXcePGdVa3Qggh2qHNXTatVit2u73J/YsXL2bq1KkAPP/88/j6\n+vLAAw+0+Drl5eXMnDmT5cuX4+/v34GShRBCuEztoL/85S/qzTffrFZVVbXYpra2Vk1JSVGXLVvW\nYpshQ4aogGyyySabbO3YhgwZ0q7M7tB++pmZmTz55JN88skn9O3bt9k2qqoyd+5cQkJCWLZsmatd\nCSGE6AQdCv2oqChqa2sJDg4GYPz48fzxj3+ksLCQefPm8dFHH7Fjxw4mTJjAiBEjUBQFgCVLlnDn\nnXd2zjsQQghxzbrNEblCCCG6nu5H5GZmZhITE0NUVBQvvPCC3uW06eGHHyY0NJSEhITG+0pKSrBa\nrURHR5OSkkJpaamOFbaupYPl3OU9VFdXM27cOEaOHElcXBzPPPMM4D71AzgcDpKSkhp3hHCn2iMj\nIxkxYgRJSUncdNNNgHvVX1paysyZM4mNjSUuLo5du3a5Tf1ff/01SUlJjZvZbGbFihXtrl/X0Hc4\nHDz22GNkZmZy6NAh1qxZ0+0P3HrooYfIzMy86r709HSsVitHjhxh0qRJpKen61Rd21o6WM5d3oOf\nnx/btm1j37597N+/n23btrFjxw63qR9g+fLlxMXFNU53ulPtiqJgs9nYu3cvu3fvBtyr/vnz5zNl\nyhQOHz7M/v37iYmJcZv6hw0bxt69e9m7dy979uyhV69ezJgxo/31t+tn3062c+dO9Y477mi8vWTJ\nEnXJkiU6VnRt8vLy1OHDhzfeHjZsmGq321VVVdWioiJ12LBhepXWbj/84Q/VrKwst3wPFRUV6pgx\nY9Qvv/zSbeovKChQJ02apG7dulW96667VFV1r38/kZGR6rlz5666z13qLy0tVW+44YYm97tL/Vf6\n+OOP1VtuuUVV1fbXr+tI//Tp00RERDTeDg8P5/Tp0zpW5Jri4mJCQ0MBCA0Npbi4WOeKrs2VB8u5\n03twOp2MHDmS0NDQxqkqd6n/iSeeYOnSpY1HsoN7/ftRFIXbb7+dMWPGsGrVKsB96s/Ly6Nfv348\n9NBDjBo1innz5lFRUeE29V/p3Xff5f777wfa//fXNfQbvt56EkVR3OJ9lZeXc88997B8+XICAgKu\neqy7vweDwcC+ffs4deoUn376Kdu2bbvq8e5a/8aNG+nfvz9JSUktLoXbXWtv8Nlnn7F37142b97M\nypUr2b59+1WPd+f66+vryc3N5dFHHyU3N5fevXs3mQrpzvU3qK2t5cMPP2TWrFlNHruW+nUN/YED\nB1JQUNB4u6CggPDwcB0rck1oaGjjUctFRUX0799f54paV1dXxz333MODDz7I9OnTAfd7DwBms5nU\n1FT27NnjFvXv3LmTDRs2cMMNN3D//fezdetWHnzwQbeovYHFYgGgX79+zJgxg927d7tN/eHh4YSH\nhzN27FgAZs6cSW5uLmFhYW5Rf4PNmzczevRo+vXrB7T//11dQ3/MmDEcPXqU/Px8amtrWbt2LdOm\nTdOzJJdMmzaN1atXA7B69erGIO2OVFXlJz/5CXFxcSxYsKDxfnd5D+fOnWvcO6GqqoqsrCySkpLc\nov7FixdTUFBAXl4e7777Lj/4wQ9466233KJ2gMrKSi5dugRARUUFW7ZsISEhwW3qDwsLIyIigiNH\njgCQnZ1NfHw8U6dOdYv6G6xZs6Zxagdc+H+3i39vaNOmTZvU6OhodciQIerixYv1LqdN9913n2qx\nWFSTyaSGh4erb7zxhnr+/Hl10qRJalRUlGq1WtULFy7oXWaLtm/friqKoiYmJqojR45UR44cqW7e\nvNlt3sP+/fvVpKQkNTExUU1ISFBffPFFVVVVt6m/gc1mU6dOnaqqqvvUfvz4cTUxMVFNTExU4+Pj\nG/9/dZf6VVVV9+3bp44ZM0YdMWKEOmPGDLW0tNSt6i8vL1dDQkLUixcvNt7X3vrl4CwhhPAiuh+c\nJYQQ4vqR0BdCCC8ioS+EEF5EQl8IIbyIhL4QQngRCX0hhPAiEvpCCOFFJPSFEMKL/H8tQsfWhVks\noAAAAABJRU5ErkJggg==\n", + "text": [ + "" + ] + } + ], "prompt_number": 10 }, { "cell_type": "code", "collapsed": false, "input": [ - "txList = DC.Examples.WennerArray.getTxList(nElecs,aSpacing)" + "problem.mapping = expmap * m2to3" + ], + "language": "python", + "metadata": {}, + "outputs": [], + "prompt_number": 11 + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "dmis = DataMisfit.l2_DataMisfit(survey)\n", + "reg = Regularization.Tikhonov(mesh,mapping=m2to3)\n", + "opt = Optimization.InexactGaussNewton(maxIter=7,tolX=1e-15)\n", + "opt.remember('xc')\n", + "invProb = InvProblem.BaseInvProblem(dmis, reg, opt)\n", + "beta = Directives.BetaEstimate_ByEig(beta0_ratio=1e1)\n", + "betaSched = Directives.BetaSchedule(coolingFactor=5, coolingRate=2)\n", + "inv = Inversion.BaseInversion(invProb, directiveList=[beta,betaSched])" ], "language": "python", "metadata": {}, @@ -99,14 +234,55 @@ "cell_type": "code", "collapsed": false, "input": [ - "fig, ax = plt.subplots(1,1, figsize = (6.5,5))\n", - "indz = 23\n", - "dat = mesh.plotSlice(sigma, grid=True, ax = ax, ind=indz)\n", - "ax.set_xlim(-300, 300)\n", - "ax.set_ylim(-300, 300)\n", - "cb = plt.colorbar(dat[0])\n", - "ax.set_title('Depth at '+str(mesh.vectorCCz[indz])+' m')\n", - "ax.plot(xyz[:,0], xyz[:,1], 'w.', ms = 3)" + "m0 = np.log(np.ones(problem.mapping.nP)*sighalf)\n", + "mopt = inv.run(m0)" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "SimPEG.InvProblem will set Regularization.mref to m0.\n", + "SimPEG.InvProblem is setting bfgsH0 to the inverse of the eval2Deriv.\n", + " ***Done using same solver as the problem***\n", + "SimPEG.l2_DataMisfit is creating default weightings for Wd." + ] + }, + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "\n", + "============================ Inexact Gauss Newton ============================" + ] + }, + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "\n", + " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", + "-----------------------------------------------------------------------------\n", + " 0 6.56e+00 1.91e+03 0.00e+00 1.91e+03 2.21e+03 0 " + ] + }, + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "\n", + " 1 6.56e+00 2.10e+02 1.61e+01 3.15e+02 1.80e+02 0 " + ] + } + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "colorbar(mesh.plotSlice(np.log10(expmap * m2to3 * mopt), normal='Y', ind=indz, grid=True)[0])" ], "language": "python", "metadata": {}, @@ -114,61 +290,61 @@ { "metadata": {}, "output_type": "pyout", - "prompt_number": 8, + "prompt_number": 31, "text": [ - "[]" + "" ] }, { "metadata": {}, "output_type": "display_data", - "png": "iVBORw0KGgoAAAANSUhEUgAAAZwAAAFRCAYAAABANxtTAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X9YVHW+B/D3KOOmpEKpgzLquMwQkDpSCrtr9eg1UmiZ\ntWxd9K5i4r08GLFubZnseqNt/dWP7cmI0vskgncXqb0ptuGklrqtG1CK1TreHAo2GH48mWJa6uDw\nvX+YEwPMMMBwzpnT+/U858E58/me7/crPt+P3/PjezRCCAEiIqIBNkjuBhAR0fcDEw4REUmCCYeI\niCTBhENERJJgwiEiIkkw4RARkSSYcCioGQwGvP3223I3g4j8wIRDAWMwGDBs2DCMGDEC4eHhmDlz\nJrZs2YJAPeq1bNkyrF271mOfRqOBRqMJyPENBgPeeeedXpVpaGhAamoqbrzxRowdOxYPPvggXC4X\nAOBPf/oThg8f7t5CQ0MxaNAgVFdXd3usWbNmYejQoe742NjYfveJSEmYcChgNBoN/vrXv+Krr77C\n559/jsceewybNm1CRkaG3E3zi0aj6XVyzMnJwahRo9DU1ITjx4/j8OHDKCgoAAD8+7//O86fP+/e\nCgoKEBUVhfj4eK/1v/jii+74kydP9rtPRErChEMDYvjw4UhNTUVpaSmKiopw4sQJAMDly5fxm9/8\nBhMnTkRERASysrJw6dIlAMChQ4eg1+uxYcMGjB49GpMmTcKf//xnAMDWrVvx5z//GU899RSGDx+O\nn/3sZ+66qqurYTabERYWhrS0NFy+fLnbNn366af4t3/7N4waNQqjR4/GL3/5S5w7dw4AsGTJEnz+\n+edITU3F8OHD8cwzz/jVzxMnTuAXv/gFhgwZAp1Oh3nz5rn72tn27duxdOlSn8fzN+Ft374dM2fO\nxEMPPYTw8HAYjUb84x//QGFhISZMmACdTofi4mK/jkUkGUEUIAaDQbz99ttd9k+YMEG8/PLLQggh\nVq1aJX72s5+Js2fPivPnz4vU1FSxZs0aIYQQBw8eFCEhIeLhhx8WTqdTHD58WISGhopTp04JIYRY\ntmyZWLt2rcexJ06cKBITE0VTU5M4c+aMiI2NddfVWU1NjThw4IBwOp3iiy++EHfccYdYtWpVj+33\n5cEHHxRLliwR33zzjWhoaBCTJ08Wu3fv7hJXV1cnBg8eLOrq6rwea9asWWL06NFi1KhRYubMmeLQ\noUNeYwsLC0VISIjYvn27aG9vF7/73e9EZGSkyM7OFk6nU+zbt08MHz5cfP31173qD9FA4gyHBty4\nceNw5swZCCHw3//93/jjH/+IsLAwXH/99VizZg127tzpEf/kk09Cq9XijjvuwN13343S0lIAV//3\nLzrNADQaDXJychAREYHw8HCkpqbi+PHj3bYjKioKc+bMgVarxahRo/DrX/8ahw8f7lff8vLy8M9/\n/hMjRozA+PHjMWPGDI/Z1zXFxcW44447MHHiRK/H2rRpE2pra9HY2Ij//M//RGpqKj777DOv8ZMm\nTUJ6ejo0Gg0WLlyIxsZG/Nd//Re0Wi2SkpIwZMgQ1NTU9Kt/RIHEhEMDrqGhATfccANOnz6Nb775\nBrfeeivCw8MRHh6O5ORknD592h0bHh6OoUOHuj9PnDgRTU1NAOD15oCIiAj3n4cOHYoLFy50G9fS\n0oK0tDTo9XqMHDkSS5YswZdfful3P5KTk90X9EtKSgAAc+fOxc9//nN88803OH36NM6cOYPVq1d3\nKVtcXIz09HSfx09ISEBoaCi0Wi2WLl2KmTNnory83Gu8Tqdz//na39no0aM99nn7uyCSAxMODaj3\n338fjY2NuO2223DjjTdi6NChsNlsOHv2LM6ePYvW1lZ89dVX7vizZ8/im2++cX/+17/+hXHjxgHw\nnnA68hWTm5uLwYMH45///CfOnTuHHTt2oL293a+yALB37173Bf1Fixbhiy++wNGjR5GdnQ2tVosb\nbrgBy5Yt65Ikjhw5gqamJtx33309tp9IzZhwKKCunfL66quv8Ne//hWLFi3CkiVLcPPNN2PQoEH4\nj//4D6xatQpffPEFAMDhcGDfvn0ex3j88cfR1taGd999F2+++SZ+/vOfA7j6P3pfp5g61t+dCxcu\nIDQ0FCNGjIDD4cDTTz/t8b1Op8Onn37qd19HjRqFsWPH4qWXXoLL5UJrayuKiopgNps94oqKinDf\nffchNDTU67HOnTuHt956C5cuXcKVK1fwpz/9Ce+++y7mzZvnd3uIlI4JhwIqNTUVI0aMwIQJE7Bh\nwwY8/PDDKCwsdH+/adMmGI1G/OhHP8LIkSORlJSEU6dOub+/di1m3LhxWLJkCbZs2YLo6GgAQEZG\nBmw2G8LDw3Hvvfd2W7+v53Ief/xxHDt2DCNHjkRqaioWLFjgEbtmzRr84Q9/QHh4OP74xz/22FeN\nRoPXX38db7zxBkaNGgWTyYQf/OAHeO6559wxly5dwmuvvdbt6bT169cjJSUFAOB0OrF27VqMGTMG\no0ePxosvvoiysjIYjUa/+xmo55GIBoxcdytcvHhRJCQkCLPZLGJjY8Vjjz0mhBDiyy+/FHfeeacw\nmUwiKSlJnD171l1m/fr1wmg0iptuukm89dZbcjWdBsjBgweFXq+XuxlEqrR3715x0003CaPRKDZu\n3NhtzIMPPiiMRqOYOnWqOHbsmHv//fffL8aMGSMmT57sEe9rvO6ObDOc6667DgcPHsTx48fx0Ucf\n4eDBg/j73/+OjRs3uv/XO2fOHGzcuBEAYLPZUFpaCpvNBqvVipUrV3qcfyciou65XC5kZ2fDarXC\nZrOhpKSky4PF5eXlqKmpgd1ux9atW5GVleX+7v7774fVau1yXG/jtTeynlIbNmwYgKunE1wuF8LD\nw7Fnzx736Yf09HTs3r0bAFBWVoZFixZBq9XCYDDAaDSiqqpKtrbTwOBpIaLAq6qqgtFohMFggFar\nRVpaGsrKyjxiOo69iYmJaG1tRXNzMwDg9ttvR3h4eJfjehuvvZE14bS3t2PatGnQ6XSYPXs2br75\nZrS0tLhv99TpdGhpaQEANDY2Qq/Xu8vq9Xo4HA5Z2k0DY9asWfj888/lbgaR6jgcDowfP979ubvx\n05+YzryN196E9LbhgTRo0CAcP34c586dw9y5c3Hw4EGP73tamJH/GyYi6pm/Y6Xo5sHq3tTRU7ys\nCeeakSNH4u6778bRo0eh0+nQ3NyMiIgINDU1YcyYMQCAyMhI1NfXu8s0NDQgMjKyy7Fu0GhwVrKW\nExF1FRUVFZBVHoZqNLjUh3LXX389zp8/7/7cefysr6/3OGPUXYy3MbYjb+O1N7IlnNOnTyMkJARh\nYWG4ePEi9u/fj8cffxwWiwVFRUVYvXo1ioqKMH/+fACAxWLB4sWL8dBDD8HhcMButyMhIaHLcc8C\nyOuw9VVvywcq/iCA2QpqT1/KBEsfvMWz/X2LD2QdwfI78FY+rxfPc/lyCcAf+lDud51WmJg+fTrs\ndjvq6uowbtw4lJaWulfLuMZisSA/Px9paWmoqKhAWFiYx2oW3fE2XnvV73vt+uijjz4S8fHxwmw2\niylTpoinnnpKCHH1Nrs5c+Z0e5vdunXrRFRUlLjpppuE1Wrt9rgAuHHjxk32LRAAiE192Lqrv7y8\nXERHR4uoqCixfv16IYQQL7/8ssditw888ICIiooSU6dOFUePHnXvT0tLE2PHjhVDhgwRer1ebNu2\nrcfxujuabzulGhqNhjMcCeJ9lQmWPniLZ/v7Fh/IOoLld+CtfB4QkBcPajQabOpDudUBqj/QFHEN\nh75jkLsBAWCQuwH9ZJC7Af1kkLsBAWCQuwEKoqZBWk19UYVJcjcgAIK9D2y//NTQh0DRyt2AAGLC\nISJSMDUN0mrqCxGR6nCGQ0REklDTIK2mvhARqQ5nOEREJAk1DdJq6gsRkepwhkNERJJgwiEiIkmo\naZBWU1+IiFRHTTMcWV/ARkRE3x+c4RARKZiaBmlVrhZNRCS3QK0W/fc+lLstQPUHmpqSp1seArfM\nOOPlq4PxwRUvRR1Ki/dWvj/H6ExNg7Sa+kJEpDpqummACYeISMHUNEirqS9ERKrDGQ4REUlCTYO0\nmvpCRKQ6nOEQEZEkmHCIiEgSahqkubQNEZGCaUN6v3XHarUiJiYGJpMJmzZt6jYmJycHJpMJZrMZ\n1dXVPZb98MMP8eMf/xhTp06FxWLB+fPnffaFCYeISOVcLheys7NhtVphs9lQUlKCkydPesSUl5ej\npqYGdrsdW7duRVZWVo9lV6xYgaeeegofffQR7rnnHjz99NM+28GEQ0SkYCEhvd86q6qqgtFohMFg\ngFarRVpaGsrKyjxi9uzZg/T0dABAYmIiWltb0dzc7LOs3W7H7bffDgC488478b//+7++OyNUBgA3\nbty4yb4Fajz7JrT3W+f6X3vtNbFixQr35x07dojs7GyPmJ/+9KfiyJEj7s9z5swRH3zwgfjLX/7i\ntexPfvITsXv3biGEEM8++6wYPny4z/6o6XqUWx4CtyYS4+Wrg/HBFS9FHUqL91a+P8forLsZS2/5\nu6ix6OWCn9u2bUNOTg6efPJJWCwWDBkyxGe8KhMOEZFaeLsJoKNDbVc3byIjI1FfX+/+XF9fD71e\n7zOmoaEBer0ebW1tXsvedNNNeOuttwAAp06dwptvvumznbyGQ0SkZIN73mZdB+QN/27rbPr06bDb\n7airq4PT6URpaSksFotHjMViQXFxMQCgoqICYWFh0Ol0Pst+8cUXAID29nb84Q9/cN9o4A1nOERE\nShaAUTokJAT5+fmYO3cuXC4XMjIyEBsbiy1btgAAMjMzkZKSgvLychiNRoSGhqKwsNBnWQAoKSnB\niy++CABYsGABli1bNtBdISKiAROgUTo5ORnJycke+zIzMz0+5+fn+10WuPrcTk5Ojt9tYMIhIlIy\nFY3SKuoKEZEKDZa7AYHDhENEpGQqGqVV1BUiIhVS0SjN26KJiEgSGtHbR0sVzt8naomIBlIghlaN\nRgMxrQ/ljgem/kBT0WTtO3kI3BIVjJevDsYHV7wUdSgt3lv5/hyjCxWN0rKdUquvr8fs2bNx8803\nY/Lkydi8eTMA4MyZM0hKSkJ0dDTuuusutLa2usts2LABJpMJMTEx2Ldvn1xNJyKSTkgfNoWSLeFo\ntVo899xzOHHiBCoqKvDiiy/i5MmT2LhxI5KSknDq1CnMmTMHGzduBADYbDaUlpbCZrPBarVi5cqV\naG9vl6v5RETS8GNpmy6bQsmWcCIiIjBt2tWTk9dffz1iY2PhcDg83smQnp6O3bt3AwDKysqwaNEi\naLVaGAwGGI1GVFVVydV8IiJpcIYTWHV1daiurkZiYiJaWlqg0+kAADqdDi0tLQCAxsZGj9VN9Xo9\nHA6HLO0lIpKMihKO7E27cOECFixYgOeffx7Dh3suc6rRaHzedebtu4MdfhoATApIS4mIulcLoO7b\nPx/0Edcnso/SgSPrDKetrQ0LFizAkiVLMH/+fABXZzXNzc0AgKamJowZMwZA9+9qiIyM7Pa4szv8\nZLIhooE2CZ7jTkDxGk7/CSGQkZGBuLg4rFq1yr3fYrGgqKgIAFBUVORORBaLBTt37oTT6URtbS3s\ndjsSEhJkaTsRkWR4Sq3/jhw5gv/5n//B1KlTER8fD+Dqbc+PPfYYFi5ciFdeeQUGgwGvvvoqACAu\nLg4LFy5EXFwcQkJCUFBQwIc8iUj9FJxAeku2rtx2221eb2s+cOBAt/tzc3ORm5s7kM0iIlIWBZ8i\n6y1F3KVGRETqx7XUiIgGQMDWUlvUh3IlXEtNMnkdtv4eg/Hy1cH44IqXog6lxXsr359jdKGiUVpF\nXSEiUiEVXcNhwiEiUjIVjdIq6goRkQqpaJTmXWpEREoWoAc/rVYrYmJiYDKZsGnTpm5jcnJyYDKZ\nYDabUV1d3WPZqqoqJCQkID4+HjNmzMD777/vsytMOEREShaApW1cLheys7NhtVphs9lQUlKCkydP\nesSUl5ejpqYGdrsdW7duRVZWVo9lH330UTz55JOorq7G73//ezz66KM+u8KEQ0SkZAGY4VRVVcFo\nNMJgMECr1SItLQ1lZWUeMR1fDZOYmIjW1lY0Nzf7LDt27FicO3cOANDa2up1fcuOXSEiIqUKwCjt\ncDgwfvx492e9Xo/KysoeYxwOBxobG72W3bhxI2677Tb85je/QXt7O9577z2f7eAMh4hIyQJwSs3f\nB+J7+7BoRkYGNm/ejM8//xzPPfccli9f7jOeMxwiIiXzY5Q+VHd186bz613q6+s9XmjZXUxDQwP0\nej3a2tq8lq2qqnKvfXnfffdhxYoVPtvJGQ4RUZCbZQDyZn23dTZ9+nTY7XbU1dXB6XSitLQUFovF\nI8ZisaC4uBgAUFFRgbCwMOh0Op9ljUYjDh8+DAB45513EB0d7bOdXEuNiGgABGwttd/3odx/da1/\n7969WLVqFVwuFzIyMrBmzRps2bIFAJCZmQkA7rvRQkNDUVhYiFtuucVrWQD44IMP8MADD+Dy5csY\nOnQoCgoK3K+b6bZdakw4eeBaagMdL0UdjA+ueCnqUFq8t/J5CGDCWd+HcrlcvJOIiHqLa6kREZEk\nVDRKq6grREQqpKJRWkVdISJSIZ5SIyIiSaholFZRV4iIVEhFo7SKukJEpEI8pUZERJJQ0Sitoq4Q\nEamQikZpVa40QEQkt4CtNPBaH8r9nCsNSCYPgVuigvHy1cH44IqXog6lxXsr359jqJkqEw4RkWrw\npgEiIpKEikZpFXWFiEiFVDRKq6grREQqxFNqREQkCRWN0irqChGRCqlolFZRV4iIVIin1IiISBIq\nGqVV1BUiIhVS0Sitoq4QEamQikZprqVGRDQAAraW2sd9KDeFa6l1sXz5crz55psYM2YMPv746t/q\nmTNn8Itf/AL/+te/YDAY8OqrryIsLAwAsGHDBmzbtg2DBw/G5s2bcdddd3V73DwEbk0kxstXB+OD\nK16KOpQW7618f44xUKxWK1atWgWXy4UVK1Zg9erVXWJycnKwd+9eDBs2DNu3b0d8fLzPsmlpafjk\nk08AAK2trQgLC0N1dbXXNgwagH757f7774fVavXYt3HjRiQlJeHUqVOYM2cONm7cCACw2WwoLS2F\nzWaD1WrFypUr0d7eLkeziYikE9KHrROXy4Xs7GxYrVbYbDaUlJTg5MmTHjHl5eWoqamB3W7H1q1b\nkZWV1WPZnTt3orq6GtXV1ViwYAEWLFjgsyuyJpzbb78d4eHhHvv27NmD9PR0AEB6ejp2794NACgr\nK8OiRYug1WphMBhgNBpRVVUleZuJiCQVgIRTVVUFo9EIg8EArVaLtLQ0lJWVecR0HHsTExPR2tqK\n5uZmv8oKIfDqq69i0aJFPrsia8LpTktLC3Q6HQBAp9OhpaUFANDY2Ai9Xu+O0+v1cDgcsrSRiEgy\ng/uwdeJwODB+/Hj35+7GT28xjY2NPZZ99913odPpEBUV5bMrir7/QaPR+LwJgDcIEJHqBWCU9nes\n7OuNBiUlJVi8eHGPcYpLODqdDs3NzYiIiEBTUxPGjBkDAIiMjER9fb07rqGhAZGRkd0e42CHnwYA\nkwa0xUT0fVcLoO7bPx/0EdcnfozSh/5xdfOm8/hZX1/vccaou5iGhgbo9Xq0tbX5LHvlyhXs2rUL\nx44d67GdijulZrFYUFRUBAAoKirC/Pnz3ft37twJp9OJ2tpa2O12JCQkdHuM2R1+MtkQ0UCbBM9x\nJ6D8uGYz6w4g77Hvts6mT58Ou92Ouro6OJ1OlJaWwmKxeMRYLBYUFxcDACoqKhAWFgadTtdj2QMH\nDiA2Nhbjxo3zqyuyWbRoEQ4fPozTp09j/Pjx+P3vf4/HHnsMCxcuxCuvvOK+LRoA4uLisHDhQsTF\nxSEkJAQFBQU8pUZEqicCsJZaSEgI8vPzMXfuXLhcLmRkZCA2NhZbtmwBAGRmZiIlJQXl5eUwGo0I\nDQ1FYWGhz7LXlJaW9nizgLsd/e9K35WUlHS7/8CBA93uz83NRW5u7kA2iYhIUVwBGqWTk5ORnJzs\nsS8zM9Pjc35+vt9lr7mWmPyhuGs4RET0nUAlHCVQUVeIiNTnyuC+XGpX5kPxXEuNiGgABGottXNX\nhvS63MgQJ9dSk0peh62/x2C8fHUwPrjipahDafHeyvfnGGqmyoRDRKQWrsHqeeUnEw4RkYK5VPSO\naSYcIiIFu8KEQ0REUnCpaJhWT0+IiFSIp9SIiEgSTDhERCQJJhwiIpIEbxogIiJJqOmmAS5tQ0Q0\nAAK1tM2HIrrX5cyaU1zaRip5CNwSFYyXrw7GB1e8FHUoLd5b+f4cozM1XcNR3Bs/iYhInVQ5wyEi\nUgveNEBERJJQ000D6ukJEZEKqekaDhMOEZGCMeEQEZEkmHCIiEgSarppgLdFExEpmAshvd66Y7Va\nERMTA5PJhE2bNnUbk5OTA5PJBLPZjOrqar/KvvDCC4iNjcXkyZOxevVqn33hDIeISMECcUrN5XIh\nOzsbBw4cQGRkJGbMmAGLxYLY2Fh3THl5OWpqamC321FZWYmsrCxUVFT4LHvw4EHs2bMHH330EbRa\nLb744guf7eAMh4hIwVwY3Outs6qqKhiNRhgMBmi1WqSlpaGsrMwjZs+ePUhPTwcAJCYmorW1Fc3N\nzT7LvvTSS1izZg20Wi0AYPTo0T77osoZTl6nn/09DuPlq4PxwRUvRR1Kiw90+c4CcQ3H4XBg/Pjx\n7s96vR6VlZU9xjgcDjQ2Nnota7fb8be//Q25ubm47rrr8Mwzz2D69Ole26HahHNt6+8xGC9PHXl9\n+i0qL35gj66seCnqUFq8t/L9OUZngXjw099FjXu74OeVK1dw9uxZVFRU4P3338fChQvx2WefeY1X\nZcIhIvo+OXHoNE4c+tLr95GRkaivr3d/rq+vh16v9xnT0NAAvV6PtrY2r2X1ej3uvfdeAMCMGTMw\naNAgfPnll7jxxhu7bQev4RARKZg/12xiZumwIC/OvXU2ffp02O121NXVwel0orS0FBaLxSPGYrGg\nuLgYAFBRUYGwsDDodDqfZefPn4933nkHAHDq1Ck4nU6vyQbgDIeISNECcZdaSEgI8vPzMXfuXLhc\nLmRkZCA2NhZbtmwBAGRmZiIlJQXl5eUwGo0IDQ1FYWGhz7IAsHz5cixfvhxTpkzBkCFD3AnLazv6\n3RMiIhowgXrwMzk5GcnJyR77MjMzPT7n5+f7XRYAtFotduzY4XcbmHCIiBSMq0UTEZEkuJYaERFJ\nggmHiIgkwYRDRESSUNNq0RrR20dLFc7fJ2qJiAZSIIZWjUaDZ8XKXpd7WFMQkPoDTZUznDwEbokK\nxstTB5e2Cb54KepQWry38v05RmdqOqUWdCsN+PNOByIitQjEatFKEVQJ59p7GaxWK2w2G0pKSnDy\n5Em5m0VERH7oMeFs3rwZZ8+elaItPfLnnQ5ERGpyBYN7vSlVjwmnpaUFM2bMwMKFC2G1WmW9EOXt\nfQ1SWX7kCJYfOfK9/V4JbThyZDmOHFketN/L/ff3ffg3Eog+KEmgXjGtBD0mnHXr1uHUqVNYvnw5\ntm/fDpPJhNzcXHz66adStM+Dv3egHezws3bAWkNEdFUtPMedQFLTNRy/UuGgQYMQEREBnU6HwYMH\n4+zZs7jvvvtw55134umnnx7oNrr5804HAJgN4PC3PwNp28yZ3+vvpW1DXrffz5y5zWd5pX8v9+9Q\nXf9GBub7vpj07XZt3DkcwGMrOYH0Vo8J5/nnn0dxcTFuvPFGrFixAs888wy0Wi3a29thMpkkTTgd\n38swbtw4lJaWoqSkRLL6iYik9r1KOGfOnMHrr7+OiRMneuwfNGgQ3njjjQFrWHd8vZeBiEiNlHwT\nQG/1mHCeeOIJr9/FxXV9s9xA8/ZeBiIiNVLyTQC9pZ6eEBGpkJpOqXEtNSKiARCotdRWimd7Xa5A\n8zDXUpNKHgK3JhLj5amDa6kFX7wUdSgt3lv5/hyjs+/VNRwiIpKPmq7hBNVaakREFLzUkzqJiFRI\nTTcNcIZDRKRggVraxp9Xu+Tk5MBkMsFsNqO6urrHsnl5edDr9YiPj0d8fDysVqvPvnCGQ0SkYIGY\n4Vx7tcuBAwcQGRmJGTNmwGKxeDw4X15ejpqaGtjtdlRWViIrKwsVFRU+y2o0Gjz00EN46KGH/GoH\nZzhERAoWiNcT+PNqlz179iA9PR0AkJiYiNbWVjQ3N/dYtje3XzPhEBEpWCBeT+DPq128xTQ2Nvos\n+8ILL8BsNiMjIwOtra0++8KEQ0SkYIG4huPvA/G9fVg0KysLtbW1OH78OMaOHYuHH37YZzyv4RAR\nKZg/13BOHzqBLw+d8Pq9P6926RzT0NAAvV6PtrY2r2XHjBnj3r9ixQqkpqb6bCdnOERECubPNZuw\nWVMRlbfIvXXW8dUuTqcTpaWlsFgsHjEWiwXFxcUAgIqKCoSFhUGn0/ks29TU5C6/a9cuTJkyxWdf\nuJYaEdEACNRaaneJsp4DO9mn+VmX+vfu3YtVq1a5X+2yZs0abNmyBQCQmZkJAMjOzobVakVoaCgK\nCwtxyy23eC0LAEuXLsXx48eh0WgwadIkbNmyBTqdznt/1Jhw8sC11AY6fqDr4FpqwRcvRR1Ki/dW\nPg+BSzhzxF97Xe5tzU+5eCcREfUOVxogIiLqJc5wiIgUjK8nICIiSajp9QTq6QkRkQqp6RoOEw4R\nkYIx4RARkSSYcIiISBK8aYCIiCShppsGVLnSABGR3AK10kC0+LDX5U5pzFxpQCp5CNwSFYyXrw7G\nB1e8FHUoLd5b+f4cozNewyEiIknwGg4REUlCTddw1NMTIiIVUtMpNS7eSUREkuAMh4hIwdQ0w2HC\nISJSMFc7Ew4REUngyhUmHCIikoDrinqGafX0hIhIhVyc4RARkRTUlHC4lhoR0QAI1Fpqg5ov9Lpc\ne8T1XEvtmtdeew15eXn4v//7P7z//vu45ZZb3N9t2LAB27Ztw+DBg7F582bcddddAICjR49i2bJl\nuHTpElJWpkhbAAAVEklEQVRSUvD88897PX4eArcmEuPlq4PxwRUvRR1Ki/dWvj/H6KzdFZhh2mq1\nYtWqVXC5XFixYgVWr17dJSYnJwd79+7FsGHDsH37dsTHx/tV9tlnn8UjjzyC06dP44YbbvDaBlke\n/JwyZQp27dqFO+64w2O/zWZDaWkpbDYbrFYrVq5c6c7SWVlZeOWVV2C322G322G1WuVoOhGRtK4M\n7v3WicvlQnZ2NqxWK2w2G0pKSnDy5EmPmPLyctTU1MBut2Pr1q3Iysryq2x9fT3279+PiRMn9tgV\nWRJOTEwMoqOju+wvKyvDokWLoNVqYTAYYDQaUVlZiaamJpw/fx4JCQkAgKVLl2L37t1SN5uISHoB\nSDhVVVUwGo0wGAzQarVIS0tDWVmZR8yePXuQnp4OAEhMTERrayuam5t7LPvQQw/hqaee8qsrilra\nprGxEXq93v1Zr9fD4XB02R8ZGQmHwyFHE4mIpHVF0/utE4fDgfHjx7s/Xxtb/YlpbGz0WrasrAx6\nvR5Tp071qysDdg0nKSkJzc3NXfavX78eqampA1UtERF14u/NVL250eDixYtYv3499u/f73f5AUs4\nHRvhr8jISNTX17s/NzQ0QK/XIzIyEg0NDR77IyMjvR7nYIefBgCTet0SIiL/1QKo+/bPB33E9ckV\nP2KqDgHvH/L6deextb6+3uOsUXcx18bftra2bst++umnqKurg9lsdsffeuutqKqqwpgxY7pth+yn\n1DpmRIvFgp07d8LpdKK2thZ2ux0JCQmIiIjAiBEjUFlZCSEEduzYgfnz53s95uwOP5lsiGigTYLn\nuBNQV/zYbpkFZOZ9t3Uyffp02O121NXVwel0orS0FBaLxSPGYrGguLgYAFBRUYGwsDDodDqvZSdP\nnoyWlhbU1taitrYWer0ex44d85psAJlui961axdycnJw+vRp3H333YiPj8fevXsRFxeHhQsXIi4u\nDiEhISgoKHBPBQsKCrBs2TJcvHgRKSkpmDdvnhxNJyKSlj8znB6EhIQgPz8fc+fOhcvlQkZGBmJj\nY7FlyxYAQGZmJlJSUlBeXg6j0YjQ0FAUFhb6LNuZP6ftZEk499xzD+65555uv8vNzUVubm6X/bfe\neis+/vjjgW4aEZGytAXmMMnJyUhOTvbYl5mZ6fE5Pz/f77KdffbZZz22gUvbEBEpmUvuBgQOEw4R\nkZIF4JSaUnAtNSKiARCotdTwdh+OM0fDtdSkkofArYnEePnqYHxwxUtRh9LivZXvzzG6UNEMR5UJ\nh4hINZhwiIhIEkw4REQkCRUlHNlXGiAiou8HznCIiJQsQA9+KgETDhGRkvHBTyIikoSKruEw4RAR\nKRkTDhERSUJFCYdL2xARDYCALW3zYh+O8wCXtpFMHgK3RAXj5auD8cEVL0UdSov3Vr4/x+hCRTMc\nVSYcIiLVYMIhIiJJ8DkcIiKSBJ/DISIiSajolBrXUiMiIklwhkNEpGQqmuEw4RARKZmKEg5PqRER\nKVlbH7ZuWK1WxMTEwGQyYdOmTd3G5OTkwGQywWw2o7q6useya9euhdlsxrRp0zBnzhzU19f77AoT\nDhGRkrn6sHU+hMuF7OxsWK1W2Gw2lJSU4OTJkx4x5eXlqKmpgd1ux9atW5GVldVj2UcffRQffvgh\njh8/jvnz5+OJJ57w2RUmHCIiJbvSh62TqqoqGI1GGAwGaLVapKWloayszCNmz549SE9PBwAkJiai\ntbUVzc3NPssOHz7cXf7ChQsYNWqUz66o8hpOXqef/T0O4+Wrg/HBFS9FHUqLD3T5LgJwDcfhcGD8\n+PHuz3q9HpWVlT3GOBwONDY2+iz729/+Fjt27MCwYcNQUVHhsx2qTTjXtv4eg/Hy1cH44IqXog6l\nxXsr359jdBGAlQb8XdS4Lwt+rlu3DuvWrcPGjRvx61//GoWFhV5jVZlwiIhUw5+VBpoOAc2HvH4d\nGRnpcUG/vr4eer3eZ0xDQwP0ej3a2tp6LAsAixcvRkpKis9m8hoOEZGS+XPNZvQsYEred1sn06dP\nh91uR11dHZxOJ0pLS2GxWDxiLBYLiouLAQAVFRUICwuDTqfzWdZut7vLl5WVIT4+3mdXOMMhIlKy\nAFzDCQkJQX5+PubOnQuXy4WMjAzExsZiy5YtAIDMzEykpKSgvLwcRqMRoaGh7lNj3soCwJo1a/DJ\nJ59g8ODBiIqKwksvveS7Hf3vChERDZgArRadnJyM5ORkj32ZmZken/Pz8/0uCwB/+ctfetUGnlIj\nIiJJcIZDRKRkfD0BERFJQkVrqTHhEBEpGRMOERFJQkWvmNaIvjxaqmD+PlFLRDSQAjG0ajQaYHYf\njnNQE5D6A02VM5w8BG6JCsbLVwfjgyteijqUFu+tfH+O0YWKTqnJclv0I488gtjYWJjNZtx77704\nd+6c+7sNGzbAZDIhJiYG+/btc+8/evQopkyZApPJhF/96ldyNJuISHoBWC1aKWRJOHfddRdOnDiB\nDz/8ENHR0diwYQMAwGazobS0FDabDVarFStXrnRPC7OysvDKK6/AbrfDbrfDarXK0XQiImkF6AVs\nSiBLwklKSsKgQVerTkxMRENDA4Cra/EsWrQIWq0WBoMBRqMRlZWVaGpqwvnz55GQkAAAWLp0KXbv\n3i1H04mIpBWAF7AphewrDWzbts29wmhjY6PHKqQd38fQcX9kZCQcDofkbSUikpyKTqkN2E0DSUlJ\naG5u7rJ//fr1SE1NBXD1PQpDhgzB4sWLA1r3wQ4/DQAmBfToRESeagHUffvngz7i+kTBCaS3Bizh\n7N+/3+f327dvR3l5Od5++233Pm/vY4iMjHSfdru2PzIy0uuxZwM4/O1PIqKBNunb7dq4c1je5iiW\nLKfUrFYrnn76aZSVleG6665z77dYLNi5cyecTidqa2tht9uRkJCAiIgIjBgxApWVlRBCYMeOHZg/\nf74cTScikpaKbhqQ5TmcBx98EE6nE0lJSQCAH//4xygoKEBcXBwWLlyIuLg4hISEoKCgwP0gZ0FB\nAZYtW4aLFy8iJSUF8+bNk6PpRETSUvBNAL0lS8Lp+Ja4znJzc5Gbm9tl/6233oqPP/54IJtFRKQ8\nvIZDRESSUFHC4VpqREQDIGBrqY3qw3FOcy01yeQhcGsiMV6+OhgfXPFS1KG0eG/l+3OMLngNh4iI\nJKGiU2pMOERESsaEQ0REklDwczW9JftaakRE5EOAFu+0Wq2IiYmByWTCpk2buo3JycmByWSC2WxG\ndXV1j2V9vWqmO0w4RERKJvqwdeJyuZCdnQ2r1QqbzYaSkhKcPHnSI6a8vBw1NTWw2+3YunUrsrKy\neizr7VUz3jDhEBGpXFVVFYxGIwwGA7RaLdLS0lBWVuYRs2fPHqSnpwO4+tqY1tZWNDc3+yzr7VUz\n3jDhEBGpnMPhwPjx492fr736xZ+YxsbGHssCnq+a8YYJh4hI5fx9IL6vD4v6+6oZ3qVGRBT0Dn27\nda/zq1/q6+s9XmrZXcy118O0tbX5LNvdq2a84QyHiEjR/HkfwUwAv+2weZo+fTrsdjvq6urgdDpR\nWloKi8XiEWOxWFBcXAwAqKioQFhYGHQ6nc+y3l4145VQGfTtng5u3LhxC+gWuPHsmz5sXesvLy8X\n0dHRIioqSqxfv14IIcTLL78sXn75ZXfMAw88IKKiosTUqVPF0aNHfZYVQgij0SgmTJggpk2bJqZN\nmyaysrJ89keVi3fmgWupDXS8FHUwPrjipahDafHeyucBgVu8E76fbeneSC7eSUREvaWetW2YcIiI\nFE09a9sw4RARKRoTDhERSYKn1IiISBLqmeHwORwiIpIEZzhERIrGU2pERCQJ9ZxSY8IhIlI09cxw\nVLnSABGR3AK30kBFH0r+iCsNSCUPgVuigvHy1cH44IqXog6lxXsr359jdKWeGY4qEw4RkXrwGg4R\nEUmCMxwiIpIEZzhERCQJznCIiEgSnOEQEZEkmHCIiEgS6jmlxsU7iYhIEpzhEBEpGk+pERGRJNRz\nSo0Jh4hI0dQzw5HlGs7atWthNpsxbdo0zJkzB/X19e7vNmzYAJPJhJiYGOzbt8+9/+jRo5gyZQpM\nJhN+9atfydFsIiIZXOnD1pXVakVMTAxMJhM2bdrUbUxOTg5MJhPMZjOqq6t7LPvaa6/h5ptvxuDB\ng3Hs2LEeeyJLwnn00Ufx4Ycf4vjx45g/fz6eeOIJAIDNZkNpaSlsNhusVitWrlzpXvE0KysLr7zy\nCux2O+x2O6xWqxxNH3C1cjcgAIK9D2y//NTQh8Bp68PmyeVyITs7G1arFTabDSUlJTh58qRHTHl5\nOWpqamC327F161ZkZWX1WHbKlCnYtWsX7rjjDr96IsspteHDh7v/fOHCBYwaNQoAUFZWhkWLFkGr\n1cJgMMBoNKKyshITJ07E+fPnkZCQAABYunQpdu/ejXnz5nV7/LxOP/uqt+W/b/FS1MH44IqXog6l\nxQe6fFf9v4ZTVVUFo9EIg8EAAEhLS0NZWRliY2PdMXv27EF6ejoAIDExEa2trWhubkZtba3XsjEx\nMb1qh2zXcH77299ix44dGDp0KKqqqgAAjY2N+NGPfuSO0ev1cDgc0Gq10Ov17v2RkZFwOBxej53X\nYeur3pYPVPxBALMV1J6+lAmWPniLZ/v7Fh/IOoLld+CtfH+O0VX/r+E4HA6MHz/e/Vmv16OysrLH\nGIfDgcbGxh7L+mvATqklJSVhypQpXbY33ngDALBu3Tp8/vnnuP/++7Fq1aqBagYRUZDr/zUcf19M\nOdAvbRuwGc7+/fv9ilu8eDFSUlIAXJ25dLyBoKGhAXq9HpGRkWhoaPDYHxkZ2e3xoqKikPfppwCC\n95Ta4QE+fqDifZUJlj54i2f7+xYfyDqC5XfQXfmoqKh+HqXzEXvn+uuv9/jceWytr6/3OGvUXcy1\n8betra3Hsn4TMjh16pT7z5s3bxa//OUvhRBCnDhxQpjNZnH58mXx2WefiR/+8Ieivb1dCCFEQkKC\nqKioEO3t7SI5OVns3btXjqYTEQWdtrY28cMf/lDU1taKy5cvC7PZLGw2m0fMm2++KZKTk4UQQrz3\n3nsiMTHR77KzZs0SH3zwQY/tkOUazpo1a/DJJ59g8ODBiIqKwksvvQQAiIuLw8KFCxEXF4eQkBAU\nFBS4p4IFBQVYtmwZLl68iJSUFK83DBARkaeQkBDk5+dj7ty5cLlcyMjIQGxsLLZs2QIAyMzMREpK\nCsrLy2E0GhEaGorCwkKfZQFg165dyMnJwenTp3H33XcjPj4ee/fu9doOjRADfNKOiIgIQbx4Z7A/\nPPrII48gNjYWZrMZ9957L86dO+f+LhjaD/h+6CtY+tCRPw/GKcHy5cuh0+kwZcoU974zZ84gKSkJ\n0dHRuOuuu9Da2ur+ztvvQi719fWYPXs2br75ZkyePBmbN28GEDx9uHTpEhITEzFt2jTExcVhzZo1\nAIKn/bIKyAlCGXz11VfuP2/evFlkZGQIIb67DuR0OkVtba2IiopyXweaMWOGqKysFEII2a8D7du3\nT7hcLiGEEKtXrxarV68WQgRP+4UQ4uTJk+KTTz4Rs2bNEkePHnXvD6Y+XHPlyhURFRUlamtrhdPp\n7PY8tVL87W9/E8eOHROTJ09273vkkUfEpk2bhBBCbNy40ee/p2v/7uTS1NQkqqurhRBCnD9/XkRH\nRwubzRZUffj666+FEFevbyQmJop33303qNovl6Cd4fT24dGmpqZuHx6VS1JSEgYNuvrXn5iY6L4L\nL1jaDwAxMTGIjo7usj+Y+nBNxwfjtFqt++E2Jbr99tsRHh7usa/jQ3vp6enuv9fufhfXnnuTS0RE\nBKZNmwbg6t1UsbGxcDgcQdWHYcOGAQCcTidcLhfCw8ODqv1yCdqEA1x9eHTChAnYvn27e1rb2Njo\ncctex4eXevPwqJS2bdvmvjU8GNvfWTD2wdtDb8GipaUFOp0OAKDT6dDS0gLA++9CKerq6lBdXY3E\nxMSg6kN7ezumTZsGnU7nPj0YTO2Xi6JXi05KSkJzc3OX/evXr0dqairWrVuHdevWYePGjVi1apX7\nrgql6Kn9wNUHYIcMGYLFixdL3Ty/+NMHNfD3wbhgoNFofPZHKX29cOECFixYgOeff97jjAWg/D4M\nGjQIx48fx7lz5zB37lwcPHjQ43ult18uik44cj08Gig9tX/79u0oLy/H22+/7d6npPYD/v8OOlJa\nH/zhz4NxSqbT6dDc3IyIiAg0NTVhzJgxALr/XSjh77ytrQ0LFizAkiVLMH/+fADB1wcAGDlyJO6+\n+24cPXo0KNsvtaA9pWa3291/LisrQ3x8PADAYrFg586dcDqdqK2thd1uR0JCAiIiIjBixAhUVlZC\nCIEdO3a4/6HLwWq14umnn0ZZWRmuu+469/5gaX9nosPd9cHYh+nTp8Nut6Ourg5OpxOlpaWwWCxy\nN8tvFosFRUVFAICioiL336u334WchBDIyMhAXFycx7JWwdKH06dPu+9Au3jxIvbv34/4+Pigab+s\nZL1loR8WLFggJk+eLMxms7j33ntFS0uL+7t169aJqKgocdNNNwmr1ere/8EHH4jJkyeLqKgo8eCD\nD8rRbDej0SgmTJggpk2bJqZNmyaysrLc3wVD+4UQ4vXXXxd6vV5cd911QqfTiXnz5rm/C5Y+dFRe\nXi6io6NFVFSUWL9+vdzN8SotLU2MHTtWaLVaodfrxbZt28SXX34p5syZI0wmk0hKShJnz551x3v7\nXcjl3XffFRqNRpjNZve//7179wZNHz766CMRHx8vzGazmDJlinjqqaeEECJo2i8nPvhJRESSCNpT\nakREFFyYcIiISBJMOEREJAkmHCIikgQTDhERSYIJh4iIJMGEQ0REkmDCISIiSTDhEAF4//33YTab\ncfnyZXz99deYPHkybDab3M0iUhWuNED0rbVr1+LSpUu4ePEixo8fj9WrV8vdJCJVYcIh+lZbWxum\nT5+OoUOH4r333vveLiFPNFB4So3oW6dPn8bXX3+NCxcu4OLFi3I3h0h1OMMh+pbFYsHixYvx2Wef\noampCS+88ILcTSJSFUW/gI1IKsXFxfjBD36AtLQ0tLe34yc/+QkOHTqEWbNmyd00ItXgDIeIiCTB\nazhERCQJJhwiIpIEEw4REUmCCYeIiCTBhENERJJgwiEiIkkw4RARkSSYcIiISBL/Dy1Jcn1wOZYq\nAAAAAElFTkSuQmCC\n", + "png": "iVBORw0KGgoAAAANSUhEUgAAAYMAAAEZCAYAAAB1mUk3AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztvXt8VNW5//9euRGEAIIGJANESUKCiYECATz4hRbDJQoq\ntrwIFkWCVjhSPVBA6U+LpdyknN9BOdHjTy6W/kqRogYFIggFtQcSQKBgqEQJkBugXCL3kGR9/5id\nnT3JzGTCkEwyPO/Xa79ea6+19l5r75nsJ7M/63kepbXWCIIgCLc0Ab6egCAIguB7xBgIgiAIYgwE\nQRAEMQaCIAgCYgwEQRAExBgIgiAIiDEQmhArV67kgQceMPfDwsI4duyY7yYkCH6EGAOhUfHll19y\n//3306ZNG9q1a8eAAQPYs2eP074XLlwgMjLypo1dWlpKWloakZGRtGrVip49e5KZmWm279q1i+Tk\nZNq1a0d4eDijR4/m5MmTN218QfAlYgyERsOPP/7Iww8/zAsvvMC5c+coLCzkd7/7Hc2aNWuQ8cvK\nyujcuTOff/45P/74I3/4wx8YPXo0x48fB+D8+fM899xzHD9+nOPHjxMWFsbTTz/dIHMThHpHC0Ij\nYffu3bpNmzYu21esWKEHDBhg7iul9Hfffae11vry5ct66tSpukuXLrp169Z6wIAB+sqVK1prrXfu\n3Kn79++v27RpoxMTE/X27ds9ntN9992nP/jgA6dte/fu1WFhYR6fSxAaM/LLQGg0dOvWjcDAQMaP\nH09mZibnzp3z+Njf/OY37Nu3j507d3L27FkWLVpEQEAAhYWFPPzww7z66qucO3eOP/7xjzz++OP8\n8MMPtZ7z1KlTHDlyhHvvvddp++eff058fLzHcxSExowYA6HREBYWxpdffolSimeeeYbw8HAeeeQR\nTp8+7fa4iooKVqxYwZIlS7jrrrsICAigX79+hISE8Oc//5mUlBSGDRsGwIMPPkjv3r3ZuHGj23Ne\nv36dJ554gvHjxxMTE1Oj/Z///Cdz5sxh0aJFN37BgtCIEGMgNCpiY2NZsWIF+fn5HDp0iKKiIl58\n8UW3x/zwww9cvXqVrl271mg7fvw4a9eu5fbbbze3f/zjH26F34qKCsaNG0doaChLly6t0f7tt9+S\nkpLCG2+8wb/927/V/SIFoREixkBotHTr1o2nnnqKQ4cOue13xx13EBoayrffflujrXPnzowbN45z\n586Z24ULF5gxY4bTc2mtSUtL4/vvv2fdunUEBgY6tB8/fpzk5GReffVVnnjiiRu/OEFoZIgxEBoN\n33zzDf/5n/9JYWEhAPn5+axevZr+/fu7PS4gIIAJEyYwdepUiouLKS8vZ+fOnZSWlvLLX/6Sjz/+\nmM2bN1NeXs7Vq1fZvn27OUZ1Jk2axL/+9S/Wr19fYxVTYWEhP/vZz3j++ed59tlnb85FC0IjQYyB\n0GgICwsjKyuLvn370rJlS/r37899993H4sWLAVBKoZQy+1vLf/zjH0lISKBPnz60a9eOl19+mYqK\nCmw2GxkZGcybN4/w8HA6d+7M4sWLqaioqDH+8ePHeeeddzhw4AAdOnQgLCyMsLAwVq9eDcC7775L\nXl4es2fPNttatWpVz3dFEBoGpXXTSm6TmZnJiy++SHl5ORMnTmTmzJm+npIgCEKTp0kZg/Lycrp1\n68Znn31GREQEffr0YfXq1cTFxfl6aoIgCE2aJvWaKDs7m6ioKCIjIwkODmbMmDFkZGT4elqCIAg3\nnbVr13LvvfcSGBjIV1995bLfkiVLSEhIID4+niVLlji0vfnmm8TFxREfH1/rW5SgmzLrBqKwsJBO\nnTqZ+zabjaysLB/OSBAEoX5ISEjgww8/5Fe/+pXLPocOHeLdd99l9+7dBAcHM2zYMB5++GG6du3K\n3//+d9avX88///lPgoOD+f77792O16R+GVgFQ0EQBH8mNjbWqcOjlX/961/07duX0NBQAgMDGThw\nIB988AEAb731Fi+//DLBwcEA3HnnnW7P1aSMQUREBPn5+eZ+fn4+NpvNoU/lihPZZJNNttq2Hj16\nePVMal7H8dq2bevVeNWJj4/niy++4OzZs1y+fJkNGzZQUFAAQG5uLp9//jn9+vVj0KBBLqP/VtKk\nXhP17t2b3Nxcjh07RseOHVmzZo257M+K9u7zNVH7QSfc4LEH4XfhMLu9B50DjbHczFvtB93PRdsu\n0IOr1W0FPcqy/wHoVGMnHNQS0C8APUA9DXoFUAxqFuh5QBmoV+F3g2D2GFDPgX4buARqGujFxnmn\ngR5plNfD5Rb28m2X4O/GcD8F9FCjz6egxxkNRcY8BwPNQG0EnQIkgFoIeibQG9QvQK8Fcqvm99Ys\nmAykG6eqXp5tlGcb29+BHdX6nDG+/e3KoNIfuf0o416NAlqAWmXMtwzUauMe9gP1AuglwAlQi0FP\nM65vcbVrHWz5PHob5T2g7zHKR0EbTszqH6CfMybSG9RE0O/CD2ktuVNdZPrvgnk97TqqM+gTwNug\n5oGeBWwzvgf9gNPGea1jdDS/Cqgi0O1wQJ2p+uyqc9slKAl13lZJ66u197H2dTVWbdx2CUpb39ix\nzgg5cMCr468Cf3DTfhTIs+z/vVq8reTkZKfe8PPmzWPEiBG1jh8bG8vMmTMZMmQILVq0oGfPnqaj\nZFlZGefOnWPXrl3s3r2b0aNHc/ToUZfnalLGICgoiKVLlzJ06FDKy8tJS0uTlUSCIPiUYDdt3Yyt\nkr9Xa9+yZYvX40+YMIEJEyYAMGvWLDp37gzYNdVRo+z/Efbp04eAgADOnDlDu3btnJ6nSRkDgOHD\nhzN8+HBfT0MQBAFomIeoOw+A06dPEx4ezokTJ/jwww/NRTWPPvoo27ZtY+DAgRw5coTS0lKXhgCa\nmGbQ1Bh0gz+FGxOD7vb1DG4Okb6ewE3g3wYF1t5JaHCa12GrCx9++CGdOnVi165dPPTQQ+Y/wUVF\nRTz00ENmv5///Ofce++9jBw5kvT0dNMrfsKECRw9epSEhARSU1P505/+5Ha8JuV05glKyYojQRA8\nx5tHoFKK/68O/Z/xcrz6pMm9JvKExiIge3xsfQnIIy3760H/wthpAWol6PFAO4sAWgLqXdATgbMW\nIbWZRTxtZjkWo2wIwqbYWlm2iqpG+gCVAnq/ZV49jP0WoKJB54JuCwHtoOIM5LcNp4s6zXEdTikh\nRKsCcrWNqOICVEfQRcZ5OoI2tDF1D+gvjfKAqvFUD8c+Z7X9f7W26golZSEAnAlsxz2qmKP6Ltpw\nnrbqCmd1c9qcvWLO6XKLAFqGVnDxagC3Xaow68E+b51rjBENeqtRHmwI8BhivPXeTDTK74L+vXFj\nehv3aiNwN6g40IeBw6BGgf4A2AzqbUN0PmgI0P9mfHaHQRtymjoMurPlnp8AHY4D6rRrYTakxI8F\n5BLvz+EvD1F/uQ5BEASf4E5AbkqIMRAEQfACf3mI+st1CIIg+AR/+WUgArIgCLc03grI6+vQf6SX\n49UnfvnLoNEIyJ7OI8jwSu3t5nx7qjxVa7T9w0MBudIjORDUWkNQDrUIv1ct9deMY0YCrS19Olq8\ngzHKRppg9byjcGsVaDP0EAAeUZtZrCcDcIZ2zFNzmKVf4TxtSFfTmKwX8wPteF+NZ7ReST6d2al+\nRn+9jSLu4riKo4s+TCnNKFb3cJehCBere2hXbnfDPxNoI1wfB+C06kLra8UAlDS7i5aX7MG6Lra4\nk+YlZwG40rqteWzzgCsUqGhsOpeOFJGtBpKkd9CeU3ysRjNCv08YF/iLSmOsXkYbzpvzBkhX05hj\nqMOvqMWs1o8CkKo+4pC252iOV985Ct87jPLAKqGdZobovBW79/NQ0J9i98J+3rjnhywCcp7h8TwU\nOOP4fVJ7QEdbvgu5jh7JYPdKdicg1yb43nbJc1G4Ln2dHdvYBOS6LhltrPilMRAEQWgo/OUh6i/X\nIQiC4BP8RTMQYyAIguAF/vIQFQFZEIRbGm8F5L116N/Ly/HqE38xag7cVAH5Bs9Vp2NvhoA8tFrd\np248kINceBQHWTyQQ0EtBf089pDXrxresXGWkNIYZYtX7wFtT8aRqI6wVKcB8Lxaxmi9EoD31Xhi\ntD1s8LenulLRoSUBJy9ScawF9FOwS8N3wBMK/n8N3wK/U/Cahn8BqxWkaigAvlDwgPGH9YWCOKN8\nWMHdRjlPQSejnF+t3tq/8jw2yxhRwBwFrxjlpxS8Z5T/TcE/NEReg4hQKLxqPz4ilK76EADfqXhG\n6PcB+FiNZrahus9WC/nEUPwfVls5brgDd1Gn0aeNz6jI4p1dZPFGPmoRkHMtocjzLR7jZ0DtAD3Q\n+Fx2OC6EUAcdPZLB8Ep2EcPMXXjrSm5lAdlfHqL+ch2CIAg+QTQDQRAEQZaWCoIgCP7zy0AEZEEQ\nbmm8FZBrJq10TQcvx6tP/PKXQZMTkEMt+WtdnW9XlSBYo22HkTvYWrfRIhhjeBZX5h5uZhGKW1hE\nyBaWnLrtLXl+7wY1EvR6IN4IEW0RjQ/rLgDEqeO8YwzyrFrFOP0OAKvUs3TRhwE4ruLsIjHAXuDf\nFfy3hkPAWwomGULx3xX8VMNB4AcFd2j4AUABGsgFYoAjxkXF4Jh1uTJIwEhgnVF+vFq94S7NAKoy\n1UZWjdEJu+jcSdtzF36m4EENCcD/q+A/NMQDaQqWGdeUpuCAUU5U2Ix41gUqmsf1nwFYp37JfP0i\nAC+r/2KTHgTAcLXdFJM7F5+uCtGda3gn7zDKRm5kjlo+r3yLl/gZqvJJY5Qt3y21y9EjGZx7JZtt\nbryTAcrKbm0BObguT9Ey78erL/zSGAiCIDQUQX5iDCTtpSAIghcEB3q+1YW1a9dy7733EhgYyFdf\nfeWy35IlS0hISCA+Pp4lS5aY9dnZ2SQlJdGzZ0/69OnD7t273Y4nxkAQBMELgoI83+pCQkICH374\nIf/n//wfl30OHTrEu+++y+7duzlw4ACffPIJ3333HQAzZsxgzpw57Nu3j9///vfMmDHD/YDaz8D+\nQlk22WSTzaPN2+eNDvd8u5HxBg0apPfu3eu0be3atTotLc3cnzNnjn799de11lqPGTNGr1mzRmut\n9V/+8hf9xBNPuB3HLzWDRiMgu/EodiDIQwF5sIu2auGqwfA4TrXsr8bMtUtrS97j1hbv4rYW79aO\nlly7cVU5eE/GtuYuVUKxtqt4d6kSBwH0FT0LgDlqHoP0JgC2q+GEnLErdaXtWsM6bZ/Hl1QJsftx\nIRpfAW4DLgOHsTv078UuII8B/mpc1BhgsVGeBsw2yrOrlf9glP8fIN0oTwYqY0dHUSVMR2OKyVHA\ntwqitF1A/lDBYxp6A79VMNe4pt8q+NQoD1UOIbOt92OangPAYvUKf9aPA/BLtc704O52+QihLeDq\nJWh2GFRv0Hvsl61SQa/GHrZ6Fuh5wAlLOOvTFm9kjPJATKp7JINzr2SzzY13MsD1MrsQ66mwW5e+\nN/NYV+fzGh8+RePj4/ntb3/L2bNnCQ0NZcOGDSQlJQGwYMECBgwYwG9+8xsqKirYuXOn23P5pTEQ\nBEFoMNw8RbdftW+uSE5O5uTJmotT582bx4gRI2odOjY2lpkzZzJkyBBatGhBz549CQy0ixNpaWm8\n8cYbPPbYY6xdu5YJEyawZcuWG7kMQRAEoVbcPEUHtbRvlbz2o2O7u4ezp0yYMIEJEyYAMGvWLDp3\ntv/Ey87O5rPPPgPg5z//ORMnTnR5DhABWRAEwTsC67DdINqNo9rp0/bohidOnODDDz9k7NixAERF\nRbFjhz2N3rZt24iJial1EL+CRiBIySabbE1n8/Z5o+M83+oy3gcffKBtNpsODQ3V7du318OGDdNa\na11YWKhTUlLMfg888IDu3r27TkxM1Nu2bTPrd+/erZOSknRiYqLu16+f/uqrr9yO55fhKJqcgNzM\nCEPtIkQ1OA9TbbZ9SlV+48q6D7CHpa7cX4ndyxjsovHvQb8KtLN4GodbxMm7QfUDbXirBrSDijOQ\n07Yr8eo7h3y+VgH0RT0fgP9SL9NfbwNgp/oZAScvAlDRoSX81fjKZeFcQN4PnFNwu4ZzrgTkY9g9\niq3exVZBeKFRnomjaGwVmVcY5aeByp/rUfaLJw8Hb+Q47KGu4zTcC/xNwc+1XUB+ScEC45peUrDV\nKA92FJAH6M0AfKmGONynZdr+n1ya+gvZOh6AbuVHaB1USklZCK0Ol6ISQB/ELiBXCvt5oKaBXow9\nzHXlooDTFm9kjLLlu6M+rfndVHtA34NT1FEwHKOdcv1a0xaQvXkE1vV5o/Z7N159IpqBIAiCN/jJ\nU9RPLkMQBMFH+MlT1E8uQxAEwUd4IQw3JvxSMxAEQfAUrzWDgXXov0M0AwfWrl3L7Nmz+de//sXu\n3bv5yU9+YrbNnz+f5cuXExgYyBtvvMGQIUMA2Lt3L+PHj+fq1aukpKQ4BGSqjt8KyCku2qqFqwYj\nZPVzlv23wUjBaxeQKz1Xwy0hkTtacu12whQtf4wLMcXMbwJjSFKHTKEzSR1yEECn6NcBeFPNcBBM\nrUKqKbDuosp7dx9Vouw3wEEFCdquGZcpCNJQ9qN98pRgF5ATgQPGRSXiOoS1tVzpaZxiTACgH5Bv\nlG2YovFdQLGCuwzR2FkI6x5U5UbGKLsIYT1c28XuTepxZurZACxUs1mpRwMwXr1fu4CcZwknnmcR\n/4tALTQ+41OWvNYY5er5sKt5u6tdoONwijqMy/DWAFyz50l256XscL4zjUtA9ho/eb/iEz8DVwGY\ncnJyWLNmDTk5OWRmZjJ58mTTik6aNIlly5aRm5tLbm4umZmZvpi6IAiCI0F12BoxPjEGsbGxTh0g\nMjIySE1NJTg4mMjISKKiosjKyqK4uJgLFy6YMTeefPJJPvroo4aetiAIQk2a1WFrxDQqW1VUVES/\nflW/X202G4WFhQQHB2Oz2cz6iIgICgsLfTFFQRAERxrVU/TGqTcB2ZMATD/96U9ZvHixqRlMmTKF\nfv368cQTTwAwceJEhg8fTmRkJC+99JIZx+OLL77g9ddf5+OPP655QSIgC4JQB7wWkH9Rez+z/9pb\nUEC+kQBMERER5Ofnm/sFBQXYbDYiIiIoKChwqI+IiHB5nkYjILsJSe1AM2OVgZtVCc7yHJttG6u8\nTc26VaCft+wvNTyOwS4gV3quhoMaB3oVdgF5MOitQGdQ0aBz4VxUc9qqK5zVzfmWrjUEZKsAahVG\nrYKpQw7kb4w/hn3AGGX3SP4amKPgFUM0rhSTvwP2Keip7Rqv03zIleVSoxwCnDHK7aqVLxvl2xyP\nbW6UI6gKVR0NbFIwXNsdk99UMEVDLFW5m+OBgQp2GMcPVDQvOQvAldZt6aXteZb3qgGM1csA+ItK\nY46eBsArajHva/s/R6PVx2YI667XvqVlaAUXrwbQIrfCuYB8whJyvMiSD/msJZw1RtnioV49pDU4\nD2tttrkJbw3AVVCn3XspO5zvtOdic41j6yBUe3o+r/GTXwY+D1RntZIjR47kr3/9K6WlpeTl5ZGb\nm0tSUhIdOnSgVatWZGVlobVm1apVPProoz6ctSAIgoEIyDfOhx9+SKdOndi1axcPPfQQw4cPB6B7\n9+6MHj2a7t27M3z4cNLT083XPunp6UycOJHo6GiioqIYNmyYL6YuCILgSANELW0IfGKrHnvsMR57\n7DGnbbNmzWLWrFk16nv16sXBgwfre2qCIAh1o5H/x+8p4oEsCMItjdcC8vO19zP7L70FBWRfclMF\nZE+9iKsfu6cOAnILI4+xixzH4DzPsdm23jFcNVQLWQ2oJdjzHINjruN2lrDVHUENBL0DuwfyPaCP\nwrm7qwTkY0TyE3WYrwx31Z+oww4CqFUYHaffAWCVepYkbU+yka0G0vpaMQAl/+oAiYbH7mGqxOTv\nqPJMzgPeVTBRQyFVgu73wB4FvY0/rD2G6At2Afhuo5ynoJNRzq/Wp2elkG14FoPd63iVgnHa7ow8\nX8HLhoCcpmCZUTZE44Bul6jo0NIhRHeMtntFH1GJDiK6NWz1Up0GwPNqmUMO6VxtX0Id+WMBwa3h\negkE5YHqAXo/dtHYmYB82hKW/Kzx0DEeUmqpo4e6Wlvzu6a2uv67Uftdh7cG4OItLiA38tc/nuKX\nxkAQBKHB8JOnqJ9chiAIgo/wk6eoz5eWCoIgNGnqaTXR9OnTiYuLIzExkVGjRlFS4jyqXmZmJrGx\nsURHR7Nw4UKz/uzZsyQnJxMTE8OQIUM4f/682/FEQBYE4ZbGawH51dr7mf1/7/l4W7ZsYfDgwQQE\nBPDSSy8BsGDBAoc+5eXldOvWjc8++4yIiAj69OnD6tWriYuLY8aMGdxxxx3MmDGDhQsXcu7cuRrH\nW/GTHziO+K2APMpF2wcuBORplv3FLgTk9qB+AXotdm/kSgG5Y5UH8o93V4Wwzg/sVCMH8ifGxB9W\nWx2EUatg+qheDcBHKpUeeicAeaWRlDS7i9bXiin5tgPcq+BrbY9O/ZCCDRoKgF8p+B/D67hSWP6B\nqjDSGOWJRvldQwCGKjG4smztM8Uov6ngNaN8B1XexTbgEQUZRrmXgr2a5lHnuNK6Lc1LztKpVT5H\nVKKDaDxYfwLAVvUwaXopAMvU88zXLwLwsvovVmu702Sq+oidxhe2v9pPgaGORnx/BhUO+jRwtCof\nNUWWHMhFdRCQUy3fhdU182k7y4tstu2x58F2ySVQRbWEubaer8jPBOR6eoomJyeb5b59+7Ju3boa\nfbKzs4mKiiIyMhKAMWPGkJGRQVxcHOvXr2fHDvvCjaeeeopBgwa5NQbymkgQBMEbGiBq6fLly0lJ\nqRmPprCwkE6dOpn7lcE9AU6dOkX79u0BaN++PadOnXI7hl/+MhAEQWgw3DxFt38L279z3e5JQM+5\nc+cSEhLC2LFja/Sr/lpca+30VblSqtZX6GIMBEEQvMHNU3RQrH2r5LXNju21BfRcuXIlGzduZOvW\nrU7bnQX3rAzi2b59e06ePEmHDh0oLi4mPNy9I4i8JhIEQfCGelpNlJmZyaJFi8jIyCA0NNRpn969\ne5Obm8uxY8coLS1lzZo1jBxp904dOXIk7733HgDvvfdercE9ZTWRIAi3NF6vJnq7Dv2f83y86Oho\nSktLadu2LQD9+/cnPT2doqIinnnmGTZs2ADApk2bePHFFykvLyctLY2XX34ZsC8tHT16NCdOnCAy\nMpL333+fNm3auJ6bPxqDRrOayE2CewdaGqs5hrruoj6tZTXRxGp17zpZTTTP2GnlIhxFuCWfQUdQ\ncaAPQ1lHzNAI+a3u4h5VzFF9FwD3qGK+1L0AGKD2OqySsa6emawXA5CupjFCvw9AER3ZqwbQS39J\nER0pVvdwlz7KqVPhZoiHilMtqkJW/AAMVrBVw3ngcQXrjK/v49WS0i8zymnVytY+1mM/tawmMlYN\n0eEaRIRC4VXCO57itOpCuD5OOKc4pJKI19l0pJjN6hGG6AwANqtHHFYQzdYzAZitFrJMjzWm8Rc2\n6wEADFFfOqzMKikLsX9EJ0rNcCDkW1Z5nbas/ioC9QLoJcCZJraayMPQFTWO9SKUhdPznbkJxuDd\nOvSfKLGJBEEQ/BOJTSQIgiDg/HV+k0OMgSAIgjf4yVPULzUDQRAET/FaM/igDv1HiWbQoPitgPwL\nF21rb5KA3BbUUNCfYg9TURlDvyNmaIQf7mzJneoi3+uWANypLppJ3BPVEYfY/FbB1JrnYIp+HbAL\nyOvUL3lc/5kztGO7Gs4gvYlTtOew+glx+ivO08YUli9cDuNiiztpeel7Lp4PM8VdwF7+rsJe7hoA\n3xh/cN2UY9nax3Jsy0vfAxB22wVzvDacN+fRjh/4Ug1hgN7MHZzhI5XKo3o1HSkiXU1zEMitovE7\nehwAz6pVDmE7rPkgKkNQ2NQZKozwCKoIVALog9iF4srP5TSocaBXGeVpoBdjF5DngZ4F/OiBgFzN\nmVVtdB0+Re0CY7rOuQTqBOjObvpYz3eicQnIXuMnT1E/uQxBEAQf4SdPUT+5DEEQBB8hq4kEQRAE\nf1lNJAKyIAi3NF4LyF/Wof8AEZAbFBGQXQjIlfkMWls8V9taBMnWlmTr4ZYY+uGY3rDX2kNoC7h6\nyX6q0BZw3FADu6jTZOt4AJLUIQfB1Comv66nAHCKcBarV5im53CGO1ipJjFev8V52pgC7QXC2Koe\nZrD+hAuEka0GkqR3cJnmphcwwCGV5JBToKs+BMB3Kt6hbO1jPTZJ2+O+h3HBHC+MC+Y82nDenF87\nfjDnfQdneFn9l4O3tfVarffAem+s9+ziVXuIsJahFXaPY7CLw5X3/7TlczlteLG+i100ngl6IXav\n44WgZwIloN4G/Zzx2b/tZwLyDR7r6nxe4ydPUT+5DEEQBB/hJ09RP7kMQRAEH+EnT1E/uQxBEAQf\n4SeriURAFgThlsZrAflwHfrHiYDcoNyo6FsdtcdPBeQWFs/V1hZBsq0l2XprSzjrcIs3bDtQHUEX\nGeftiBl2uXVQqUNoa6tgavVM/rN+HIDztOF5tYylOo1ztOEVtZg5ehrnaWMKtBcI4x31Is/q/+IK\nt7FKPcs4/Q7XCOF9NZ7ReiUA76vxPKpXA/CRSjXDZH+sRjuUrX2sx47T7wDQnMvmeGFcMOfRhvPm\n/G7nvDnvNpznl2qdeU2/VOscrtV6D6z3xnrPrPdS7zc+o1MWr+OzFi/xMxbv8R9BzTI8y88Yn/E0\n7IKuCMgen89rvMht3JjwS2MgCILQYPjJU9RPLkMQBMFH+MlTVHIgC4IgeENQHbY6MH36dOLi4khM\nTGTUqFGUlJQ47ZeZmUlsbCzR0dEsXLiwzsebaD8DkE022WTzePP2eVNxxvOtLuNt3rxZl5eXa621\nnjlzpp45c2aNPmVlZbpr1646Ly9Pl5aW6sTERJ2Tk+Px8VZ88gNn+vTpfPLJJ4SEhNC1a1dWrFhB\n69atAZg/fz7Lly8nMDCQN954gyFDhgCwd+9exo8fz9WrV0lJSWHJkiUuz39TBWQXolqtx+7C9wLy\nC5b9Jdh9wtzdAAAgAElEQVTz44JdQK70XG1hESRbW7yR24FKAb0Re8jrAdjd7ttV5UYGo2yIcCoc\nzurmALRVVxwEU2vI5m26PwDnaMPjahPr9HAuEMZ49T4r9WguEmYKtJe5jRnqTV7XU7jMbcxWC5mt\nZ1JKCPPUHGbpVwCYp+YwU88GYKGazTQ9B8AUgCvL1j7WYytDbDej1BwvjAtMUit5S48njAumUBzG\nBR5Rm8nQQ2jDeQaqbHboJAAGqmyHkN5WT2OXovFBo5wAhiO0XRCuFPPPWkT+Eov4X2LJe1xifMYv\nYBd03636Tqh3Hb87am3N75r61PV3Xe0CneC8DYzxjoK+x00f6/mO4nG+5BrHepE/2en5boKAXF5P\nT9Hk5GSz3LdvX9atW1ejT3Z2NlFRUURGRgIwZswYMjIyiIuL8+h4Kz55TTRkyBC+/vprDhw4QExM\nDPPnzwcgJyeHNWvWkJOTQ2ZmJpMnTzaXYU2aNIlly5aRm5tLbm4umZmZvpi6IAiCA+VBnm83yvLl\ny0lJSalRX1hYSKdOncx9m81GYWGhx8db8ckvA1cWKyMjg9TUVIKDg4mMjCQqKoqsrCy6dOnChQsX\nSEqy//f15JNP8tFHHzFs2DBfTF8QBMHkWrMQl21fbK/gyx0VlpoKh/bk5GROnjxZ47h58+YxYsQI\nAObOnUtISAhjx46t0c8Tvyp3x1vxuQ6+fPlyUlPti6CLioro16/qt2qllQsODsZms5n1ERERTq2f\nIAhCQ1Me6NoF+f7Bgdw/uGp/wZwrDu1btmxxe+6VK1eyceNGtm7d6rQ9IiKC/Px8cz8/P9/hWVnb\n8Q54rGbUkQcffFDHx8fX2NavX2/2+cMf/qBHjRpl7j///PP6z3/+s7mflpam//a3v+k9e/boBx98\n0Kz//PPP9cMPP+x0XBqBICWbbLI1nc0bAP29bunxVpfxNm3apLt3766///57l32uX7+u77nnHp2X\nl6evXbvmICB7cryVevtlcCMWr7qVKygowGazERERQUFBgUN9RESEy3PfkgLy+Gp1K6kpIM8ydkJB\nvYrdI9kazrqFRahsDeoXoNdiF5ArvWFbW0IrY5RzjXI0DmKyVTC15vk9rLsAcIEwktQhsnU8Fwnj\nZ2on23R/LnMbD6utfKIHc5nbGK0+5n09gss0N0XmcoJIU39xCBf9lnETKkVfZ2VrTmLrsZUexCGU\nmuPdxmVzHmFcMIXi5lw25x3GBeLUcfOa4tRxirV9McRdqsQxPLUr0diIh68GYA9TDXZBuFLML3Hh\ndVxiCVt9yZL3+KLx+Y+3fBdEQHZ+vpsgIJfVU3CiKVOmUFpaar5W79+/P+np6RQVFfHMM8+wYcMG\ngoKCWLp0KUOHDqW8vJy0tDTi4uLcHu8Kn7wmyszMZNGiRezYsYPQ0Ko0QSNHjmTs2LFMnTqVwsJC\ncnNzSUpKQilFq1atyMrKIikpiVWrVvHrX//aF1MXBEFwoLyeHqO5ublO6zt27MiGDRvM/eHDhzN8\n+HCPj3eFT4yBK4vVvXt3Ro8eTffu3QkKCiI9Pd0USNLT0xk/fjxXrlwhJSVFxGNBEBoF5X4SttQn\nxsCdxZo1axazZs2qUd+rVy8OHjxYn9MSBEGoM/5iDCSEtSAItzTePAKVUqZe5Alx6riEsG5I/FZA\nHuWi7QMXAvLzlv2lVIW0bmHxXHXljdwC1NOgV2AXjZ2FtsYoW8Vkq2eyRTC15kz+XrcE4DK30UWd\n5rgO5zLNTSH2CrfxE3WYr3Qc12hGf7WfnboHpYSYIu41QhiivmSzHgDAEPWlQ+hoT8qVntA/UztN\nD+IQSs3xbuMyieoIB3QMt3GZaFVArrbRnCvY1BkKdDvCyi/QOqjUQSy/boSACW4N+oRxDzq7EY03\nGuUUMCJs28Xh50C/bZQrP6NLHngdXwO1Cgyt3F62fHfUByIgm+e7GR7IfvIY9Y+rEARB8BH+8ppI\njIEgCIIXiDEQBEEQ6s3PoKERAVkQhFsabwXknbqHx/37q/0iIDckt6SAPK5a3SrMHLhg5MGt9EgO\ntXiutrB4I1tzI7ewCJgtLN6wLSxiMkb5U6M8tJqYbBVMjxrle6q8lK+1sAvKVy/B5dua01Zd4axu\nzjVCuEuVUKxbU0ozU2QuJcQUccsIrOH5e0h3BSBefedR2XpsZbjtEEpNcbgZ17hTXeR73ZKQ8lJT\nKL7tUinBreF6CQRdokZOaAeP7D1GuTdmeGo1sMrTWI2sEo1VqnG/wS7KVn4Wlyyf0SUPvI6vYc9z\nbOQ9VqtBj7R8F9aDtsTLAVBbRUC+UeQ1kSAIgkAprqOWNiXEGAiCIHiBv2gGYgwEQRC8wF/8DERA\nFgThlsZbAfmT6gKMGx5WW0VAbkj8VkAe6aJtfZVYaNatxiEvsnrX4pHcDNRiwyO5Gah5Rnjrlpbw\nyM08EJMxymuN8i+qedNaBVOrkFoZYqqFISgfxS5qG0KsbgYB7aDiDJSGVonMZYEBtAyt4OLVAMqD\ngmp4/lrzL3tSth5b6SEdWEaVOHzNHopbnwauGl7EJ7ALppV5oK+B6gF6v3F9Pap5F1vFdavobr1/\nhmisnjPuN3gmGl914XV8yVhUYCw4UB+AtmQ8VBtBD6z2fdnh+u9G7REB2R0iIAuCIAiiGQiCIAhQ\nSjNfT+GmEODrCQiCIDRlygn0eKsL06dPJy4ujsTEREaNGkVJSYnTfpmZmcTGxhIdHc3ChQtrtC9e\nvJiAgADOnj3rdjwRkAVBuKXxVkBeqUd73H+8et/j8bZs2cLgwYMJCAjgpZdeAmDBggUOfcrLy+nW\nrRufffYZERER9OnTh9WrV5upL/Pz83nmmWf45ptv2Lt3L23btnU5nl++JmpyAnKo4QHqZlGC2lqL\ngFwtP7Ja6+iVrFZZBOVQiwhZVzG5GVVhrjHKK4zy09W8aa1etlYh1RCWaUZVPuVQixDbzCLQBhme\nvLlAoEVwDrIIuhhlqxewq7K1v7VseEhzzYU4XGaI33uwC7cDDKH4ElX5oSuvz5V3sfU+We+f8c+c\nmmkIxuAoGl+l9lDVV43P/BdG/fqq74ta77g4QX16AwJynPM28541pIB8g8e6Op+31NfS0spMkAB9\n+/Zl3bp1NfpkZ2cTFRVFZGQkAGPGjCEjI8M0BlOnTuX111/nkUceqXU8eU0kCILgBfX1msjK8uXL\nSUlJqVFfWFhIp06dzH2bzUZhYSEAGRkZ2Gw27rvvPo/G8MtfBoIgCA2FNw/55ORkTp48WaN+3rx5\njBgxAoC5c+cSEhLC2LFja/Rz9Vr8ypUrzJs3jy1btph1tb2eEmMgCILgBe6MwZHtxeRuL3bZbn1Y\nO2PlypVs3LiRrVu3Om2PiIggPz/f3M/Pz8dms/Hdd99x7NgxEhMTASgoKKBXr15kZ2cTHu7cUUME\nZEEQbmm8FZAX68ke95+m0j0eLzMzk2nTprFjxw7uuOMOp33Kysro1q0bW7dupWPHjiQlJTkIyJXc\nfffdIiB7Q6MTkGu+MrS3baRGeGv1gaNXslpNVZ7kIIsI6YmYHGrxhm1GVT5ejPISo/xCNW9aV8Ky\n4bFMM0NkXW+UK4XYIKryLAcZHsw7sAvIlcJtkEV85gbLVk/hSlG7zDJ2mWVO1yxzvWZ4W6/FLtxW\n98h+1yhPrCYUG97FapohzGOI9K8a5d9Tlaf6mgeexq5E4zLjO2F8X9RGx++W2lrzu6n+0cACcmfP\n+tY49kRjFJDrx+lsypQplJaWmkJy//79SU9Pp6ioiGeeeYYNGzYQFBTE0qVLGTp0KOXl5aSlpdUw\nBODZP8l+aQwEQRAaivoyBrm5uU7rO3bsyIYNG8z94cOHM3z4cLfnOnr0qNt2EGMgCILgFf4SjqLW\npaU/+9nPHKwQwLPPPltvExIEQWhKlBPk8daYqdUY5OXlsXDhQl577TWzbvfu3fU6KUEQhKZCQ/gZ\nNAS1GoM2bdqwbds2Tp06xYgRIzh//nxDzEsQBKFJ4C/GoNalpT179mTfvn2Afc3r4sWLOXfuHAUF\nBQ0ywboiS0sFQagL3i4tnabneNx/sXql6Sa3ee6558zy+PHjSUhI4L//+7/rdVLeclOXlt7gudSe\nOiwtDTJiwwx03UXtwGXyG2eJb6rHK1JrLUtNgyzLE5tZli0GgXob9HNG2dmS0yBLzByMsnVppHXJ\npKvlp0st86hMmhNkLMV8F/sS0qeNpalBlqWbQcbSVOsyVSdLVs1ln87Kro41loZSZhm73DKnMstc\nrxnXscSor4zbhJNlo0asIfWqsUwXY8musYRULcZMOqSWUhU/qgzUSmM58DVjaXCqUa5MXFNmWU56\njaoESWWOS5XVVsfvltpRc8m02gW6B05R+2VpqTsauxbgKbVexa9+9SuH/V69erF8+fJ6m5AgCEJT\norG//vEU/zBpgiAIPsJfjIFPopa+8sorJCYm0qNHDwYPHuwQW2P+/PlER0cTGxvL5s2bzfq9e/eS\nkJBAdHQ0L7zwgi+mLQiCUIMyAj3eGjM+MQYzZszgwIED7N+/n0cffdRctpqTk8OaNWvIyckhMzOT\nyZMnm2LLpEmTWLZsGbm5ueTm5pKZmemLqQuCIDhwy/gZ1AdhYWFm+eLFi2YQpoyMDFJTUwkODiYy\nMpKoqCiysrIoLi7mwoULJCUlAfDkk0/y0Ucf+WLqgiAIDvjL0lKfmarf/va3rFq1iubNm5OdnQ1A\nUVER/fpVLXOoTNQQHByMzWYz6yMiIswEDoIgCL6klBBfT+GmUG+/DJKTk0lISKixffzxx4A9YcOJ\nEyd4+umnefHFF+trGoIgCPWKv2gGPs9ncOLECVJSUjh06JCZ7Lky+fOwYcN47bXX6NKlCz/96U85\nfPgwAKtXr2bHjh28/fbbNc4nTmeCINQFb53OHtd/9rj/OvXLput0Vh/k5uYSHR0N2HWCnj17AjBy\n5EjGjh3L1KlTKSwsJDc3l6SkJJRStGrViqysLJKSkli1ahW//vWvXZ6/0TideZoLoZmHTmcu8h2o\nrTUd0tSnjvkP1EaLY1qQxXEp0BITP8ji3GR1TAu0OECFWhzTMMpWpyljoZdaUs2xyuqkNssyj99X\nc1p71aivzJ+ApRxkOHJZnbpuVrkyEX2ZZR5lVOV0wFIux3my+srrtt4P630yHMrUu8Z9xbjHhjOg\nWm1xFCy3fEZXLfkJyl04l5VZvkNlRn4Cw+lR/cPxu+jMwUztB52AU9RB0NHO28y53tJOZ437P35P\n8YkxePnll/nmm28IDAyka9euvPXWWwB0796d0aNH0717d4KCgkhPTzf/009PT2f8+PFcuXKFlJQU\nhg0b5oupC4IgOCDGwAv+9re/uWybNWsWs2bNqlHfq1cvDh48WJ/TEgRBqDONXQvwFJ8sLRUEQfAX\n6svPYPr06cTFxZGYmMioUaMoKSlx2i8zM5PY2Fiio6NZuHChQ9ubb75JXFwc8fHxzJw50+14YgwE\nQRC8oJQQj7e6MGTIEL7++msOHDhATEwM8+fPr9GnvLyc559/nszMTHJycli9erW50Obvf/8769ev\n55///CeHDh3iN7/5jdvxxBgIgiB4QX0tLU1OTiYgwP6I7tu3r9O0AdnZ2URFRREZGUlwcDBjxowh\nIyMDgLfeeouXX36Z4OBgAO68806344kxEARB8IKGCEexfPlyUlJSatQXFhbSqVMnc7/SURfsqzY/\n//xz+vXrx6BBg9izZ4/bMRp3sAxBEIRGjjeriZKTkzl58mSN+nnz5jFixAjA7qAbEhLC2LFja/Rz\n51dVVlbGuXPn2LVrF7t372b06NEcPXrUZX8xBoIgCF7gzhhc2P4VF7bvc9m+ZcsWt+deuXIlGzdu\nZOvWrU7bIyIiHKI+5+fnm6F7bDYbo0aNAqBPnz4EBARw5swZ2rVr5/RcYgwEQRC8wJ0xuG1QH24b\n1MfcL35thcfnzczMZNGiRezYsYPQ0FCnfXr37k1ubi7Hjh2jY8eOrFmzhtWr7Sn8Hn30UbZt28bA\ngQM5cuQIpaWlLg0BNIJwFDcbCUchCEJd8DYcRZz+yuP+h9VPPB4vOjqa0tJS2rZtC0D//v1JT0+n\nqKiIZ555hg0bNgCwadMmXnzxRcrLy0lLS+Pll18G4Pr160yYMIH9+/cTEhLC4sWLGTRokOtr8Udj\n0OTCUQQ6hg9wej437eofNUNZVA9fUT1khRnSINAS6gBLOciSXxdL2RrKAqNshFAww1pUlq1hFqzl\ncZZ5OAt5gYtyEFX5mrnBsrvzgz20wiocw0WMwx7iofI6ypxcq7VsvTfW+1f9HleWh1o+k8rPzFV4\niXIjnITx3TLLZZbva5kRXsIIOWEtm/vVQk+og67zHKvD/h2OwltjEKMPeNz/iEqU2ESCIAj+iISj\nEARBEPwmHIUYA0EQBC9o7OksPcU/rkIQBMFHyGsiQRAEQYyBIAiCANdK/SMHshgDQRAELygv84/H\nqH9chSAIgo8oL5PXRIIgCLc8YgwEQRAEyq6LMRAEQbjlqSj3j8eof1yFIAiCr/CT10R+GahOEATB\nU7wNVMc3dTi+m5JAdQ1Jo4laWodja+vvLgqqNYqlqzqX+9UiptZaDrJE08S7MkFGNFUjUudNLd/o\nnMo8uAeuIoc6KxufqfXzrVF2Fl203IgkakQXdVouMyKKGtFGzXIZqNyqSKPWsrN9s85F1NFaI5KW\nGdFEPYxEWpe+N/NYV+fzmrKbcI5GgF8aA0EQhAZDjIEgCIIgxkAQBEGA676ewM0hwNcTEARBaNKU\n12GrA9OnTycuLo7ExERGjRpFSUmJ036ZmZnExsYSHR3NwoULzfrs7GySkpLo2bMnffr0Yffu3W7H\nE2MgCILgDWV12OrAkCFD+Prrrzlw4AAxMTHMnz+/Rp/y8nKef/55MjMzycnJYfXq1Rw+fBiAGTNm\nMGfOHPbt28fvf/97ZsyY4XY8MQaCIAjecLUOWx1ITk4mIMD+iO7bty8FBQU1+mRnZxMVFUVkZCTB\nwcGMGTOGjIwMAO666y7z18T58+eJiIhwO55oBoIgCN7QAALy8uXLSU1NrVFfWFhIp06dzH2bzUZW\nVhYACxYsYMCAAfzmN7+hoqKCnTt3uh1DjIEgCII3eGEMkpOTOXnyZI36efPmMWLECADmzp1LSEgI\nY8eOrdHPnZNtWloab7zxBo899hhr165lwoQJbNmyxWV/MQaCIAje4M4YHNoOX2932ezu4QywcuVK\nNm7cyNatW522R0REkJ+fb+7n5+djs9kA+yukzz77DICf//znTJw40e1YPtUMFi9eTEBAAGfPnjXr\n5s+fT3R0NLGxsWzevNms37t3LwkJCURHR/PCCy/4YrqCIAg1ue5m6zYIRs2u2upAZmYmixYtIiMj\ng9DQUKd9evfuTW5uLseOHaO0tJQ1a9YwcuRIAKKiotixYwcA27ZtIyYmxu14PvtlkJ+fz5YtW+jS\npYtZl5OTw5o1a8jJyaGwsJAHH3yQ3NxclFJMmjSJZcuWkZSUREpKCpmZmQwbNsxX0xcEQbBTxyWj\nnjJlyhRKS0tJTk4GoH///qSnp1NUVMQzzzzDhg0bCAoKYunSpQwdOpTy8nLS0tKIi7PHJ3nnnXf4\n93//d65du0bz5s1555133I7nM2MwdepUXn/9dR555BGzLiMjg9TUVIKDg4mMjCQqKoqsrCy6dOnC\nhQsXSEpKAuDJJ5/ko48+EmMgCILvqScBOTc312l9x44d2bBhg7k/fPhwhg8fXqNf7969TTHZE3xi\nDDIyMrDZbNx3330O9UVFRfTrVxVNzWazUVhYSHBwsPkeDOzvyQoLCxtsvoIgCC6p45LRxkq9GQNX\nKvncuXOZP3++gx5ws0O6zrbYiUFhMKjVTT29IAhNlO1X7dtNRWITuceVSn7o0CHy8vJITEwEoKCg\ngF69epGVlVVDGS8oKMBmsxEREeHgcFFQUODWgWK2e98KQRBuUQaF2rdKXvvxJpzUT4xBg68mio+P\n59SpU+Tl5ZGXl4fNZuOrr76iffv2jBw5kr/+9a+UlpaSl5dHbm4uSUlJdOjQgVatWpGVlYXWmlWr\nVvHoo4829NQFQRBqUk/hKBoan/sZWJ0munfvzujRo+nevTtBQUGkp6eb7enp6YwfP54rV66QkpIi\n4rEgCI0DP4laKmkvBUG4pfE67eXcOhz/W0l72aCYqQO9xCENYT0fW1t/d+1qf1VKRLPuoGOdu/06\nlV2lWryBsss2N2kbb6hspGy0pm90Vfa0n9uykZbRmqKxrmVvj/f0vO7qPGmrS58b6Xszj3V1Pq+R\n1USCIAhCY9cCPEWMgSAIgjf4iWYgxkAQBMEb6ikcRUMjxkAQBMEb5DWRIAiCIMZAEARBEM1AEARB\nAK75egI3BzEGgiAI3iCviQRBEAR5TSQIgiD4zdJSn+ZAFgRBaPLUU9TS6dOnExcXR2JiIqNGjaKk\npMRpvwkTJtC+fXsSEhxj0pw9e5bk5GRiYmIYMmQI58+fdzueGANBEARvqCdjMGTIEL7++msOHDhA\nTEwM8+fPd9rv6aefJjMzs0b9ggULSE5O5siRIwwePJgFCxa4HU+MgSAIgjdcr8NWB5KTkwkIsD+i\n+/bt65Dgy8oDDzzA7bffXqN+/fr1PPXUUwA89dRTfPTRR27HE81AEATBGxpgaeny5ctJTU2t0zGn\nTp2iffv2ALRv355Tp0657S/GQBAEwRu8WFrqKlf8vHnzGDFiBGDPGx8SEsLYsWNveBylVK25XsQY\nCIIgeIO71z8/bocL2102u8oVX8nKlSvZuHEjW7durfO02rdvz8mTJ+nQoQPFxcWEh4e77S+agSAI\ngjeUu9laDIIOs6u2OpCZmcmiRYvIyMggNDS0ztMaOXIk7733HgDvvfderXnjxRgIgiB4Qz2tJpoy\nZQoXL14kOTmZnj17MnnyZACKiop46KGHzH6pqancf//9HDlyhE6dOrFixQoAXnrpJbZs2UJMTAzb\ntm3jpZdecjuevCYSBEHwhnoKR5Gbm+u0vmPHjmzYsMHcX716tdN+bdu25bPPPvN4PDEGgiAI3iDh\nKARBEASJWioIgiBI1FJBEAQBv3lNpLTW2teTuJnU5lghCIJgxZtHoFIKbq/D8eeUV+PVJ375y0D3\nuDnnUftv/FxqP+iE2vuZ/Q+6768Ogo5z0Xa4Zps6DDrasp/rYr8c1FHQ9xj1zsploE6A7mzUuyt3\nNMpFbsoW3xd1umrfabkc1BnQ7Yz6GyiXtraXQ0pqL3vaz9tybfOu87UGenAvndxzd3UObR2dt5l9\nimrv49C3s2d9axx74saPdXU+r5HXRIIgCIIYA0EQBMFvNAMxBoIgCN7gJ78MJByFIAiCIMZAEARB\nEGMgCIIg4CNjMHv2bGw2Gz179qRnz55s2rTJbJs/fz7R0dHExsayefNms37v3r0kJCQQHR3NCy+8\n4ItpC4IgOKGe8l42MD4xBkoppk6dyr59+9i3bx/Dhw8HICcnhzVr1pCTk0NmZiaTJ082HTQmTZrE\nsmXLyM3NJTc312kCaEEQhIannmJYNzA+e03kzAsvIyOD1NRUgoODiYyMJCoqiqysLIqLi7lw4QJJ\nSUkAPPnkk7UmdxYEQWgY5JeBV7z55pskJiaSlpbG+fPnAXvSBpvNZvax2WwUFhbWqI+IiKCwsLDB\n5ywIglCTK3XYGi/15mfgKtHz3LlzmTRpEq+++ioAr7zyCtOmTWPZsmU3bezZxVXlQS1hUNhNO7Ug\nCE2Y7Vft282lcf/H7yn1ZgxqS/RcycSJExkxYgRg/48/Pz/fbCsoKMBmsxEREUFBQYFDfUREhMtz\nzr7rBictCIJfMyjUvlXy2o8346z1owVMnz6dTz75hJCQELp27cqKFSto3bp1jX4TJkxgw4YNhIeH\nc/DgwTofX4lPXhMVF1f96/7hhx+SkGCP0DZy5Ej++te/UlpaSl5eHrm5uSQlJdGhQwdatWpFVlYW\nWmtWrVpVa3JnQRCEhqF+NIMhQ4bw9ddfc+DAAWJiYpg/f77Tfk8//bTTBTWeHl+JT8JRzJw5k/37\n96OU4u677+Z//ud/AOjevTujR4+me/fuBAUFkZ6eboakTk9PZ/z48Vy5coWUlBSGDRvmi6kLgiBU\no35+GSQnJ5vlvn37sm7dOqf9HnjgAY4dO3bDx1fiE2Pwpz/9yWXbrFmzmDVrVo36Xr16OfwEEgRB\naBzUv2awfPlyUlNT6/V4CVQnCILgFe5WCe0FvnLZ6mqhzbx580wtde7cuYSEhDB27Ngbmp2nx4sx\nEARB8Ap3r4kSja0Sx1WTtS20WblyJRs3bmTr1q03NLO6HC/GQBAEwSvq5zVRZmYmixYtYseOHYSG\nhtZ+gJfHS6A6QRAEr6ifcBRTpkzh4sWLJCcn07NnTyZPngzYnXMfeughs19qair3338/R44coVOn\nTqxYscLt8a6QXwaCIAheUT+/DHJzc53Wd+zYkQ0bNpj7q1evrtPxrhBjIAiC4BWNOwCdp4gxEARB\n8AoJRyEIgiA08gB0niLGQBAEwSvkl4EgCIIgmoEgCILgL78MxM+gHtl+0dcz8J7t/vE6lB1+8M/b\ndv945vghkvZSqIXtl3w9A+8RY9B4EGPQWPGPtJfymkgQBMEr/OA/DcQYCIIgeIl//HxWWmvt60nc\nTAYNGsSOHTt8PQ1BEJoAAwcOZPv27Td8fGXyLU+5/fbbOXv27A2PV5/4nTEQBEEQ6o4IyIIgCIIY\nA0EQBEGMwU1l8eLFBAQEOLwTnD9/PtHR0cTGxrJ582azfu/evSQkJBAdHc0LL7zgi+nWYPr06cTF\nxZGYmMioUaMoKSkx25rSdVjJzMwkNjaW6OhoFi5c6OvpuCQ/P5+f/vSn3HvvvcTHx/PGG28AcPbs\nWZKTk4mJiWHIkCGcP3/ePMbVZ9IYKC8vp2fPnmbqxqZ6HbcUWrgpnDhxQg8dOlRHRkbqM2fOaK21\n/grupXAAAATlSURBVPrrr3ViYqIuLS3VeXl5umvXrrqiokJrrXWfPn10VlaW1lrr4cOH602bNvls\n7pVs3rxZl5eXa621njlzpp45c6bWuuldRyVlZWW6a9euOi8vT5eWlurExESdk5Pj62k5pbi4WO/b\nt09rrfWFCxd0TEyMzsnJ0dOnT9cLFy7UWmu9YMECt59J5WfXGFi8eLEeO3asHjFihNZaN9nruJWQ\nXwY3ialTp/L666871GVkZJCamkpwcDCRkZFERUWRlZVFcXExFy5cICkpCYAnn3ySjz76yBfTdiA5\nOZmAAPtXom/fvhQUFABN7zoqyc7OJioqisjISIKDgxkzZgwZGRm+npZTOnToQI8ePQBo2bIlcXFx\nFBYWsn79ep566ikAnnrqKfP+OvtMsrOzfTZ/KwUFBWzcuJGJEyeijfUpTfE6bjXEGNwEMjIysNls\n3HfffQ71RUVF2Gw2c99ms1FYWFijPiIigsLCwgabrycsX76clJQUoOleR2FhIZ06dTL3K+fd2Dl2\n7Bj79u2jb9++nDp1ivbt2wPQvn17Tp06Bbj+TBoD//Ef/8GiRYvMfyyAJnkdtxridOYhycnJnDx5\nskb93LlzmT9/vsO7Tt2IV+u6uo558+aZ73fnzp1LSEgIY8eObejp3VTquga8MXDx4kUef/xxlixZ\nQlhYmEObUsrtNTWG6/3kk08IDw+nZ8+eLtfvN4XruBURY+AhW7ZscVp/6NAh8vLySExMBOw/kXv1\n6kVWVhYRERHk5+ebfQsKCrDZbERERJivYCrrIyIi6vcCDFxdRyUrV65k48aNbN261axrjNfhCdXn\nnZ+f7/BfaGPj+vXrPP7444wbN45HH30UsP8XffLkSTp06EBxcTHh4eGA88+kMdz7//3f/2X9+vVs\n3LiRq1ev8uOPPzJu3Lgmdx23JL4WLfwNZwLytWvX9NGjR/U999xjCq9JSUl6165duqKiotEIr5s2\nbdLdu3fX33//vUN9U7uOSq5fv67vuecenZeXp69du9aoBeSKigo9btw4/eKLLzrUT58+XS9YsEBr\nrfX8+fNrCK/OPpPGwvbt2/XDDz+stW7a13GrIMbgJnP33XebxkBrrefOnau7du2qu3XrpjMzM836\nPXv26Pj4eN21a1c9ZcoUX0y1BlFRUbpz5866R48eukePHnrSpElmW1O6DisbN27UMTExumvXrnre\nvHm+no5LvvjiC62U0omJieb937Rpkz5z5owePHiwjo6O1snJyfrcuXPmMa4+k8bC9u3bzdVETfk6\nbhUkHIUgCIIgq4kEQRAEMQaCIAgCYgwEQRAExBgIgiAIiDEQBEEQEGMgCIIgIMZAEARBQIyBIAiC\ngBgDwY/ZvXs3iYmJXLt2jUuXLhEfH09OTo6vpyUIjRLxQBb8mldeeYWrV69y5coVOnXqxMyZM309\nJUFolIgxEPya69ev07t3b5o3b87OnTslPLIguEBeEwl+zQ8//MClS5e4ePEiV65c8fV0BKHRIr8M\nBL9m5MiRjB07lqNHj1JcXMybb77p6ykJQqNEktsIfsuf/vQnmjVrxpgxY6ioqOD+++9n+/btDBo0\nyNdTE4RGh/wyEARBEEQzEARBEMQYCIIgCIgxEARBEBBjIAiCICDGQBAEQUCMgSAIgoAYA0EQBAEx\nBoIgCALwfwEwulPiyuCzdQAAAABJRU5ErkJggg==\n", "text": [ - "" + "" ] } ], - "prompt_number": 8 + "prompt_number": 31 }, { "cell_type": "code", "collapsed": false, "input": [ - "survey = DC.SurveyDC(txList)\n", - "problem = DC.ProblemDC(mesh)\n", - "problem.pair(survey)" + "dpred = survey.dpred(mopt)" ], "language": "python", "metadata": {}, "outputs": [], - "prompt_number": 13 + "prompt_number": 22 }, { "cell_type": "code", "collapsed": false, "input": [ - "datasyn = survey.dpred(sigma)" + "plot(np.c_[dpred,survey.dobs])" ], "language": "python", "metadata": {}, "outputs": [ { - "ename": "KeyboardInterrupt", - "evalue": "", - "output_type": "pyerr", - "traceback": [ - "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m\n\u001b[1;31mKeyboardInterrupt\u001b[0m Traceback (most recent call last)", - "\u001b[1;32m\u001b[0m in \u001b[0;36m\u001b[1;34m()\u001b[0m\n\u001b[1;32m----> 1\u001b[1;33m \u001b[0mdatasyn\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0msurvey\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mdpred\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0msigma\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m", - "\u001b[1;32mC:\\Users\\SEOGI\\simpeg\\SimPEG\\Utils\\CounterUtils.pyc\u001b[0m in \u001b[0;36mwrapper\u001b[1;34m(self, *args, **kwargs)\u001b[0m\n\u001b[0;32m 81\u001b[0m \u001b[0mcounter\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mgetattr\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;34m'counter'\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mNone\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 82\u001b[0m \u001b[1;32mif\u001b[0m \u001b[0mtype\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mcounter\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mis\u001b[0m \u001b[0mCounter\u001b[0m\u001b[1;33m:\u001b[0m \u001b[0mcounter\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mcount\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0m__class__\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0m__name__\u001b[0m\u001b[1;33m+\u001b[0m\u001b[1;34m'.'\u001b[0m\u001b[1;33m+\u001b[0m\u001b[0mf\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0m__name__\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m---> 83\u001b[1;33m \u001b[0mout\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mf\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;33m*\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;33m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 84\u001b[0m \u001b[1;32mreturn\u001b[0m \u001b[0mout\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 85\u001b[0m \u001b[1;32mreturn\u001b[0m \u001b[0mwrapper\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", - "\u001b[1;32mC:\\Users\\SEOGI\\simpeg\\SimPEG\\Utils\\codeutils.pyc\u001b[0m in \u001b[0;36mrequiresVarWrapper\u001b[1;34m(self, *args, **kwargs)\u001b[0m\n\u001b[0;32m 222\u001b[0m \u001b[1;32mif\u001b[0m \u001b[0mgetattr\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mvar\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mNone\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mis\u001b[0m \u001b[0mNone\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 223\u001b[0m \u001b[1;32mraise\u001b[0m \u001b[0mException\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mextra\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 224\u001b[1;33m \u001b[1;32mreturn\u001b[0m \u001b[0mf\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;33m*\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;33m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 225\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 226\u001b[0m \u001b[0mdoc\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mrequiresVarWrapper\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0m__doc__\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", - "\u001b[1;32mC:\\Users\\SEOGI\\simpeg\\SimPEG\\Survey.pyc\u001b[0m in \u001b[0;36mdpred\u001b[1;34m(self, m, u)\u001b[0m\n\u001b[0;32m 305\u001b[0m \u001b[0mWhere\u001b[0m \u001b[0mP\u001b[0m \u001b[1;32mis\u001b[0m \u001b[0ma\u001b[0m \u001b[0mprojection\u001b[0m \u001b[0mof\u001b[0m \u001b[0mthe\u001b[0m \u001b[0mfields\u001b[0m \u001b[0monto\u001b[0m \u001b[0mthe\u001b[0m \u001b[0mdata\u001b[0m \u001b[0mspace\u001b[0m\u001b[1;33m.\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 306\u001b[0m \"\"\"\n\u001b[1;32m--> 307\u001b[1;33m \u001b[1;32mif\u001b[0m \u001b[0mu\u001b[0m \u001b[1;32mis\u001b[0m \u001b[0mNone\u001b[0m\u001b[1;33m:\u001b[0m \u001b[0mu\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mprob\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mfields\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mm\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 308\u001b[0m \u001b[1;32mreturn\u001b[0m \u001b[0mUtils\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mmkvc\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mprojectFields\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mu\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 309\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n", - "\u001b[1;32mC:\\Users\\SEOGI\\Documents\\simpegdc\\simpegDC\\BaseDC.pyc\u001b[0m in \u001b[0;36mfields\u001b[1;34m(self, m)\u001b[0m\n\u001b[0;32m 145\u001b[0m \u001b[0mAinv\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mSolver\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mA\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 146\u001b[0m \u001b[0mQ\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0msurvey\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mgetRhs\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mmesh\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 147\u001b[1;33m \u001b[0mPhi\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mAinv\u001b[0m \u001b[1;33m*\u001b[0m \u001b[0mQ\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 148\u001b[0m \u001b[1;32mreturn\u001b[0m \u001b[0mPhi\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 149\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n", - "\u001b[1;32mC:\\Users\\SEOGI\\simpeg\\SimPEG\\Utils\\SolverUtils.pyc\u001b[0m in \u001b[0;36m__mul__\u001b[1;34m(self, b)\u001b[0m\n\u001b[0;32m 48\u001b[0m \u001b[0mX\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mi\u001b[0m\u001b[1;33m]\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0msolver\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0msolve\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mb\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mi\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 49\u001b[0m \u001b[1;32melse\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m---> 50\u001b[1;33m \u001b[0mX\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mi\u001b[0m\u001b[1;33m]\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mfun\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mA\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mb\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mi\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m**\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 51\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 52\u001b[0m \u001b[1;32mif\u001b[0m \u001b[0mcheckAccuracy\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", - "\u001b[1;31mKeyboardInterrupt\u001b[0m: " + "metadata": {}, + "output_type": "pyout", + "prompt_number": 24, + "text": [ + "[,\n", + " ]" + ] + }, + { + "metadata": {}, + "output_type": "display_data", + "png": "iVBORw0KGgoAAAANSUhEUgAAAXUAAAEACAYAAABMEua6AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl80/Xhx/FX2qRHUii0QIGm2oPSlqsUy6HiLF6ACiLH\nLDplypDpcHNu6tCpuJ9TOvW3qajDazqdBafOIkIRxAqOKSKXjsNySahcLS29jyTf3x+dnf31gEpL\n+Kbv5+ORB03ySfJOH+27Hz75HhbDMAxERMQvBPg6gIiItB+VuoiIH1Gpi4j4EZW6iIgfUamLiPgR\nlbqIiB85Yann5uaSnJxMYmIiWVlZTe5/7LHHSEtLIy0tjcGDB2O1WikpKemQsCIi0jpLa9upezwe\nkpKSWLVqFdHR0QwfPpzs7GxSUlKaHb906VL+9Kc/sWrVqg4LLCIiLWt1pr5+/Xr69etHbGwsNpuN\nzMxMcnJyWhz/+uuvM3369HYPKSIiJ6fVUi8oKCAmJqbhutPppKCgoNmxlZWVrFixgilTprRvQhER\nOWmtlrrFYjnpJ3r33XcZPXo03bp1O+VQIiLy/VhbuzM6OhqXy9Vw3eVy4XQ6mx27aNGiVpde+vXr\nx+7du79nTBGRzikhIYFdu3ad/AOMVtTV1Rnx8fHG3r17jZqaGiM1NdXYtm1bk3ElJSVGRESEUVlZ\n2eJzneClzngPPPCAryOcEjPnN3N2w1B+XzN7/rZ2Z6szdavVyoIFCxg7diwej4eZM2eSkpLCwoUL\nAZg9ezYA77zzDmPHjiU0NPT7/jESEZF20GqpA4wfP57x48c3uu3bMv/WjBkzmDFjRvsmExGRNtMe\npScpIyPD1xFOiZnzmzk7KL+vmT1/W7W681G7vpDFwml6KRERv9HW7tRMXUTEj6jURUT8iEpdRMSP\nqNRFRPyISl1ExI+o1EVE/IhKXUTEj6jURUT8iEpdRMSPqNRFRPyI35b61q3gdvs6hYjI6eWXpb56\nNaSlwa9/7eskIiKnl9+VussF110Hb7wBubnwn0O/N+v5z58n881M9hbvbfPreL31/xOorYXqaqis\nhLq6UwguItIO/KrUa2pg6lT45S9hyhRYuhTuv79+5v5dtZ5abll6C3/85I8kRSYx4oURPPPZM3gN\nb7PPaxiQnw+vvQa33QbnjKwieNAygrsXEhYG4eEQGQm9e8Orr9aPFxHxBb869O7s2VBUBFnP7eGH\nb05jYM+BjAq8lXmzRvLxWgv9+8ORiiNMeWMK3UO689rk1+ga3JUdhTu4MedGQqwhvDDhBRIiEjAM\n2LSpvsgXLwarFVLPO0z14GfYYPkz/SLj+OrYDi6Ku4gfD/0x4/uN54stNm64Afr3hz//GXr16tC3\nKyKdQFu703Sl7va6qa6rJSzY3uj2l16CRx+Fxe/vYcJbY/jFyF9gGAbPbniWuvJwqtfeysuPJfPT\nldfyo0EzmBI5jy+2BlBQAA4H2B0ePqr9E0uOPUKC9wr2b+sN5b0YPawX4y7sxmfl/+AfO/5B5sBM\nfjHqFyT3SOZ49XH+vu3vvLz5ZfKP5XNr+q3cPep+5s2z8PLL8MwzMHnyKb9lEenE/LrUPV4PSfdP\nZnfNp5y9/Y+c3y2TQQMt9OoFc+fCa8t2M+vji5g7ei4/Tf8pAF7Dy/u73+dnLz/DPmMtzo3PcyRv\nKvHxMHQonHVW/Xp4eXn95VDtbgLi80gedpSQyCMcrTzC0cqj/OCsHzA7fTY97D2azZZflM+0v0/j\n5nNu5tbht7JuHcyYAUlJMGsWXH452Gyn9PZFpBPy61Kf+MwveX/LFl77yYPcs/Y2gup6MbLwaQ5v\nT2TSjbv5H1fjQv8ujwdycgzOOsvCwIHQEefI3nVsF+e9eB5Lr13KiOgRVFZCdjb85S+wa1f9B7g/\n/jEMHlw/vqK2gmc+e4ZnNzzLHefewZwRc1p8bo/XQ1ltGd1CurV/cBE5Y/ltqT/0/gLmLXuadyet\nY3xGd9xeN0988gSPfPwIs4bN4vUvX+ee0fcwO332iZ+sA729/W3uWHEHG2dvJCI0ouH2/Hx4+WX4\n61+hW68KnJOeYUPQ44yJ/wE3pd3EnGVzuHHojdxzwT1YLJaGx7nd8PXRo9z03nUYAW7W3LS6mVcV\nEX/ll6W+ZMd7TPvrT/hZ6D/53/viG923//h+7lp5F5fEX8JPhv2kPaKesl+t+BU7i3ayZPoSAiz/\n3cCotKaUp9c/y6Nr/0hE2YUUv3MfcY5BTJ4MFQEH+XPFZXQrHI9jXRbfFFgoK4O6qE9g2g+xHkkn\nOGovZY9t8uE7E5HTzVSlXuup5VD5IQ6WHeRg+UGKq4qJ7x7P4KjBDbPczYc2c/7CS0na+C6fvT2K\nwMDTkfbU1HnqyHglgysTr2TuBXMprCzkiU+e4NkNzzK231jmjp7LoF6DcLthzRrIyYGgIIh0HuOl\n2vEM7jmUBeOeYdHuZ3lk3e94YeILuA8nce2SiVQ/ttPXb09ETqMzutSLq4pZtWcVubtyWblnJd+U\nfUMvRy/6dulLn7A+dAvpxu7i3Xx55EscNgeDeg1iU8G/qVvyJ/7992lER5+OpO3jQOkBhj8/nMsS\nLuPdne8ybcA07jr/LhIiElp9XFlNGRMXTcR13EVYUBhv/fAtEiIS+GjTAS5ZNIq6rAOn6R2IyJmg\nraVu7cAsTUTNjyHGuIB+jOMK404CSvpR+O9ADh+G3UeguBi6d4dhUQZdol1Yor7A8s8gXrv3UlMV\nOoCzq5NFUxaxLH8ZX9zyBdFdT+4NdAnuwrJrl/Hq1le5fsj1hNrqP9GN6GLHE1DZkZFFxA+c1pn6\nfQ9WYiUUw6jf67JbN4iKqt9JJyqq/npxMRw+/N9Lnz4wffrpSHhmO3i0hr5PdsX4nxpfRxGR0+iM\nXn45TS/ll6qrDUIfsVL3QA3WgNP6HywR8aG2dqdfHfvFnwUHW6DOwfFKLcGISMtOWOq5ubkkJyeT\nmJhIVlZWs2Py8vJIS0tj0KBBZGRktHdGASwWsNQ5OHq8wtdRROQM1ur/4z0eD3PmzGHVqlVER0cz\nfPhwJk6cSEpKSsOYkpISfvazn7FixQqcTieFhYUdHrqzCvDYKSqtBJN9aCwip0+rM/X169fTr18/\nYmNjsdlsZGZmkpOT02jM66+/zpQpU3A6nQD06NH8sVHk1AV6HRSVaqYuIi1rtdQLCgqIiYlpuO50\nOikoKGg0Jj8/n2PHjjFmzBjS09N59dVXOyapYDXsHCtXqYtIy1pdfvnuMUhaUldXx8aNG/nggw+o\nrKzk3HPPZdSoUSQmJrZbSKlnM/RBqYi0rtVSj46OxuVyNVx3uVwNyyzfiomJoUePHoSGhhIaGsoP\nfvADtmzZ0mypz5s3r+HrjIwMfajaRjYclFRopi7iz/Ly8sjLy/vej291O3W3201SUhIffPABffv2\nZcSIEWRnZzf6oHTHjh3MmTOHFStWUFNTw8iRI1m8eDEDBgxo/ELaTv2UOW/P5JrUSTx+Y6avo4jI\nadKuhwmwWq0sWLCAsWPH4vF4mDlzJikpKSz8z9mcZ8+eTXJyMuPGjWPIkCEEBAQwa9asJoUu7SM4\nwMHxKs3URaRl2qPURAbeeRuD+iSy+I6f+zqKiJwm2qPUj9mtDspr9EGpiLRMpW4idpuDilotv4hI\ny1TqJuIIslNZp5m6iLRMpW4iYcEOKt2aqYtIy1TqJtIl2E61V6UuIi1TqZtI1xAHNR4tv4hIy1Tq\nJtI11E6NoZm6iLRMpW4i3RwOag3N1EWkZSp1E+nucFBn0UxdRFqmUjeR7mF23Cp1EWmFSt1EIrs6\n8ARo+UVEWqZSN5GILna8Vs3URaRlKnUT6RnuwFCpi0grVOomEtE1BAJrqXN7fB1FRM5QKnUTCQy0\nQJ2dY+VaVxeR5qnUTcbidlB4XKUuIs1TqZtMgMdOUanW1UWkeSp1k7F6HRSVqdRFpHkqdZOxeh0U\na01dRFqgUjcZK3aKyzVTF5HmqdRNJsji4HilZuoi0jyVuskEW+wcr9JMXUSap1I3mZAAB6UqdRFp\ngUrdZEKtDsprtPwiIs1TqZtMSKCd8lrN1EWkeSp1k3EEOaio1UxdRJqnUjcZh81OZZ1m6iLSPJW6\nyYQFO6hyq9RFpHkqdZPpGuKg2qPlFxFp3glLPTc3l+TkZBITE8nKympyf15eHuHh4aSlpZGWlsZD\nDz3UIUGlXpdQO9VezdRFpHnW1u70eDzMmTOHVatWER0dzfDhw5k4cSIpKSmNxl144YUsWbKkQ4NK\nvXC7g1pDM3URaV6rM/X169fTr18/YmNjsdlsZGZmkpOT02ScYRgdFlAa62a3U4tm6iLSvFZLvaCg\ngJiYmIbrTqeTgoKCRmMsFgvr1q0jNTWVyy+/nG3btnVMUgGge5iDOpW6iLSg1eUXi8VywicYNmwY\nLpcLu93O8uXLmTRpEl999VWzY+fNm9fwdUZGBhkZGW0KK9Dd4cAToOUXEX+Vl5dHXl7e9358q6Ue\nHR2Ny+VquO5yuXA6nY3GdOnSpeHr8ePHc+utt3Ls2DEiIiKaPN93S12+n4gudjwBmqmL+Kv/P+F9\n8MEH2/T4Vpdf0tPTyc/PZ9++fdTW1rJ48WImTpzYaMzhw4cb1tTXr1+PYRjNFrq0jx7hDryBmqmL\nSPNanalbrVYWLFjA2LFj8Xg8zJw5k5SUFBYuXAjA7NmzefPNN3n22WexWq3Y7XYWLVp0WoJ3Vj26\n2jFsmqmLSPMsxmnadMVisWgrmXbg8Xqx/s5K9T1ugoO075iIv2trd6oVTCYwIADcoRQer/J1FBE5\nA6nUTSjAbafwuJZgRKQplboJBXgcFJXqw1IRaUqlbkKBXjvHyjRTF5GmVOomZDUcFFeo1EWkKZW6\nCdlwUFKh5RcRaUqlbkJB2Cmp1ExdRJpSqZtQcICD0irN1EWkKZW6CYUE2imt0kxdRJpSqZtQaKCD\nshqVuog0pVI3IbvVQUWtll9EpCmVugnZbXYqajVTF5GmVOomFBbsoLJOM3URaUqlbkKOIDtVHs3U\nRaQplboJdQlxUK1SF5FmqNRNKDzUQY1Xyy8i0pRK3YTC7XZqDM3URaQplboJdbM7qEMzdRFpSqVu\nQt3C7LgtmqmLSFMqdROKCHOo1EWkWSp1E4rs6sATqOUXEWlKpW5CPbra8QZqpi4iTanUTSiyqwNs\nFXi9vk4iImcalboJhQWHgq2KqirD11FE5AyjUjehwIBA8ARTeLzK11FE5AyjUjepALeDolJ9WCoi\njanUTSrAY6eoVB+WikhjKnWTsnodFJWp1EWksROWem5uLsnJySQmJpKVldXiuM8++wyr1crbb7/d\nrgGleVbslFRo+UVEGmu11D0eD3PmzCE3N5dt27aRnZ3N9u3bmx139913M27cOAxDW2ScDjbDQXGF\nZuoi0lirpb5+/Xr69etHbGwsNpuNzMxMcnJymox76qmnmDp1Kj179uywoNJYsMXB8UrN1EWksVZL\nvaCggJiYmIbrTqeTgoKCJmNycnK45ZZbALBYLB0QU/6/oAA7pVWaqYtIY9bW7jyZgr799tuZP38+\nFosFwzBaXX6ZN29ew9cZGRlkZGScdFBpLCTAQVm1Sl3E3+Tl5ZGXl/e9H99qqUdHR+NyuRquu1wu\nnE5nozGff/45mZmZABQWFrJ8+XJsNhsTJ05s8nzfLXU5NaFWO+U1Wn4R8Tf/f8L74IMPtunxrZZ6\neno6+fn57Nu3j759+7J48WKys7MbjdmzZ0/D1zfeeCMTJkxottClfdltDsprNVMXkcZaLXWr1cqC\nBQsYO3YsHo+HmTNnkpKSwsKFCwGYPXv2aQkpTTlsDirKNFMXkcYsxmnaBvHbNXdpH5Mem8++w8Vs\nfrTlfQdExPza2p3ao9SkugQ7qHJr+UVEGlOpm1SXEDvVXi2/iEhjKnWT6mp3UOPVTF1EGlOpm1R4\nqINaQzN1EWlMpW5S3Rx26tBMXUQaU6mbVHeHA7dFpS4ijanUTSqiix13gJZfRKQxlbpJRXZ14A3U\nTF1EGlOpm1SPrg68gZqpi0hjKnWT6uawg62CujpfJxGRM4lK3aQcQXawVVJerkMviMh/tXpALzlz\nWQOsYFg5VlpD9+4hp/x8a79ey7MbniWlRwpDooYwJGoIZ3c7mwCL/u6LmIlK3cQC3A4KSytI4NRK\nfU/xHqb9fRp3nncnhZWFPLfxObYe3srx6uOMib2ImcNu4vLEy+v/kJyEQ8eLeOVfSym1HKDEfZCD\n5d/wTdk3RHeN5rcX/Ja0PmmnlFdEWqZSN7FAj4NjZZVA5Pd+jvLacq5adBU3JvyWIzlz2LcPir4G\n99dQXVXEykHv8PG5Wbi7zuaapBu4Y8xNJPdMavI8hw/Da++6+POW/2V32CuEHrwEb2E/aouS6GoZ\nQ0RQH47HbWD0l1cQHzyCmQnzuGzIUPr3B6t+CkXajQ69a2KhdyXx4iU5XHtZ8vd6vNfwctXfprLn\n35Ecfek5br3FQv/+EBsLZ58NvXvD0aPwwQfwZt4OVha9RFXiqwRarDhK0+ovZWl4jvdmf8/nsCTn\nMCbiRuZP/CVDE6IBqKuDwkI4cgRcLti+q4p3v3mO9bYsrIdHYtk4i8v6X8iEcQ7Gjat/TRH5r7Z2\np0rdxLreOYxHzn2en00+p82PNQyY/OSDvLfjfW7wribr4WAiTzDhNwz46iuDbQf3suP4Jr4q3cTO\n0k0cqt7LjGHXctuoW4kIjTip16+qq+L5jc+TvfktNh3+nK5lwynbcilx3su4/yfncM01FnQOcxGV\neqcSeedo7kh9hHt/dEGbHvfFF/DDB/7Bnv6/4J0r1jP+At9Oj8try/lo30fk5r/PW1+8x/HD3eiz\n+zc8d8fVXJQR6NNsIr6mk2R0IkEWB8crT36v0ooKuOsu+MG0LRwYejNrfva2zwsdICwojCv6X8FT\nVzzBgbu/4rWb74PzH2PskhRSb3yBTVtrfB1RxDRU6iYWEuCgtOrk9ipdsgQGDIDdRwpwzJrAi1Oe\nZmRMegcnbLsASwBXp1xF/l3/4r2fPkdNwpukv5LMrDv3Ul7u63QiZz6VuomFBNoprW55pl5YCK++\nCuPHw513wjMvlLP33An8bMQt/HDgD09j0razWCxc1j+DHb/N5X8uv4PFQZeRNOww77xTv7YvIs3T\nxmQmFhrooLymvtS9Xjh0CPbtg48+gqVL4csv4aKLYNo0yJzu4Zp3ppPWO43fjP6Nb4O30T0X30Zt\nYBF/ixjPXfd/yEsvhTN3Lng8UFb230tGBiQk+DqtiG+p1E3MHmTno48riYuDb76B7t0hJgbOOw8e\neAAuvBCCg+vH/mL5HVTVVfHnK/+MxYSblTxw4QMUVhay9faJnHMwl1tuCcXhgLAuBjVRa9hi/yMj\nvrqaFVkzfB1VxKe09YuJ/eq931J8NJh7L7iPqD51HKlxsa9kH1V1VVgsFgIsAViw8MmBT1j878Ws\nm7mObiHdfB37e/MaXq57+zoq6yp5Y+ob5OzM4bF1j1FSXUJE3WAOfh3O10+85OuYIu2qrd2pmbqJ\n9QwPY9HOp/jg3Rc4VH6I3mG9ie0Wi8PmwMDAMAwMDIIDg3nv2vdMXehQ/yHqK5NeYWL2RHo+2pMh\nUUO494J7mZA0gT+v+JC79zzk64giPqeZuokVVhay9fBW4rrF4ezqxBZo83Wk06LaXc3e4r2k9Exp\nuO3zPXsY/vQYvI9/7cNkIu1POx9Jp1TrriP4wTCO3F5Oz8jO8cdNOgftfCSdUpDVhq2mL5/s2O/r\nKCI+pVIXv9HVE8fGPXt9HUPEp05Y6rm5uSQnJ5OYmEhWVlaT+3NyckhNTSUtLY1zzjmH1atXd0hQ\nkROJCorj39/s8XUMEZ9qdU3d4/GQlJTEqlWriI6OZvjw4WRnZ5OS8t8PqCoqKnA4HAB88cUXXH31\n1ezatavpC2lNXTrY5fMf4khJBRvmP+LrKCLtpl3X1NevX0+/fv2IjY3FZrORmZlJTk5OozHfFjpA\neXk5PXr0aGNkkfaR1Cueg1VafpHOrdVSLygoICYmpuG60+mkoKCgybh33nmHlJQUxo8fz5NPPtn+\nKUVOwtDYOIpRqUvn1urORye7O/mkSZOYNGkSa9eu5frrr2fnzp3Njps3b17D1xkZGWRkZJx0UJET\nGZUUR3XIHgwDnWBDTCsvL4+8vLzv/fhWSz06OhqXy9Vw3eVy4XQ6Wxx/wQUX4Ha7KSoqIrKZ0+h8\nt9RF2lv/vlEYtgr2Hyrn7D5hvo4j8r38/wnvgw8+2KbHt7r8kp6eTn5+Pvv27aO2tpbFixczceLE\nRmN2797dsIi/ceNGgGYLXaSjWSwWQqri+GSHlmCk82p1pm61WlmwYAFjx47F4/Ewc+ZMUlJSWLhw\nIQCzZ8/mrbfe4q9//Ss2m42wsDAWLVp0WoKLNCfciGPj3r1cM2awr6OI+IQOEyB+Je2e2zirSwI5\nc2/3dRSRdqHDBEinFtstjn0lWn6RzkulLn4lpXc8h2pU6tJ5qdTFrwyLi+O4RaUunZdKXfzKuSlx\n1Nj34vXq8xvpnFTq4leie3TF4glm+/6jvo4i4hMqdfE7oTXxfPqVlmCkc1Kpi9/pThxbvlapS+ek\nUhe/09cex84jOq66dE4qdfE78d3j+LpUM3XpnFTq4ncGRcdzpFalLp2TSl38zjnxcZQFqtSlc1Kp\ni98ZmXwWdaEHqHW7fR1F5LRTqYvfiQgPJqC6F5v3HPB1FJHTTqUufsmhbdWlk2r1eOoiZhUZGMcX\nrr3AGF9HaaSitoJl+csoqS5hpHMkA3sOJDAg0NexxI+o1MUvRTviyC88uW3VD5YdZEfhDvKP5ZNf\nlE/+sXy6BHfhjlF3kNYn7ZSzlNWUsfSrpby5/U1W7VnFKOcoeof15vF/PU5BWQHpfdMZFT2K6YOn\nMyRqyCm/nnRuKnXxS/0i4/j40IoTjnts3WP8fu3vGdxrMAndE4kOSSQjYhT7SvYx/tWJJHRN4eaB\nd3FhzMX06WMhOLj+cYZhUFxdjMPmINga3Og5vYaXLYe2sHLPSlbtWcUnBz4hvdcFDLRM5eaq59j5\nRiSlgXDbJTD8ymMUhaxnneufjHttHOf0PYd7L7iXUc5RHfFtkU5AZz4Sv/T439fyP5/cTcnj61oc\n87u837Ng7Ss43lzNsX1OKiogIgIiI8Fuhxp3LUXRr1OU9Aeos+M9NIiw6P1YurmoDDxAsDWIak8V\nDmsXIoL70C2wDwEeO9vL1xHkjqR7ySUE7ruEo5+NIYRw0tJouNTUwKpVsHIlBAbCpZfCtTdUszP0\nL2T9M4v47vHce8G9XBx/8Wn8rsmZqK3dqVIXv7RmUwEXLUrHnXWwyX2GYXDdC/N4a/sbnLdrNY/e\n34eEBAgPh4BmNh3wGl6W5y/n66LDlBfEULA9hu2fxrDxEwcGXrr0KiK010GCIg8S0rWUAeEjGeg8\ni7POgpgYiI+HqKjmcxoG7NgBubnw1FNw9tkw9946DvbI5oG8+/n5yJ9zx7l3tPN3R8xEpS4CVFZ5\ncfzeQelvi+gSYm+4fds2gyv/dA8FjqW8lLGK665qoW19oK4OXn8dHn4YevSA2Xe6uH/fBdz3g/uY\nOWymr+OJj7S1O7WmLn7JHhpAYMVZLNv6KTZ3BO+tLuKjzwo5ELySHoM3sPe2D+nbrYevYzZis8GM\nGfCjH8Hf/w6PPBBDnfE+t5dkUFrYlV9eNs3XEcUENFMXv9Xr5hkUdVtJYE0kvbr0IKFPD1Ljopk3\n5n4iQiN8He+EDAO2boXH/7aFv1kuI37rK/w2cxwzZvg6mZxOWn4R+Y8NG8BigWHD6v81s7X71jHh\nb1dR++o/OPjpaMLDfZ1ITheVuoifen/3+0x46QZeG3aAaVO0ctpZtLU7dZgAEZO4LOEyokJiePnD\nNb6OImcwlbqIiUwZOJk1R99G/+mVlqjURUxk9gWTqTr7H/x7m9fXUeQMpVIXMZHknkl0DerO88s+\n9XUUOUOdVKnn5uaSnJxMYmIiWVlZTe7/29/+RmpqKkOGDOH8889n69at7R5UROpdEj2ZnPy3fR1D\nzlAnLHWPx8OcOXPIzc1l27ZtZGdns3379kZj4uPjWbNmDVu3buW+++7j5ptv7rDAIp3dzy+dwn7H\nW1RUaGFdmjphqa9fv55+/foRGxuLzWYjMzOTnJycRmPOPfdcwv+z4ezIkSM5cEBnnBHpKOcnDCE4\n2MJflm/xdRQ5A52w1AsKCoiJiWm47nQ6KSgoaHH8iy++yOWXX94+6USkCYvFwnDHFP664S1fR5Ez\n0AlL3dKGXfE+/PBDXnrppWbX3UWk/cw8bzJbarSuLk2dcLe06OhoXC5Xw3WXy4XT6WwybuvWrcya\nNYvc3Fy6d+/e7HPNmzev4euMjAwyMjLanlhEuC5jBDfllrBq8w4uGZrs6zjSjvLy8sjLy/vejz/h\nYQLcbjdJSUl88MEH9O3blxEjRpCdnU1KSkrDmP3793PRRRfx2muvMWpU82ds0WECRNpX8h23kRTd\nh5xf3ePrKNKB2v0wAVarlQULFjB27FgGDBjANddcQ0pKCgsXLmThwoUA/O53v6O4uJhbbrmFtLQ0\nRowY8f3fgYiclCkDJvPRES3BSGM6oJeISR087Cb6f/uy41fr6d8rFoA6Tx2bD20mPCSc/pH9fRtQ\n2oWO0ijSifSc+RNGjQwktV9PPt7/MRu+2UB893gKKwvp5ejFDwf+kGkDppEYmdjhWbyGl69LvmZH\n4Q62F26n2l3NPRdoaehUqdRFOpGb7v8XOWX3Ex80Cus351OxcxTf7OlGTZ0HW/w/cSe/QVXsmwTV\n9MFxPB1rxdlYy2MJKDsbW8VZ9AiLoE9kGFG9LPTqVX+e1qoqg6Lqoxyu3k+Rez/H6r6hpKaIMk8R\nFUYh1QFFWGxVWIPrsAa5CbS5sdhqKAv4mrDA7pztSCGxezIfHX2TVTe8z+Cowb7+NpmaSl2kE9mz\nB559FpzO+pNcf3uy69BQqKyEigooq/DwacE69pZv43DNPg7VfM3Byn0cqnJRWltMjaeaEEtXgrzh\n4LVSHlDZlo5+AAAJmUlEQVRAkMVO94CziLSeRY/gvvR0RNI7vAd9u0XijIwkCAfHi62UllgpLbFR\nUmTDXXQWpUe7cuQIHD0KX/a5i1tmW/njlQ/7+ttkaip1EWmTOk8dpTWlHK85Tp2nDmdXJ44gxyk/\n7+U3bebT2KspvG9Pm/Z3kcZ0kgwRaRNboI1IeyTx3eNJ6pHULoUOMPfGVMpLQviX65N2eT45OSp1\nEekQo0dbCHddy2MrXvd1lE5FpS4iHcJigVtGT2fZ/jdwe92+jtNpqNRFpMPc8eN+uAvP5q2Nq30d\npdNQqYtIhwkPh/Tga3l0ebavo3QaKnUR6VDzpl7Dpqp3qKip8nWUTkGlLiIdatzoPjhKh/Hwm8t8\nHaVTUKmLSIeb3P9a/vKZlmBOB5W6iHS4h66bzCH7Sv69+3iHvYZhGNrBEe1RKiKnSfw9k+hdfDU3\nj5yBYYDXa1Bn1HKoysXu8s18XbOZAs9mjli+INyIJdk6jmFdxzG4x1AiIizY7fWHPwgNhZAQg+Pe\nb/iyeANfFG5gS+FnbDm6gfNizmPJ9By/2oNVhwkQkTPSC2tz+OkH12AxrHgttXgtdQQYNuze3kR5\n04i2DiU2eCixjkEUVOWzpSKXfCOXGsroVpJBnbeGats31AYdpC7kEAG14diOnoP1aDqBh4bDoaHU\nZY5l8czHmZDkP+dJVqmLyBmruKqYwIBAggODsQXaCLCceAV497HdfLz/Y8KCwujTpQ99u/Sld1hv\nQqwhjcZ5vTBk6rsUpv2GA/duwRpwwrN1moJKXUQ6rfx8gwF/uIj7J2dy3/jZvo7TLnRALxHptBIT\nLfw86XEeWjuP41Vlvo7jE5qpi4hf8Xig9y03kJ54FsvvfMjXcU6ZZuoi0qkFBsLi2b9nxbFn+WTb\nAV/HOe1U6iLidy46J4YLHT9l8tP30tkWCLT8IiJ+qbiijF4P9afHrl8QavTAatixGnZshBJiC8ER\nHExYSP2/XewhLJwfx5m4ebu2fhER+Y8lX67mtQ3/oMpdSZW7khpPFVXuSqrd1dS4a6jx1FDrrcHt\nraPw/l2+jtsslbqIiB/RB6UiIp2YSl1ExI+o1EVE/IhKXUTEj5xUqefm5pKcnExiYiJZWVlN7t+x\nYwfnnnsuISEhPP744+0eUkRETs4JS93j8TBnzhxyc3PZtm0b2dnZbN++vdGYyMhInnrqKX796193\nWFBfy8vL83WEU2Lm/GbODsrva2bP31YnLPX169fTr18/YmNjsdlsZGZmkpOT02hMz549SU9Px2az\ndVhQXzP7D4aZ85s5Oyi/r5k9f1udsNQLCgqIiYlpuO50OikoKOjQUCIi8v2csNT96bRQIiJ+zziB\nf/3rX8bYsWMbrj/88MPG/Pnzmx07b94847HHHmv2voSEBAPQRRdddNGlDZeEhIQT1XQjJzzfU3p6\nOvn5+ezbt4++ffuyePFisrOzmx3b2q6su3admcdVEBHxJyd17Jfly5dz++234/F4mDlzJnPnzmXh\nwoUAzJ49m0OHDjF8+HBKS0sJCAigS5cubNu2jbCwsA5/AyIi8l+n7YBeIiLS8Tp8j9IT7bh0prnp\nppuIiopi8ODBDbcdO3aMSy+9lP79+3PZZZdRUlLiw4Stc7lcjBkzhoEDBzJo0CCefPJJwDzvobq6\nmpEjRzJ06FAGDBjA3LlzAfPkh/p9O9LS0pgwYQJgruyxsbEMGTKEtLQ0RowYAZgrf0lJCVOnTiUl\nJYUBAwbw6aefmib/zp07SUtLa7iEh4fz5JNPtj1/m1bg28jtdhsJCQnG3r17jdraWiM1NdXYtm1b\nR77kKVuzZo2xceNGY9CgQQ233XnnnUZWVpZhGIYxf/584+677/ZVvBM6ePCgsWnTJsMwDKOsrMzo\n37+/sW3bNlO9h4qKCsMwDKOurs4YOXKksXbtWlPlf/zxx41rr73WmDBhgmEY5vr5iY2NNYqKihrd\nZqb8N9xwg/Hiiy8ahlH/81NSUmKq/N/yeDxG7969jf3797c5f4eW+rp16xptOfPII48YjzzySEe+\nZLvYu3dvo1JPSkoyDh06ZBhGfWkmJSX5KlqbXXXVVcbKlStN+R4qKiqM9PR048svvzRNfpfLZVx8\n8cXG6tWrjSuvvNIwDHP9/MTGxhqFhYWNbjNL/pKSEiMuLq7J7WbJ/10rVqwwRo8ebRhG2/N36PKL\nv+y4dPjwYaKiogCIiori8OHDPk50cvbt28emTZsYOXKkqd6D1+tl6NChREVFNSwlmSX/L3/5Sx59\n9FECAv77q2WW7FC/X8oll1xCeno6zz//PGCe/Hv37qVnz57ceOONDBs2jFmzZlFRUWGa/N+1aNEi\npk+fDrT9+9+hpe6POy5ZLBZTvK/y8nKmTJnCE088QZcuXRrdd6a/h4CAADZv3syBAwdYs2YNH374\nYaP7z9T8S5cupVevXqSlpbW4ee+Zmv1b//znP9m0aRPLly/n6aefZu3atY3uP5Pzu91uNm7cyK23\n3srGjRtxOBzMnz+/0ZgzOf+3amtreffdd5k2bVqT+04mf4eWenR0NC6Xq+G6y+XC6XR25Et2iKio\nKA4dOgTAwYMH6dWrl48Tta6uro4pU6Zw/fXXM2nSJMB87wEgPDycK664gs8//9wU+detW8eSJUuI\ni4tj+vTprF69muuvv94U2b/Vp08foP54TldffTXr1683TX6n04nT6WT48OEATJ06lY0bN9K7d29T\n5P/W8uXLOeecc+jZsyfQ9t/dDi317+64VFtby+LFi5k4cWJHvmSHmDhxIq+88goAr7zySkNRnokM\nw2DmzJkMGDCA22+/veF2s7yHwsLChk/3q6qqWLlyJWlpaabI//DDD+Nyudi7dy+LFi3ioosu4tVX\nXzVFdoDKykrKysoAqKio4P3332fw4MGmyd+7d29iYmL46quvAFi1ahUDBw5kwoQJpsj/rezs7Ial\nF/gev7sdvN5vLFu2zOjfv7+RkJBgPPzwwx39cqcsMzPT6NOnj2Gz2Qyn02m89NJLRlFRkXHxxRcb\niYmJxqWXXmoUFxf7OmaL1q5da1gsFiM1NdUYOnSoMXToUGP58uWmeQ9bt2410tLSjNTUVGPw4MHG\nH/7wB8MwDNPk/1ZeXl7D1i9myb5nzx4jNTXVSE1NNQYOHNjw+2qW/IZhGJs3bzbS09ONIUOGGFdf\nfbVRUlJiqvzl5eVGZGSkUVpa2nBbW/Nr5yMRET+i09mJiPgRlbqIiB9RqYuI+BGVuoiIH1Gpi4j4\nEZW6iIgfUamLiPgRlbqIiB/5PypJbjVFwhbIAAAAAElFTkSuQmCC\n", + "text": [ + "" ] } ], - "prompt_number": 14 + "prompt_number": 24 }, { "cell_type": "code",