An Explicit Formula for the Matrix Logarithm

In various areas of applied numerics, the problem of calculating the logarithm of a matrix A emerges. Since series expansions of the logarithm usually do not converge well for matrices far away from the identity, the standard numerical method calculates successive square roots. In this article, a new algorithm is presented that relies on the computation of successive matrix exponentials. Convergence of the method is demonstrated for a large class of initial matrices and favor...

In this work, new closed-form formulas for the matrix exponential are provided. Our method is direct and elementary, it gives tractable and manageable formulas not current in the extensive literature on this essential subject. Moreover, others are recuperated and generalized. As a consequence, we easily obtain the Chevalley{Jordan decomposition and the spectral projections of any matrix. In addition, closed-form expressions for the arbitrary positive powers of matrices and th...

We show an explicit formula, with a quite easy deduction, for the exponential matrix $e^{tA}$ of a real square matrix $A$ of order $n\times n$. The elementary method developed requires neither Jordan canonical form, nor eigenvectors, nor resolution of linear systems of differential equations, nor resolution of linear systems with constant coefficients, nor matrix inversion, nor complex integration, nor functional analysis. The basic tools are power series and the method of pa...

How to calculate the exponential of matrices in an explicit manner is one of fundamental problems in almost all subjects in Science. Especially in Mathematical Physics or Quantum Optics many problems are reduced to this calculation by making use of some approximations whether they are appropriate or not. However, it is in general not easy. In this paper we give a very useful formula which is both elementary and getting on with computer.

In these notes, we consider the problem of finding the logarithm or the square root of a real matrix. It is known that for every real n x n matrix, A, if no real eigenvalue of A is negative or zero, then A has a real logarithm, that is, there is a real matrix, X, such that e^X = A. Furthermore, if the eigenvalues, xi, of X satisfy the property -pi < Im(xi) < pi, then X is unique. It is also known that under the same condition every real n x n matrix, A, has a real square root...

The most popular method for computing the matrix logarithm is a combination of the inverse scaling and squaring method in conjunction with a Pad\'e approximation, sometimes accompanied by the Schur decomposition. The main computational effort lies in matrix-matrix multiplications and left matrix division. In this work we illustrate that the number of such operations can be substantially reduced, by using a graph based representation of an efficient polynomial evaluation schem...

We consider in this work quantities that can be obtained as limits of powers of parametrized matrices, for instance the inverse matrix or the logarithm of the determinant. Under the assumption of affine dependence in the parameters, we use the Empirical Interpolation Method (EIM) to derive an approximation for powers of these matrices, from which we derive a nonintrusive approximation for the aforementioned limits. We derive upper bounds of the error made by the obtained form...

In this article, we focus on the error that is committed when computing the matrix logarithm using the Gauss--Legendre quadrature rules. These formulas can be interpreted as Pad\'e approximants of a suitable Gauss hypergeometric function. Empirical observation tells us that the convergence of these quadratures becomes slow when the matrix is not close to the identity matrix, thus suggesting the usage of an inverse scaling and squaring approach for obtaining a matrix with this...

An algorithm for numerically computing the exponential of a matrix is presented. We have derived a polynomial expansion of $e^x$ by computing it as an initial value problem using a symbolic programming language. This algorithm is shown to be comparable in operation count and convergence with the state--of--the--art method which is based on a Pade approximation of the exponential matrix function. The present polynomial form, however, is more reliable because the evaluation req...

We consider the computation of the matrix logarithm by using numerical quadrature. The efficiency of numerical quadrature depends on the integrand and the choice of quadrature formula. The Gauss--Legendre quadrature has been conventionally employed; however, the convergence could be slow for ill-conditioned matrices. This effect may stem from the rapid change of the integrand values. To avoid such situations, we focus on the double exponential formula, which has been develope...