An Arnoldi code for computing selected eigenvalues of sparse, real, unsymmetric matrices

  • Authors:
  • Jennifer A. Scott

  • Affiliations:
  • Rutherford Appleton Lab., Oxon, UK

  • Venue:
  • ACM Transactions on Mathematical Software (TOMS)
  • Year:
  • 1995

Quantified Score

Hi-index 0.00

Visualization

Abstract

Arnoldi methods can be more effective than subspace iterationmethods for computing the dominant eigenvalues of a large, sparse, real,unsymmetric matrix. A code, EB12, for thesparse, unsymmetric eigenvalue problem based on a subspace iterationalgorithm, optionally combined with Chebychev acceleration, has recentlybeen described by Duff and Scott and is included in the HarwellSubroutine Library. In this article we consider variants of the methodof Arnoldi and discuss the design and development of a code to implementthese methods. The new code, which is calledEB13, offers the user the choice of abasic Arnoldi algorithm, an Arnoldi algorithm with Chebychevacceleration, and a Chebychev preconditioned Arnoldi algorithm. Eachmethod is available in blocked and unblocked form. The code may be usedto compute either the rightmost eigenvalues, the eigenvalues of largestabsolute value, or the eigenvalues of largest imaginary part. Theperformance of each option in the EB13package is compared with that of subspace iteration on a range of testproblems, and on the basis of the results, advice is offered to the useron the appropriate choice of method.—Author's Abstract