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

PDF
paper.pdf Download (289kB) |

## Abstract

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 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: | 16 Feb 2006 |

Last Modified: | 20 Oct 2017 14:12 |

URI: | http://eprints.maths.manchester.ac.uk/id/eprint/156 |

### Actions (login required)

View Item |