You are here: MIMS > EPrints
MIMS EPrints

2006.13: A Schur--Parlett Algorithm for Computing Matrix Functions

2006.13: Philip I. Davies and Nicholas J. Higham (2003) A Schur--Parlett Algorithm for Computing Matrix Functions. SIAM Journal On Matrix Analysis and Applications, 25 (2). pp. 464-485. ISSN 1095-7162

Full text available as:

PDF - Requires a PDF viewer such as GSview, Xpdf or Adobe Acrobat Reader
283 Kb

DOI: 10.1137/S0895479802410815


An algorithm for computing matrix functions is presented. It employs a Schur decomposition with reordering and blocking followed by the block form of a recurrence of Parlett, with functions of the nontrivial diagonal blocks evaluated via a Taylor series. A parameter is used to balance the conflicting requirements of producing small diagonal blocks and keeping the separations of the blocks large. The algorithm is intended primarily for functions having a Taylor series with an infinite radius of convergence, but it can be adapted for certain other functions, such as the logarithm. Novel features introduced here include a convergence test that avoids premature termination of the Taylor series evaluation and an algorithm for reordering and blocking the Schur form. Numerical experiments show that the algorithm is competitive with existing special-purpose algorithms for the matrix exponential, logarithm, and cosine. Nevertheless, the algorithm can be numerically unstable with the default choice of its blocking parameter (or in certain cases for all choices), and we explain why determining the optimal parameter appears to be a very difficult problem. A MATLAB implementation is available that is much more reliable than the function \texttt{funm} in MATLAB~6.5 (R13).

Item Type:Article
Uncontrolled Keywords:matrix function, matrix exponential, matrix logarithm, matrix cosine, Taylor series, Schur decomposition, Parlett recurrence, sep function, LAPACK, MATLAB
Subjects:MSC 2000 > 15 Linear and multilinear algebra; matrix theory
MSC 2000 > 65 Numerical analysis
MIMS number:2006.13
Deposited By:Nick Higham
Deposited On:16 February 2006

Download Statistics: last 4 weeks
Repository Staff Only: edit this item