mirror of
https://github.com/wassname/simpeg.git
synced 2026-10-28 13:30:34 +08:00
update the inverse notebook
This commit is contained in:
rowanc1 committed
1 parent
b55e960e0d
commit
f60a415b59
1 file changed
+226
-50
+226
-50
@@ -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": [
|
||||
"[<matplotlib.lines.Line2D at 0x109dd95d0>]"
|
||||
]
|
||||
},
|
||||
{
|
||||
"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+zYse7B57XXXkN9fb3fvnDAISJSsGAc4Wg0GknbEkIE1LbNmzdjw4Line truncated
|
||||
"text": [
|
||||
"<matplotlib.figure.Figure at 0x1095a2d90>"
|
||||
]
|
||||
}
|
||||
],
|
||||
"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": [
|
||||
"[<matplotlib.lines.Line2D at 0x102073390>,\n",
|
||||
" <matplotlib.lines.Line2D at 0x102073610>]"
|
||||
]
|
||||
},
|
||||
{
|
||||
"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+6IiOHLine truncated
|
||||
"text": [
|
||||
"<matplotlib.figure.Figure at 0x109dc9d90>"
|
||||
]
|
||||
}
|
||||
],
|
||||
"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": [
|
||||
"[<matplotlib.lines.Line2D at 0xab90c50>]"
|
||||
"<matplotlib.colorbar.Colorbar instance at 0x10bbbbd40>"
|
||||
]
|
||||
},
|
||||
{
|
||||
"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+fmYO3cuXC4XMjIyEBsbiy1btgAAMjMzkZKSgvLychiNRoSGhqKwsNBnWQLine truncated
|
||||
"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+1YW1a9dy7733EhLine truncated
|
||||
"text": [
|
||||
"<matplotlib.figure.Figure at 0xa5a5550>"
|
||||
"<matplotlib.figure.Figure at 0x10cd55b90>"
|
||||
]
|
||||
}
|
||||
],
|
||||
"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<ipython-input-14-b30b5a6cd927>\u001b[0m in \u001b[0;36m<module>\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": [
|
||||
"[<matplotlib.lines.Line2D at 0x10b6164d0>,\n",
|
||||
" <matplotlib.lines.Line2D at 0x10b616750>]"
|
||||
]
|
||||
},
|
||||
{
|
||||
"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\nA6nUTSjAbafwuJZgRKQplboJBXgcFJXqw1IRaUqlbkKBXjvHyjRTF5GmVOomZDUcFFeo1ELine truncated
|
||||
"text": [
|
||||
"<matplotlib.figure.Figure at 0x10c3c3fd0>"
|
||||
]
|
||||
}
|
||||
],
|
||||
"prompt_number": 14
|
||||
"prompt_number": 24
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
|
||||
Reference in new issue
Block a user