Shift-Invert Arnoldi Approximation to the Toeplitz Matrix Exponential

  • Authors:
  • Spike T. Lee;Hong-Kui Pang;Hai-Wei Sun

  • Affiliations:
  • ma76522@umac.mo and ya87402@umac.mo and HSun@umac.mo;-;-

  • Venue:
  • SIAM Journal on Scientific Computing
  • Year:
  • 2010

Quantified Score

Hi-index 0.00

Visualization

Abstract

The shift-invert Arnoldi method is employed to generate an orthonormal basis from the Krylov subspace corresponding to a real Toeplitz matrix and an initial vector. The vectors and recurrence coefficients produced by this method are exploited to approximate the Toeplitz matrix exponential. Toeplitz matrix inversion formula and rapid Toeplitz matrix-vector multiplications are utilized to lower the computational costs. For convergence analysis, a sufficient condition is established to guarantee that the error bound is independent of the norm of the matrix. Numerical results are given to demonstrate the efficiency of the method.