A parallel sweeping preconditioner for heterogeneous 3D Helmholtz equations
arXiv:1204.0111 · doi:10.1137/120871985
Abstract
A parallelization of a sweeping preconditioner for 3D 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(γ^2 N^{4/3}) and O(γ N log N), where γ(ω) 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".
Final version of document
References in corpus (2)
Cited by in corpus (20)
- Detection of lensing substructure using ALMA observations of the dusty galaxy SDP.81
- A dispersion minimizing scheme for the 3-D Helmholtz equation based on ray theory
- L-Sweeps: A scalable, parallel preconditioner for the high-frequency Helmholtz equation
- An improved sweeping domain decomposition preconditioner for the Helmholtz equation
- Schwarz methods by domain truncation
- Recursive Sweeping Preconditioner for the 3D Helmholtz Equation
- A time-domain preconditioner for the Helmholtz equation
- Distributed-memory Hierarchical Interpolative Factorization
- Multigrid-Augmented Deep Learning Preconditioners for the Helmholtz Equation using Compact Implicit Layers
- Parallel accelerated cyclic reduction preconditioner for three-dimensional elliptic PDEs with variable coefficients
- An iterative domain decomposition, spectral finite element method on non-conforming meshes suitable for high frequency Helmholtz problems
- Domain Decomposition with local impedance conditions for the Helmholtz equation with absorption
- A Class of Iterative Solvers for the Helmholtz Equation: Factorizations, Sweeping Preconditioners, Source Transfer, Single Layer Potentials, Polarized Traces, and Optimized Schwarz Methods
- A hybrid shifted Laplacian multigrid and domain decomposition preconditioner for the elastic Helmholtz equations
- Sparsify and sweep: an efficient preconditioner for the Lippmann-Schwinger equation
- Additive Sweeping Preconditioner for the Helmholtz Equation
- A short note on the nested-sweep polarized traces method for the 2D Helmholtz equation
- Exponentially convergent multiscale methods for high frequency heterogeneous Helmholtz equations
- Domain decomposition preconditioners for high-order discretisations of the heterogeneous Helmholtz equation
- The method of polarized traces for the 2D Helmholtz equation