A fixed point iteration for computing the matrix logarithm
Numerical Analysis
2007-07-19 v1 Data Structures and Algorithms
Abstract
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 favorable choices of the initial matrix are discussed.
Cite
@article{arxiv.0707.2701,
title = {A fixed point iteration for computing the matrix logarithm},
author = {Gernot Schaller},
journal= {arXiv preprint arXiv:0707.2701},
year = {2007}
}
Comments
4 pages, 3 figures, comments welcome