A Novel Deflation Technique for Solving Quadratic Eigenvalue Problems
by
Tsung-Ming Huang
Wen-Wei Lin
Vol. 9 No. 1 (2014) P.57~P.84
ABSTRACT
In this paper we propose numerical algorithms for solving large-scale quadratic eigenvalue problems for which a set of eigenvalues closest to a fixed target and the associated
eigenvectors are of interest. The desired eigenvalues are usually with smallest modulo in the spectrum. The algorithm based on the
quadratic Jacobi-Davidson (QJD) algorithm is proposed to find the first smallest eigenvalue closest to the target. To find the
successive eigenvalues closest to the target, we propose a novel explicit non-equivalence low-rank deflation technique. The technique transforms the smallest eigenvalue to infinity, while all other eigenvalues remain unchanged. Thus, the original second
smallest eigenvalue becomes the smallest of the new quadratic eigenvalue problem, which can then be solved by the $QJD$ algorithm. To compare with locking and restarting quadratic
eigensolver, our numerical experience shows that the QJD method combined with our explicit non-equivalence deflation is robust and efficient.
KEYWORDS
Quadratic eigenvalue problems, quadratic Jacobi-Davidson method, non-euqivalence low-rank deflation
MATHEMATICAL SUBJECT CLASSIFICATION 2010
Primary: 15A18, 47A15, 47J10
MILESTONES
Received: 2013-06-24
Revised : 2013-10-19
Accepted: 2013-10-28
Download Full Content