wassname/simpeg
1
0
Fork 0
mirror of https://github.com/wassname/simpeg.git synced 2026-08-03 13:10:16 +08:00
Code Issues Packages Projects Releases Wiki Activity

Get DC inversion working with the changes in the data class.

Browse Source
This commit is contained in:
rowanc1
2014-01-14 12:22:58 -08:00
parent 45988c2094
commit 5797fa4110
1 changed files with 4 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files
SimPEG/examples/DC.py
+4 -4
View File
@@ -203,7 +203,7 @@ if __name__ == '__main__':
problem = DCProblem(mesh)
problem.P = P
problem.RHS = q
dobs, Wd = problem.createSyntheticData(mSynth, std=0.05)
data = problem.createSyntheticData(mSynth, std=0.05)
u = problem.field(mSynth)
u = problem.reshapeFields(u)
@@ -211,13 +211,13 @@ if __name__ == '__main__':
# plt.show()
# Now set up the problem to do some minimization
problem.dobs = dobs
problem.std = dobs*0 + 0.05
# problem.dobs = dobs
# problem.std = dobs*0 + 0.05
m0 = mesh.gridCC[:,0]*0+sig2
opt = inverse.InexactGaussNewton(maxIterLS=20, maxIter=10, tolF=1e-6, tolX=1e-6, tolG=1e-6, maxIterCG=6)
reg = inverse.Regularization(mesh)
inv = inverse.Inversion(problem, reg, opt, beta0=1e4)
inv = inverse.Inversion(problem, reg, opt, data, beta0=1e4)
# Check Derivative
derChk = lambda m: [inv.dataObj(m), inv.dataObjDeriv(m)]