From 8aa823ba5d2b0292b0a9eeabf550219d49839f1a Mon Sep 17 00:00:00 2001 From: Gudni Karl Date: Thu, 12 Feb 2015 19:39:26 -0800 Subject: [PATCH] Fixing up ProblemMT, SurveyMT and adding Sources folder --- MT Script-3D_layerTest.ipynb | 49 ++++++++++++++++++++++-------------- simpegMT/ProblemMT.py | 24 +++++------------- simpegMT/__init__.py | 4 ++- 3 files changed, 39 insertions(+), 38 deletions(-) diff --git a/MT Script-3D_layerTest.ipynb b/MT Script-3D_layerTest.ipynb index 16066797..fd96492a 100644 --- a/MT Script-3D_layerTest.ipynb +++ b/MT Script-3D_layerTest.ipynb @@ -1,7 +1,7 @@ { "metadata": { "name": "", - "signature": "sha256:2a604d362593cf35c24bee160845949e5d870c81c4db8253338e6e041bc0e5dd" + "signature": "sha256:090d72f749f721ee8a9ec2a32ec92230d1a8d709fd7a2d443382e8716798661a" }, "nbformat": 3, "nbformat_minor": 0, @@ -21,7 +21,7 @@ "language": "python", "metadata": {}, "outputs": [], - "prompt_number": 58 + "prompt_number": 1 }, { "cell_type": "code", @@ -40,7 +40,7 @@ ] } ], - "prompt_number": 59 + "prompt_number": 2 }, { "cell_type": "code", @@ -56,13 +56,13 @@ { "metadata": {}, "output_type": "pyout", - "prompt_number": 60, + "prompt_number": 3, "text": [ "2529.536000000001" ] } ], - "prompt_number": 60 + "prompt_number": 3 }, { "cell_type": "code", @@ -74,7 +74,7 @@ "language": "python", "metadata": {}, "outputs": [], - "prompt_number": 61 + "prompt_number": 4 }, { "cell_type": "code", @@ -97,7 +97,7 @@ ] } ], - "prompt_number": 62 + "prompt_number": 5 }, { "cell_type": "code", @@ -118,9 +118,9 @@ { "metadata": {}, "output_type": "pyout", - "prompt_number": 63, + "prompt_number": 6, "text": [ - "" + "" ] }, { @@ -128,11 +128,11 @@ "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAXIAAAEKCAYAAAAPVd6lAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl4FeX1wPHvyUIEWQRU9kUptuyLCCq1pCoWtaJUqmgr\nVnFFBbeiSJVgf7XiXtdKhVYUpGgVRURBa1C0gCAg+1IIJpCAENayhITz+2MmIQlzb+5N7jbhfJ5n\nHubOvDP3HG5yMvedd2ZEVTHGGONfSfEOwBhjTNVYITfGGJ+zQm6MMT5nhdwYY3zOCrkxxvicFXJj\njPE5K+TGGONzVshNwhORa0VkoYjsFZEtIvKRiPR2150hIm+LyA8isktElorIPSJiP9vmuGE/7Cah\nici9wLPA/wGnAi2Al4D+ItIGmA9sAjqq6knAr4EzgTrxidiY2BO7stMkKhGpB+QAv1PVf3msfxOo\np6qXxTw4YxKIHZGbRHYOcALwXoD1FwDvxC4cYxKTFXKTyBoC21X1SJD1uTGMx5iEZIXcJLIdwMlB\nTlzuAJrGMB5jEpIVcpPI/gMcAgYEWP8pcGXswjEmMVkhNwlLVXcDjwAvicjlIlJLRFJF5GIRGQuM\nBs4VkSdEpBGAiPxIRN5wT5Qac1ywQm4Smqo+A9wL/AHYBnwPDAXeU9UNOCdEWwMrRGQXzsnPb4C9\ncQnYmDiw4YfGGBNhItIPeA5IBl5T1bEebZ4HLgb24wyxXVzZ97MjcmOMiSARSQZeBPoB7YFrRKRd\nuTaXAD9S1bbALcArVXlPK+TGGBNZPYH1qpqlqoeBKcDl5dr0B14HUNX5wEnF53kqwwq5McZEVjMg\nu9TrHHdZRW2aV/YNrZAbY0xkhXriUSq53TFSKrthLImInZE1xoRMVcsXybCEW3PKvd9mnJu7FWuB\nc8RNkDbN3WWV4otC7siIdwARpTqajIwMxoyp0s9bQlEdDYBIOvDzuMYSScV5ZWRkkJGREd9goqC6\n5SUSmd+p/wux3R+OXbQQaCsirYEtwNXANeXafADcCUwRkbOBXaq6tZKh+qmQG2NM7KRWcjtVLRSR\nO4FPcIYfjlfVVSJyq7v+VVX9SEQuEZH1wP+AG6oSqxVyY4zxUJXiqKozgZnllr1a7vWdVXiLMuxk\nZxylp6fHO4QoaR3vAKKiun5e1TWvqqoZ4pQIrJDHUfX9BTot3gFERXX9vKprXlWVGuKUCKxrxRhj\nPPipOPopVmOMiZlEOdoOhRVyY4zx4Kfi6KdYjTEmZuyI3BhjfM4KuTHG+FyiDC0MhRVyY4zx4Kfi\n6KdYjTEmZqxrxRhjfM5PxdFPsRpjTMzYEbkxxvicn4qjn2I1xpiY8dMRud00yxhjPETj7oci0kBE\nZovIWhGZJSInBWg3UkRWiMgyEZksImnB9muF3BhjPETp7ocPArNV9QzgM/d1Ge6ThW4GuqtqJ5yH\nUwwKtlMr5MYY4yElxClM/YHX3fnXgSs82uwBDgO1RCQFqEUFz/O0PnJjjPGQGmp1LAxrt41KPZtz\nK9CofANVzReRp4HvgQPAJ6r6abCdWiE3xhgPKQGq45dFMPdI4O1EZDbQ2GPVqNIvVFVFRD22bwPc\njfOord3A2yLyG1WdFDDWwOEYY8zxKzXZe/n5yXB+qdeP/6/selXtG2ifIrJVRBqrap6INAG2eTTr\nAXytqjvcbd4FzgUCFnLrIzfGGA8pKaFNYfoAuN6dvx6Y5tFmNXC2iNQUEQEuBFYG26kVcmOM8ZCa\nFtoUpseBviKyFufA/nEAEWkqIjMAVHUpMBFYCHznbjcu2E6ta8UYY7xEoTqqaj7OEXb55VuAS0u9\nfgJ4ItT9WiE3xhgvPqqOPgrVGGNiyEfV0UehGmNMDAUYtZKIrJAbY4wXH1VHH4VqjDExFP6IlLix\nQm6MMV58VB19FKoxxsSQj6qjj0I1xpgYspOdxhjjcz6qjj4K1RhjYshH1dFHoRpjTAz5qDr6KFRj\njIkhG35ojDE+56PqaLexNcYYL8khTmEQkV+LyAoRKRKR7kHanSQi74jIKhFZKSJnB9uvFXIPeXn3\n0b17EwDmzPkdgwZ1LFnXvv0pTJ06kDVr7qSw8GHGjbvsmO379GlFUdEjx0w33NA1ZjmUV9WcAJKS\nhAce6M3q1Xdw4MAo8vLu46WXLolJ/IFUNa+///1yz8+qsPBhGjasGbM8TAKKztOXlwEDgC8qaPcX\n4CNVbQd0BlZVFKoppU2b+tSqlcrixbmkpibRo0dT5s79vmR9zZopZGXt5v3313Dvveegeswj90p0\n6/Yqubl7S17v2XMoqrEHEqmc/vGPy+nVqzkjRsxmyZI86tRJo3Xrk2KVxjEikdewYTMZMWJ2yWsR\nYdq0q9m3r4AdOw7EJA+ToKJzP/LV4PycBSIi9YDzVPV6d5tCnGd3BmSFvJzevVsyf/5mVOGss5qx\nY8d+cnL2lKxftCiXRYtyARgypFvQfW3fvp8fftgf1XhDEYmc0tNbM2hQRzp3/iurV28vWb58udcj\nB2MjEnnt3VvA3r0FJa/btm1Ar17N+fWv345u8Cbxxe+CoNOAH0Tk70AXYBEwXFUDFhMr5K6dOx9A\nVUlLSyEpScjPH0FqajJpacnk549AFRo2DPmBHQDMnXsDtWqlsn59Pq++uog33viu4o0iKJI5XXll\nOzZs2MlFF7Vh+vRrqFEjmf/8J5v7759dpnjGQjQ+q2K33tqDvLx9TJu2OsJRG98JUB0zt0BmbuDN\nRGQ20Nhj1UOqOj3Ed+4O3Kmq34jIc8CDwCNhhnr86dz5FUSEefOGcNttM1iyJI8pU65k8uTlvP9+\neL/UW7bsZejQGSxcuIUjR5RLLmnLuHGX8aMfNWD06MzoJOAhkjm1aVOfli3r8dvfdmLIkA8oKCji\nT386n3//ezAdO75CQUFRlLI4ViTzKq1GjWSuv74Lr766iCNHAneZmePECd6L0093pmJjFpddr6p9\nq/jOOUCOqn7jvn4Hp5AHZIXclZ29h06dTiU1NZnp09dQu3YNunZtTP/+U9i+PbzukXXr8lm3Lr/k\n9eLFeSQnJ3H//ecwZsycmBWJSOaUlCSkpaUwePC0kq6Vq69+h9zc+7jkkrYxPYKNZF6lDRzYnvr1\nT2DcuEURjNb4VvS7Vjw7ylU1T0SyReQMVV2L84zPFcF2ZIUcWL78dlq2rEdKShKpqcns3v1gSeHa\nsGEYAO3avcTmzXsr2FNg8+fncOKJNTjllFps3fq/SIUeUKRzys3dh6qW6R/fvn0/27fvp2XLelHJ\nwUs0P6vbbjuTTz75L99/H/S8kjleRKE6isgA4HngZGCGiCxW1YtFpCnwN1UtfgDzXcAkEakB/Be4\nIcah+k+/fpOoUSOZCRP6M3PmeqZOXcHo0X04dKiIxx+fCziFrCq6d2/C/v2Hq3TEGI5I5/TFF5sY\nPLgLZ5zRkLVrdwDQoEFNTj65FllZu6KSg5dofVbt2p1M794tGTDgn5EO2fhVdEatvAe857F8C3Bp\nqddLgbNC3a8VciAnZw9JSULnzo245ZYP2bhxF506NSIjI5ONG8sWqZSUJDp0OAWAOnXSaNiwJl26\nNKKgoIhVq5yj1bvvPptNm3axcuUPqMIvftGGUaPO48UXv6GoKDbdKpHO6a23ljNq1HlMmNCfYcM+\n5vDhIsaOvZB163Ywc+a6mOQUjbyK3XprD7Zs2cv06WtilotJcNXlNrYiMgHnr8Q2Ve3kLssAbgJ+\ncJs9pKoz3XUjgRuBImCYqs5yl58J/APn9MFHqjrcXZ4GTMQ5Q7sDuFpVN0Uwv5B169aYQ4eKWLt2\nB3XrptGhwyl88cWxoTRrVodvv70VAFWle/cmDBjQjqysXbRp8zwAycnCY49dQIsWdTl8+Ajr1u1g\n2LCPmTBh8TH780tOBw8WcuGFb/Dcc78gM/N69u8/TGZmFhde+AaHDx/xbV4AJ5yQwnXXdeaFFxYQ\n5LIAc7zx0WGuBLugRUTOA/YBE0sV8tHAXlV9plzb9sBknK8DzYBPgbaqqiKyAGcozQIR+Qh4XlU/\nFpGhQEdVHSoiVwMDVHWQRxwKGZHIN2GojgZAZEycI4mc6pgTHM3L+IOIoKqBr7gJbR+qj4bY9hGq\n/H5VFfQSfVX9Etjpscor6MuBt1T1sKpmAeuBXiLSBKijqgvcdhOBK9z5/sDr7vy/gAvCC98YY6Ik\nLcQpAVT2Xit3ichSERkvIsXXaDfFGf9YLAfnyLz88s3uctx/s+HoZagi0qCSMRljTORE514rUVGZ\nMF4Bir90/BF4GhgSsYgCGD36aBdQeno66enp0X7LmKiOX9urY04mcWVmZpKZmRn5HSdIkQ5F2KGq\nasnNNUTkNaD4ktPNQItSTZvjHIlvdufLLy/epiWwRURSgHqqmo+HjIyMcEM1xhwHyh/YjRkToXM0\n1WXUihcRaaKqxXcaGIBzW0aAD4DJIvIMTpdJW2CBe7Jzj4j0AhYA1+EMiC/e5npgHjAQ+CzQ+44J\ncrcwPxrtnmSuTnlVx5ygeuZVHXOCo3lFRHU5IheRt4A+wMkikg2MBtJFpCugwEbgVgBVXSkiU4GV\nQCEwVI8OiRmKM/ywJs7ww4/d5eOBN0RkHc7ww2NGrBhjTFxUl0Kuqtd4LJ4QpP1jwGMeyxcBnTyW\nHwKuqjhMY4yJserctWKMMceFAHc/TERWyI0xxouPqqOPQjXGmBjyUdeKPXzZGGO8ROGCIBF5UkRW\nuRdUvus+nzNQ22QRWSwiFT5VyAq5McZ4ic6VnbOADqraBVgLjAzSdjjOKMAKx1RaITfGGC/JIU5h\nUNXZqlp8u9D5lL1YsoSINAcuAV4jwJOESrM+cmOM8RL9USs3Am8FWPcs8Hugbig7skJujDFeKnmy\nU0RmA409Vj2kqtPdNqOAAlWd7LH9L3GeAbFYRNJDeU8r5MYY4yVAdcxcBJnfBt5MVfsG262I/A6n\n2yTQbbvPBfqLyCU43wvqishEVR0cZqjGGHOcC1Ad03s5U7Ex40PfpYj0w+ky6aOqB73aqOpDwENu\n+z7A/cGKONjJTmOM8RadUSsvALWB2e7QwpcBRKSpiMwIsE2Fo1bsiNwYY7xE4YIgVW0bYPkWnOcj\nl18+B5hT0X6tkBtjjBcfVUcfhWqMMTGUIM/jDIUVcmOM8eKj6uijUI2X0S9Gbl9j7ozcvqqiOuZk\nfMhH1dFHoRpjTAz5qDr6KFRjjIkd9dFtbK2QG2OMhyIfVUcfhWqMMbFjhdwYY3zuUFqNEFsWRDWO\nUFghN8YYD0XJ/ukkt0JujDEeinz00E4r5MYY46HQCrkxxvhbkY/Ko93G1hhjPBSRHNIUDhH5o4gs\nFZElIvKZiLTwaNNCRD4XkRUislxEhlW0XyvkxhjjIRqFHHhCVbuoaldgGjDao81h4B5V7QCcDdwh\nIu2C7dQ/3x2MMSaGDhHq8MPQqereUi9rA9s92uQBee78PhFZBTQFVgXarxVyY4zxEK0+chH5E3Ad\nsB/niDtY29ZAN2B+sHbWtWKMMR4q27UiIrNFZJnHdBmAqo5S1ZbAP4BnA72/iNQG3gGGq+q+YLHa\nEbkxxngI1P+9MPN/LMzcH3A7Ve0b4ltMBj7yWiEiqcC/gDdVdVpFO7JCbowxHgKNI++aXpeu6XVL\nXo8bc0w3d0Ai0lZV17kvLwcWe7QRYDywUlWfC2W/VsiNMcZDlPrI/ywiPwaKgP8CtwOISFPgb6p6\nKdAb+C3wnYgUF/qRqvpxoJ1aITfGGA/RuERfVQcGWL4FuNSdn0uY5y+tkBtjjIeCKAw/jBYr5D5X\nHZ9JWR1zMv5j91oxxhif89O9VvwTqTHGxJCfbmNrFwSVc19eHk26dwfgd3Pm0HHQoJJ1p7Rvz8Cp\nU7lzzRoeLizksnHjPPfRoG1bfvPxx4zct4/7t23j0pdfJqVmzZjEH0hV8zqxUSMGvPkmty9bxh8K\nCvjtrFkxiz2Yqub14/79uXbGDO7dsoWR+/Zx+7Jl9LzrrpjF76WqOTXu2pXrP/+c+3JzGXXgAMOz\nsrj4+edJq1v3mLaxEonfq2InNmrEfbm5PFJURO0mTaIWc5TutRIVVshLqd+mDam1apG7eDFJqak0\n7dGD7+fOLVmfUrMmu7OymPPoo2xduhRVPWYfqSeeyODPPqOooIDx55zDO1ddRZt+/eg/fnwsUykj\nEnmlpKVxYMcO/vP002z49FPwaBNrkcirVZ8+ZH/1Ff+84gpe7tCBr598kgv+/GfO/f3vY5lKiUjk\nVHjwIIsnTOCNvn15/kc/4oMhQzj9oou44vXXY5lKiUjkVEKEX02aRM78oFesR0QhySFNicC6Vkpp\n2bs3m+fPB1WanXUW+3fsYE9OTsn63EWLyF20CIBuQ4Z47qPTtddSq2FD3r32Wgr2OVfVfnTHHVz7\n4Yd8NnIkuzdtin4i5UQir93ff8/Hw4cDTvGr06xZ9AOvQCTymnXffWVeL504kSbdu9Phqqv4+skn\noxd8AJHIafvq1Wxfvbrk9d7Nm1n48sv0Ge11o73oi0ROxfo8/DCFBw8y79ln+fFll0U17gLSorr/\nSLJCDjywcyeqSkpaGpKUxIj8fJJTU0lOS2NEfj6o8kTDhiHtq0Xv3mR//XVJEQfYMHs2euQILc49\nN6aFPJJ5JZJo53VC/fplPr9YiGZOdZs3p93AgaybOTPCUQcX6Zxap6fT7aabeLVbN07t2DGKkTsS\npdskFFbIgVc6d0ZEGDJvHjNuu428JUu4csoUlk+ezOr33w9rX3WaNGFfXl6ZZUcKCzmQn0+dKPbn\neYlkXokkmnm16tOHjoMG8c8BAyIUbWiikdONX31F465dSTnhBNZ/8gkfVHC0G2mRzOnEU09lwBtv\n8N7gwRzYsSNKEZeVKN0mobA+cmBPdjZp9eqRnJrKmunTObBzJ427dmX5lCnsyc5mT3Z2yPsK2r8X\nY5HMK5FEK69mvXpx9XvvkTl6NOs+8ryXUdREI6e3r7qKV7t1Y+qVV3JSq1YM/Oc/oxB5YJHM6VeT\nJrF04kSyPv+8zHLntiTRUURKSFMiSIwo4uj25cup17IlSSkpJKem8uDu3UhSEilpaQzbsAGAl9q1\nY+/mzSHtb19uLnVblH16U1JKCjUbNGBvbm7E4w8k0nklimjl1apPH6754AO+fOwxvho7NhqhBxSt\nnIrb71i7lr25uQz5+mtO/slPyvSfR0ukczrt/PNp1adPyUno4gI+PCuLxa+9xoyhQyOeg3Wt+Mik\nfv1IrlGD/hMmsH7mTFZMnUqf0aMpOnSIuY8/DjjFOVTZX31Fv7/8hRq1a5f0s57ety+SlET2V19F\nJQcvkc7rGHH65hGNvNpecgkDp07l8z/8gXnPhXSzuYiK+mcFJCU7RSkpJTa/8pHO6eVyfeLNevbk\n8gkTePOii/hhVcAH51SJFXIf2ZOTgyQl0ahzZz685RZ2bdxIo06dyMzIYNfGjWXaJqWkcEqHDgCk\n1alDzYYNadSlC0UFBWx3f5iWTZ7Mzx5+mF9Nnsy/R42iVsOGXPLSSyyfMoXd33/v27wAGnXpAkDN\nBg2oUacOjTp3BhG2Ll3q27zaDxzIryZN4svHHmPZ5Mmc2KgRAFpUxP7tod+eNJFy6jZkCAd37uSH\nlSspPHiQUzt25MKxY9myaBHbli/3ZU7byxXrE0891Vm+Zg3/27o1KjlUm0IuIhNw7si1TVU7ucsa\nAP8EWgFZwFWqustdNxK4EecWjcNUdZa7/Eycp2GcAHykqsPd5WnARKA7sAO4WlVjPj6vcbduFB06\nxI61a0mrW5dTOnRg0xdfHNOuTrNm3Prtt4DTF96ke3faDRjArqwsnm/TBoDD+/cz8cILufiFFxjy\nn/9QeOAAK99+m0/uvTemOUFk8wJK2hS3u3XxYlSVP8boKK9YJPPqMXQokpxMn0ceoc8jj5RsWz73\naItkTkcKCzlv1Cjqt2lDUkoKe7KzWfXuuzEfThnpn7/yon0+6pCPhh9KsP8METkP2AdMLFXInwC2\nq+oTIvIAUF9VHxSR9jhPvDgLaAZ8CrRVVRWRBcCdqrpARD4CnlfVj0VkKNBRVYeKyNXAAFUd5BGH\nZkQ07fgb7f6/j4niyZpYq445QfXMqzrmBE5eIoKqVikxEdEnNLQrfEfIC2G/n4jcBzwJnKyq+R7r\n+wHPAcnAa6oa9MRN0FErqvolsLPc4v5A8SVirwNXuPOXA2+p6mFVzQLWA71EpAlQR1UXuO0mltqm\n9L7+BVwQLB5jjImVaF2iLyItgL6AZ++DiCQDLwL9gPbANSLSLtg+KzP8sJGqFndKbQUaufNNgZxS\n7XJwjszLL9/sLsf9NxtAVQuB3W7XjTHGxFUUL9F/BhgRZH1PYL2qZqnqYWAKzoFyQFXq3HS7TWIy\nfEFLXV6cnp5Oenp6LN426kYn0LjzSKmOOUH1zKs65JSZmUlmZiYAGRkZEdtvNMaIi8jlQI6qfhdk\nDHzJAa4rB+gVbL+ViXSriDRW1Ty322Sbu3wzUHoAdXM3gM3ufPnlxdu0BLaISApQz6u/CCL7ARlj\nqo/yB3ZjxoyJyH4DdZtkZW5iU2bgMRkiMhto7LFqFDASuKh0c492Yf91rUwh/wC4Hhjr/jut1PLJ\nIvIMzl+UtsAC96h9j4j0AhYA1wHPl9vXPGAg8FnAd21VvU7KsMn9rKpTXtUxJ6ieeW2qvic7IyVQ\nIW+Rfjot0k8vef3FmLll1qtqX6/tRKQjcBqw1D0abw4sEpGeqrqtVNPyB8UtKNs9fYyKhh++BfQB\nThaRbOAR4HFgqogMwR1+6Aa/UkSmAiuBQmCoHh0SMxRn+GFNnOGHxU+DHg+8ISLrcIYfHjNixRhj\n4uFQhJ/ZqarLOXpOERHZCJzp0QuxEGgrIq2BLcDVwDXB9h20kKtqoI0vDND+MeAxj+WLgE4eyw/h\n/iEwxphEEoP7qJR8fRCRpsDfVPVSVS0UkTuBT3CGH45X1aCXrx73V3YaY4yXaF/Zqaqnl5rfgnPx\nZfHrmUDI9x22Qm6MMR6qzSX6xhhzvPLT/citkBtjjIdEudd4KPwTqTHGxJB1rRhjjM8VRHj4YTRZ\nITfGGA/WR26MMT5nfeTGGONz1kdujDE+Z4XcGGN8zvrIjTHG56yP3BhjfM6GHxpjjM/5qWulMs/s\nNMaYaq+IlJCmyhCR+0TkSKBnFIvISBFZISLLRGSyiKQF258VcmOM8VBEckhTuESkBdAX8HxenPtA\niZuB7qraCeee5EEfumOF3BhjPESrkAPPACOCrN8DHAZquc8yroXz+LeArI/cGGM8RGMcuYhcDuSo\n6ncS4HmpqpovIk8D3wMHgE9U9dNg+7VCbowxHg4RtFs6IBGZDTT2WDUKGAlcVLq5x/ZtgLuB1sBu\n4G0R+Y2qTgr0nlbIjTHGQ6Aj8v2Z37A/c2HA7VS1r9dyEekInAYsdY/GmwOLRKSnqm4r1bQH8LWq\n7nC3exc4F7BCbowx4QhUyNPSzyYt/eyS1/lj/hrS/lR1OdCo+LWIbATOVNX8ck1XAw+LSE3gIM7D\n7hcE27ed7DTGGA+FJIc0VYEWz4hIUxGZAaCqS4GJwELgO7fJuGA7siNyY4zxEO1L9FX19FLzW4BL\nS71+Angi1H1ZITfGGA9290NjjPE5K+TGGONzhwrsplnGGONrRYX+KY/+idQYY2KoqNC6Vowxxtes\nkBtjjM8VHrZCbowxvnakyD/l0T+RGmNMLFnXijHG+NxB/5RH/0RqjDGxVBjvAEJnhdwYY7xYITfG\nGJ/zUSG329gaY4yXwyFOYRCRDBHJEZHF7tQvQLuTROQdEVklIitF5GyvdsWskFfWwjzo2N2ZnzoH\n+pd6yPXA62Fj0bHTuT+PT6yhCpYTwAk14YE/w9wNsPYgzMuGYX+IfZzhCpbXlM+9P6uVe+MTa6gq\n+qx+dyd8ugJW7YMFm+Gpv0PDU2IfZ5juy8ujSXcnr9/NmUPHQUfzkuRkzv3977lj1Soe2r+fO9es\nocftt0cvmKIQp/Ao8IyqdnOnjwO0+wvwkaq2AzoDq4Lt1LpWKqNVG6hZC1YshtRU6NwDvplbtk1R\nEfRsCqUfsLp7Z2zjDEdFOSUlwT9mQK3a8OAtsGEN1G8I9U+OX8yhqCivWwY4y4slJcEH38CcQL9f\nCaCinPoPglFPw6jbYO6n0LQF/Omv8MxEuP7i+MVdgfpt2pBaqxa5ixeTlJpK0x49+H7u0bx+PmYM\n3W++mek330ze0qW0OPdcLhs3jqKCAhaPHx/5gKLXteL91OXilSL1gPNU9XoAVS3EeXZnQFbIK6NH\nb1gyH1Shy1mwcwfk5hzbLn977GOrrIpyunIwdOgOP2vjrAPYkh2fWMNRUV57dpVt/9MLoXEzeDO0\nx3fFRUU5de0Fq7+DqX93Xm/JhrfGwT1j4hNviFr27s3m+U5ezc46i/07drAn52heXa6/nq+feoo1\nH3wAwO5Nm2jWsyfnjRoVnUJ+MPK7dN0lIoNxngB0n6qW+yHkNOAHEfk70AVYBAxX1f2BdmiFPBzf\n7QQUaqSBJMF3+ZCS6rz+Lt/9xWrotE1Ohi/WO90RG9bAuKfg3x/FNXxPoeZ08ZWwdAHcdA8MuA4K\nD8NXn8HjDybmN41wPqvSfnsbLP/WmRJNqDllzoSrh0Cvn8H8L+CURnDJr+GzD+OdgacHdu5EVUlJ\nS0OSkhiRn09yairJaWmMyHfyeqJhQ5LT0ig6dKjMtoUHD3JSq1bUbd68TNGPiEoekYvIbKCxx6pR\nwCvAo+7rPwJPA0PKtUsBugN3quo3IvIc8CDwSKD3tEIejn6dna6SafPgodtg5RJ4cQq8PxlmvX+0\n3X9Xw/03wKqlzi/ZL6+C8dPhgZuOHiUlilBzatUGmrd2uoxuHwgn1oaHn4W/TYOr+sQt/IBCzau0\nUxvDBZfBw3fENtZQhZrTF7Pg0bth4idOV1FKilPEH7gpfrEH8UrnzogIQ+bNY8Ztt5G3ZAlXTpnC\n8smTWf3+0bzWz5xJz2HD2PDZZ/ywYgXNevak2403oqrUado0doV8WSYszwy4mar2DWX3IvIaMN1j\nVQ6Qo6qOsBpvAAASaElEQVTfuK/fwSnkAVkhD8eWbPhJJ+co6NPpTjFr3xWG9C/bjbJ4vjMVW7IA\n6jWA2x5IvEIeak7inhe/axDscbvrfn8jTP8G2neBlUtjH3swoeZV2lU3wsEDTmFMRKHmdOFlzh/Z\nP94DC76EJs3hoSfhyQlw93Xxiz+APdnZnNqpE8mpqayZPp0atWvTuGtXpvTvz/7tR/P6ePhwLv3r\nX7ltyRJUlb2bN/Pta6/x0wcfRI8ciXxggQp5u3RnKjYl9C4rEWmiqrnuywHAsvJtVDVPRLJF5AxV\nXQtcCKwItl8r5KGavRyatnSOblJSYflu52inRhp8ucFpc0E7yNvsvf2S+XD5tbGLNxTh5LQt1zmx\ntqfUOZd1K51/m7VKrEJemc9KBAbdDO9PggMBuyLjJ5yc7ngIpr15tJ9/7Qr43z54+wt4+hHI3hi/\nPMq5ffly6rVsSVJKCsmpqTy4ezeSlERKWhrDNjh5vdSuHXs3b+bgrl38a9Ag3k1O5sRTT2Vfbm7J\nqJWdbtuICnNoYYjGikhXnNErG4FbAUSkKfA3VS1+APNdwCQRqQH8F7gh2E6tkIdqcD9IreEc1WTO\nhA+nwt2joeAQvPy402ZbbuDtO3aHLd/HJtZQhZPTgi/g1hFQuw7sc4fmtfmx829OVsxDD6oyn1V6\nP2jWEia9Gvt4QxFOTiJOF1hpeuTougQyqV8/kmvUoP+ECayfOZMVU6fSZ/Roig4dYu7jTl77cst+\nVlpUVLKs4zXXkDVnDgfy8yMfXPhDCyukqoMDLN8CXFrq9VLgrFD3a+PIQ5Wb4xSsdp3hk/eco5qf\ndHL6HrM3OlPx17u7RzuFoVUbaNsehj/ifG1/7Zm4pnCMcHJ642U4uN8Zwta2vTNa4vG/wbxMWPVd\nPLM4Vjh5Fbv2VqcLLNFyKRZOTh+/6/y8/eo6aNEazvopjHnBOWfzfRSOXKtgT04Ou7KyaNS5M6vf\ne49dGzfSqFMn1n74Ibs2bmTXxo0l3SZNzjyT9gMHUv/002l+9tn8+u23adS5Mx8PGxad4ApDnBKA\nHZGHo0M3OHQINqyFOnWhbQfnSLW82nXgjy/BKY2dPtf1q2Dor+GTabGPuSKh5vTDVrjmfHj4Gadf\nfFc+/HsGPP5A7GMORah5ATRqCj+/BEbeEtsYwxVqTn99whnBcsdIaPqKM8Ty63/D2JGxjzkEjbt1\no+jQIXasXUta3bqc0qEDm744Nq+UtDR+9sgjNGjThqKCArLmzGHCuefyw8qV0QksesMPI05UtXIb\nimQBe3C+gBxW1Z4i0gD4J9AKyAKuKh4jKSIjgRvd9sNUdZa7/EzgH8AJOFcyDfd4L9WWlQozcW1y\n/99bJdZX3SqpjjlB9czLzWlMgnW1VNVoVUQEVa1SYiKivBRibbyj6u9XVVXpWlEg3b3MtKe77EFg\ntqqeAXzmvkZE2gNXA+2BfsDLIiU/Qa8AQ1S1LdA20L0HjDEmpnzUtVLVPvLyf4X6A6+7868DV7jz\nlwNvqephVc0C1gO9RKQJUEdVF7jtJpbaxhhj4uc4KeQKfCoiC0XkZndZI1Xd6s5vBRq5801xBrkX\nywGaeSzf7C43xpj4isLdD6OlKic7e6tqroicAswWkdWlV6qqikjlOuC9bIrcrhJKdcyrOuYE1TKv\n0ZU8R3ZciMLww2ipdCEvvjpJVX8QkfeAnsBWEWnsXpnUBNjmNt8MtCi1eXOcI/HN7nzp5Z5X1GRk\nZJTMp6enk56eXtnQjTHVSGZmJpmZmZHfcXUftSIitYBkVd0rIicCs4AxOJeS7lDVsSLyIHCSqj7o\nnuycjFPsmwGfAj9yj9rnA8OABcAM4Pny9+i1USs+UR1zguqZV3XMCWBTBEetjAyxNv45/qNWKntE\n3gh4zx14kgJMUtVZIrIQmCoiQ3CHHwKo6koRmQqsxDk9MFSP/gUZijP8sCbO8MMEvhG0Mea4kSD9\n36GoVCFX1Y1AV4/l+ThH5V7bPAY85rF8EdCpMnEYY0zUHA995MYYU60lyNDCUFghN8YYL1bIjTHG\n53zUR253PzTGGC+HQpzCJCJ3icgqEVkuImODtEsWkcUi4vUUoTLsiNwYY7xEoWtFRH6OcyuTzqp6\n2L2gMpDhOCP96lS0XzsiN8YYL9G5RP924M+qehicCyq9GolIc+AS4DWOvafVMayQG2OMl6IQp/C0\nBX4mIvNEJFNEegRo9yzweyCkh5Fa14oxxngJ1LWyPRN2ZAbcTERmA409Vo3Cqbn1VfVsETkLmAqc\nXm77XwLbVHWxiKSHEqoVcmOM8RKokJ+U7kzF1o4ps1pV+wbapYjcDrzrtvtGRI6ISENV3VGq2blA\nfxG5BOeBO3VFZGKg532Cda0YY4y36PSRTwPOBxCRM4Aa5Yo4qvqQqrZQ1dOAQcC/gxVxsEJujDHe\nojP8cAJwuogsA94CBgOISFMRmRFgmwrv3mVdK8YY4yUKww/d0SrXeSzfAlzqsXwOMKei/VohN8YY\nLz66stMKuTHGeLG7HxpjjM/ZTbOMMcbnrJAbY4zPWR+5Mcb4XCXubBgvVsiNMcaLda0YY4zPWdeK\nMcb4nA0/NMYYn7OuFWOM8Tkr5MYY43PWR26MMT5nR+TGGGPKE5EpwI/dlycBu1S1W7k2LYCJwKk4\nt7Adp6rPB9uvFXJjjIkRVR1UPC8iTwG7PJodBu5R1SUiUhtYJCKzVXVVoP1aITfGmBgTEQGuAn5e\nfp2q5gF57vw+EVkFNAWskBtjTHiierbzPGCrqv43WCMRaQ10A+YHa2eF3BhjPAU62/mFO3kTkdlA\nY49VD6nqdHf+GmBysHd3u1XeAYar6r5gba2QG2OMp0BH5Oe4U7HHyqxV1b7B9ioiKcAAoHuQNqnA\nv4A3VXVaRZFaITfGGE8HorXjC4FV7nM6j+H2n48HVqrqc6HsMCmCwRljTDVyOMQpbFcDb5VeICJN\nRWSG+7I38Fvg5yKy2J36BduhHZEbY4yn6FwRpKo3eCzbAlzqzs8lzINsK+TGGOPJP9foWyE3xhhP\n/rlG3wq5McZ48s8RuZ3srIyFedDRHTk0dQ70H1R2fdee8O5XsGY/LNgMv/8TiMQ+znAFy6tte3h5\nKny+BjYUwuPj4hNjZQTL66obYMq/4dttsHw3TP8GLr8mPnGGI1hOP7sI3vvayWnNfpizDu57FFJ8\ncNxW0e9WsbbtYNU+WF8QxWAOhDjFnw8+2QTTqg3UrAUrFkNqKnTuAd/MPbq+SXN4czZ89DaMGAKn\nnQFPTnAK+RMPxS/uilSU1wk1IScLZr8PN90LqnELNSwV5XXOz+Hj9+D/7ofd+fCLAfDMRCgshBlv\nxy/uYCrKae9ueO1ZWLsc9u11CuOfx8GJdeDRe+IXd0UqyqvYCTXhpanw1WfQJ+hgjiqyrpXqq0dv\nWDLfKWRdzoKdOyA35+j6394Oe3bBiJuc1+tXw9MPw8gn4C+PwqGD8Ym7IhXltWyRMwFcPSQ+MVZG\nRXndM7hs+9eehV594JdXJW4hryinxfOdqVhuDpydDmf3iXmoYakor2J/fAkWfOHkmH5xFAPyT9eK\nFfJQfbcTUKiRBpIE3+VDSqrz+rt894evofPD+OWsstvO+QQefRE6doNF/4lL+AGFmpffVCWvevXh\n+w0xDTcklc2pzY8hvR98/G7MQw5JOHn96jrodCb0Pwv6R7sLzI7Iq59+nZ3ukWnz4KHbYOUSeHEK\nvD8ZZr1/tN0pjeGbL8tu+0Oe8++pTWIXb6hCzctvKpvXgN9A116QMSx2sYYq3JzmZUP9k6FGDXj7\n7/DkH2IfcyhCzetHP4FRT8GgdCiIZt94Mf8ckdvJzlBtyYY69ZwjhU+nw+6d0L4rfDDFWbclO94R\nVo7ldVTf/k5f8ogbYeXS2MdckXBzurI3XNoN7rkOfvYLyPhLfOKuSCh51agBL78NT/0B1gW8m2uE\nFYY4xZ8dkYdi9nJo2tI565+S6oxuSEpyvvp96X4Fv6Ad5G2GbbnHHnmf3Mj5d1tubOOuSDh5+Ull\n8rrsanjq7/DATTAt6E3p4qMyOW3+3vl3/WooKoK/TIKxI+HA/tjHH0ioeaWkOCOn/viSM4FzFJ+U\n5IxcefpheGVshIPzzxG5FfJQDO4HqTWc0SeZM+HDqXD3aCg4BC8/7rQpLtKLvoIB15XdPr0f7P8f\nLF8c27grEk5efhJuXoNugjHPOyc+P3onPjFXpKqfVXJy2X8TRah5icBFHctue9EVcM8YuLgLbN8W\nheASY2hhKKyQhyI3x/nL364zjLwFsjfCTzrBsxnOfGlvvAKD74Sxf3NGQLRqA/c+Cv94IfFGrIST\nV0oKnNHBmT+xDtRvCO27wOGCGH7VDVE4eQ252xlR9PAdzrmNU9xvTwUFzlf8RBFOTjffC+tXwcZ1\nzonCzj3gwbEwa5ozHDGRhJNX+Z+zLj29l0eMHZFXPx26waFDsGEt1KkLbTs4Q6DKy9sM110EDz8D\nHy50hiJOfjVxTzSFmlfjZjDjW2de1Rmb/IsBztjy89rENOSQhJrXDcOcQvLYX52p2LxMuOaCmIUb\nklBzSk5x/jg1bw1Hjjif0esvwoSQ7ogae6Hm5SWq1zMkRv93KEQT4MIO9xaNzwHJwGuqOrbcetWW\ncQkteja5/++tfHDFZ6iqY05QPfOqjjkBbFJEBFWtUmIiovByiK2Hhvx+ItITeBFIxflLMVRVv/Fo\nF7Qmlhf3USsikoyTWD+gPXCNiLSLb1SxkZmZGe8QoiIzwXqQIsXyOt5EZdTKE8DDqtoNeMR9XUZl\namLcCznQE1ivqlmqehiYAlwe55hiwgq5v1hex5uoPFgiF6jnzp8EeA0JC7smJkIfeTOg9ADYHKBX\nnGIxxhhXVPrIHwTmishTOAfS53i0CbsmJkIhD62TflP8+/Kjojrmdc9oyMiIdxSRVx3z2qROTtUt\nr4io3PBDEZkNNPZYNQoYBgxT1fdE5NfABKD8w5rDLgpxP9kpImcDGaraz309EjhSunPfOfFgjDGh\niczJzsi/n4jsUdW67rwAu1S1Xrk2FdbE8hLhiHwh0FZEWgNbcB5MWuZuOFX9UIwxJhxRrDnrRaSP\nqs4BzgfWerSpsCaWF/dCrqqFInIn8AnOUJvxqppgV5gYY0xE3AK8JCJpOH03twCISFPgb6p6aWVq\nYty7VowxxlRNIgw/DEpE+onIahFZJyIPxDueiohIloh8JyKLRWSBu6yBiMwWkbUiMktETirVfqSb\n22oRuajU8jNFZJm7Lua3rRORCSKyVUSWlVoWsTxEJE1E/ukunycireKYV4aI5Lif2WIRubjUOr/k\n1UJEPheRFSKyXESGuct9/5mZEKhqwk44XyvWA61xroRaArSLd1wVxLwRaFBu2RPACHf+AeBxd769\nm1Oqm+N6jn5LWgD0dOc/AvrFOI/zgG7AsmjkAQwFXnbnrwamxDGv0cC9Hm39lFdjoKs7XxtYA7Sr\nDp+ZTRVPiX5E7teLhcqfKOkPvO7Ovw5c4c5fDrylqodVNQvnl6mXiDQB6qjqArfdxFLbxISqfgmU\nv2tUJPMova9/ATG5sUmAvODYzwz8lVeeqi5x5/cBq3DGI/v+MzMVS/RC7jUwvlmcYgmVAp+KyEIR\nudld1khVt7rzWwH3Fns0xcmpWHF+5ZdvJjHyjmQeJZ+tqhYCu0WkQZTiDsVdIrJURMaX6n7wZV7u\naIduwHyq92dmXIleyP14Jra3OvdRuBi4Q0TOK71Sne+lfsyrjOqSh+sV4DSgK84l1E/HN5zKE5Ha\nOEfLw1W1zD1rq9lnZkpJ9EK+GWhR6nULyh4tJBxVzXX//QF4D6d7aKuINAZwv7oW3wW/fH7NcfLb\n7M6XXp4Ij+mJRB45pbZp6e4rBainqvnRCz0wVd2mLuA1nM+sOEbf5CUiqThF/A1VneYurpafmSkr\n0Qt5ycB4EamBc4LlgzjHFJCI1BKROu78icBFwDKcmK93m10PFP+SfQAMEpEaInIa0BZYoKp5wB4R\n6SUiAlxXapt4ikQe73vsayDwWSwS8OIWuGIDcD4z8FFebhzjgZWqWvrG49XyMzPlxPtsa0UTThfF\nGpyTMSPjHU8FsZ6GMxJgCbC8OF6gAfApzlVcs4CTSm3zkJvbauAXpZafiVNQ1gPPxyGXt3CuKivA\n6Re9IZJ5AGnAVGAdMA9oHae8bsQ5ofcdsBSn0DXyYV4/BY64P3uL3alfdfjMbKp4sguCjDHG5xK9\na8UYY0wFrJAbY4zPWSE3xhifs0JujDE+Z4XcGGN8zgq5Mcb4nBVyY4zxOSvkxhjjc/8PdCg9jetU\nmeQAAAAASUVORK5CYII=\n", "text": [ - "" + "" ] } ], - "prompt_number": 63 + "prompt_number": 6 }, { "cell_type": "code", @@ -147,7 +147,7 @@ "language": "python", "metadata": {}, "outputs": [], - "prompt_number": 64 + "prompt_number": 7 }, { "cell_type": "code", @@ -161,7 +161,7 @@ "language": "python", "metadata": {}, "outputs": [], - "prompt_number": 65 + "prompt_number": 8 }, { "cell_type": "code", @@ -178,7 +178,7 @@ "ez_x = np.zeros((M.nEz,1),dtype=complex)\n", "# Assign the source to ex_x\n", "for i in arange(M.vnEx[0]):\n", - " for j in arange(M.vnEx[2]):\n", + " for j in arange(M.vnEx[1]):\n", " ex_x[i,j,:] = e0_1d\n", "eBG_x = np.vstack((simpeg.Utils.mkvc(M.r(ex_x,'Ex','Ex','V'),2),ey_x,ez_x))\n", "rhs_x = ABG.dot(eBG_x)" @@ -186,27 +186,38 @@ "language": "python", "metadata": {}, "outputs": [], - "prompt_number": 66 + "prompt_number": 10 }, { "cell_type": "code", "collapsed": false, "input": [ - "M.vnEy" + "for i in arange(M.vnEx[0]):\n", + " for j in arange(M.vnEx[1]):\n", + " ex_x[i,j,:] = e0_1d" ], "language": "python", "metadata": {}, + "outputs": [], + "prompt_number": 15 + }, + { + "cell_type": "code", + "collapsed": false, + "input": [], + "language": "python", + "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "pyout", - "prompt_number": 67, + "prompt_number": 17, "text": [ - "array([21, 20, 19])" + "(19,)" ] } ], - "prompt_number": 67 + "prompt_number": 17 }, { "cell_type": "code", diff --git a/simpegMT/ProblemMT.py b/simpegMT/ProblemMT.py index a6669fa7..7762086d 100644 --- a/simpegMT/ProblemMT.py +++ b/simpegMT/ProblemMT.py @@ -85,14 +85,17 @@ class MTProblem(Problem.BaseProblem): e = Ainv * rhs # is this e? Src = self.survey.getSources(freq) + # Stroe the fields F[Src, 'e'] = e - F[Src, 'b'] = self.mesh.edgeCurl * e # ??? + F[Src, 'b'] = self.mesh.edgeCurl * e return F def getA(self, freq): """ + Function to get the A matrix. + :param float freq: Frequency :rtype: scipy.sparse.csr_matrix :return: A @@ -104,6 +107,8 @@ class MTProblem(Problem.BaseProblem): return C.T*mui*C + 1j*omega(freq)*sig def getAbg(self, freq): """ + Function to get the A matrix for the background model. + :param float freq: Frequency :rtype: scipy.sparse.csr_matrix :return: A @@ -132,26 +137,9 @@ class MTProblem(Problem.BaseProblem): """ raise NotImplementedError() - getAbg(freq) - - """ Put this in MT.Sources.EldadsSource - from simpegMT.Utils import get1DEfields - # Get a 1d solution for a halfspace background - mesh1d = simpeg.Mesh.TensorMesh([M.hz],np.array([M.x0[2]])) - e0_1d = get1DEfields(mesh1d,M.r(sigBG,'CC','CC','M')[0,0,:],freq) - # Setup x (east) polarization (_x) - ex_x = np.zeros(M.vnEx,dtype=complex) - ey_x = np.zeros((M.nEy,1),dtype=complex) - ez_x = np.zeros((M.nEz,1),dtype=complex) - # Assign the source to ex_x - for i in arange(M.vnEx[0]): - for j in arange(M.vnEx[2]): - ex_x[i,j,:] = e0_1d - eBG_x = np.vstack((simpeg.Utils.mkvc(M.r(ex_x,'Ex','Ex','V'),2),ey_x,ez_x)) - rhs_x = ABG.dot(eBG_x) """ Txs = self.survey.getTransmitters(freq) diff --git a/simpegMT/__init__.py b/simpegMT/__init__.py index 1fc4b26e..e75af9a5 100644 --- a/simpegMT/__init__.py +++ b/simpegMT/__init__.py @@ -1,4 +1,6 @@ # from EM import * import Utils # import FDEM -import Base +# import Sources +# import PromblemMT +# import SurveyMT