A Two-Level Block Preconditioned Jacobi–Davidson Method for Multiple and Clustered Eigenvalues of Elliptic Operators

Qigang Liang,Wei Wang,Xuejun Xu
DOI: https://doi.org/10.1137/23m1580711
IF: 3.039
2024-04-23
SIAM Journal on Numerical Analysis
Abstract:SIAM Journal on Numerical Analysis, Volume 62, Issue 2, Page 998-1019, April 2024. In this paper, we propose a two-level block preconditioned Jacobi–Davidson (BPJD) method for efficiently solving discrete eigenvalue problems resulting from finite element approximations of [math]th ([math]) order symmetric elliptic eigenvalue problems. Our method works effectively to compute the first several eigenpairs, including both multiple and clustered eigenvalues with corresponding eigenfunctions, particularly. The method is highly parallelizable by constructing a new and efficient preconditioner using an overlapping domain decomposition (DD). It only requires computing a couple of small scale parallel subproblems and a quite small scale eigenvalue problem per iteration. Our theoretical analysis reveals that the convergence rate of the method is bounded by [math], where [math] is the diameter of subdomains and [math] is the overlapping size among subdomains. The constant [math] is independent of the mesh size [math] and the internal gaps among the target eigenvalues, demonstrating that our method is optimal and cluster robust. Meanwhile, the [math]-dependent constant [math] decreases monotonically to 1, as [math], which means that more subdomains lead to the better convergence rate. Numerical results supporting our theory are given.
mathematics, applied
What problem does this paper attempt to address?