## 2010.91: A Schur--Pad\'e Algorithm for Fractional Powers of a Matrix

2010.91:
Nicholas J. Higham and Lijing Lin
(2010)
*A Schur--Pad\'e Algorithm for Fractional Powers of a Matrix.*

There is a more recent version of this eprint available. Click here to view it.

Full text available as:

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

Other (Zip file containing MATLAB codes) 4 Kb |

## Abstract

A new algorithm is developed for computing arbitrary real powers $A^p$ of a matrix $A\in\mathbb{C}^{n\times n}$. The algorithm starts with a Schur decomposition, takes $k$ square roots of the triangular factor $T$, evaluates an $[m/m]$ Pad\'e approximant of $(1-x)^p$ at $I - T^{1/2^k}$, and squares the result $k$ times. The parameters $k$ and $m$ are chosen to minimize the cost subject to achieving double precision accuracy in the evaluation of the Pad\'e approximant, making use of a result that bounds the error in the matrix Pad\'e approximant by the error in the scalar Pad\'e approximant with argument the norm of the matrix. The Pad\'e approximant is evaluated from the continued fraction representation in bottom-up fashion, which is shown to be numerically stable. In the squaring phase the diagonal and first superdiagonal are computed from explicit formulae for $T^{p/2^j}$, yielding increased accuracy. Since the basic algorithm is designed for $p\in(-1,1)$, a criterion for reducing an arbitrary real $p$ to this range is developed, making use of bounds for the condition number of the $A^p$ problem. How best to compute $A^k$ for a negative integer $k$ is also investigated. In numerical experiments the new algorithm is found to be superior in accuracy and stability to several alternatives, including the use of an eigendecomposition and approaches based on the formula $A^p = \exp(p\log(A))$.

Item Type: | MIMS Preprint |
---|---|

Uncontrolled Keywords: | matrix power, matrix root, fractional power, primary matrix function, Schur decomposition, Pad\'e approximation, Pad\'e approximant, matrix logarithm, matrix exponential, MATLAB, CICADA |

Subjects: | MSC 2000 > 15 Linear and multilinear algebra; matrix theory MSC 2000 > 65 Numerical analysis |

MIMS number: | 2010.91 |

Deposited By: | Nick Higham |

Deposited On: | 24 February 2011 |

### Available Versions of this Item

- A Schur--Pad\'e Algorithm for Fractional Powers of a Matrix (deposited 05 October 2011)
- A Schur--Pad\'e Algorithm for Fractional Powers of a Matrix (deposited 01 June 2011)
- A Schur--Pad\'e Algorithm for Fractional Powers of a Matrix (deposited 24 February 2011)
**[Currently Displayed]** - A Schur--Pad\'e Algorithm for Fractional Powers of a Matrix (deposited 19 October 2010)

Download Statistics: last 4 weeks

Repository Staff Only: edit this item