A Parallel Sweeping Preconditioner for Heterogeneous 3D Helmholtz Equations.

Jack Poulson,Bjoern Engquist,Siwei Li,Lexing Ying
DOI: https://doi.org/10.1137/120871985
IF: 2.968
2013-01-01
SIAM Journal on Scientific Computing
Abstract:A parallelization of a sweeping preconditioner for three-dimensional Helmholtz equations without large cavities is introduced and benchmarked for several challenging velocity models. The setup and application costs of the sequential preconditioner are shown to be $O(\gamma^2 N^{4/3})$ and $O(\gamma N \log N)$, where $\gamma(\omega)$ denotes the modestly frequency-dependent number of grid points per perfectly matched layer. Several computational and memory improvements are introduced relative to using black-box sparse-direct solvers for the auxiliary problems, and competitive runtimes and iteration counts are reported for high-frequency problems distributed over thousands of cores. Two open-source packages are released along with this paper: Parallel Sweeping Preconditioner (PSP) and the underlying distributed multifrontal solver, Clique.
What problem does this paper attempt to address?