An efficient boundary integral equation method applicable to 3D dynamic rupture simulations

This paper presents a fast H-matrices-based boundary integral equation method that can be applied to various wave-related problems. We propose an efficient algorithm for convolution in time direction for the intermediate domain between P- and S-waves using the exponentiation notation of the integral kernel. Furthermore, lattice H-matrices are used to address computational inefficiencies in parallel environments due to the hierarchical structure of H-matrices. 3D stress wave propagation simulations confirm the accuracy of the method.