maxIterCG=20, tolCG=1e-3, LSalwaysPass=True, stepOffBoundsFact=1e-8) invProb = InvProblem.BaseInvProblem(dmis, reg, opt) # LIST OF DIRECTIVES # betaest = Directives.BetaEstimate_ByEig() IRLS = Directives.Update_IRLS(f_min_change=1e-6, minGNiter=2, beta_tol=1e-2, coolingRate=2) update_SensWeight = Directives.UpdateSensWeighting() update_Jacobi = Directives.UpdatePreCond() ProjSpherical = Directives.ProjSpherical() JointAmpMVI = Directives.JointAmpMVI() betaest = Directives.BetaEstimate_ByEig(beta0_ratio = 1e+3) saveModel = Directives.SaveUBCModelEveryIteration(mapping=actvMap, saveComp=True) saveModel.fileName = work_dir+out_dir + 'JOINT_MVIC_A' inv = Inversion.BaseInversion(invProb, directiveList=[betaest, IRLS, update_SensWeight, update_Jacobi, saveModel]) # Run JOINT mrec = inv.run(mstart) #NOTE - Would like to have dpred working on both surveys dpred = invProb.getFields(mrec)
tolCG=1e-3, LSalwaysPass=True, stepOffBoundsFact=1e-8) invProb = InvProblem.BaseInvProblem(dmis, reg, opt) # LIST OF DIRECTIVES # betaest = Directives.BetaEstimate_ByEig() IRLS = Directives.Update_IRLS(f_min_change=1e-6, minGNiter=2, beta_tol=1e-2, coolingRate=2) update_SensWeight = Directives.UpdateSensWeighting() update_Jacobi = Directives.UpdatePreCond() ProjSpherical = Directives.ProjSpherical() JointAmpMVI = Directives.JointAmpMVI() betaest = Directives.BetaEstimate_ByEig(beta0_ratio=1e+3) saveModel = Directives.SaveUBCModelEveryIteration(mapping=actvMap, saveComp=True) saveModel.fileName = work_dir + out_dir + 'JOINT_MVIC_A' inv = Inversion.BaseInversion(invProb, directiveList=[ betaest, JointAmpMVI, IRLS, update_SensWeight, update_Jacobi, saveModel ]) # Run JOINT mrec = inv.run(mstart) #%% Switch to JOINT MVIS-MAI