'rho_min': -1e3,
'rho_max': 1e3,
}
inv_r2_in_dict = {}
inv_r2_in_dict.update(r2_in_dict)
inv_r2_in_dict.update({
'reg_elems': [[1, len(mr2.mesh_dict['elements']), main_res]],
'res_matrix':
1,
'inv_dict':
inv_options
})
# Load forward model data
meas_fwd_data = pyres_utils.load_fwd_output(work_dir=work_dir)
inv_protocol_dict = {'meas_data': meas_fwd_data}
run_inv_r2_dict = {
'protocol_dict': inv_protocol_dict,
'r2_in_dict': inv_r2_in_dict,
'run_bool': True
}
iR2_inv.run_all(**run_inv_r2_dict)
#%% ######### Plot gmsh mesh ################
#m_fig,m_ax,m_tri_obj = plot_utils.plot_tri_mesh(mr2.mesh_dict)
#%%
inv_r2_in_dict = {}
inv_r2_in_dict.update(r2_in_dict)
inv_r2_in_dict.update({
'reg_elems':
[[1,
int((xx_mesh.shape[0] - 1) * (yy_mesh.shape[0] - 1)), main_res * 10.]],
'res_matrix':
1,
'inv_dict':
inv_options,
'output_domain': []
})
# Load forward model data
meas_fwd_data = pyres_utils.load_fwd_output(work_dir=out_fwd_dirs[0])
inv_protocol_dict = {'meas_data': meas_fwd_data}
run_inv_r2_dict = {
'protocol_dict': inv_protocol_dict,
'r2_in_dict': inv_r2_in_dict,
'run_bool': True
}
iR2_inv.run_all(**run_inv_r2_dict)
#%% Run difference inversion
if not os.path.isdir(inv_work_dir):
os.makedirs(inv_work_dir)
inv_r2_in_dict = {}
inv_r2_in_dict.update(r2_in_dict)
inv_r2_in_dict.update({
'reg_elems':
[[1,
int((xx_mesh.shape[0] - 1) * (yy_mesh.shape[0] - 1)), main_res * 10.]],
'res_matrix':
1,
'inv_dict':
inv_options,
'output_domain': []
})
# Load forward model data
meas_fwd_data = pyres_utils.load_fwd_output(work_dir=out_fwd_dirs[0])
inv_protocol_dict = {'meas_data': meas_fwd_data}
run_inv_r2_dict = {
'protocol_dict': inv_protocol_dict,
'r2_in_dict': inv_r2_in_dict,
'run_bool': True
}
iR2_inv.run_all(**run_inv_r2_dict)
#%% Run difference inversion
if not os.path.isdir(inv_work_dir):
os.makedirs(inv_work_dir)
