A positivity preserving inexact Noda iteration for computing the smallest eigenpair of a large irreducible M-matrix

Zhongxiao Jia*, Wen-Wei Lin, Ching Sung Liu

*Corresponding author for this work

Research output: Contribution to journalArticle

7 Scopus citations

Abstract

In this paper, based on the Noda iteration, we present inexact Noda iterations (INI), to find the smallest eigenvalue and the associated positive eigenvector of a large irreducible nonsingular $$M$$M-matrix. The positivity of approximations is critical in applications, and if the approximations lose the positivity then they may be meaningless and could not be interpreted. We propose two different inner tolerance strategies for solving the inner linear systems involved, and prove that the convergence of resulting INI algorithms is globally linear and superlinear with the convergence order $$\frac{1+\sqrt{5}}{2}$$1+52, respectively. The proposed INI algorithms are structure preserving and maintains the positivity of approximate eigenvectors. We also revisit the exact Noda iteration and establish a new quadratic convergence result. All the above is first done for the problem of computing the Perron root and the positive Perron vector of an irreducible nonnegative matrix and is then adapted to computing the smallest eigenpair of the irreducible nonsingular $$M$$M-matrix. Numerical examples illustrate that the proposed INI algorithms are practical, and they always preserve the positivity of approximate eigenvectors. We compare them with the Jacobi–Davidson method, the implicitly restarted Arnoldi method and the explicitly restarted Krylov–Schur method, all of which cannot guarantee the positivity of approximate eigenvectors, and illustrate that the overall efficiency of the INI algorithms is competitive with and can be considerably higher than the latter three methods.

Original languageEnglish
Pages (from-to)645-679
Number of pages35
JournalNumerische Mathematik
Volume130
Issue number4
DOIs
StatePublished - 19 Aug 2015

Keywords

  • 15B99
  • 65F15
  • 65F50
  • Primary 15A18
  • Secondary 15B48

Fingerprint Dive into the research topics of 'A positivity preserving inexact Noda iteration for computing the smallest eigenpair of a large irreducible M-matrix'. Together they form a unique fingerprint.

  • Cite this