define a new window_type, synthetic, for synthetic test. - Githubissues

uafgeotools / mtuq

moment tensor uncertainty quantification

BSD 2-Clause "Simplified" License

65 stars 22 forks source link

define a new window_type, synthetic, for synthetic test. #269

Open Liang-Ding opened 4 weeks ago

Liang-Ding commented 4 weeks ago

We define a new 'window_type' named 'synthetic' for synthetic examples in our HMC package, seisHMC, since it relies on the functions in MTUQ. Here are the examples using synthetic waveforms:

GridSearch + DoubleCouple: https://github.com/Liang-Ding/seishmc/blob/main/examples/GridSearch.SYN.DoubleCouple.py
HMC + DoubleCouple: https://github.com/Liang-Ding/seishmc/blob/main/examples/HMC.SYN.DoubleCouple.py
GridSearch + FMT: https://github.com/Liang-Ding/seishmc/blob/main/examples/GridSearch.SYN.FullMomentTensor.py
HMC + FMT: https://github.com/Liang-Ding/seishmc/blob/main/examples/HMC.SYN.FullMomentTensor.py We achieve lower waveform misfits in our examples using HMC algorithm. Check the figure: misfit comparison

Additional results are avaliable at: https://github.com/Liang-Ding/seishmc/blob/main/doc/gallery.md.