Fasi, Massimiliano and Higham, Nicholas J. (2018) Multiprecision Algorithms for Computing the Matrix Logarithm. SIAM Journal On Matrix Analysis And Applications, 39 (1). pp. 472491. ISSN 10957162
This is the latest version of this item.
Text
fahi18.pdf  Accepted Version Download (528kB) 
Abstract
Two algorithms are developed for computing the matrix logarithm in floating point arithmetic of any specified precision. The backward errorbased approach used in the state of the art inverse scaling and squaring algorithms does not conveniently extend to a multiprecision environment, so instead we choose algorithmic parameters based on a forward error bound. We derive a new forward error bound for Pad\'{e} approximants that for highly nonnormal matrices can be much smaller than the classical bound of Kenney and Laub. One of our algorithms exploits a Schur decomposition while the other is transformationfree and uses only the computational kernels of matrix multiplication and the solution of multiple righthand side linear systems. For double precision computations the algorithms are competitive with the state of the art algorithm of AlMohy, Higham, and Relton implemented in \texttt{logm} in MATLAB\@. They are intended for computing environments providing multiprecision floating point arithmetic, such as Julia, MATLAB via the Symbolic Math Toolbox or the Multiprecision Computing Toolbox, or Python with the mpmath or SymPy packages. We show experimentally that the algorithms behave in a forward stable manner over a wide range of precisions, unlike existing alternatives.
Item Type:  Article 

Uncontrolled Keywords:  multiprecision arithmetic, matrix logarithm, principal logarithm, inverse scaling and squaring method, Fr\'{e}chet derivative, Pad\'{e} approximation, Taylor approximation, forward error analysis, MATLAB, logm. 
Subjects:  MSC 2010, the AMS's Mathematics Subject Classification > 15 Linear and multilinear algebra; matrix theory MSC 2010, the AMS's Mathematics Subject Classification > 65 Numerical analysis 
Depositing User:  Nick Higham 
Date Deposited:  17 Mar 2018 09:00 
Last Modified:  17 Mar 2018 09:00 
URI:  http://eprints.maths.manchester.ac.uk/id/eprint/2628 
Available Versions of this Item

Multiprecision Algorithms for Computing the Matrix Logarithm. (deposited 11 May 2017)

Multiprecision Algorithms for Computing the Matrix Logarithm. (deposited 25 Oct 2017)
 Multiprecision Algorithms for Computing the Matrix Logarithm. (deposited 17 Mar 2018 09:00) [Currently Displayed]

Multiprecision Algorithms for Computing the Matrix Logarithm. (deposited 25 Oct 2017)
Actions (login required)
View Item 