期刊文献+

On convergence of the inexact Rayleigh quotient iteration with the Lanczos method used for solving linear systems 被引量:2

On convergence of the inexact Rayleigh quotient iteration with the Lanczos method used for solving linear systems
原文传递
导出
摘要 For the Hermitian inexact Rayleigh quotient iteration (RQI), we consider the local convergence of the inexact RQI with the Lanczos method for the linear systems involved. Some attractive properties are derived for the residual, whose norm is ξk, of the linear system obtained by the Lanczos method at outer iteration k + 1. Based on them, we make a refined analysis and establish new local convergence results. It is proved that (i) the inexact RQI with Lanezos converges quadratically provided that ξk ≤ξ with a constant ξ≥) 1 and (ii) the method converges linearly provided that ξk is bounded by some multiple of 1/‖τk‖ with ‖τk‖ the residual norm of the approximate eigenpair at outer iteration k. The results are fundamentally different from the existing ones that always require ξk 〈 1, and they have implications on effective implementations of the method. Based on the new theory, we can design practical criteria to control ξk to achieve quadratic convergence and implement the method more effectively than ever before. Numerical experiments confirm our theory and demonstrate that the inexact RQI with Lanczos is competitive to the inexact RQI with MINRES. For the Hermitian inexact Rayleigh quotient iteration(RQI),we consider the local convergence of the in exact RQI with the Lanczos method for the linear systems involved.Some attractive properties are derived for the residual,whose norm is ξk,of the linear system obtained by the Lanczos method at outer iteration k+1.Based on them,we make a refned analysis and establish new local convergence results.It is proved that(i) the inexact RQI with Lanczos converges quadratically provided that ξk≤ξ with a constant ξ1 and (ii) the method converges linearly provided that ξk is bounded by some multiple of1/||rk|| with rkthe residual norm of the approximate eigenpair at outer iteration k.The results are fundamentally diferent from the existing ones that always require ξk<1,and they have implications on efective implementations of the method.Based on the new theory,we can design practical criteria to control ξkto achieve quadratic convergence and implement the method more efectively than ever before.Numerical experiments confrm our theory and demonstrate that the inexact RQI with Lanczos is competitive to the inexact RQI with MINRES.
作者 JIA ZhongXiao
出处 《Science China Mathematics》 SCIE 2013年第10期2145-2160,共16页 中国科学:数学(英文版)
基金 supported by National Basic Research Program of China(Grant No.2011CB302400) National Natural Science Foundation of China(Grant No.11071140)
关键词 HERMITIAN inexact RQI CONVERGENCE inner iteration outer iteration LANCZOS 迭代收敛 线性系统 sz法 瑞利商 Lanczos法 局部收敛性 求解 埃尔米特
  • 相关文献

参考文献26

  • 1Berns-Müller J, Graham I G, Spence A. Inexact inverse iteration for symmetric matrices. Linear Algebra Appl, 2006,416: 389-413.
  • 2Berns-Müller J, Spence A. Inexact inverse iteration with variable shift for nonsymmetric generalized eigenvalue problems. SIAM J Matrix Anal Appl, 2006, 28: 1069-1082.
  • 3Daniel G W, Gragg W B, Kaufmann L, et al. Reorthogonalization and stable algorithms for updating the Gram-Schmidt QR factorization. Math Comput, 1976, 30: 772-795.
  • 4Duff I S, Grimes R G, Lewis J G. User's guide for the Harwell-Boeing sparse matrix collection (Release 1). Technical Report, RAL-92-086. Rutherford Appleton Laboratory, UK, 1992. Data available at http://math.nist.gov/MarketMatrix.
  • 5Freitag M A, Spence A. Convergence of inexact inverse iteration with application to preconditioned iterative solves. BIT, 2007, 47: 27-44.
  • 6Freitag M A, Spence A. Convergence theory for inexact inverse iteration applied to the generalised nonsymmetric eigenproblem. Electron Trans Numer Anal, 2007, 28: 40-64.
  • 7Freitag M A, Spence A. Rayleigh quotient iteration and simplified Jacobi-Davidson method with preconditioned iterative solves. Linear Algebra Appl, 2008, 428: 2049-2060.
  • 8Freitag M A, Spence A. A tuned preconditioner for inexact inverse iteration applied to Hermitian eigenvalue problems. IMA J Numer Anal, 2008, 28: 522-551.
  • 9Golub G H, van Loan C F. Matrix Computations, 3rd Ed. Baltimore: The John Hopkins Univ Press, 1996.
  • 10Hochestenbach M E, Notay Y. Controlling inner iterations in the Jacobi-Davidson method. SIAM J Matrix Anal Appl,2009, 31: 460-477.

