부산대학교 도서관

To link to full-text access for this article, visit this link: http://dx.doi.org/10.1016/j.jcp.2007.03.033 Byline: C.D. Riyanti (a), A. Kononov (b), Y.A. Erlangga (c), C. Vuik (a), C.W. Oosterlee (a), R.-E. Plessix (d), W.A. Mulder (d) Keywords: Helmholtz equation; Krylov subspace method; Preconditioner; Multigrid method Abstract: We investigate the parallel performance of an iterative solver for 3D heterogeneous Helmholtz problems related to applications in seismic wave propagation. For large 3D problems, the computation is no longer feasible on a single processor, and the memory requirements increase rapidly. Therefore, parallelization of the solver is needed. We employ a complex shifted-Laplace preconditioner combined with the Bi-CGSTAB iterative method and use a multigrid method to approximate the inverse of the resulting preconditioning operator. A 3D multigrid method with 2D semi-coarsening is employed. We show numerical results for large problems arising in geophysical applications. Author Affiliation: (a) Delft Institute of Applied Mathematics, Delft University of Technology, Mekelweg 4, 2628 CD Delft, The Netherlands (b) Computational Physics Group, Delft University of Technology, The Netherlands (c) Scientific Computing, Technical University, Berlin, Germany (d) Shell International E&P, P.O. Box 60, 2280 AB Rijswijk, The Netherlands Article History: Received 30 September 2006; Revised 27 March 2007; Accepted 29 March 2007