TY - JOUR

T1 - A parallel two-level polynomial Jacobi–Davidson algorithm for large sparse PDE eigenvalue problems

AU - Cheng, Yu Fen

AU - Hwang, Feng Nan

N1 - Publisher Copyright:
© 2017 Elsevier Ltd

PY - 2017/10

Y1 - 2017/10

N2 - Many scientific and engineering applications require accurate, fast, robust, and scalable numerical solution of large sparse algebraic polynomial eigenvalue problems (PEVP's) that arise from some appropriate discretization of partial differential equations. The polynomial Jacobi-Davidson (PJD) algorithm has been numerically shown as a promising approach for the PEVP's to finding the interior spectrum. The PJD algorithm is a subspace method, which extracts the candidate eigenpair from a search space and the space updated by embedding the solution of the correction equation at the JD iteration. In this research, we develop and study the two-level PJD algorithm for PEVP's with emphasis on the application of the dissipative acoustic cubic eigenvalue problem. The proposed two-level PJD algorithm consists of two important ingredients: A good initial basis for the search space is constructed on the fine-level by using the interpolation of the coarse solution of the same eigenvalue problem in order to enhance the robustness of the algorithm. Also, an efficient and scalable two-level preconditioner based on the Schwarz framework is used for the correction equation. Some numerical examples obtained on a parallel cluster of computers are given in order to demonstrate the robustness and scalability of our PJD algorithm.

AB - Many scientific and engineering applications require accurate, fast, robust, and scalable numerical solution of large sparse algebraic polynomial eigenvalue problems (PEVP's) that arise from some appropriate discretization of partial differential equations. The polynomial Jacobi-Davidson (PJD) algorithm has been numerically shown as a promising approach for the PEVP's to finding the interior spectrum. The PJD algorithm is a subspace method, which extracts the candidate eigenpair from a search space and the space updated by embedding the solution of the correction equation at the JD iteration. In this research, we develop and study the two-level PJD algorithm for PEVP's with emphasis on the application of the dissipative acoustic cubic eigenvalue problem. The proposed two-level PJD algorithm consists of two important ingredients: A good initial basis for the search space is constructed on the fine-level by using the interpolation of the coarse solution of the same eigenvalue problem in order to enhance the robustness of the algorithm. Also, an efficient and scalable two-level preconditioner based on the Schwarz framework is used for the correction equation. Some numerical examples obtained on a parallel cluster of computers are given in order to demonstrate the robustness and scalability of our PJD algorithm.

KW - Acoustic wave equation

KW - Cubic eigenvalue problems

KW - Domain decomposition

KW - Jacobi–Davidson methods

KW - Parallel computing

KW - Two-level Schwarz preconditioner

UR - http://www.scopus.com/inward/record.url?scp=85019579960&partnerID=8YFLogxK

U2 - 10.1016/j.advengsoft.2017.05.011

DO - 10.1016/j.advengsoft.2017.05.011

M3 - 期刊論文

AN - SCOPUS:85019579960

SN - 0965-9978

VL - 112

SP - 222

EP - 230

JO - Advances in Engineering Software

JF - Advances in Engineering Software

ER -