二级参考文献26

  • 1Berns-Müller J, Graham I G, Spence A. Inexact inverse iteration for symmetric matrices. Linear Algebra Appl, 2006,416: 389-413.
  • 2Berns-Müller J, Spence A. Inexact inverse iteration with variable shift for nonsymmetric generalized eigenvalue problems. SIAM J Matrix Anal Appl, 2006, 28: 1069-1082.
  • 3Daniel G W, Gragg W B, Kaufmann L, et al. Reorthogonalization and stable algorithms for updating the Gram-Schmidt QR factorization. Math Comput, 1976, 30: 772-795.
  • 4Duff I S, Grimes R G, Lewis J G. User's guide for the Harwell-Boeing sparse matrix collection (Release 1). Technical Report, RAL-92-086. Rutherford Appleton Laboratory, UK, 1992. Data available at http://math.nist.gov/MarketMatrix.
  • 5Freitag M A, Spence A. Convergence of inexact inverse iteration with application to preconditioned iterative solves. BIT, 2007, 47: 27-44.
  • 6Freitag M A, Spence A. Convergence theory for inexact inverse iteration applied to the generalised nonsymmetric eigenproblem. Electron Trans Numer Anal, 2007, 28: 40-64.
  • 7Freitag M A, Spence A. Rayleigh quotient iteration and simplified Jacobi-Davidson method with preconditioned iterative solves. Linear Algebra Appl, 2008, 428: 2049-2060.
  • 8Freitag M A, Spence A. A tuned preconditioner for inexact inverse iteration applied to Hermitian eigenvalue problems. IMA J Numer Anal, 2008, 28: 522-551.
  • 9Golub G H, van Loan C F. Matrix Computations, 3rd Ed. Baltimore: The John Hopkins Univ Press, 1996.
  • 10Hochestenbach M E, Notay Y. Controlling inner iterations in the Jacobi-Davidson method. SIAM J Matrix Anal Appl,2009, 31: 460-477.

共引文献1

同被引文献27

  • 1Berns-Müller J, Graham I G, Spence A. Inexact inverse iteration for symmetric matrices. Linear Algebra Appl, 2006,416: 389-413.
  • 2Berns-Müller J, Spence A. Inexact inverse iteration with variable shift for nonsymmetric generalized eigenvalue problems. SIAM J Matrix Anal Appl, 2006, 28: 1069-1082.
  • 3Daniel G W, Gragg W B, Kaufmann L, et al. Reorthogonalization and stable algorithms for updating the Gram-Schmidt QR factorization. Math Comput, 1976, 30: 772-795.
  • 4Duff I S, Grimes R G, Lewis J G. User's guide for the Harwell-Boeing sparse matrix collection (Release 1). Technical Report, RAL-92-086. Rutherford Appleton Laboratory, UK, 1992. Data available at http://math.nist.gov/MarketMatrix.
  • 5Freitag M A, Spence A. Convergence of inexact inverse iteration with application to preconditioned iterative solves. BIT, 2007, 47: 27-44.
  • 6Freitag M A, Spence A. Convergence theory for inexact inverse iteration applied to the generalised nonsymmetric eigenproblem. Electron Trans Numer Anal, 2007, 28: 40-64.
  • 7Freitag M A, Spence A. Rayleigh quotient iteration and simplified Jacobi-Davidson method with preconditioned iterative solves. Linear Algebra Appl, 2008, 428: 2049-2060.
  • 8Freitag M A, Spence A. A tuned preconditioner for inexact inverse iteration applied to Hermitian eigenvalue problems. IMA J Numer Anal, 2008, 28: 522-551.
  • 9Golub G H, van Loan C F. Matrix Computations, 3rd Ed. Baltimore: The John Hopkins Univ Press, 1996.
  • 10Hochestenbach M E, Notay Y. Controlling inner iterations in the Jacobi-Davidson method. SIAM J Matrix Anal Appl,2009, 31: 460-477.

引证文献2

二级引证文献3

相关作者

内容加载中请稍等...

相关机构

内容加载中请稍等...

相关主题

内容加载中请稍等...

浏览历史

内容加载中请稍等...
;
使用帮助 返回顶部