Program codes for Bayesian simultaneous fault and slip inversion in heterogeneous crustal structures

This repository provides programs for estimating fault geometry and slip distribution of an earthquake from crustal deformation data. The Hamiltonian Monte Carlo method is used by integrating the formulation of dislocation potentials and differential ensemble modeling using physics-informed neural networks. The programs are written in Python using PyTorch and Pyro.