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


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.

Quadratic eigenvalue problems, quadratic Jacobi-Davidson method, non-euqivalence low-rank deflation

Primary: 15A18, 47A15, 47J10


Received: 2013-06-24
Revised : 2013-10-19
Accepted: 2013-10-28

Download Full Content