A new method for accelerating Arnoldi algorithms for large scale Eigenproblems

We propose a new method for accelerating the convergence of the implicitly restarted Arnoldi (IRA) algorithm for the solution of large sparse nonsymmetric eigenvalue problems. A new relationship between the residual of the current step and the residual in the previous step is derived and we use this relationship to develop a technique for dynamically switching the Krylov subspace dimension at successive cycles. We give numerical results for various difficult nonsymmetric eigenvalue problems that demonstrate the capability of the dynamic switching strategy for significantly accelerating the convergence of Arnoldi algorithms. For some large scale difficult eigenvalue problems that arise in the fields of computational fluid dynamics, electrical engineering and materials science, our strategy leads to significant reductions in the number of matrix-vector products, orthogonalization costs and computational time.

[1]  Zhongxiao Jia,et al.  Polynomial characterizations of the approximate eigenvectors by the refined Arnoldi method and an implicitly restarted refined Arnoldi algorithm , 1999 .

[2]  W. Arnoldi The principle of minimized iterations in the solution of the matrix eigenvalue problem , 1951 .

[3]  Takashi Nodera,et al.  The DEFLATED-GMRES(m, k) method with switching the restart frequency dynamically , 2000 .

[4]  Richard B. Lehoucq,et al.  Implicitly Restarted Arnoldi Methods and Subspace Iteration , 2001, SIAM J. Matrix Anal. Appl..

[5]  Chao Yang,et al.  ARPACK users' guide - solution of large-scale eigenvalue problems with implicitly restarted Arnoldi methods , 1998, Software, environments, tools.

[6]  Y. Saad Variations on Arnoldi's method for computing eigenelements of large unsymmetric matrices , 1980 .

[7]  John G. Lewis,et al.  Sparse matrix test problems , 1982, SGNM.

[8]  Z. Jia,et al.  Refined iterative algorithms based on Arnoldi's process for large unsymmetric eigenproblems , 1997 .

[9]  Y. Saad,et al.  Numerical Methods for Large Eigenvalue Problems , 2011 .

[10]  Zhongxiao Jia,et al.  The Convergence of Generalized Lanczos Methods for Large Unsymmetric Eigenproblems , 1995, SIAM J. Matrix Anal. Appl..

[11]  Ronald B. Morgan,et al.  GMRES with Deflated Restarting , 2002, SIAM J. Sci. Comput..

[13]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[14]  R. Morgan,et al.  Harmonic projection methods for large non-symmetric eigenvalue problems , 1998 .

[15]  Ronald B. Morgan,et al.  On restarting the Arnoldi method for large nonsymmetric eigenvalue problems , 1996, Math. Comput..

[16]  Danny C. Sorensen,et al.  Implicit Application of Polynomial Filters in a k-Step Arnoldi Method , 1992, SIAM J. Matrix Anal. Appl..

[17]  Kesheng Wu,et al.  Dynamic Thick Restarting of the Davidson, and the Implicitly Restarted Arnoldi Methods , 1998, SIAM J. Sci. Comput..