Abstract:We revisit the shift-and-invert Arnoldi method proposed in [S. Lee, H. Pang, and H. Sun. {\it Shift-invert Arnoldi approximation to the Toeplitz matrix exponential}, SIAM J. Sci. Comput., 32: 774--792, 2010] for numerical approximation to the product of Toeplitz matrix exponential with a vector. In this approach, one has to solve two large scale Toeplitz linear systems in advance. However, if the desired accuracy is high, the cost will be prohibitive. Therefore, it is interesting to investigate how to solve the Toeplitz systems inexactly in this method. The contribution of this paper is in three regards. First, we give a new stability analysis on the Gohberg-Semencul formula (GSF) and define the GSF condition number of a Toeplitz matrix. It is shown that, when the size of the Toeplitz matrix is large, our result is sharper than the one given in [M. Gutknecht and M. Hochbruck. {\it The stability of inversion formulas for Toeplitz matrices}, Linear Algebra Appl., 223/224: 307--324, 1995]. Second, we establish a relation between the error of Toeplitz systems and the residual of Toeplitz matrix exponential. We show that if the GSF condition number of the Toeplitz matrix is medium sized, then the Toeplitz systems can be solved in a low accuracy. Third, based on this relationship, we present a practical stopping criterion for relaxing the accuracy of the Toeplitz systems, and propose an inexact shift-and-invert Arnoldi algorithm for the Toeplitz matrix exponential problem. Numerical experiments illustrate the numerical behavior of the new algorithm, and show the effectiveness of our theoretical results.

The Inverses of Block Toeplitz Matrices

On Group Inverse of Singular Toeplitz Matrices

Explicit formulas for the inverses of Toeplitz matrices, with applications

On Drazin inverse of singular Toeplitz matrix

Explicit inverse of symmetric, tridiagonal near Toeplitz matrices Part II: with weakly diagonally dominant Toeplitz

Modified Newton'S Algorithm for Computing the Group Inverses of Singular Toeplitz Matrices

Structured Matrices And Their Generalized Inverses

On the invertibility of matrices with a double saddle-point structure

The Intrinsic Toeplitz Structure and Its Applications in Algebraic Riccati Equations

Computing Moore-Penrose Inverses of Toeplitz Matrices by Newton's Iteration

Inequalities on $ 2\times 2 $ block accretive partial transpose matrices

New Representations of the Group Inverse of 2×2 Block Matrices

Integer Levinson Algorithm for the Inversion of any Nonsingular Hermitian Toeplitz Matrix

Explicit inverse of symmetric, tridiagonal near Toeplitz matrices with strictly diagonally dominant Toeplitz part

Generalized Inverses of Tensors Via a General Product of Tensors

Generalized Inverses and a Block-Rank Equation.

Cayley Transform for Toeplitz and Dual Matrices

Representations of Generalized Inverses and Drazin Inverse of Partitioned Matrix with Banachiewicz-Schur Forms

Block Toeplitz Matrices: Multiplicative Properties

Efficient Computation for Invertibility Sequence of Banded Toeplitz Matrices

Inexact Shift-and-Invert Arnoldi for Toeplitz Matrix Exponential