Skip to main content Accessibility help

Sweeping preconditioners for elastic wave propagation with spectral element methods

  • Paul Tsuji (a1), Jack Poulson (a2), Björn Engquist (a3) and Lexing Ying (a4)


We present a parallel preconditioning method for the iterative solution of the time-harmonic elastic wave equation which makes use of higher-order spectral elements to reduce pollution error. In particular, the method leverages perfectly matched layer boundary conditions to efficiently approximate the Schur complement matrices of a block LDLT factorization. Both sequential and parallel versions of the algorithm are discussed and results for large-scale problems from exploration geophysics are presented.



Hide All
[1] F. Aminzadeh, J. Brac and T. Kunz, 3-D Salt and Overthrust models. In SEG/EAGE 3-D Modeling Series 1. Tulsa, OK (1997).
[2] Berenger, J.P., A perfectly matched layer for the absorption of electromagnetic waves. J. Comput. Phys. 114 (1994) 185200.
[3] Bramble, J. and Pasciak, J., A note on the existence and uniqueness of solutions of frequency domain elastic wave problems: a priori estimates in H 1. J. Math Anal. Appl. 345 (2008) 396404.
[4] Chew, W.C. and Weedon, W.H., A 3D perfectly matched medium from modified Maxwell’s equations with stretched coordinates. Microwave Optical Tech. Lett. 7 (1994) 599604.
[5] J. Choi, J.J. Dongarra, R. Pozo and D.W. Walker, ScaLAPACK: A scalable linear algebra library for distributed memory concurrent computers, in Proc. of the Fourth Symposium on the Frontiers of Massively Parallel Computation, IEEE Comput. Soc. Press(1992) 120–127.
[6] Engquist, B. and Ying, L., Sweeping preconditioner for the Helmholtz equation: hierarchical matrix representation. Commun. Pure Appl. Math. 64 (2011) 697735.
[7] Engquist, B. and Ying, L., Sweeping preconditioner for the Helmholtz equation: moving perfectly matched layers. Multiscale Model. Simul. 9 (2011) 686710.
[8] Erlangga, Y.A., Vuik, C. and Oosterlee, C.W., On a class of preconditioners for solving the Helmholtz equation. Appl. Numer. Math. 50 (2004) 409425.
[9] O.G. Ernst and M.J. Gander, Why it is difficult to solve Helmholtz problems with classical iterative methods, in vol. 83 of Numerical Analysis of Multiscale Problems. Edited by I. Graham, T. Hou, O. Lakkis and R. Scheichl. Springer-Verlag (2011) 325–361.
[10] Grasedyck, L. and Hackbusch, W., Construction and arithmetics of ℋ-matrices. Computing 70 (2003) 295334.
[11] Gupta, A., Karypis, G. and Kumar, V., A highly scalable parallel algorithm for sparse matrix factorization. IEEE Transactions on Parallel and Distributed Systems 8 (1997) 502520.
[12] A. Gupta, S. Koric and T. George, Sparse matrix factorization on massively parallel computers, in Proc. of the Conference on High Performance Computing, Networking, Storage and Analysis. Portland, OR (2009).
[13] Harari, I. and Albocher, U., Studies of FE/PML for exterior problems of time-harmonic elastic waves. Comput. Methods Appl. Mech. Eng. 195 (2006) 38543879.
[14] T. Hughes, The finite element method: linear static and dynamic finite element analysis. Prentice-Hall, Inc. (1987).
[15] G. Karniadakis, Spectral/hp element methods for CFD. Oxford University Press (1999).
[16] Komatitsch, D. and Tromp, J., Introduction to the spectral element method for three-dimensional seismic wave propagation. Geophys. J. Int. 139 (1999) 806822.
[17] Liu, J.. The multifrontal method for sparse matrix solution: theory and practice. SIAM Review 34 (1992) 82109.
[18] Patera, A., A spectral element method for fluid dynamics: laminar flow in a channel expansion. J. Comput. Phys. 54 (1984) 468488.
[19] J. Poulson, B. Engquist, S. Li and L. Ying, A parallel sweeping preconditioner for heterogeneous 3D Helmholtz equations e-prints ArXiv (2012).
[20] Poulson, J., Marker, B., van de Geijn, R.A., Hammond, J.R. and Romero, N.A., Elemental: A new framework for distributed memory dense matrix computations. ACM Trans. Math. Software 39.
[21] Raghavan, P., Efficient parallel sparse triangular solution with selective inversion. Parallel Proc. Lett. 8 (1998) 2940.
[22] Saad, Y. and Schultz, M.H., A generalized minimum residual algorithm for solving nonsymmetric linear systems. SIAM J. Sci. Stat. Comput. 7 (1986) 856869.
[23] Schreiber, R., A new implementation of sparse Gaussian elimination. ACM Trans. Math. Software 8 (1982) 256276.
[24] Tsuji, P., Engquist, B. and Ying, L., A sweeping preconditioner for time-harmonic Maxwell’s equations with finite elements. J. Comput. Phys. 231 (2012) 37703783.
[25] Tsuji, P. and Ying, L., A sweeping preconditioner for Yee’s finite difference approximation of time-harmonic Maxwell’s equations. Frontiers of Mathematics in China 7 (2012) 347363.


Sweeping preconditioners for elastic wave propagation with spectral element methods

  • Paul Tsuji (a1), Jack Poulson (a2), Björn Engquist (a3) and Lexing Ying (a4)


Full text views

Total number of HTML views: 0
Total number of PDF views: 0 *
Loading metrics...

Abstract views

Total abstract views: 0 *
Loading metrics...

* Views captured on Cambridge Core between <date>. This data will be updated every 24 hours.

Usage data cannot currently be displayed.