Approximating the Logarithm of a Matrix to Specified Accuracy

The standard inverse scaling and squaring algorithm for computing the matrix logarithm begins by transforming the matrix to Schur triangular form in order to facilitate subsequent matrix square root and Pade approximation computations. A transformation-free form of this method that exploits incomplete Denman--Beavers square root iterations and aims for a specified accuracy (ignoring roundoff) is presented. The error introduced by using approximate square roots is accounted for by a novel splitting lemma for logarithms of matrix products. The number of square root stages and the degree of the final Pade approximation are chosen to minimize the computational work. This new method is attractive for high-performance computation since it uses only the basic building blocks of matrix multiplication, LU factorization and matrix inversion.

[1]  N. Higham.  Computing real square roots of a real matrix , 1987 .

[2]  A. Laub,et al.  Rational iterative methods for the matrix sign function , 1991 .

[3]  N. Higham.  The matrix sign decomposition and its relation to the polar decomposition , 1994 .

[4]  A. Laub,et al.  The matrix sign function , 1995, IEEE Trans. Autom. Control..

[5]  G. A. Baker.  Essentials of Padé approximants , 1975 .

[6]  Nicholas J. Higham,et al.  Evaluating Padé Approximants of the Matrix Logarithm , 2000, SIAM J. Matrix Anal. Appl..

[7]  R. Ward.  Numerical Computation of the Matrix Exponential with Accuracy Estimate , 1977 .

[8]  A. Laub,et al.  Condition Estimates for Matrix Functions , 1989 .

[9]  Nicholas J. Higham,et al.  Stable iterations for the matrix square root , 1997, Numerical Algorithms.

[10]  Shirley Dex,et al.  JR 旅客販売総合システム（マルス）における運用及び管理について , 1991 .

[11]  A. Laub,et al.  Padé error estimates for the logarithm of a matrix , 1989 .

[12]  C. Loan,et al.  Nineteen Dubious Ways to Compute the Exponential of a Matrix , 1978 .

[13]  N. Higham.  Newton's method for the matrix square root , 1986 .

[14]  Gene H. Golub,et al.  Matrix Computations, Third Edition , 1996 .

[15]  Luca Dieci,et al.  Conditioning and Padé Approximation of the Logarithm of a Matrix , 2000, SIAM J. Matrix Anal. Appl..

[16]  Burton H. Singer,et al.  The Representation of Social Processes by Markov Models , 1976, American Journal of Sociology.

[17]  Å. Björck,et al.  A Schur method for the square root of a matrix , 1983 .

[18]  Benedetta Morini,et al.  Computational Techniques for Real Logarithms of Matrices , 1996, SIAM J. Matrix Anal. Appl..

[19]  E. Denman,et al.  The matrix sign function and computations in systems , 1976 .

[20]  Charles R. Johnson,et al.  Topics in Matrix Analysis , 1991 .

[21]  Richard C. Allen,et al.  An Analysis of an Inverse Problem in Ordinary Differential Equations , 1981 .