from viscojapan.inversion.regularization import Roughening, Composite
from viscojapan.inversion.basis_function import BasisMatrix
epochs = [0]
fault_file = '../../fault_model/fault_bott60km.h5'
basis = BasisMatrix.create_from_fault_file(fault_file)
rough = Roughening.create_from_fault_file(fault_file)
inv = OccamDeconvolution(
file_G0 = '../../green_function/G0_He63km_VisM1.0E19_Rake90.h5',
files_Gs = ['../../green_function/G1_He63km_VisM1.0E19_Rake80.h5'],
nlin_par_names = ['rake'],
file_d = '../../cumu_post_with_seafloor.h5',
file_sd = '../sd/sd_seafloor_inf.h5',
file_incr_slip0 = 'slip0/slip0.h5',
filter_sites_file = 'sites_with_seafloor',
epochs = epochs,
regularization = rough,
basis = basis,
)
inv.set_data_except_L()
for nth, alpha in enumerate(logspace(-4, 0, 30)):
reg = Composite().add_component(component = rough,
arg = alpha,
arg_name = 'roughening')
inv.regularization = reg
inv.set_data_L()
inv.invert()
inv.predict()
from alphas import alphas
epochs = [0]
fault_file = '../../greens_function/fault_He50km.h5'
basis = BasisMatrix.create_from_fault_file(fault_file)
rough = Roughening.create_from_fault_file(fault_file)
inv = OccamDeconvolution(
file_G0='../../greens_function/G_Rake80.h5',
files_Gs=['../../greens_function/G_Rake90.h5'],
nlin_par_names=['rake'],
file_d='../../true_model/d_simu.h5',
file_sd='../../sites_sd/sites_sd.h5',
file_incr_slip0='slip0/m_true.h5',
filter_sites_file='sites_with_seafloor',
epochs=epochs,
regularization=rough,
basis=basis,
)
inv.set_data_except_L()
for nth, alpha in enumerate(alphas):
reg = Composite().add_component(component=rough,
arg=alpha,
arg_name='roughening')
inv.regularization = reg
inv.set_data_L()
inv.invert()
inv.predict()
from viscojapan.inversion.regularization import Roughening, Composite
from viscojapan.inversion.basis_function import BasisMatrix
epochs = [0]
fault_file = '../../fault_model/fault_bott60km.h5'
basis = BasisMatrix.create_from_fault_file(fault_file)
rough = Roughening.create_from_fault_file(fault_file)
inv = OccamDeconvolution(
file_G0='../../green_function/G0_He63km_VisM1.0E19_Rake90.h5',
files_Gs=['../../green_function/G1_He63km_VisM1.0E19_Rake80.h5'],
nlin_par_names=['rake'],
file_d='../../cumu_post_with_seafloor.h5',
file_sd='../sd/sd_seafloor_inf.h5',
file_incr_slip0='slip0/slip0.h5',
filter_sites_file='sites_with_seafloor',
epochs=epochs,
regularization=rough,
basis=basis,
)
inv.set_data_except_L()
for nth, alpha in enumerate(logspace(-4, 0, 30)):
reg = Composite().add_component(component=rough,
arg=alpha,
arg_name='roughening')
inv.regularization = reg
inv.set_data_L()
inv.invert()
inv.predict()
fault_file = '../fault_model/fault_bott40km.h5'
epochs = epochs[0:4]
basis = BasisMatrix.create_from_fault_file(fault_file, num_epochs = len(epochs))
inv = OccamDeconvolution(
file_G0 = '../green_function/G0_He40km_Vis5.8E18_Rake81.h5',
files_Gs = ['../green_function/G1_He40km_Vis1.1E19_Rake81.h5',
'../green_function/G2_He45km_Vis5.8E18_Rake81.h5',
'../green_function/G3_He40km_Vis5.8E18_Rake90.h5'
],
nlin_par_names = ['log10(visM)','log10(He)','rake'],
file_d = '../obs/cumu_post_with_seafloor.h5',
file_sd = 'sd/sd_uniform_seafloor_inf.h5',
file_incr_slip0 = 'slip0/incr_slip0.h5',
filter_sites_file = 'sites_2EXPs',
epochs = epochs,
regularization = None,
basis = basis,
)
inv.set_data_except(excepts=['L'])
reg_temp = 0
reg_edge = 0.06
for nrough, reg_rough in enumerate(reg_roughes):
outfname = 'outs/nrough_%02d.h5'%(nrough)