diff --git a/notebooks/DCinverse.ipynb b/notebooks/DCinverse.ipynb new file mode 100644 index 00000000..2cf8b20b --- /dev/null +++ b/notebooks/DCinverse.ipynb @@ -0,0 +1,185 @@ +{ + "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": 2 + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "mesh = Mesh.TensorMesh([hx, hy, hz], 'CCN')" + ], + "language": "python", + "metadata": {}, + "outputs": [], + "prompt_number": 3 + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "sighalf = 1e-2\n", + "sigma = np.ones(mesh.nC)*sighalf\n", + "p0 = np.r_[-50., 50., -50.]\n", + "p1 = np.r_[ 50.,-50., -150.]\n", + "blk_ind = Utils.ModelBuilder.getIndecesBlock(p0, p1, mesh.gridCC)\n", + "sigma[blk_ind] = 1e-3" + ], + "language": "python", + "metadata": {}, + "outputs": [], + "prompt_number": 4 + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "nElecs = 31 \n", + "x_temp = np.linspace(-250, 250, nElecs)\n", + "aSpacing = x_temp[1]-x_temp[0]\n", + "y_temp = 0.\n", + "xyz = Utils.ndgrid(x_temp, np.r_[y_temp], np.r_[0.])" + ], + "language": "python", + "metadata": {}, + "outputs": [], + "prompt_number": 10 + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "txList = DC.Examples.WennerArray.getTxList(nElecs,aSpacing)" + ], + "language": "python", + "metadata": {}, + "outputs": [], + "prompt_number": 12 + }, + { + "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)" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "metadata": {}, + "output_type": "pyout", + "prompt_number": 8, + "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", + "text": [ + "" + ] + } + ], + "prompt_number": 8 + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "survey = DC.SurveyDC(txList)\n", + "problem = DC.ProblemDC(mesh)\n", + "problem.pair(survey)" + ], + "language": "python", + "metadata": {}, + "outputs": [], + "prompt_number": 13 + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "datasyn = survey.dpred(sigma)" + ], + "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: " + ] + } + ], + "prompt_number": 14 + }, + { + "cell_type": "code", + "collapsed": false, + "input": [], + "language": "python", + "metadata": {}, + "outputs": [] + } + ], + "metadata": {} + } + ] +} \ No newline at end of file