diff --git a/.gitignore b/.gitignore index 2c242ac7..1d24ca05 100644 --- a/.gitignore +++ b/.gitignore @@ -280,3 +280,11 @@ notebooks/MT3D_dobs.npy notebooks/MT3D_dtrue.npy notebooks/trueModel.vtr notebooks/018-InversionModel-2015-11-13-10-29.npy +notebooks/scipy2015/SeogiModelMT3Dfor21Dinv/xyzmod* +notebooks/scipy2015/SeogiModelMT3Dfor21Dinv/seogiModel_dobs.npy +notebooks/MT Script-3D_InversionTest-Copy1.ipynb +notebooks/scipy2015/SeogiModelMT3Dfor21Dinv/model3D_zdet1Dinv_paddedMesh.vtr +notebooks/scipy2015/SeogiModelMT3Dfor21Dinv/model3D_zxy1Dinv_paddedMesh.vtr +notebooks/scipy2015/SeogiModelMT3Dfor21Dinv/model3D_zyx1Dinv_paddedMesh.vtr +notebooks/scipy2015/SeogiModelMT3Dfor21Dinv/modelTDpaddedMesh.vtr +notebooks/scipy2015/SeogiModelMT3Dfor21Dinv/modelTDuniMesh.vtr diff --git a/notebooks/scipy2015/SeogiModelMT3Dfor21Dinv/MT3DforData1Dinv.ipynb b/notebooks/scipy2015/SeogiModelMT3Dfor21Dinv/MT3DforData1Dinv.ipynb index 46cafc81..f5de911c 100644 --- a/notebooks/scipy2015/SeogiModelMT3Dfor21Dinv/MT3DforData1Dinv.ipynb +++ b/notebooks/scipy2015/SeogiModelMT3Dfor21Dinv/MT3DforData1Dinv.ipynb @@ -27,7 +27,7 @@ }, { "cell_type": "code", - "execution_count": 19, + "execution_count": 3, "metadata": { "collapsed": false }, @@ -55,7 +55,7 @@ }, { "cell_type": "code", - "execution_count": 20, + "execution_count": 4, "metadata": { "collapsed": true }, @@ -79,7 +79,7 @@ }, { "cell_type": "code", - "execution_count": 21, + "execution_count": 5, "metadata": { "collapsed": false }, @@ -108,7 +108,7 @@ }, { "cell_type": "code", - "execution_count": 22, + "execution_count": 7, "metadata": { "collapsed": false }, @@ -122,7 +122,7 @@ }, { "cell_type": "code", - "execution_count": 23, + "execution_count": 8, "metadata": { "collapsed": false }, @@ -208,7 +208,7 @@ }, { "cell_type": "code", - "execution_count": 13, + "execution_count": 12, "metadata": { "collapsed": false }, @@ -219,7 +219,7 @@ "(1088,)" ] }, - "execution_count": 13, + "execution_count": 12, "metadata": {}, "output_type": "execute_result" } @@ -230,7 +230,7 @@ }, { "cell_type": "code", - "execution_count": 30, + "execution_count": 13, "metadata": { "collapsed": false }, @@ -248,14 +248,14 @@ "# Assign the dobs\n", "survey.dtrue = d_true\n", "survey.dobs = d_obs\n", - "survey.std = np.abs(survey.dobs*std) + 0.01*np.linalg.norm(survey.dobs) #survey.dobs*0 + std\n", + "survey.std = np.abs(survey.dobs*std) #+ 0.01*np.linalg.norm(survey.dobs) #survey.dobs*0 + std\n", "# Assign the data weight\n", - "survey.Wd = 1/survey.std #(abs(survey.dobs)*survey.std)" + "Wd = 1/survey.std #(abs(survey.dobs)*survey.std)" ] }, { "cell_type": "code", - "execution_count": 31, + "execution_count": 14, "metadata": { "collapsed": false }, @@ -267,7 +267,7 @@ }, { "cell_type": "code", - "execution_count": 32, + "execution_count": 15, "metadata": { "collapsed": false }, @@ -275,12 +275,13 @@ "source": [ "# Run the setup\n", "mt1DdataZyxList = simpegmt.Utils.dataUtils.convert3Dto1Dobject(mtSeogiDobs,'zyx')\n", - "mt1DdataZxyList = simpegmt.Utils.dataUtils.convert3Dto1Dobject(mtSeogiDobs,'zxy')" + "mt1DdataZxyList = simpegmt.Utils.dataUtils.convert3Dto1Dobject(mtSeogiDobs,'zxy')\n", + "mt1DdataZdetList = simpegmt.Utils.dataUtils.convert3Dto1Dobject(mtSeogiDobs,'zdet')" ] }, { "cell_type": "code", - "execution_count": 33, + "execution_count": 16, "metadata": { "collapsed": false }, @@ -299,7 +300,7 @@ }, { "cell_type": "code", - "execution_count": 49, + "execution_count": 17, "metadata": { "collapsed": false }, @@ -325,7 +326,7 @@ " 1.00000000e-08, 1.00000000e-08, 1.00000000e-08])" ] }, - "execution_count": 49, + "execution_count": 17, "metadata": {}, "output_type": "execute_result" } @@ -336,7 +337,7 @@ }, { "cell_type": "code", - "execution_count": 66, + "execution_count": 22, "metadata": { "collapsed": false }, @@ -356,6 +357,7 @@ " opt.remember('xc')\n", " # Data misfit\n", " dmis = simpeg.DataMisfit.l2_DataMisfit(data.survey)\n", + " dmis.Wd = 1./data.survey.std\n", " # Regularization\n", " regMesh = simpeg.Mesh.TensorMesh([m1d.hx[problem.mapping.sigmaMap.maps[-1].indActive]],m1d.x0)\n", " reg = simpeg.Regularization.Tikhonov(regMesh)\n", @@ -394,7 +396,7 @@ }, { "cell_type": "code", - "execution_count": 67, + "execution_count": 23, "metadata": { "collapsed": false }, @@ -407,28 +409,51 @@ "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.\n", "============================ Inexact Gauss Newton ============================\n", " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", "-----------------------------------------------------------------------------\n", - " 0 9.42e+01 1.17e+02 0.00e+00 1.17e+02 3.97e+01 0 \n", + " 0 1.37e+03 1.07e+03 0.00e+00 1.07e+03 3.77e+02 0 \n", + " 1 1.37e+03 4.40e+01 1.93e-02 7.05e+01 7.43e+01 0 \n", "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.1766e+01\n", - "0 : |xc-x_last| = 5.8426e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.9728e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.9728e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", + "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.0682e+02\n", + "1 : |xc-x_last| = 1.4393e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", + "0 : |proj(x-g)-x| = 7.4277e+01 <= tolG = 1.0000e-01\n", + "0 : |proj(x-g)-x| = 7.4277e+01 <= 1e3*eps = 1.0000e-02\n", + "0 : maxIter = 10 <= iter = 2\n", "------------------------- DONE! -------------------------\n" ] } ], "source": [ + "dat = mt1DdataZyxList[0]\n", "mopt = runInversionModel(dat,problem,m1d,'zyx')\n" ] }, { "cell_type": "code", - "execution_count": 69, + "execution_count": 25, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 25, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "dat.survey" + ] + }, + { + "cell_type": "code", + "execution_count": 31, "metadata": { "collapsed": false }, @@ -441,2076 +466,101 @@ "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.\n", "============================ Inexact Gauss Newton ============================\n", " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", "-----------------------------------------------------------------------------\n", - " 0 1.43e+02 1.17e+02 0.00e+00 1.17e+02 3.97e+01 0 \n", + " 0 1.69e+03 1.07e+03 0.00e+00 1.07e+03 3.77e+02 0 \n", + " 1 1.69e+03 3.88e+01 1.72e-02 6.80e+01 6.98e+01 0 \n", "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.1766e+01\n", - "0 : |xc-x_last| = 5.8497e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.9728e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.9728e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", + "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.0682e+02\n", + "1 : |xc-x_last| = 1.3066e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", + "0 : |proj(x-g)-x| = 6.9756e+01 <= tolG = 1.0000e-01\n", + "0 : |proj(x-g)-x| = 6.9756e+01 <= 1e3*eps = 1.0000e-02\n", + "0 : maxIter = 10 <= iter = 2\n", "------------------------- DONE! -------------------------\n", "Inversion_modAt-700_-500\n", "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.\n", "============================ Inexact Gauss Newton ============================\n", " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", "-----------------------------------------------------------------------------\n", - " 0 2.42e+02 1.03e+02 0.00e+00 1.03e+02 3.61e+01 0 \n", + " 0 1.12e+03 1.10e+03 0.00e+00 1.10e+03 3.79e+02 0 \n", + " 1 1.12e+03 5.89e+01 2.11e-02 8.26e+01 9.39e+01 0 \n", "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.0355e+01\n", - "0 : |xc-x_last| = 5.5520e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.6125e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.6125e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", + "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.1005e+02\n", + "1 : |xc-x_last| = 1.8888e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", + "0 : |proj(x-g)-x| = 9.3881e+01 <= tolG = 1.0000e-01\n", + "0 : |proj(x-g)-x| = 9.3881e+01 <= 1e3*eps = 1.0000e-02\n", + "0 : maxIter = 10 <= iter = 2\n", "------------------------- DONE! -------------------------\n", "Inversion_modAt-700_-300\n", "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.\n", "============================ Inexact Gauss Newton ============================\n", " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", "-----------------------------------------------------------------------------\n", - " 0 1.51e+02 1.14e+02 0.00e+00 1.14e+02 3.79e+01 0 \n", + " 0 7.28e+02 1.14e+03 0.00e+00 1.14e+03 3.81e+02 0 \n", + " 1 7.28e+02 6.54e+01 2.62e-02 8.44e+01 8.78e+01 0 \n", "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.1504e+01\n", - "0 : |xc-x_last| = 5.8931e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.7875e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.7875e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", + "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.1411e+02\n", + "1 : |xc-x_last| = 1.8407e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", + "0 : |proj(x-g)-x| = 8.7820e+01 <= tolG = 1.0000e-01\n", + "0 : |proj(x-g)-x| = 8.7820e+01 <= 1e3*eps = 1.0000e-02\n", + "0 : maxIter = 10 <= iter = 2\n", "------------------------- DONE! -------------------------\n", "Inversion_modAt-700_-100\n", "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.\n", "============================ Inexact Gauss Newton ============================\n", " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.20e+02 1.11e+02 0.00e+00 1.11e+02 3.65e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.1190e+01\n", - "0 : |xc-x_last| = 6.1753e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.6489e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.6489e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-700_100\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 7.99e+01 1.32e+02 0.00e+00 1.32e+02 3.96e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.3305e+01\n", - "0 : |xc-x_last| = 7.0601e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.9626e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.9626e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-700_300\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.96e+02 3.22e+01 0.00e+00 3.22e+01 2.29e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 3.3205e+00\n", - "1 : |xc-x_last| = 2.2199e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 2.2906e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 2.2906e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-700_500\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 4.09e+02 3.27e+01 0.00e+00 3.27e+01 2.17e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 3.3715e+00\n", - "1 : |xc-x_last| = 2.0546e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 2.1726e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 2.1726e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-700_700\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.07e+02 1.36e+02 0.00e+00 1.36e+02 3.94e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.3697e+01\n", - "0 : |xc-x_last| = 7.3447e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.9383e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.9383e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-500_-700\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 9.15e+01 1.07e+02 0.00e+00 1.07e+02 3.74e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.0761e+01\n", - "0 : |xc-x_last| = 5.6974e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.7352e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.7352e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-500_-500\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 8.51e+01 1.09e+02 0.00e+00 1.09e+02 3.78e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.0984e+01\n", - "0 : |xc-x_last| = 5.7587e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.7827e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.7827e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-500_-300\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 7.56e+01 1.02e+02 0.00e+00 1.02e+02 3.72e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.0297e+01\n", - "0 : |xc-x_last| = 5.3214e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.7174e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.7174e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-500_-100\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 9.83e+01 1.20e+02 0.00e+00 1.20e+02 3.90e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.2079e+01\n", - "0 : |xc-x_last| = 6.3216e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.9040e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.9040e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-500_100\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.58e+02 1.22e+02 0.00e+00 1.22e+02 3.66e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.2345e+01\n", - "0 : |xc-x_last| = 7.1814e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.6612e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.6612e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-500_300\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 2.00e+02 2.94e+01 0.00e+00 2.94e+01 2.19e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 3.0361e+00\n", - "1 : |xc-x_last| = 2.1725e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 2.1854e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 2.1854e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-500_500\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 2.77e+02 3.00e+01 0.00e+00 3.00e+01 2.16e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 3.1010e+00\n", - "1 : |xc-x_last| = 2.0539e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 2.1590e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 2.1590e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-500_700\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 3.99e+01 1.26e+02 0.00e+00 1.26e+02 3.80e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.2719e+01\n", - "0 : |xc-x_last| = 6.9852e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.7970e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.7970e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-300_-700\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 2.51e+02 1.11e+02 0.00e+00 1.11e+02 3.83e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.1223e+01\n", - "0 : |xc-x_last| = 5.7989e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.8255e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.8255e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-300_-500\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 8.77e+01 9.40e+01 0.00e+00 9.40e+01 3.37e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 9.5027e+00\n", - "0 : |xc-x_last| = 5.4759e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.3691e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.3691e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-300_-300\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.17e+02 1.11e+02 0.00e+00 1.11e+02 3.87e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.1162e+01\n", - "0 : |xc-x_last| = 5.8012e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.8683e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.8683e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-300_-100\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 9.88e+01 1.28e+02 0.00e+00 1.28e+02 4.05e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.2864e+01\n", - "0 : |xc-x_last| = 6.4960e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 4.0539e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 4.0539e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-300_100\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 7.13e+01 1.28e+02 0.00e+00 1.28e+02 3.79e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.2900e+01\n", - "0 : |xc-x_last| = 7.2851e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.7902e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.7902e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-300_300\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 2.24e+02 3.04e+01 0.00e+00 3.04e+01 2.02e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 3.1373e+00\n", - "1 : |xc-x_last| = 2.2801e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 2.0213e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 2.0213e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-300_500\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 5.66e+02 2.80e+01 0.00e+00 2.80e+01 2.08e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 2.8952e+00\n", - "1 : |xc-x_last| = 1.6263e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 2.0780e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 2.0780e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-300_700\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 7.19e+01 1.28e+02 0.00e+00 1.28e+02 3.79e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.2919e+01\n", - "0 : |xc-x_last| = 7.0614e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.7861e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.7861e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-100_-700\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 2.07e+02 1.15e+02 0.00e+00 1.15e+02 3.87e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.1552e+01\n", - "0 : |xc-x_last| = 6.0418e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.8734e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.8734e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-100_-500\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.09e+02 9.16e+01 0.00e+00 9.16e+01 3.48e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 9.2615e+00\n", - "0 : |xc-x_last| = 5.1793e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.4758e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.4758e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-100_-300\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 2.29e+02 8.86e+01 0.00e+00 8.86e+01 3.67e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 8.9632e+00\n", - "0 : |xc-x_last| = 4.5904e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.6724e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.6724e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-100_-100\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 9.95e+01 1.15e+02 0.00e+00 1.15e+02 3.85e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.1620e+01\n", - "0 : |xc-x_last| = 6.1282e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.8471e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.8471e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-100_100\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 8.89e+01 1.41e+02 0.00e+00 1.41e+02 4.18e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.4166e+01\n", - "0 : |xc-x_last| = 8.0841e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 4.1759e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 4.1759e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-100_300\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 2.30e+02 3.16e+01 0.00e+00 3.16e+01 2.35e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 3.2630e+00\n", - "1 : |xc-x_last| = 2.1631e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 2.3472e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 2.3472e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-100_500\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 2.64e+02 2.85e+01 0.00e+00 2.85e+01 2.10e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 2.9527e+00\n", - "1 : |xc-x_last| = 2.1866e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 2.1030e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 2.1030e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-100_700\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 7.26e+01 1.24e+02 0.00e+00 1.24e+02 3.75e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.2467e+01\n", - "0 : |xc-x_last| = 6.9961e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.7496e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.7496e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt100_-700\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.81e+02 1.22e+02 0.00e+00 1.22e+02 3.78e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.2313e+01\n", - "0 : |xc-x_last| = 6.8371e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.7832e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.7832e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt100_-500\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.97e+02 5.86e+01 0.00e+00 5.86e+01 3.10e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 5.9638e+00\n", - "0 : |xc-x_last| = 3.7925e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.0994e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.0994e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt100_-300\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 2.97e+02 3.01e+01 0.00e+00 3.01e+01 1.75e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 3.1064e+00\n", - "1 : |xc-x_last| = 2.4332e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 1.7523e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 1.7523e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt100_-100\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.53e+02 1.27e+02 0.00e+00 1.27e+02 3.85e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.2840e+01\n", - "0 : |xc-x_last| = 7.2442e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.8528e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.8528e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt100_100\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.59e+02 1.48e+02 0.00e+00 1.48e+02 3.99e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.4939e+01\n", - "0 : |xc-x_last| = 8.1940e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.9944e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.9944e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt100_300\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 3.97e+02 3.06e+01 0.00e+00 3.06e+01 2.35e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 3.1638e+00\n", - "1 : |xc-x_last| = 1.9566e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 2.3465e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 2.3465e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt100_500\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 3.34e+02 2.94e+01 0.00e+00 2.94e+01 1.96e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 3.0363e+00\n", - "1 : |xc-x_last| = 1.8749e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 1.9639e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 1.9639e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt100_700\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.41e+02 1.33e+02 0.00e+00 1.33e+02 3.99e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.3369e+01\n", - "0 : |xc-x_last| = 7.1080e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.9943e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.9943e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt300_-700\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.07e+02 1.44e+02 0.00e+00 1.44e+02 3.82e+01 0 \n", - " 1 1.07e+02 3.54e+01 2.89e-03 3.57e+01 3.11e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.4521e+01\n", - "1 : |xc-x_last| = 2.9250e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.1116e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.1116e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 2\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt300_-500\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 4.99e+02 3.43e+01 0.00e+00 3.43e+01 2.40e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 3.5267e+00\n", - "1 : |xc-x_last| = 2.0777e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 2.3969e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 2.3969e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt300_-300\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 2.66e+03 1.22e+03 0.00e+00 1.22e+03 4.79e+02 0 \n", - " 1 2.66e+03 1.28e+02 1.59e-03 1.32e+02 6.85e+01 0 \n", - " 2 2.66e+03 5.99e+01 3.95e-03 7.04e+01 5.02e+00 0 Skip BFGS \n", - " 3 3.32e+02 5.91e+01 4.02e-03 6.04e+01 1.54e+01 0 Skip BFGS \n", - " 4 3.32e+02 4.14e+01 2.81e-02 5.07e+01 1.42e+01 0 \n", - " 5 3.32e+02 3.98e+01 2.39e-02 4.78e+01 6.27e+00 0 \n", - " 6 4.15e+01 3.88e+01 2.61e-02 3.99e+01 1.51e+01 0 Skip BFGS \n", - " 7 4.15e+01 3.50e+01 5.58e-02 3.73e+01 1.40e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.2199e+02\n", - "0 : |xc-x_last| = 8.5871e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 1.3985e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 1.3985e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 8\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt300_-100\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 6.86e+01 1.67e+02 0.00e+00 1.67e+02 4.17e+01 0 \n", - " 1 6.86e+01 5.21e+01 3.65e-03 5.24e+01 4.46e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.6779e+01\n", - "0 : |xc-x_last| = 3.1821e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 4.4603e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 4.4603e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 2\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt300_100\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 8.73e+01 1.66e+02 0.00e+00 1.66e+02 3.96e+01 0 \n", - " 1 8.73e+01 3.65e+01 6.23e-03 3.70e+01 3.67e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.6655e+01\n", - "1 : |xc-x_last| = 2.4568e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.6714e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.6714e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 2\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt300_300\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 3.99e+02 3.13e+01 0.00e+00 3.13e+01 2.30e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 3.2312e+00\n", - "1 : |xc-x_last| = 1.8739e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 2.3039e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 2.3039e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt300_500\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 2.46e+02 2.54e+01 0.00e+00 2.54e+01 1.98e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 2.6358e+00\n", - "1 : |xc-x_last| = 2.1834e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 1.9823e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 1.9823e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt300_700\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 5.99e+01 1.26e+02 0.00e+00 1.26e+02 3.65e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.2712e+01\n", - "0 : |xc-x_last| = 7.1861e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.6454e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.6454e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt500_-700\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 6.37e+01 1.48e+02 0.00e+00 1.48e+02 4.04e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.4888e+01\n", - "0 : |xc-x_last| = 8.8279e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 4.0428e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 4.0428e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt500_-500\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 2.47e+02 4.52e+01 0.00e+00 4.52e+01 2.99e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 4.6219e+00\n", - "0 : |xc-x_last| = 3.0713e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 2.9892e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 2.9892e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt500_-300\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 7.76e+02 3.70e+02 0.00e+00 3.70e+02 1.80e+02 0 \n", - " 1 7.76e+02 4.80e+01 1.94e-03 4.95e+01 2.63e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 3.7135e+01\n", - "1 : |xc-x_last| = 1.4930e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 2.6307e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 2.6307e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 2\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt500_-100\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 5.30e+01 1.56e+02 0.00e+00 1.56e+02 4.11e+01 0 \n", - " 1 5.30e+01 4.53e+01 3.32e-03 4.54e+01 3.88e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.5691e+01\n", - "0 : |xc-x_last| = 3.1356e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.8848e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.8848e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 2\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt500_100\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.42e+02 1.42e+02 0.00e+00 1.42e+02 3.78e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.4327e+01\n", - "0 : |xc-x_last| = 9.1426e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.7802e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.7802e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt500_300\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.39e+02 3.72e+01 0.00e+00 3.72e+01 2.63e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 3.8246e+00\n", - "1 : |xc-x_last| = 2.8630e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 2.6294e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 2.6294e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt500_500\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 4.86e+02 3.13e+01 0.00e+00 3.13e+01 2.17e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 3.2257e+00\n", - "1 : |xc-x_last| = 1.8364e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 2.1663e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 2.1663e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt500_700\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.42e+02 1.27e+02 0.00e+00 1.27e+02 3.80e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.2782e+01\n", - "0 : |xc-x_last| = 7.0677e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.7996e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.7996e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt700_-700\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 9.44e+01 1.07e+02 0.00e+00 1.07e+02 3.61e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.0777e+01\n", - "0 : |xc-x_last| = 6.0624e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.6147e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.6147e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt700_-500\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.22e+02 7.34e+01 0.00e+00 7.34e+01 3.25e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 7.4419e+00\n", - "0 : |xc-x_last| = 4.7542e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.2461e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.2461e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt700_-300\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 9.49e+01 5.54e+01 0.00e+00 5.54e+01 3.00e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 5.6400e+00\n", - "0 : |xc-x_last| = 3.5995e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 2.9982e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 2.9982e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt700_-100\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 9.95e+01 1.16e+02 0.00e+00 1.16e+02 3.79e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.1677e+01\n", - "0 : |xc-x_last| = 6.2824e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.7893e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.7893e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt700_100\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.57e+02 1.37e+02 0.00e+00 1.37e+02 3.98e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.3786e+01\n", - "0 : |xc-x_last| = 8.1967e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.9796e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.9796e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt700_300\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 3.73e+02 3.31e+01 0.00e+00 3.31e+01 2.39e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 3.4053e+00\n", - "1 : |xc-x_last| = 1.8793e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 2.3851e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 2.3851e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt700_500\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 3.36e+02 3.40e+01 0.00e+00 3.40e+01 2.37e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 3.5047e+00\n", - "1 : |xc-x_last| = 1.9647e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 2.3722e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 2.3722e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt700_700\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 2.19e+02 1.37e+02 0.00e+00 1.37e+02 3.91e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.3813e+01\n", - "0 : |xc-x_last| = 7.5331e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.9057e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.9057e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-700_-700\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.07e+02 1.01e+02 0.00e+00 1.01e+02 3.61e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.0201e+01\n", - "0 : |xc-x_last| = 5.3873e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.6096e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.6096e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-700_-500\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.21e+02 9.66e+01 0.00e+00 9.66e+01 3.53e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 9.7594e+00\n", - "0 : |xc-x_last| = 5.3855e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.5277e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.5277e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-700_-300\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.11e+02 9.29e+01 0.00e+00 9.29e+01 3.48e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 9.3881e+00\n", - "0 : |xc-x_last| = 5.1323e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.4781e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.4781e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-700_-100\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.02e+02 8.11e+01 0.00e+00 8.11e+01 3.47e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 8.2061e+00\n", - "0 : |xc-x_last| = 4.3819e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.4684e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.4684e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-700_100\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 2.25e+02 4.34e+01 0.00e+00 4.34e+01 2.54e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 4.4386e+00\n", - "0 : |xc-x_last| = 3.3773e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 2.5426e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 2.5426e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-700_300\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 6.59e+02 1.21e+02 0.00e+00 1.21e+02 6.78e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.2221e+01\n", - "0 : |xc-x_last| = 3.3931e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 6.7813e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 6.7813e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-700_500\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.20e+03 1.20e+02 0.00e+00 1.20e+02 6.91e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.2144e+01\n", - "0 : |xc-x_last| = 3.1214e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 6.9149e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 6.9149e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-700_700\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.79e+02 4.43e+01 0.00e+00 4.43e+01 2.53e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 4.5340e+00\n", - "0 : |xc-x_last| = 3.6308e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 2.5344e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 2.5344e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-500_-700\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.83e+02 1.09e+02 0.00e+00 1.09e+02 3.93e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.1040e+01\n", - "0 : |xc-x_last| = 5.4541e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.9340e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.9340e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-500_-500\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.44e+02 9.56e+01 0.00e+00 9.56e+01 3.63e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 9.6589e+00\n", - "0 : |xc-x_last| = 4.9909e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.6339e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.6339e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-500_-300\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.22e+02 8.76e+01 0.00e+00 8.76e+01 3.28e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 8.8584e+00\n", - "0 : |xc-x_last| = 5.0487e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.2789e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.2789e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-500_-100\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.27e+02 7.38e+01 0.00e+00 7.38e+01 3.27e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 7.4766e+00\n", - "0 : |xc-x_last| = 4.6538e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.2669e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.2669e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-500_100\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.92e+02 4.61e+01 0.00e+00 4.61e+01 2.55e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 4.7117e+00\n", - "0 : |xc-x_last| = 4.1157e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 2.5464e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 2.5464e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-500_300\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.12e+03 9.17e+01 0.00e+00 9.17e+01 5.80e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 9.2749e+00\n", - "1 : |xc-x_last| = 2.6911e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 5.8036e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 5.8036e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-500_500\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 6.66e+02 9.90e+01 0.00e+00 9.90e+01 5.87e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.0000e+01\n", - "0 : |xc-x_last| = 3.2072e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 5.8652e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 5.8652e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-500_700\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.57e+02 3.95e+01 0.00e+00 3.95e+01 2.20e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 4.0479e+00\n", - "0 : |xc-x_last| = 3.9463e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 2.2028e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 2.2028e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-300_-700\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 2.28e+02 1.00e+02 0.00e+00 1.00e+02 3.69e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.0138e+01\n", - "0 : |xc-x_last| = 5.2547e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.6862e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.6862e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-300_-500\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 7.29e+01 9.71e+01 0.00e+00 9.71e+01 3.54e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 9.8059e+00\n", - "0 : |xc-x_last| = 5.2673e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.5373e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.5373e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-300_-300\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 9.56e+01 9.27e+01 0.00e+00 9.27e+01 3.44e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 9.3670e+00\n", - "0 : |xc-x_last| = 5.0924e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.4424e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.4424e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-300_-100\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.39e+02 6.65e+01 0.00e+00 6.65e+01 3.07e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 6.7482e+00\n", - "0 : |xc-x_last| = 4.2930e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.0690e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.0690e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-300_100\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.79e+02 3.95e+01 0.00e+00 3.95e+01 2.24e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 4.0532e+00\n", - "0 : |xc-x_last| = 3.5885e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 2.2417e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 2.2417e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-300_300\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 4.89e+02 9.68e+01 0.00e+00 9.68e+01 5.87e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 9.7833e+00\n", - "0 : |xc-x_last| = 3.4074e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 5.8749e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 5.8749e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-300_500\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 8.08e+02 8.66e+01 0.00e+00 8.66e+01 5.25e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 8.7551e+00\n", - "1 : |xc-x_last| = 2.9192e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 5.2476e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 5.2476e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-300_700\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 2.95e+02 3.47e+01 0.00e+00 3.47e+01 1.99e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 3.5731e+00\n", - "0 : |xc-x_last| = 3.2502e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 1.9861e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 1.9861e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-100_-700\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 2.57e+02 1.01e+02 0.00e+00 1.01e+02 3.68e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.0234e+01\n", - "0 : |xc-x_last| = 5.5230e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.6814e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.6814e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-100_-500\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.23e+02 1.06e+02 0.00e+00 1.06e+02 3.72e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.0747e+01\n", - "0 : |xc-x_last| = 5.8517e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.7156e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.7156e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-100_-300\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.34e+02 1.05e+02 0.00e+00 1.05e+02 3.67e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.0599e+01\n", - "0 : |xc-x_last| = 5.6348e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.6674e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.6674e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-100_-100\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.37e+02 8.15e+01 0.00e+00 8.15e+01 3.41e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 8.2455e+00\n", - "0 : |xc-x_last| = 5.0343e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.4144e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.4144e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-100_100\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 3.25e+02 3.83e+01 0.00e+00 3.83e+01 2.23e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 3.9293e+00\n", - "0 : |xc-x_last| = 3.1900e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 2.2253e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 2.2253e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-100_300\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 5.63e+02 8.86e+01 0.00e+00 8.86e+01 5.13e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 8.9615e+00\n", - "0 : |xc-x_last| = 3.0613e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 5.1257e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 5.1257e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-100_500\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 3.67e+02 9.28e+01 0.00e+00 9.28e+01 5.40e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 9.3781e+00\n", - "0 : |xc-x_last| = 3.5506e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 5.4019e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 5.4019e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt-100_700\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 5.02e+02 4.30e+01 0.00e+00 4.30e+01 2.17e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 4.3998e+00\n", - "0 : |xc-x_last| = 3.1262e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 2.1679e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 2.1679e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt100_-700\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.44e+02 8.98e+01 0.00e+00 8.98e+01 3.62e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 9.0810e+00\n", - "0 : |xc-x_last| = 4.7947e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.6169e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.6169e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt100_-500\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.01e+02 1.17e+02 0.00e+00 1.17e+02 3.97e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.1822e+01\n", - "0 : |xc-x_last| = 6.9117e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.9668e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.9668e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt100_-300\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.64e+02 1.17e+02 0.00e+00 1.17e+02 3.93e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.1801e+01\n", - "0 : |xc-x_last| = 6.3223e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.9263e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.9263e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt100_-100\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.40e+02 6.70e+01 0.00e+00 6.70e+01 3.18e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 6.8005e+00\n", - "0 : |xc-x_last| = 4.3751e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.1786e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.1786e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt100_100\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.65e+02 3.76e+01 0.00e+00 3.76e+01 1.86e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 3.8645e+00\n", - "0 : |xc-x_last| = 4.1042e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 1.8617e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 1.8617e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt100_300\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 3.24e+02 1.14e+02 0.00e+00 1.14e+02 6.64e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.1451e+01\n", - "0 : |xc-x_last| = 4.2288e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 6.6371e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 6.6371e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt100_500\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 5.36e+02 8.90e+01 0.00e+00 8.90e+01 5.48e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 8.9953e+00\n", - "0 : |xc-x_last| = 3.1393e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 5.4782e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 5.4782e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt100_700\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 2.70e+02 3.77e+01 0.00e+00 3.77e+01 2.04e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 3.8660e+00\n", - "0 : |xc-x_last| = 3.3812e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 2.0401e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 2.0401e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt300_-700\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.99e+02 4.73e+01 0.00e+00 4.73e+01 2.65e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 4.8348e+00\n", - "0 : |xc-x_last| = 3.3721e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 2.6479e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 2.6479e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt300_-500\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 6.62e+02 4.67e+02 0.00e+00 4.67e+02 2.26e+02 0 \n", - " 1 6.62e+02 5.78e+01 1.42e-03 5.87e+01 3.06e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 4.6812e+01\n", - "1 : |xc-x_last| = 2.0385e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.0562e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.0562e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 2\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt300_-300\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.84e+03 1.82e+03 0.00e+00 1.82e+03 6.43e+02 0 \n", - " 1 1.84e+03 1.86e+02 2.02e-03 1.89e+02 9.10e+01 0 \n", - " 2 1.84e+03 6.11e+01 6.12e-03 7.24e+01 1.05e+01 0 Skip BFGS \n", - " 3 2.31e+02 5.60e+01 7.68e-03 5.78e+01 1.65e+01 0 Skip BFGS \n", - " 4 2.31e+02 3.65e+01 3.73e-02 4.51e+01 1.20e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.8186e+02\n", - "1 : |xc-x_last| = 1.6068e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 1.2001e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 1.2001e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 5\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt300_-100\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 5.65e+02 7.36e+01 0.00e+00 7.36e+01 3.86e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 7.4581e+00\n", - "0 : |xc-x_last| = 3.7379e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.8553e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.8553e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt300_100\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.89e+02 3.55e+01 0.00e+00 3.55e+01 1.55e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 3.6522e+00\n", - "0 : |xc-x_last| = 3.6420e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 1.5467e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 1.5467e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt300_300\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 5.81e+02 1.11e+02 0.00e+00 1.11e+02 7.01e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.1172e+01\n", - "0 : |xc-x_last| = 3.2476e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 7.0108e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 7.0108e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt300_500\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 9.73e+02 8.52e+01 0.00e+00 8.52e+01 5.21e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 8.6151e+00\n", - "1 : |xc-x_last| = 2.6690e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 5.2083e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 5.2083e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt300_700\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 3.30e+02 3.68e+01 0.00e+00 3.68e+01 2.06e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 3.7823e+00\n", - "0 : |xc-x_last| = 2.9969e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 2.0632e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 2.0632e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt500_-700\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.21e+02 6.93e+01 0.00e+00 6.93e+01 3.20e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 7.0267e+00\n", - "0 : |xc-x_last| = 3.9405e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.2005e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.2005e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt500_-500\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 2.39e+02 3.47e+01 0.00e+00 3.47e+01 2.65e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 3.5688e+00\n", - "1 : |xc-x_last| = 2.3554e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 2.6546e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 2.6546e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt500_-300\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 3.16e+02 2.65e+01 0.00e+00 2.65e+01 2.05e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 2.7466e+00\n", - "1 : |xc-x_last| = 1.6641e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 2.0456e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 2.0456e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt500_-100\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 2.29e+02 3.58e+01 0.00e+00 3.58e+01 2.11e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 3.6774e+00\n", - "0 : |xc-x_last| = 3.0232e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 2.1091e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 2.1091e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt500_100\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.80e+02 3.31e+01 0.00e+00 3.31e+01 1.67e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 3.4065e+00\n", - "0 : |xc-x_last| = 3.5591e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 1.6716e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 1.6716e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt500_300\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 7.24e+02 1.59e+02 0.00e+00 1.59e+02 9.55e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.5965e+01\n", - "0 : |xc-x_last| = 3.4100e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 9.5514e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 9.5514e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt500_500\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 5.66e+02 9.71e+01 0.00e+00 9.71e+01 5.76e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 9.8086e+00\n", - "0 : |xc-x_last| = 3.1907e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 5.7591e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 5.7591e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt500_700\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 2.28e+02 3.86e+01 0.00e+00 3.86e+01 2.28e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 3.9599e+00\n", - "0 : |xc-x_last| = 3.3778e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 2.2819e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 2.2819e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt700_-700\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 9.55e+01 9.82e+01 0.00e+00 9.82e+01 3.59e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 9.9157e+00\n", - "0 : |xc-x_last| = 5.3459e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.5874e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.5874e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt700_-500\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 8.43e+01 1.29e+02 0.00e+00 1.29e+02 3.82e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.2966e+01\n", - "0 : |xc-x_last| = 7.1708e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.8151e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.8151e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt700_-300\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.41e+02 1.27e+02 0.00e+00 1.27e+02 3.70e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.2768e+01\n", - "0 : |xc-x_last| = 7.2504e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.6990e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.6990e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt700_-100\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 1.33e+02 9.86e+01 0.00e+00 9.86e+01 3.72e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 9.9560e+00\n", - "0 : |xc-x_last| = 5.2234e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 3.7166e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 3.7166e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt700_100\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 2.31e+02 4.70e+01 0.00e+00 4.70e+01 2.66e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 4.7956e+00\n", - "0 : |xc-x_last| = 3.4674e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 2.6619e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 2.6619e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt700_300\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 5.97e+02 1.29e+02 0.00e+00 1.29e+02 7.34e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.3048e+01\n", - "0 : |xc-x_last| = 3.6080e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 7.3417e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 7.3417e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt700_500\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 5.51e+02 1.15e+02 0.00e+00 1.15e+02 6.64e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 1.1598e+01\n", - "0 : |xc-x_last| = 3.3062e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 6.6369e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 6.6369e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n", - "Inversion_modAt700_700\n", - "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.\n", - "============================ Inexact Gauss Newton ============================\n", - " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", - "-----------------------------------------------------------------------------\n", - " 0 2.81e+02 3.93e+01 0.00e+00 3.93e+01 2.33e+01 0 \n", - "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 4.0343e+00\n", - "0 : |xc-x_last| = 2.9948e+00 <= tolX*(1+|x0|) = 2.9759e+00\n", - "0 : |proj(x-g)-x| = 2.3254e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 2.3254e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 1\n", - "------------------------- DONE! -------------------------\n" + "-----------------------------------------------------------------------------\n" + ] + }, + { + "ename": "KeyboardInterrupt", + "evalue": "", + "output_type": "error", + "traceback": [ + "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m", + "\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[0;32m 1\u001b[0m \u001b[1;32mif\u001b[0m \u001b[0mTrue\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 2\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0mdat\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mmt1DdataZyxList\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m----> 3\u001b[1;33m \u001b[0mrunInversionModel\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mdat\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mproblem\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mm1d\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;34m'zyx'\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 4\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 5\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0mdat\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mmt1DdataZxyList\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", + "\u001b[1;32m\u001b[0m in \u001b[0;36mrunInversionModel\u001b[1;34m(data, problem, m1d, nameflag)\u001b[0m\n\u001b[0;32m 39\u001b[0m \u001b[0mm_0\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mnp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mlog\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;36m1e-2\u001b[0m\u001b[1;33m+\u001b[0m\u001b[1;36m0\u001b[0m\u001b[1;33m*\u001b[0m\u001b[0mproblem\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0msigmaPrimary\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mproblem\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mmapping\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0msigmaMap\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mmaps\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;33m-\u001b[0m\u001b[1;36m1\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mindActive\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 40\u001b[0m \u001b[0mproblem\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0msurvey\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mmtrue\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mm_0\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m---> 41\u001b[1;33m \u001b[0mmopt\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0minv\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mrun\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mm_0\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 42\u001b[0m \u001b[1;31m# Save the model\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 43\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n", + "\u001b[1;32m/media/gudni/ExtraDrive1/Codes/python/simpeg/SimPEG/Utils/CounterUtils.pyc\u001b[0m in \u001b[0;36mwrapper\u001b[1;34m(self, *args, **kwargs)\u001b[0m\n\u001b[0;32m 90\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 91\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[0mcountTic\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---> 92\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 93\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[0mcountToc\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[0;32m 94\u001b[0m \u001b[1;32mreturn\u001b[0m \u001b[0mout\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", + "\u001b[1;32m/media/gudni/ExtraDrive1/Codes/python/simpeg/SimPEG/Inversion.pyc\u001b[0m in \u001b[0;36mrun\u001b[1;34m(self, m0)\u001b[0m\n\u001b[0;32m 60\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0minvProb\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mstartup\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mm0\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 61\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mdirectiveList\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mcall\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;34m'initialize'\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m---> 62\u001b[1;33m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mm\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mopt\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mminimize\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0minvProb\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mevalFunction\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0minvProb\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mcurModel\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 63\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mdirectiveList\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mcall\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;34m'finish'\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 64\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n", + "\u001b[1;32m/media/gudni/ExtraDrive1/Codes/python/simpeg/SimPEG/Utils/CounterUtils.pyc\u001b[0m in \u001b[0;36mwrapper\u001b[1;34m(self, *args, **kwargs)\u001b[0m\n\u001b[0;32m 90\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 91\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[0mcountTic\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---> 92\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 93\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[0mcountToc\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[0;32m 94\u001b[0m \u001b[1;32mreturn\u001b[0m \u001b[0mout\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", + "\u001b[1;32m/media/gudni/ExtraDrive1/Codes/python/simpeg/SimPEG/Optimization.pyc\u001b[0m in \u001b[0;36mminimize\u001b[1;34m(self, evalFunction, x0)\u001b[0m\n\u001b[0;32m 206\u001b[0m \u001b[1;32mwhile\u001b[0m \u001b[0mTrue\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 207\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mdoStartIteration\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 208\u001b[1;33m \u001b[0mself\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[0mg\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mH\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mevalFunction\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mxc\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mreturn_g\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mTrue\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mreturn_H\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mTrue\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 209\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mprintIter\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 210\u001b[0m \u001b[1;32mif\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mstoppingCriteria\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m \u001b[1;32mbreak\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", + "\u001b[1;32m/media/gudni/ExtraDrive1/Codes/python/simpeg/SimPEG/Utils/CounterUtils.pyc\u001b[0m in \u001b[0;36mwrapper\u001b[1;34m(self, *args, **kwargs)\u001b[0m\n\u001b[0;32m 90\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 91\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[0mcountTic\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---> 92\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 93\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[0mcountToc\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[0;32m 94\u001b[0m \u001b[1;32mreturn\u001b[0m \u001b[0mout\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", + "\u001b[1;32m/media/gudni/ExtraDrive1/Codes/python/simpeg/SimPEG/InvProblem.pyc\u001b[0m in \u001b[0;36mevalFunction\u001b[1;34m(self, m, return_g, return_H)\u001b[0m\n\u001b[0;32m 124\u001b[0m \u001b[0mout\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;33m(\u001b[0m\u001b[0mf\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 125\u001b[0m \u001b[1;32mif\u001b[0m \u001b[0mreturn_g\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 126\u001b[1;33m \u001b[0mphi_dDeriv\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mdmisfit\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mevalDeriv\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mm\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mu\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mu\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 127\u001b[0m \u001b[0mphi_mDeriv\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mreg\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mevalDeriv\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[0;32m 128\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n", + "\u001b[1;32m/media/gudni/ExtraDrive1/Codes/python/simpeg/SimPEG/Utils/CounterUtils.pyc\u001b[0m in \u001b[0;36mwrapper\u001b[1;34m(self, *args, **kwargs)\u001b[0m\n\u001b[0;32m 90\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 91\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[0mcountTic\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---> 92\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 93\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[0mcountToc\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[0;32m 94\u001b[0m \u001b[1;32mreturn\u001b[0m \u001b[0mout\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", + "\u001b[1;32m/media/gudni/ExtraDrive1/Codes/python/simpeg/SimPEG/DataMisfit.pyc\u001b[0m in \u001b[0;36mevalDeriv\u001b[1;34m(self, m, u)\u001b[0m\n\u001b[0;32m 156\u001b[0m \u001b[0msurvey\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[0m\n\u001b[0;32m 157\u001b[0m \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[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[1;32m--> 158\u001b[1;33m \u001b[1;32mreturn\u001b[0m \u001b[0mprob\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mJtvec\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mm\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mWd\u001b[0m \u001b[1;33m*\u001b[0m \u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mWd\u001b[0m \u001b[1;33m*\u001b[0m \u001b[0msurvey\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mresidual\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mm\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mu\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[0mu\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mu\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 159\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 160\u001b[0m \u001b[1;33m@\u001b[0m\u001b[0mUtils\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mtimeIt\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", + "\u001b[1;32m/media/gudni/ExtraDrive1/Codes/python/simpegmt/simpegMT/BaseMT.py\u001b[0m in \u001b[0;36mJtvec\u001b[1;34m(self, m, v, u)\u001b[0m\n\u001b[0;32m 125\u001b[0m \u001b[1;31m# Get the\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 126\u001b[0m \u001b[0mdA_duIT\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mATinv\u001b[0m \u001b[1;33m*\u001b[0m \u001b[0mPTv\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 127\u001b[1;33m \u001b[0mdA_dmT\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mgetADeriv_m\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mfreq\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mu_src\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mmkvc\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mdA_duIT\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0madjoint\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mTrue\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 128\u001b[0m \u001b[0mdRHS_dmT\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mgetRHSDeriv_m\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mfreq\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mmkvc\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mdA_duIT\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0madjoint\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mTrue\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 129\u001b[0m \u001b[1;31m# Make du_dmT\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", + "\u001b[1;32m/media/gudni/ExtraDrive1/Codes/python/simpegmt/simpegMT/ProblemMT1D/Problems.pyc\u001b[0m in \u001b[0;36mgetADeriv_m\u001b[1;34m(self, freq, u, v, adjoint)\u001b[0m\n\u001b[0;32m 68\u001b[0m \u001b[1;31m#\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 69\u001b[0m \u001b[0mu_src\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mu\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;34m'e_1dSolution'\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m---> 70\u001b[1;33m \u001b[0mdMf_dsig\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[0mgetFaceInnerProductDeriv\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mcurModel\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0msigma\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mu_src\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;33m*\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mcurModel\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0msigmaDeriv\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 71\u001b[0m \u001b[1;32mif\u001b[0m \u001b[0madjoint\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 72\u001b[0m \u001b[1;32mreturn\u001b[0m \u001b[1;36m1j\u001b[0m \u001b[1;33m*\u001b[0m \u001b[0momega\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mfreq\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;33m*\u001b[0m \u001b[1;33m(\u001b[0m \u001b[0mdMf_dsig\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mT\u001b[0m \u001b[1;33m*\u001b[0m \u001b[0mv\u001b[0m \u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", + "\u001b[1;32m/media/gudni/ExtraDrive1/Codes/python/simpeg/SimPEG/PropMaps.pyc\u001b[0m in \u001b[0;36mfget\u001b[1;34m(self)\u001b[0m\n\u001b[0;32m 107\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 108\u001b[0m \u001b[0mm\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'%sModel'\u001b[0m\u001b[1;33m%\u001b[0m\u001b[0mprop\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mname\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 109\u001b[1;33m \u001b[1;32mreturn\u001b[0m \u001b[0mmapping\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mderiv\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 110\u001b[0m \u001b[1;32mreturn\u001b[0m \u001b[0mproperty\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mfget\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mfget\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 111\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n", + "\u001b[1;32m/media/gudni/ExtraDrive1/Codes/python/simpeg/SimPEG/Maps.pyc\u001b[0m in \u001b[0;36mderiv\u001b[1;34m(self, m)\u001b[0m\n\u001b[0;32m 181\u001b[0m \u001b[0mmi\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mm\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 182\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0mmap_i\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mreversed\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mmaps\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 183\u001b[1;33m \u001b[0mderiv\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mmap_i\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mderiv\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mmi\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;33m*\u001b[0m \u001b[0mderiv\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 184\u001b[0m \u001b[0mmi\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mmap_i\u001b[0m \u001b[1;33m*\u001b[0m \u001b[0mmi\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 185\u001b[0m \u001b[1;32mreturn\u001b[0m \u001b[0mderiv\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", + "\u001b[1;32m/media/gudni/ExtraDrive1/Codes/python/simpeg/SimPEG/Maps.pyc\u001b[0m in \u001b[0;36mderiv\u001b[1;34m(self, m)\u001b[0m\n\u001b[0;32m 249\u001b[0m \u001b[0;31m\\\u001b[0m\u001b[0;31m\\\u001b[0m\u001b[0mfrac\u001b[0m\u001b[1;33m{\u001b[0m\u001b[0;31m\\\u001b[0m\u001b[0mpartial\u001b[0m\u001b[0;31m \u001b[0m\u001b[0;31m\\\u001b[0m\u001b[0mexp\u001b[0m\u001b[1;33m{\u001b[0m\u001b[0mm\u001b[0m\u001b[1;33m}\u001b[0m\u001b[1;33m}\u001b[0m\u001b[1;33m{\u001b[0m\u001b[0;31m\\\u001b[0m\u001b[0mpartial\u001b[0m \u001b[0mm\u001b[0m\u001b[1;33m}\u001b[0m \u001b[1;33m=\u001b[0m\u001b[0;31m \u001b[0m\u001b[0;31m\\\u001b[0m\u001b[0;31m\\\u001b[0m\u001b[0mtext\u001b[0m\u001b[1;33m{\u001b[0m\u001b[0msdiag\u001b[0m\u001b[1;33m}\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0;31m\\\u001b[0m\u001b[0mexp\u001b[0m\u001b[1;33m{\u001b[0m\u001b[0mm\u001b[0m\u001b[1;33m}\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 250\u001b[0m \"\"\"\n\u001b[1;32m--> 251\u001b[1;33m \u001b[1;32mreturn\u001b[0m \u001b[0mUtils\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0msdiag\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mnp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mexp\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mUtils\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mmkvc\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mm\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 252\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 253\u001b[0m \u001b[1;32mclass\u001b[0m \u001b[0mReciprocalMap\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mIdentityMap\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", + "\u001b[1;32m/media/gudni/ExtraDrive1/Codes/python/simpeg/SimPEG/Utils/matutils.pyc\u001b[0m in \u001b[0;36msdiag\u001b[1;34m(h)\u001b[0m\n\u001b[0;32m 38\u001b[0m \u001b[1;32mdef\u001b[0m \u001b[0msdiag\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mh\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 39\u001b[0m \u001b[1;34m\"\"\"Sparse diagonal matrix\"\"\"\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m---> 40\u001b[1;33m \u001b[1;32mreturn\u001b[0m \u001b[0msp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mspdiags\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mmkvc\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mh\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;36m0\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mh\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0msize\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mh\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0msize\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mformat\u001b[0m\u001b[1;33m=\u001b[0m\u001b[1;34m\"csr\"\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 41\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 42\u001b[0m \u001b[1;32mdef\u001b[0m \u001b[0msdInv\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mM\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", + "\u001b[1;32m/home/gudni/anaconda/lib/python2.7/site-packages/scipy/sparse/construct.pyc\u001b[0m in \u001b[0;36mspdiags\u001b[1;34m(data, diags, m, n, format)\u001b[0m\n\u001b[0;32m 58\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 59\u001b[0m \"\"\"\n\u001b[1;32m---> 60\u001b[1;33m \u001b[1;32mreturn\u001b[0m \u001b[0mdia_matrix\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mdata\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mdiags\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mshape\u001b[0m\u001b[1;33m=\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mm\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mn\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0masformat\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mformat\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 61\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 62\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n", + "\u001b[1;32m/home/gudni/anaconda/lib/python2.7/site-packages/scipy/sparse/base.pyc\u001b[0m in \u001b[0;36masformat\u001b[1;34m(self, format)\u001b[0m\n\u001b[0;32m 211\u001b[0m \u001b[1;32mreturn\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 212\u001b[0m \u001b[1;32melse\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 213\u001b[1;33m \u001b[1;32mreturn\u001b[0m \u001b[0mgetattr\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;34m'to'\u001b[0m \u001b[1;33m+\u001b[0m \u001b[0mformat\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 214\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 215\u001b[0m \u001b[1;31m###################################################################\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", + "\u001b[1;32m/home/gudni/anaconda/lib/python2.7/site-packages/scipy/sparse/dia.pyc\u001b[0m in \u001b[0;36mtocsr\u001b[1;34m(self)\u001b[0m\n\u001b[0;32m 225\u001b[0m \u001b[1;32mdef\u001b[0m \u001b[0mtocsr\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 226\u001b[0m \u001b[1;31m#this could be faster\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 227\u001b[1;33m \u001b[1;32mreturn\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mtocoo\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mtocsr\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 228\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 229\u001b[0m \u001b[1;32mdef\u001b[0m \u001b[0mtocsc\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", + "\u001b[1;32m/home/gudni/anaconda/lib/python2.7/site-packages/scipy/sparse/dia.pyc\u001b[0m in \u001b[0;36mtocoo\u001b[1;34m(self)\u001b[0m\n\u001b[0;32m 250\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 251\u001b[0m \u001b[1;32mfrom\u001b[0m \u001b[1;33m.\u001b[0m\u001b[0mcoo\u001b[0m \u001b[1;32mimport\u001b[0m \u001b[0mcoo_matrix\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 252\u001b[1;33m \u001b[1;32mreturn\u001b[0m \u001b[0mcoo_matrix\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mdata\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mrow\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mcol\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mshape\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mshape\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 253\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 254\u001b[0m \u001b[1;31m# needed by _data_matrix\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", + "\u001b[1;32m/home/gudni/anaconda/lib/python2.7/site-packages/scipy/sparse/coo.pyc\u001b[0m in \u001b[0;36m__init__\u001b[1;34m(self, arg1, shape, dtype, copy)\u001b[0m\n\u001b[0;32m 204\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mdata\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mdata\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mastype\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mdtype\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 205\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 206\u001b[1;33m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0m_check\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 207\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 208\u001b[0m \u001b[1;32mdef\u001b[0m \u001b[0mgetnnz\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0maxis\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mNone\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", + "\u001b[1;32m/home/gudni/anaconda/lib/python2.7/site-packages/scipy/sparse/coo.pyc\u001b[0m in \u001b[0;36m_check\u001b[1;34m(self)\u001b[0m\n\u001b[0;32m 257\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 258\u001b[0m \u001b[1;32mif\u001b[0m \u001b[0mnnz\u001b[0m \u001b[1;33m>\u001b[0m \u001b[1;36m0\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 259\u001b[1;33m \u001b[1;32mif\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mrow\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mmax\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[0mshape\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;36m0\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 260\u001b[0m \u001b[1;32mraise\u001b[0m \u001b[0mValueError\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;34m'row index exceeds matrix dimensions'\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 261\u001b[0m \u001b[1;32mif\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mcol\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mmax\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[0mshape\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;36m1\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", + "\u001b[1;31mKeyboardInterrupt\u001b[0m: " ] } ], "source": [ - "for dat in mt1DdataZyxList:\n", - " runInversionModel(dat,problem,m1d,'zyx')\n", + "if True:\n", + " for dat in mt1DdataZyxList:\n", + " runInversionModel(dat,problem,m1d,'zyx')\n", "\n", - "for dat in mt1DdataZxyList:\n", - " runInversionModel(dat,problem,m1d,'zxy')" + " for dat in mt1DdataZxyList:\n", + " runInversionModel(dat,problem,m1d,'zxy')\n", + " for dat in mt1DdataZdetList:\n", + " runInversionModel(dat,problem,m1d,'zdet') \n" ] }, { @@ -2524,14 +574,39 @@ }, { "cell_type": "code", - "execution_count": 24, + "execution_count": 32, + "metadata": { + "collapsed": false + }, + "outputs": [], + "source": [ + "# import multiprocessing, os, numpy as np\n", + "# from glob import glob\n", + "# import itertools as it\n", + "# import subprocess\n", + "\n", + "\n", + "# def runInv_star(args):\n", + "# return runInversionModel(*args)\n", + "\n", + "\n", + "# if __name__ == '__main__':\n", + "# pool = multiprocessing.Pool(processes=6)\n", + "# pool.map(runInversionZxy,it.izip(mt1DdataZxyList,it.repeat( (problem,m1d,'zxy') )))\n", + "# pool.join()" + ] + }, + { + "cell_type": "code", + "execution_count": 34, "metadata": { "collapsed": false }, "outputs": [], "source": [ "seogizxy1Dmods = np.vstack(([np.load(f) for f in glob('xyzmod_zxy*.npy')]))\n", - "seogizyx1Dmods = np.vstack(([np.load(f) for f in glob('xyzmod_zyx*.npy')]))" + "seogizyx1Dmods = np.vstack(([np.load(f) for f in glob('xyzmod_zyx*.npy')]))\n", + "seogizdet1Dmods = np.vstack(([np.load(f) for f in glob('xyzmod_zdet*.npy')]))" ] }, { @@ -2545,7 +620,7 @@ }, { "cell_type": "code", - "execution_count": 29, + "execution_count": 35, "metadata": { "collapsed": false, "scrolled": true @@ -2558,24 +633,27 @@ "mod3d_zxy1dMod = griddata(seogizxy1Dmods[:,0:3],seogizxy1Dmods[:,3],mesh3d.gridCC,'linear')\n", "mod3d_zxy1dMod[airInd] = np.log(1e-8)\n", "mod3d_zyx1dMod = griddata(seogizyx1Dmods[:,0:3],seogizyx1Dmods[:,3],mesh3d.gridCC,'linear')\n", - "mod3d_zyx1dMod[airInd] = np.log(1e-8)" + "mod3d_zyx1dMod[airInd] = np.log(1e-8)\n", + "mod3d_zdet1dMod = griddata(seogizdet1Dmods[:,0:3],seogizdet1Dmods[:,3],mesh3d.gridCC,'linear')\n", + "mod3d_zdet1dMod[airInd] = np.log(1e-8)" ] }, { "cell_type": "code", - "execution_count": 31, + "execution_count": 36, "metadata": { "collapsed": false }, "outputs": [], "source": [ "simpeg.Utils.meshutils.writeVTRFile('model3D_zxy1Dinv_paddedMesh.vtr',mesh3d,{'S/m':np.exp(mod3d_zxy1dMod)})\n", - "simpeg.Utils.meshutils.writeVTRFile('model3D_zyx1Dinv_paddedMesh.vtr',mesh3d,{'S/m':np.exp(mod3d_zyx1dMod)})" + "simpeg.Utils.meshutils.writeVTRFile('model3D_zyx1Dinv_paddedMesh.vtr',mesh3d,{'S/m':np.exp(mod3d_zyx1dMod)})\n", + "simpeg.Utils.meshutils.writeVTRFile('model3D_zdet1Dinv_paddedMesh.vtr',mesh3d,{'S/m':np.exp(mod3d_zdet1dMod)})" ] }, { "cell_type": "code", - "execution_count": 68, + "execution_count": 52, "metadata": { "collapsed": false }, @@ -2590,13 +668,160 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 50, "metadata": { "collapsed": false }, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "[(794.79531859664246, 56.319860569098076),\n", + " (674.66872912184408, 57.635385492661392),\n", + " (559.52741397235366, 58.423832448812171),\n", + " (459.43398962262597, 58.410853192444343),\n", + " (380.28287396790699, 57.529223784875711),\n", + " (322.8521167724395, 56.001248998922648),\n", + " (283.43537890603488, 54.287024515800717),\n", + " (255.71182522194459, 52.81526241219477),\n", + " (233.73967183794295, 51.557597236493933),\n", + " (216.97743323521141, 50.281238562685871),\n", + " (204.41678397612844, 49.05493425669367),\n", + " (196.11922783991341, 47.732975834266547),\n", + " (192.81620390650565, 46.626437679858043),\n", + " (191.83932042041988, 45.994944912824344),\n", + " (191.13999703648773, 45.673623213889364),\n", + " (190.40201746939533, 45.487392219519187),\n", + " (189.72482164986508, 45.362566052988335)]" + ] + }, + "execution_count": 50, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ - "np.exp(m_0)" + "simpegmt.Utils.dataUtils.getAppRes(problem.dataPair(problem.survey,problem.survey.dpred(mopt)))" + ] + }, + { + "cell_type": "code", + "execution_count": 53, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "model3D_zdet1Dinv_paddedMesh.vtr xyzmod_zxy_300_500.npy\r\n", + "model3D_zxy1Dinv_paddedMesh.vtr xyzmod_zxy_-300_-700.npy\r\n", + "model3D_zyx1Dinv_paddedMesh.vtr xyzmod_zxy_-300_700.npy\r\n", + "modelTDpaddedMesh.vtr xyzmod_zxy_300_-700.npy\r\n", + "modelTDuniMesh.vtr xyzmod_zxy_300_700.npy\r\n", + "MT3DforData1Dinv. xyzmod_zxy_-500_-100.npy\r\n", + "MT3DforData1Dinv.ipynb xyzmod_zxy_-500_100.npy\r\n", + "MT3DforData1Dinv.py xyzmod_zxy_500_-100.npy\r\n", + "seogiModel_dobs.npy xyzmod_zxy_500_100.npy\r\n", + "seogiModel_MTdata.npy xyzmod_zxy_-500_-300.npy\r\n", + "simpegTDmodel.con xyzmod_zxy_-500_300.npy\r\n", + "xyzmod_zdet_-100_-100.npy xyzmod_zxy_500_-300.npy\r\n", + "xyzmod_zdet_-100_100.npy xyzmod_zxy_500_300.npy\r\n", + "xyzmod_zdet_100_-100.npy xyzmod_zxy_-500_-500.npy\r\n", + "xyzmod_zdet_100_100.npy xyzmod_zxy_-500_500.npy\r\n", + "xyzmod_zdet_-100_-300.npy xyzmod_zxy_500_-500.npy\r\n", + "xyzmod_zdet_-100_300.npy xyzmod_zxy_500_500.npy\r\n", + "xyzmod_zdet_100_-300.npy xyzmod_zxy_-500_-700.npy\r\n", + "xyzmod_zdet_100_300.npy xyzmod_zxy_-500_700.npy\r\n", + "xyzmod_zdet_-100_-500.npy xyzmod_zxy_500_-700.npy\r\n", + "xyzmod_zdet_-100_500.npy xyzmod_zxy_500_700.npy\r\n", + "xyzmod_zdet_100_-500.npy xyzmod_zxy_-700_-100.npy\r\n", + "xyzmod_zdet_100_500.npy xyzmod_zxy_-700_100.npy\r\n", + "xyzmod_zdet_-100_-700.npy xyzmod_zxy_700_-100.npy\r\n", + "xyzmod_zdet_-100_700.npy xyzmod_zxy_700_100.npy\r\n", + "xyzmod_zdet_100_-700.npy xyzmod_zxy_-700_-300.npy\r\n", + "xyzmod_zdet_100_700.npy xyzmod_zxy_-700_300.npy\r\n", + "xyzmod_zdet_-300_-100.npy xyzmod_zxy_700_-300.npy\r\n", + "xyzmod_zdet_-300_100.npy xyzmod_zxy_700_300.npy\r\n", + "xyzmod_zdet_300_-100.npy xyzmod_zxy_-700_-500.npy\r\n", + "xyzmod_zdet_300_100.npy xyzmod_zxy_-700_500.npy\r\n", + "xyzmod_zdet_-300_-300.npy xyzmod_zxy_700_-500.npy\r\n", + "xyzmod_zdet_-300_300.npy xyzmod_zxy_700_500.npy\r\n", + "xyzmod_zdet_300_-300.npy xyzmod_zxy_-700_-700.npy\r\n", + "xyzmod_zdet_300_300.npy xyzmod_zxy_-700_700.npy\r\n", + "xyzmod_zdet_-300_-500.npy xyzmod_zxy_700_-700.npy\r\n", + "xyzmod_zdet_-300_500.npy xyzmod_zxy_700_700.npy\r\n", + "xyzmod_zdet_300_-500.npy xyzmod_zyx_-100_-100.npy\r\n", + "xyzmod_zdet_300_500.npy xyzmod_zyx_-100_100.npy\r\n", + "xyzmod_zdet_-300_-700.npy xyzmod_zyx_100_-100.npy\r\n", + "xyzmod_zdet_-300_700.npy xyzmod_zyx_100_100.npy\r\n", + "xyzmod_zdet_300_-700.npy xyzmod_zyx_-100_-300.npy\r\n", + "xyzmod_zdet_300_700.npy xyzmod_zyx_-100_300.npy\r\n", + "xyzmod_zdet_-500_-100.npy xyzmod_zyx_100_-300.npy\r\n", + "xyzmod_zdet_-500_100.npy xyzmod_zyx_100_300.npy\r\n", + "xyzmod_zdet_500_-100.npy xyzmod_zyx_-100_-500.npy\r\n", + "xyzmod_zdet_500_100.npy xyzmod_zyx_-100_500.npy\r\n", + "xyzmod_zdet_-500_-300.npy xyzmod_zyx_100_-500.npy\r\n", + "xyzmod_zdet_-500_300.npy xyzmod_zyx_100_500.npy\r\n", + "xyzmod_zdet_500_-300.npy xyzmod_zyx_-100_-700.npy\r\n", + "xyzmod_zdet_500_300.npy xyzmod_zyx_-100_700.npy\r\n", + "xyzmod_zdet_-500_-500.npy xyzmod_zyx_100_-700.npy\r\n", + "xyzmod_zdet_-500_500.npy xyzmod_zyx_100_700.npy\r\n", + "xyzmod_zdet_500_-500.npy xyzmod_zyx_-300_-100.npy\r\n", + "xyzmod_zdet_500_500.npy xyzmod_zyx_-300_100.npy\r\n", + "xyzmod_zdet_-500_-700.npy xyzmod_zyx_300_-100.npy\r\n", + "xyzmod_zdet_-500_700.npy xyzmod_zyx_300_100.npy\r\n", + "xyzmod_zdet_500_-700.npy xyzmod_zyx_-300_-300.npy\r\n", + "xyzmod_zdet_500_700.npy xyzmod_zyx_-300_300.npy\r\n", + "xyzmod_zdet_-700_-100.npy xyzmod_zyx_300_-300.npy\r\n", + "xyzmod_zdet_-700_100.npy xyzmod_zyx_300_300.npy\r\n", + "xyzmod_zdet_700_-100.npy xyzmod_zyx_-300_-500.npy\r\n", + "xyzmod_zdet_700_100.npy xyzmod_zyx_-300_500.npy\r\n", + "xyzmod_zdet_-700_-300.npy xyzmod_zyx_300_-500.npy\r\n", + "xyzmod_zdet_-700_300.npy xyzmod_zyx_300_500.npy\r\n", + "xyzmod_zdet_700_-300.npy xyzmod_zyx_-300_-700.npy\r\n", + "xyzmod_zdet_700_300.npy xyzmod_zyx_-300_700.npy\r\n", + "xyzmod_zdet_-700_-500.npy xyzmod_zyx_300_-700.npy\r\n", + "xyzmod_zdet_-700_500.npy xyzmod_zyx_300_700.npy\r\n", + "xyzmod_zdet_700_-500.npy xyzmod_zyx_-500_-100.npy\r\n", + "xyzmod_zdet_700_500.npy xyzmod_zyx_-500_100.npy\r\n", + "xyzmod_zdet_-700_-700.npy xyzmod_zyx_500_-100.npy\r\n", + "xyzmod_zdet_-700_700.npy xyzmod_zyx_500_100.npy\r\n", + "xyzmod_zdet_700_-700.npy xyzmod_zyx_-500_-300.npy\r\n", + "xyzmod_zdet_700_700.npy xyzmod_zyx_-500_300.npy\r\n", + "xyzmod_zxy_-100_-100.npy xyzmod_zyx_500_-300.npy\r\n", + "xyzmod_zxy_-100_100.npy xyzmod_zyx_500_300.npy\r\n", + "xyzmod_zxy_100_-100.npy xyzmod_zyx_-500_-500.npy\r\n", + "xyzmod_zxy_100_100.npy xyzmod_zyx_-500_500.npy\r\n", + "xyzmod_zxy_-100_-300.npy xyzmod_zyx_500_-500.npy\r\n", + "xyzmod_zxy_-100_300.npy xyzmod_zyx_500_500.npy\r\n", + "xyzmod_zxy_100_-300.npy xyzmod_zyx_-500_-700.npy\r\n", + "xyzmod_zxy_100_300.npy xyzmod_zyx_-500_700.npy\r\n", + "xyzmod_zxy_-100_-500.npy xyzmod_zyx_500_-700.npy\r\n", + "xyzmod_zxy_-100_500.npy xyzmod_zyx_500_700.npy\r\n", + "xyzmod_zxy_100_-500.npy xyzmod_zyx_-700_-100.npy\r\n", + "xyzmod_zxy_100_500.npy xyzmod_zyx_-700_100.npy\r\n", + "xyzmod_zxy_-100_-700.npy xyzmod_zyx_700_-100.npy\r\n", + "xyzmod_zxy_-100_700.npy xyzmod_zyx_700_100.npy\r\n", + "xyzmod_zxy_100_-700.npy xyzmod_zyx_-700_-300.npy\r\n", + "xyzmod_zxy_100_700.npy xyzmod_zyx_-700_300.npy\r\n", + "xyzmod_zxy_-300_-100.npy xyzmod_zyx_700_-300.npy\r\n", + "xyzmod_zxy_-300_100.npy xyzmod_zyx_700_300.npy\r\n", + "xyzmod_zxy_300_-100.npy xyzmod_zyx_-700_-500.npy\r\n", + "xyzmod_zxy_300_100.npy xyzmod_zyx_-700_500.npy\r\n", + "xyzmod_zxy_-300_-300.npy xyzmod_zyx_700_-500.npy\r\n", + "xyzmod_zxy_-300_300.npy xyzmod_zyx_700_500.npy\r\n", + "xyzmod_zxy_300_-300.npy xyzmod_zyx_-700_-700.npy\r\n", + "xyzmod_zxy_300_300.npy xyzmod_zyx_-700_700.npy\r\n", + "xyzmod_zxy_-300_-500.npy xyzmod_zyx_700_-700.npy\r\n", + "xyzmod_zxy_-300_500.npy xyzmod_zyx_700_700.npy\r\n", + "xyzmod_zxy_300_-500.npy\r\n" + ] + } + ], + "source": [ + "ls" ] }, { @@ -2694,13 +919,40 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 42, "metadata": { "collapsed": false }, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "[(979.21054686559228, 60.117584762903384),\n", + " (742.46348350358323, 60.178574585305192),\n", + " (575.52692078307155, 59.504482832217768),\n", + " (448.01468147227973, 57.739339607718549),\n", + " (380.47444523645635, 58.103945329921899),\n", + " (329.83191047669334, 55.337196869237211),\n", + " (289.04291928642527, 54.655448967096305),\n", + " (266.75073534275322, 53.22914208405799),\n", + " (233.25010687113831, 51.415310785084358),\n", + " (235.86124553266205, 49.932520459157615),\n", + " (204.17536410665218, 49.937485696638511),\n", + " (216.58363015060209, 48.71243562166476),\n", + " (213.79479472770714, 45.446839687730105),\n", + " (194.30255638458473, 46.234641994868568),\n", + " (198.73561291019692, 47.04819962721875),\n", + " (206.89598365667399, 44.89666240461478),\n", + " (213.17403683446534, 43.140057028057178)]" + ] + }, + "execution_count": 42, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ - "simpegmt.Utils.dataUtils.getAppRes(mt1DdataList[0])" + "simpegmt.Utils.dataUtils.getAppRes(mt1DdataZdetList[0])" ] }, { @@ -2720,18 +972,6 @@ "display_name": "Python 2", "language": "python", "name": "python2" - }, - "language_info": { - "codemirror_mode": { - "name": "ipython", - "version": 2 - }, - "file_extension": ".py", - "mimetype": "text/x-python", - "name": "python", - "nbconvert_exporter": "python", - "pygments_lexer": "ipython2", - "version": "2.7.10" } }, "nbformat": 4, diff --git a/simpegMT/BaseMT.py b/simpegMT/BaseMT.py index 10f7360b..6105d106 100644 --- a/simpegMT/BaseMT.py +++ b/simpegMT/BaseMT.py @@ -26,7 +26,6 @@ class BaseMTProblem(BaseFDEMProblem): # Remove fields that are not needed del odict['hook'] del odict['setKwargs'] - del odict['PropMap'] # Return the dict return odict @@ -123,8 +122,6 @@ class BaseMTProblem(BaseFDEMProblem): # Get the adjoint projectFieldsDeriv # PTv needs to be nE, PTv = rx.projectFieldsDeriv(src, self.mesh, u, mkvc(v[src, rx],2), adjoint=True) # wrt u, need possibility wrt m - # Need to reshape PTv - PTv = np.hstack((mkvc(PTv[:len(PTv)/2],2),mkvc(PTv[len(PTv)/2::],2))) # Get the dA_duIT = ATinv * PTv dA_dmT = self.getADeriv_m(freq, u_src, mkvc(dA_duIT), adjoint=True) diff --git a/simpegMT/SurveyMT.py b/simpegMT/SurveyMT.py index 4fe906b9..6c08264c 100644 --- a/simpegMT/SurveyMT.py +++ b/simpegMT/SurveyMT.py @@ -228,7 +228,6 @@ class RxMT(Survey.BaseRx): # Calculate the complex derivative PDeriv_complex = Hd * (ZijN_uV - Zij * Hd_uV ) - # ero # Extract the real number for the real/imag components. Pv = np.array(getattr(PDeriv_complex, real_or_imag)) elif adjoint: @@ -299,7 +298,8 @@ class RxMT(Survey.BaseRx): # Calculate the complex derivative PDeriv_real = ZijN_uV(aHd*v) - aHd_uV(Zij.T*aHd*v)# - # NOTE: .toarray() is to return a non-sparse array which is needed for for Ainv* operation. Might want to take care of this elsewhere. + # NOTE: Need to reshape the output to go from 2*nU array to a (nU,2) matrix for each polarization + PDeriv_real = PDeriv_real.reshape((mesh.nE,2)) # Extract the data if real_or_imag == 'imag': Pv = 1j*PDeriv_real