1887

Abstract

Summary

Frequency-domain elastic wave modeling relies on an efficient linear solver for the large, sparse, ill-conditioned linear system derived from the discretization of the elastic wave equation. Direct solvers are mostly based on LU decomposition. These methods are efficient for multiple right-hand sides problems, but they require significant memory resources. Conversely, iterative solvers benefit from the sparsity of the system, but they require sophisticated preconditioners to converge due to the system ill-conditioning. In this study, we investigate the performance of an iterative method named CARP-CG for frequency-domain elastic wave modeling. The CARP-CG method transforms the original system into a symmetric positive semi-definite system by cyclic row-projections. This system is efficiently solved with the conjugate gradient (CG) method. The cyclic row-projection transformation can be seen as a purely algebraic preconditioning technique which is easy to implement. The algorithm can be parallelized through a row-block decomposition combined with component-averaging operations. Numerical experiments on the 2D frequency-domain elastic problem with high Poisson’s ratio exhibit a good scalability of CARP-CG. Comparisons for different frequencies between CARP-CG and standard Krylov iterative solvers (GMRES and CG on the normal equations) emphasize the robustness and the fast convergence of the method.

Loading

Article metrics loading...

/content/papers/10.3997/2214-4609.20141633
2014-06-16
2024-04-28
Loading full text...

Full text loading...

References

  1. Amestoy, P., Duff, I.S. and L’Excellent, J.Y.
    [2000] Multifrontal parallel distributed symmetric and unsymmetric solvers. Computer Methods in Applied Mechanics and Engineering, 184(2–4), 501–520.
    [Google Scholar]
  2. Erlangga, Y.A. and Nabben, R.
    [2008] On a multilevel Krylov method for the Helmholtz equation preconditioned by shifted Laplacian. Electronic Transactions on Numerical Analysis, 31, 403–424.
    [Google Scholar]
  3. Frayssé, V., Giraud, L., Gratton, S. and Langou, J.
    [2003] A set of GMRES routines for real and complex arithmetics on high performance computers. Tech. Rep. 3, CeRFACS, Technical Report TR/PA/03/3.
    [Google Scholar]
  4. Gordon, D. and Gordon, R.
    [2008] CGMN revisited: Robust and efficient solution of stiff linear systems derived from elliptic partial differential equations. ACM Trans. on Mathematical Software, 35(3), 18:1–18:27.
    [Google Scholar]
  5. [2010] CARP-CG: A robust and efficient parallel solver for linear systems, applied to strongly convection dominated PDEs. Parallel Computing, 36, 495–515.
    [Google Scholar]
  6. [2013] Robust and highly scalable parallel solution of the Helmholtz equation with large wave numbers. Journal of Computational and Applied Mathematics, 237(1), 182–196.
    [Google Scholar]
  7. Kaczmarz, S.
    [1937] Angenäherte Auflösung von Systemen linearer Gleichungen (English translation by Jason Stockmann). Bulletin International de l’Académie Polonaise des Sciences et des Lettres, 35, 355–357.
    [Google Scholar]
  8. Saad, Y.
    [2003] Iterative methods for sparse linear systems. SIAM, Philadelphia.
    [Google Scholar]
  9. van Leeuwen, T., Gordon, D., Gordon, R. and Herrmann, F.
    [2012] Preconditioning the Helmholtz equations via row projections. Expanded Abstracts, EAGE, A002.
    [Google Scholar]
  10. Virieux, J.
    [1986] P-SV wave propagation in heterogeneous media: velocity stress finite difference method. Geophysics, 51, 889–901.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/papers/10.3997/2214-4609.20141633
Loading
/content/papers/10.3997/2214-4609.20141633
Loading

Data & Media loading...

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error