A Refined Second-order Arnoldi (RSOAR) Method for the Quadratic Eigenvalue Problem and Implicit Restarting

Zhongxiao Jia,Yuquan Sun
2010-01-01
Abstract:Based on the orthogonal projection principle and a second-order Arnoldi procedure that generates an orthonormal basis of the second-order Krylov subspace, the second-order Arnoldi (SOAR) method has been proposed for the quadratic eigenvalue problem. It can be used to compute some eigenpairs of the quadratic eigenvalue problem and preserve physical structures of the original problem. This paper proposes a refined second-order Arnoldi (RSOAR) method that is based on the refined projection principle and has the same structure preserving properties. To make the methods practical, restarting is generally necessary. However, the SOAR based methods cannot be restarted effectively. As the first step, to correct this serious deficiency, we instead use a modified SOAR procedure starting with a general vector. It is well known that if deflation occurs then implicit restarting cannot be applied to the MSOAR procedure. To this end, we propose an effective strategy to treat the critical deflation issue in the MSOAR procedure so that implicit restarting is applicable. We then develop implicitly restarted MSOAR and RSOAR algorithms, in which we propose the exact and refined shifts for the two algorithms, respectively. Numerical examples illustrate the efficiency of the restarted algorithms and the superiority of the restarted RSOAR algorithm.
What problem does this paper attempt to address?