A parallel sweeping preconditioner for frequency-domain seismic wave propagation
MetadataShow full item record
AbstractWe present a parallel implementation of Engquist and Ying's sweeping preconditioner, which exploits radiation boundary conditions in order to form an approximate block LDLT factorization of the Helmholtz operator with only O(N4/3) work and an application (and memory) cost of only O(N logN). The approximate factorization is then used as a preconditioner for GMRES, and we show that essentially O(1) iterations are required for convergence, even for the full SEG/EAGE over-thrust model at 30 Hz. In particular, we demonstrate the solution of said problem in a mere 15 minutes on 8192 cores of TACC's Lonestar, which may be the largest-scale 3D heterogeneous Helmholtz calculation to date. Generalizations of our parallel strategy are also briefly discussed for time-harmonic linear elasticity and Maxwell's equations.
CitationPoulson J, Engquist B, Li S, Ying L (2012) A parallel sweeping preconditioner for frequency-domain seismic wave propagation. SEG Technical Program Expanded Abstracts 2012. Available: http://dx.doi.org/10.1190/segam2012-0765.1.
SponsorsThe authors thank TACC for access to Lonestar, and particularly Bill Barth for suggesting that we enforce local memory allocations and Tommy Minyard for organizing the large runs. We also thank Sergey Fomel, Anshul Gupta, Edgar Solomonik, and Laura Grigori for helpful discussions. The work is financially supported by KAUST and the TCCS sponsors.
PublisherSociety of Exploration Geophysicists