mirror of
https://github.com/wassname/simpeg.git
synced 2026-10-29 13:40:11 +08:00
Verification for DC
This commit is contained in:
SEOGI KANG committed
1 parent
1af495b910
commit
4b2800ccef
1 file changed
+257
@@ -0,0 +1,257 @@
|
||||
{
|
||||
"metadata": {
|
||||
"name": ""
|
||||
},
|
||||
"nbformat": 3,
|
||||
"nbformat_minor": 0,
|
||||
"worksheets": [
|
||||
{
|
||||
"cells": [
|
||||
{
|
||||
"cell_type": "code",
|
||||
"collapsed": false,
|
||||
"input": [
|
||||
"from SimPEG import *\n",
|
||||
"import simpegDC as DC\n",
|
||||
"%pylab inline"
|
||||
],
|
||||
"language": "python",
|
||||
"metadata": {},
|
||||
"outputs": [
|
||||
{
|
||||
"output_type": "stream",
|
||||
"stream": "stdout",
|
||||
"text": [
|
||||
"Populating the interactive namespace from numpy and matplotlib\n"
|
||||
]
|
||||
}
|
||||
],
|
||||
"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)]"
|
||||
],
|
||||
"language": "python",
|
||||
"metadata": {},
|
||||
"outputs": [],
|
||||
"prompt_number": 30
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"collapsed": false,
|
||||
"input": [
|
||||
"mesh = Mesh.TensorMesh([hx, hy, hz], 'CCN')"
|
||||
],
|
||||
"language": "python",
|
||||
"metadata": {},
|
||||
"outputs": [],
|
||||
"prompt_number": 31
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"collapsed": false,
|
||||
"input": [
|
||||
"sighalf = 1e-2\n",
|
||||
"sigma = np.ones(mesh.nC)*sighalf"
|
||||
],
|
||||
"language": "python",
|
||||
"metadata": {},
|
||||
"outputs": [],
|
||||
"prompt_number": 32
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"collapsed": false,
|
||||
"input": [
|
||||
"xtemp = np.linspace(-150, 150, 21)\n",
|
||||
"ytemp = np.linspace(-150, 150, 21)\n",
|
||||
"xyz_rxP = Utils.ndgrid(xtemp-10., ytemp, np.r_[0.])\n",
|
||||
"xyz_rxN = Utils.ndgrid(xtemp+10., ytemp, np.r_[0.])\n",
|
||||
"xyz_rxM = Utils.ndgrid(xtemp, ytemp, np.r_[0.])"
|
||||
],
|
||||
"language": "python",
|
||||
"metadata": {},
|
||||
"outputs": [],
|
||||
"prompt_number": 45
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"collapsed": false,
|
||||
"input": [
|
||||
"fig, ax = plt.subplots(1,1, figsize = (5,5))\n",
|
||||
"mesh.plotSlice(sigma, grid=True, ax = ax)\n",
|
||||
"ax.plot(xyz_rxP[:,0],xyz_rxP[:,1], 'w.')\n",
|
||||
"ax.plot(xyz_rxN[:,0],xyz_rxN[:,1], 'r.', ms = 3)"
|
||||
],
|
||||
"language": "python",
|
||||
"metadata": {},
|
||||
"outputs": [
|
||||
{
|
||||
"metadata": {},
|
||||
"output_type": "pyout",
|
||||
"prompt_number": 46,
|
||||
"text": [
|
||||
"[<matplotlib.lines.Line2D at 0xca18f60>]"
|
||||
]
|
||||
},
|
||||
{
|
||||
"metadata": {},
|
||||
"output_type": "display_data",
|
||||
"png": "iVBORw0KGgoAAAANSUhEUgAAAVIAAAFRCAYAAAAmQSVBAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XtwVPUZ//HP5kbUYoIXQpqIyyVcMo0UG1CxjvQnQcho\nJFNLoVahYv/Q6fzqaFtaRjtxOiqOta3KMEMtRdROQX9TLo6FRqbgVKVJ60hHm7ZJa5AQQhC5GiCb\nbM7vD8KykbNnz/I9e/bC+zWzMznP9/l+n8PO+rjnnN2zAcuyLAEAzltOqncAADIdjRQADNFIAcAQ\njRQADNFIAcAQjRQADNFIkdFefPFF3XTTTZHt4cOHa/fu3anbIVyQaKRIe2+//bZmzJih4uJiXX75\n5frqV7+qv//977a5x48fVzAY9LT+ihUrVF1drcLCQn3nO98ZMtbS0qLq6mpddtllKi4u1o033qi3\n337b0/pIf3mp3gHAybFjx3Tbbbdp1apVmj9/vnp7e/WXv/xFw4YN820fysrK9Oijj+pPf/qTTp48\nec7Ya6+9FmneK1as0J133qn9+/f7tn9IPd6RIq21trYqEAjom9/8pgKBgAoLC1VTU6Oqqirb/Jyc\nHH300UeSpJMnT+rhhx9WMBhUcXGxbrrpJp06dUqS9Ne//lUzZszQiBEj9OUvf1lvvfVWzH2or6/X\nHXfcocsvv/ycsaKiIo0ZM0aBQEDhcFg5OTkqLS314F+OTMI7UqS1iRMnKjc3V4sXL9aCBQt03XXX\nacSIEa7m/uAHP9C//vUv7dy5UyUlJWpublZOTo46Ozt122236ZVXXtGcOXO0bds2ff3rX9e///1v\nXXHFFTHXc/o2dXFxsXp6evTFL35Rf/7znxP+dyKz8Y4UaW348OF6++23FQgE9N3vflcjR47UHXfc\noQMHDjjOGxgY0Jo1a/Tss8+qtLRUOTk5uv7661VQUKBXXnlFtbW1mjNnjiRp1qxZqq6u1h//+EfH\nNQOBQMyxI0eO6OjRo1qwYIG+8Y1vODZdZB8aKdLepEmTtGbNGnV0dOjDDz/Uvn379OCDDzrOOXjw\noE6dOqVx48adM/bxxx/rtdde04gRIyKPd955J+55zXjN8eKLL9by5cvV2tqqDz74IP4/DFmDRoqM\nMnHiRC1atEgffvihY94VV1yhwsJC/fe//z1nbPTo0br77rt1+PDhyOP48eP60Y9+5Lim0zvSM8Lh\nsAYGBnTxxRfHzUX2oJEirf3nP//RL37xC3V2dkqSOjo69Pvf/1433HCD47ycnBzde++9euihh9TV\n1aVwOKydO3cqFArp29/+tl5//XU1NjYqHA7r1KlT2rFjR6TG553J6e/vVzgcVm9vr8LhsCRp27Zt\n2rVrl8LhsI4dO6aHHnpIEydO1Pjx4719IpDWaKRIa8OHD1dTU5Ouu+46feELX9ANN9yga665Rs88\n84yk0+8So98pRv/985//XFVVVZo2bZouv/xy/eQnP9HAwIDKy8u1adMmPfHEExo5cqRGjx6tZ555\nRgMDA7b78LOf/UwXX3yxnnrqKb3yyiu66KKL9Pjjj0s6fW504cKFKi4u1sSJE/XJJ59o8+bNSXxG\nkI4C3NgZAMzwjhQADNFIAcAQjRQADNFIAcBQ1n1F1M1n/QAgUTfffLN27NhhO5Z1V+1PN9IGHyo1\neFzH7XrbJX3NYH4i+U45dmNuYqbbXq555rlMZo1kr+kUjzeWSI6b/FivTZNabiRjTfs6sdolh/YA\nYIhGCgCGaKQZJ5jqHcgiwVTvQJYJpnoHUoZGmnHGpHoHsgjPpbcu3OeTRgoAhmikAGCIRgoAhmik\nAGCIRgoAhmikAGCIRgoAhmikAGCIRgoAhmikAGCIRgoAhmikAGCIRgoAhmikAGCIRgoAhmikAGCI\nRgoAhmikAGCIRgoAhmikAGAoYMX6xfsMFQgEUr0LALJUrHaZ5/N++KTBpxpe1jFdL9H5bvKdcuzG\n3MRMt5Oxph81krWmUzzeWCI5JvlezfVzzVh17HFoDwCGaKQAYIhGCgCGaKQAYIhGCgCGaKQAYIhG\nCgCGaKQAYIhGCgCGaKQAYIhGCgCGaKQAYIhGCgCGaKQAYIhGCgCGaKQAYIhGCgCGaKQAYIjfbAIA\nl/jNpqTU8LKO6XqJzneT75RjN+YmZrqdjDX9qJGsNZ3i8cYSyTHJ92qun2vGqmOPQ3sAMEQjBQBD\nNFIAMJTSRhoMBnXNNddo6tSpmj59uiTp0KFDqqmp0YQJEzR79mwdOXIkkv/kk0+qoqJCkyZNUmNj\nY6p2GwCGSGkjDQQC2rFjh95//301NzdLkpYvX66amhq1trbqlltu0fLlyyVJLS0tWr9+vVpaWrR1\n61Y98MADGhgYSOXuA4CkNDi0//zHCTZv3qxFixZJkhYtWqSNGzdKkjZt2qSFCxcqPz9fwWBQ48eP\njzRfAEillL8jnTVrlqqrq/XCCy9Ikrq7u1VSUiJJKikpUXd3tyRp3759Ki8vj8wtLy9XZ2en/zsN\nAJ+T0s+RvvPOOyotLdUnn3yimpoaTZo0ach4IBBw/IA9H74HkA5S2khLS0slSVdeeaXq6+vV3Nys\nkpIS7d+/X6NGjVJXV5dGjhwpSSorK1NHR0dk7t69e1VWVhZj5e1RfwcljUnK/gPIZu2SdrvKTNmh\n/YkTJ3T8+HFJUk9PjxobG1VVVaW6ujqtXbtWkrR27VrNmzdPklRXV6d169YpFAqpvb1dbW1tkSv9\n5/pa1IMmCuB8jNHQXhJbyt6Rdnd3q76+XpLU39+vu+66S7Nnz1Z1dbXmz5+v1atXKxgM6tVXX5Uk\nVVZWav78+aqsrFReXp5WrlzJoT2AtJCyRjpmzBjt2rXrnPhll12mbdu22c5ZtmyZli1bluxdA4CE\npPzjTwCQ6WikAGCI+5ECgEvcjzQpNbysY7peovPd5Dvl2I25iZluJ2NNP2oka02neLyxRHJM8r2a\n6+easerY49AeAAzRSAHAEI0UAAzRSAHAEI0UAAzRSAHAEI0UAAzRSAHAEI0UAAzRSAHAEN+1BwCX\n+K59Ump4Wcd0vUTnu8l3yrEbcxMz3U7Gmn7USNaaTvF4Y4nkmOR7NdfPNWPVscehPQAYopECgCEa\nKQAYopECgCEaKQAYopECgCEaKQAY4gP5AOASH8hPSg0v65iul+h8N/lOOXZjbmKm28lY048ayVrT\nKR5vLJEck3yv5vq5Zqw69ji0BwBDNFIAMEQjBQBDNFIAMEQjBQBDfPwJAFzi409JqeFlHdP1Ep3v\nJt8px27MTcx0Oxlr+lEjWWs6xeONJZJjku/VXD/XjFXHHof2AGCIRgoAhmikAGCIi00A4BIXm5JS\nw8s6puslOt9NvlOO3ZibmOl2Mtb0o0ay1nSKxxtLJMck36u5fq4Zq449Du0BwBCNFAAM0UgBwBAX\nmwDAJS42JaWGl3VM10t0vpt8pxy7MTcx0+1krOlHjWSt6RSPN5ZIjkm+V3P9XDNWHXsc2gOAIRop\nABjiHCkAuMQ50qTU8LKO6XqJzneT75RjN+YmZrqdjDX9qJGsNZ3i8cYSyTHJ92qun2vGqmOPQ3sA\nMEQjBQBDnCMFAJc4R5qUGl7WMV0v0flu8p1y7MbcxEy3k7GmHzWStaZTPN5YIjkm+V7N9XPNWHXs\ncWgPAIZopABgiHOkAOAS50iTUsPLOqbrJTrfTb5Tjt2Ym5jz9qpVpZowYbbGjh2h4tr/o/7wCfX3\nFWjYsDwd3fLnGLEb1R8uUH9feDC2QcW19eoPW4Oxfh3d8u7gXEubl/5f1T/9vAoLc9XbG1Z/+IQ2\nL/296p9e4hA7of5wgTYv/e1gTOrt1eB6v1X903epsPCSwblna4QHLLW1fapweEA3am9CzwXnSFO5\nZqw69rK0kSJTTZgwQTNnBk9vXDpM0jCF+gZUkJ+jotFFMWKXSlJUbPRgnmzmSrW1tSr69a8lScOG\n5UkaptraChX9ephD7NLBudExRcWKouYOrTFu7Ai1tn2ajKcLaYJGirRy4sQJSdKRI6dU/O67am5u\n1tGjl6umZlycWKeOHj01GDtiE4vK+8c/VPPuu6djxYWn1/vHfhexzqjYERUXF38uFp13ukZzc6dm\nz35ZR9Urf941IRU4R4q0smbNGtXW1ioUCqmkvl6f9fSov79fRUVFOrBhgyexLY88ojt/+UtZlqWc\nQCCpsZ6eHrW2tiocDuvGVD+5MBarXWZpI23woVKDx3VM10t0vpt8pxy7MTcx5+3t22/WzJkzT2/M\nmCHp7OG53n3Xk9iBjRs1ct68IXt5YOObGjmvxsPY2RoHD55Qa9unnCM1nuvnmvZ1YrVLPv6EtDJ2\n7FhJUjg8IEnq7++PvHidYwNRsbBN7GxeX1+fJOnMfxP9/f3q6wu7iA1ExSybWHReX+Tv3R8fNXxW\nkO5opEgre/bskSTl5p5+aebl5Sln8HSNcywnKpZrEzubV1BQIEk6cxYoLy9PBQW5LmI5UbGATSw6\nryDyd/DqIsNnBWnPyjKSeGTw44033rAsy7IOHz5sWZZlNTU1WY2NjRkba2pqsoqKilL+vPLw5hFL\nll61b/Cphpd1TNdLdL6bfKccuzE3MeftAwcO6MCBHoVCAV1Sc70qeo6pv79AoaJCfbbhTzFi01TR\nE1Z//8Bg7P/pkpo5qugJDcby9dmGbYNzQ9ryyA8Umna9LrIG1BfIUUXPMW15ZH2cmKWKnrC2PLJ2\nMNavvkDe4HprFZo2TRdZgcG5Z2tU9IT0p8+KFNZwzpEaz/VzzVh17HFoj7QSDAY1cuQlKi8vUn5+\njkYUF6uoqFAF+TkOsXyNKC6MipUP5p2J5UfNLdSsWbNUkJ+jYQV5kfVmzRobJ5Y/OPdMrCLine truncated
|
||||
"text": [
|
||||
"<matplotlib.figure.Figure at 0xca186d8>"
|
||||
]
|
||||
}
|
||||
],
|
||||
"prompt_number": 46
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"collapsed": false,
|
||||
"input": [
|
||||
"rx = DC.DipoleRx(xyz_rxP, xyz_rxN)\n",
|
||||
"tx = DC.DipoleTx([-200, 0, -12.5],[+200, 0, -12.5], [rx])"
|
||||
],
|
||||
"language": "python",
|
||||
"metadata": {},
|
||||
"outputs": [],
|
||||
"prompt_number": 47
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"collapsed": false,
|
||||
"input": [
|
||||
"survey = DC.SurveyDC([tx])\n",
|
||||
"problem = DC.ProblemDC(mesh)\n",
|
||||
"problem.pair(survey)"
|
||||
],
|
||||
"language": "python",
|
||||
"metadata": {},
|
||||
"outputs": [],
|
||||
"prompt_number": 48
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"collapsed": false,
|
||||
"input": [
|
||||
"data = survey.dpred(sigma)"
|
||||
],
|
||||
"language": "python",
|
||||
"metadata": {},
|
||||
"outputs": [],
|
||||
"prompt_number": 49
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"collapsed": false,
|
||||
"input": [
|
||||
"a = np.random.randn(3)"
|
||||
],
|
||||
"language": "python",
|
||||
"metadata": {},
|
||||
"outputs": [],
|
||||
"prompt_number": 50
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"collapsed": false,
|
||||
"input": [
|
||||
"print (a.reshape([1,-1])).repeat(3, axis = 0)\n",
|
||||
"print (a.reshape([1,-1])).repeat(3, axis = 0).sum(axis=1)"
|
||||
],
|
||||
"language": "python",
|
||||
"metadata": {},
|
||||
"outputs": [
|
||||
{
|
||||
"output_type": "stream",
|
||||
"stream": "stdout",
|
||||
"text": [
|
||||
"[[ 0.46013066 1.09015012 1.85863975]\n",
|
||||
" [ 0.46013066 1.09015012 1.85863975]\n",
|
||||
" [ 0.46013066 1.09015012 1.85863975]]\n",
|
||||
"[ 3.40892053 3.40892053 3.40892053]\n"
|
||||
]
|
||||
}
|
||||
],
|
||||
"prompt_number": 51
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"collapsed": false,
|
||||
"input": [
|
||||
"def DChalf(txlocP, txlocN, rxloc, sigma, I=1.):\n",
|
||||
" rp = (txlocP.reshape([1,-1])).repeat(rxloc.shape[0], axis = 0)\n",
|
||||
" rn = (txlocN.reshape([1,-1])).repeat(rxloc.shape[0], axis = 0)\n",
|
||||
" rP = np.sqrt(((rxloc-rp)**2).sum(axis=1))\n",
|
||||
" rN = np.sqrt(((rxloc-rn)**2).sum(axis=1))\n",
|
||||
" return I/(sigma*2.*np.pi)*(1/rP-1/rN)"
|
||||
],
|
||||
"language": "python",
|
||||
"metadata": {},
|
||||
"outputs": [],
|
||||
"prompt_number": 52
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"collapsed": false,
|
||||
"input": [
|
||||
"data_analP = DChalf(np.r_[-200, 0, 0.],np.r_[+200, 0, 0.], xyz_rxP, sighalf)\n",
|
||||
"data_analN = DChalf(np.r_[-200, 0, 0.],np.r_[+200, 0, 0.], xyz_rxN, sighalf)\n",
|
||||
"data_anal = data_analP-data_analN"
|
||||
],
|
||||
"language": "python",
|
||||
"metadata": {},
|
||||
"outputs": [],
|
||||
"prompt_number": 53
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"collapsed": false,
|
||||
"input": [
|
||||
"Data_anal = data_anal.reshape((21, 21), order = 'F')\n",
|
||||
"Data = data.reshape((21, 21), order = 'F')\n",
|
||||
"X = xyz_rxM[:,0].reshape((21, 21), order = 'F')\n",
|
||||
"Y = xyz_rxM[:,1].reshape((21, 21), order = 'F')"
|
||||
],
|
||||
"language": "python",
|
||||
"metadata": {},
|
||||
"outputs": [],
|
||||
"prompt_number": 54
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"collapsed": false,
|
||||
"input": [
|
||||
"fig, ax = plt.subplots(1,2, figsize = (12, 5))\n",
|
||||
"vmin = np.r_[data, data_anal].min()\n",
|
||||
"vmax = np.r_[data, data_anal].max()\n",
|
||||
"dat0 = ax[0].contourf(X, Y, Data, 60, vmin = vmin, vmax = vmax)\n",
|
||||
"dat1 = ax[1].contourf(X, Y, Data_anal, 60, vmin = vmin, vmax = vmax)\n",
|
||||
"cb0 = plt.colorbar(dat1, orientation = 'horizontal', ax = ax[0])\n",
|
||||
"cb1 = plt.colorbar(dat1, orientation = 'horizontal', ax = ax[1])"
|
||||
],
|
||||
"language": "python",
|
||||
"metadata": {},
|
||||
"outputs": [
|
||||
{
|
||||
"metadata": {},
|
||||
"output_type": "display_data",
|
||||
"png": "iVBORw0KGgoAAAANSUhEUgAAAtIAAAEdCAYAAAAhPIBWAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztvX9wHdV9//2+sZWhxBIS/iEbXaUythwhW8ji6xjSZ0id\ngik0XxxoUxfTBhIg04G0SZ8wCc1kkjqdCdj9lmkpjSdMxpk4mbaQdGLMt4k9hm8e8yMtGDB+mid2\nigA5kWXZxJYFMiY1dvb542qv9u49Z/fs2bO7Z3ffrxmNpN295549e/Z93/dzPudsxXEcB4QQQggh\nhJBIvCvrChBCCCGEEJJHaKQJIYQQQgjRgEaaEEIIIYQQDWikCSGEEEII0YBGmhBCCCGEEA1opAkh\nhBBCCNEgtpG+7bbb0NnZiYGBgfq2jRs3olqtYmhoCENDQ9i5c2d933333Yfe3l709fVh9+7dcd+e\nEEJIBKjZhBBijkrcdaSffvppzJkzB7fccgt+8pOfAAC+8pWvoLW1FZ/97Gcbjj1w4ABuvvlmPP/8\n8xgbG8PVV1+Nl19+Ge96FwPjhBCSBtRsQggxR2w1vPLKK9HR0dG0XeTPd+zYgQ0bNqClpQU9PT1Y\nunQp9u7dG7cKhBBCFKFmE0KIORILKzz44IMYHBzE7bffjsnJSQDAkSNHUK1W68dUq1WMjY0lVQVC\nCCGKULMJISQ6s5Mo9M4778SXv/xlAMCXvvQl3H333di6davw2EqlIth2IYCTSVSNEEISZcmSJXjl\nlVeyrkYkqNmEkLISV7MTMdILFiyo/33HHXfg+uuvBwB0dXVhdHS0vu/w4cPo6uoSlHASwEYTNTFQ\nRlz+DcD/zLoSBijKeQA8Fxux4TxeN1LKq69uNFJOmtij2dIaJli2Hxv6oimKci5FOQ+A56KDGW2W\nEVezE0ntGB8fr/+9ffv2+uzwdevW4eGHH8aZM2cwMjKC4eFhrF692uA7L/D9EEKIKuXVj+w0O4hy\nXgtCiB+7tTl2RHrDhg148skncfz4cXR3d+MrX/kK9uzZg/3796NSqWDx4sV46KGHAAD9/f1Yv349\n+vv7MXv2bGzZskU4TBgN+xqVEFIEvNqSbEQkTbLXbBnUckKICn6tyFafYy9/lwQ1od4YcESeBPdl\nAMuyroQBinIeAM/FRvJ0HmGivVG4AkaRCddsGTZqeZ76YhhFOZeinAfAc0kaHVMdT7MTyZFOBhsF\nVwXbOpkuRTkPgOdiI3k6j2JGqtPDdi3PU18MoyjnUpTzAHguSePqS3rabLmRtl1wCSHlhqZaHeo5\nISQt0tNmi400RZcQkieoWc2wTQghWZNslNpiI00IISSf0EATQmwjmSg1jTQhhBAD0DwTQvKCuSh1\nYo8IJ4QQQgghxF7iBwBopAkhhBBCCNGAqR2RuTDrCgiYyLoCpNCwzxNCCCEiaKStNAlRMXEONCbF\npAj9W0TU82L/JoA99wP7I0kaG/p6Ofp5iYy0DZ3KZlTapxw3RX5gn1ZH1lbs0/knj/eBbp3ZX8tH\nHvu3SzkCHgU00nnudLZTjpsiO9h304cGOz/w/ojWBuzD9sK+LCafelwAI80OaS8618buG0Yd9st8\n479+RemXeYD3jhlU25F92wzst8khalt7+m1OjTQ7bHHhtSU24u2X9gh4/uH9nj2McgfDPmon9pjr\nHBlpdmZCiA3QVMeDWp5f8jzKyH5XfLLRZsuNNDs+IcRmaKrVoJaXF157kgXppedZbKR58xFC8gQ1\nawa2BSHEJpIz1hYbaUIIIfmB5pkQkhfM6RUfEU4IIYQQQogGNNKEEEIIIYRoUPLUjk7D5R0zXB4h\nZYb3JyGEELspkJE2/aGrg0od+GFOSDb3q8578n7Nniz6Cq87IXKy8lt23pc5NdI2mGZdwupuZ0ch\nJBp5vke9iM6D96g5bO0nUerF/kCKgq33o4udwcqcGGnbL65Jgs6Vgk1sokz3pRf/efO+VKOo/UV2\nXuwXxFaKei8CWeizxUa6yBdaF5pskja8D8NhG4kpe7vQYJOsKfs9CKRhrGOv2nHbbbehs7MTAwMD\n9W0TExNYu3Ytli1bhmuuuQaTk5P1fffddx96e3vR19eH3bt3x317UqdT8kOICuw/ZSE5zWbfUYP3\nGTEN9Vsd8+0U20h/4hOfwK5duxq2bdq0CWvXrsXLL7+Mq666Cps2bQIAHDhwAI888ggOHDiAXbt2\n4a677sKvf/3ruFUggfAGI17YH8oONdtGeE8SVdhXbCO2kb7yyivR0dHRsO2xxx7DrbfeCgC49dZb\n8eijjwIAduzYgQ0bNqClpQU9PT1YunQp9u7dG7cKRAvejMWH15g0Q83OC7x/CftAHkgkR/rYsWPo\n7Kxd8M7OThw7VstJOXLkCK644or6cdVqFWNjYzHfLc5jHs09a70YcBJV/qHQyqFWyEhXs4k+XEWm\n2FC/80jikw0rlQoqlUrgfjHf8vx9JYD3G6yVS9xnrRf7w7XxpqZY20mZhDfu/ZrUe/8IwP60KpI4\nZjR75fSPiCyuY5G1mgGQ/EL9ThbZfb8fJjU7ESPd2dmJo0ePYuHChRgfH8eCBQsAAF1dXRgdHa0f\nd/jwYXR1dUlK+WwSVTNMUMcomnBTrO2g6MKbpVnW5Xemf1y2ZVURbcxo9senf9t4DcPqVCS9plbb\nS1H128Z7HpDX63fQ+EU/nmbHzpEWsW7dOmzbVqvYtm3bcMMNN9S3P/zwwzhz5gxGRkYwPDyM1atX\nJ1EFC7jQ91M0mLuVPGVo46LfJ/nAjGbn+Rr6+2Gez8VP0TXEZorY9kW5T8zVP3ZEesOGDXjyySdx\n/PhxdHd346//+q/xl3/5l1i/fj22bt2Knp4efPe73wUA9Pf3Y/369ejv78fs2bOxZcuWwCHEYuG/\nWEWKgLgwf0+fooisCnkV3mJAzY6CqK/mXbsZsU6GImo4tVqFiuM4TtaV8FMT6v8362qkSN6FWYey\nincRxTaMsonxICyU1UQpn2Z7KZp+l1Wbo1I0LS+bTnuJp9kWP9lQRpzOa6tAFDHqEUYZotdFE1oV\n8iTGRdQSkj5FG20sgzbrUDQ9z5NW243FRjqJThtWpk1iUTRxViHvQ45FE1oVbBbjJK9HGa+1Dmm1\nk01aUUTtzrs2R6Wo93dZ9dqL+b5rcWrH0ayr4cNG4SiCQEfBtmtQVLGVYaMI23gNFpY0tcM2zXax\nSTeKrNk2tXMUbNQQU9im2ba2dTzNppE2gi0CUmSRJuljiwjbKr4yaKTzQda6XRa9zrqdveRNS3TI\nWrfz2MbxNNvi1I48Ycuwl/cGKotIE7NkLcJAPoWY5I+sdbsseq16P8dp/7JrRpa6Xfa2tzoibbJa\nWRlbW76JF1mkSXxsMM9AdoJs+n0rJY1IJ3HOWWpoVu9NvSaqZKXdRdFql3iaXRIjLSJtkaSpJjZR\nVvOcxvvRSKdLGtpKU01sIQvtTkM3s4xs00gbpMiCLIIiXS5sMc9AsYWZRjpbiqjj1OpyU0TzbFNK\nCI10giQtljTVJGnKZJ5tEWYaabsomo5Tq8tD2vpd1iVDaaSnX+T5O9EzSkI0aaiJaYpuoBMU5dha\nUlIjXTF4zrnT8CTLDYJ6XTyKYp4TNs6V8EOUcYpqpEsryjTVRJcim+eUjLMJYopyHjGu2TISeYsi\nBEeo1fmG5lmIaW2WQSOtAY11RCjSdmODgS6xcfZDI50+xt467/pNrc4XRTDQhspMyzSLKKyRnjdd\nrRMpvanxVihClMMPRdo+irZ+aE4jGnM9fx8vqZGeF/OcTWm9kaanoSZJkncDbaA8U9o8N/yQUGJq\ntv1GOg5xhdlKQU663DAo0HZQtPVDLRBmE4IM0EibpJAanqZ2U6/touQG2haN9kMjHYFMRTlJ8eQk\nl/JQtGWQYpRdUFHOI5VKBbhE4Zx/aegNdbXcKg2noS4PeTbQGZlnE/o8X/G4g0U10iqiLENVrDMx\n1jTURAca6NSMs6r4yogpynkktma76BptHS23xlQz7aO4lNRA62i1rnGOq9cAjXQoUYQ5VWNdtLVN\nAQq0aYpmniOWn2bE2YQYu5TVSH9AcM6TMQuOaqxTj1Sb0lnqdbHIq4HOgXmOo9Xtku3/EU+zZ2u/\nMmlkJ+yiKtBuo6sIsvei6giy25GUrofbYZMS0KTLF+EVD4q0HkXLfdYoP46BjirMUUU5TJdII3F1\n3Ht9ktTwSNrtxZTOZqnX1Gpz0EArkVSQIyN9ttdIhxHUYCJxTkuQARpqADTVUSjayhua5aeVSxcm\nzDTLeszz/H1c8TX+tg4y1v7rFqbjbr+goQ6BhjoeeR45TNFAJ2GedbV6XvghUbDXSItOVEecTZtq\nQF2YtQw1kIyIZmWoAQq1jKIa6IhlJynILkmYZ8NiXChkbROm4brGOkjHdYIiNNQklLyPHqYwTyWJ\neSlRtToFnbbXSIsIahCZQJs21UB0YY4sykmaahsMNVBesc764SkWLV8XRZSzMs80y+qoRKSjGuww\n/XZJylSX2lAD5dXpIGigQzFpoFWNc1StNqjt9hpp/0mGRTJUjncviEyQ45hqVVHOTJDTKleVskU/\naKAbUBXmtA10hmJcKKIaZhUTHtVUm0z9iGWoTWisqXKiciHKo9FhlNBAA8kEPEyNEKroLyPSHnSF\nOchQA+Gm2rSh1jLTQLKGOssnJhbdUGdtoAHrnkCoIsxp5z2Hia2qGJc5v1p07n59VdFxU6Y6iqGO\nkvKRaXSaqXnZUEITbauBlmlIhhqdqJHu6elBW1sbZs2ahZaWFuzduxcTExP4oz/6I/z85z9HT08P\nvvvd76K9XXBmKqIMhOdSBxlq7/uYilKriLK2mQaSMdRZm2mgeGJtg4EGcvlkqzSHBUX6UWLine truncated
|
||||
"text": [
|
||||
"<matplotlib.figure.Figure at 0xc854160>"
|
||||
]
|
||||
}
|
||||
],
|
||||
"prompt_number": 55
|
||||
}
|
||||
],
|
||||
"metadata": {}
|
||||
}
|
||||
]
|
||||
}
|
||||
Reference in new issue
Block a user