示例#1
0
    def test_StaticInversion(self):

        fault_file = join(self.share_dir, 'fault_He50km_east.h5')

        rough = Roughening.create_from_fault_file(fault_file)
        reg = Composite().add_component(rough, 1., 'roughening')

        basis = BasisMatrix.create_from_fault_file(fault_file)

        inv = StaticInversion(
            file_G=join(self.share_dir, 'G.h5'),
            file_d=join(self.share_dir, 'cumu_post_with_seafloor.h5'),
            file_sd=join(self.share_dir, 'sites_sd.h5'),
            file_sites_filter=join(self.share_dir, 'sites_with_seafloor'),
            regularization=reg,
            basis=basis,
        )

        inv.set_data_all()
        inv.run()
        inv.save(join(self.outs_dir, 'out.h5'), overwrite=True)
    def test_StaticInversion(self):        

        fault_file = join(self.share_dir, 'fault_He50km_east.h5')

        rough = Roughening.create_from_fault_file(fault_file)
        reg = Composite().add_component(rough, 1., 'roughening')

        basis = BasisMatrix.create_from_fault_file(fault_file)

        inv = StaticInversion(
            file_G = join(self.share_dir, 'G.h5'),
            file_d = join(self.share_dir,'cumu_post_with_seafloor.h5'),
            file_sd = join(self.share_dir, 'sites_sd.h5'),
            file_sites_filter = join(self.share_dir, 'sites_with_seafloor'),
            regularization = reg,
            basis = basis,
        )

        inv.set_data_all()
        inv.run()
        inv.save(join(self.outs_dir,'out.h5'), overwrite=True)
示例#3
0
from numpy import logspace

from viscojapan.inversion import OccamDeconvolution
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')
示例#4
0
from numpy import logspace

from viscojapan.inversion import OccamDeconvolution
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')