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
