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