Higham, Nicholas J. and Lin, Lijing (2013) An Improved SchurPade Algorithm for Fractional Powers of a Matrix and their Frechet Derivatives. SIAM. J. Matrix Anal. & Appl., 34 (3). pp. 13411360. ISSN 10957162
This is the latest version of this item.
PDF
hili13.pdf Download (357kB) 
Abstract
The SchurPadÃ�Æ�Ã�Â© algorithm [N. J. Higham and L. Lin, A SchurPadÃ�Æ�Ã�Â© algorithm for fractional powers of a matrix, SIAM J. Matrix Anal. Appl., 32(3):10561078, 2011] computes arbitrary real powers $A^t$ of a matrix $A\in\mathbb{C}^{n\times n}$ using the building blocks of Schur decomposition, matrix square roots, and PadÃ�Æ�Ã�Â© approximants. We improve the algorithm by basing the underlying error analysis on the quantities $\(I A)^k\^{1/k}$, for several small $k$, instead of $\IA\$. We extend the algorithm so that it computes along with $A^t$ one or more FrÃ�Æ�Ã�Â©chet derivatives, with reuse of information when more than one FrÃ�Æ�Ã�Â©chet derivative is required, as is the case in condition number estimation. We also derive a version of the extended algorithm that works entirely in real arithmetic when the data is real. Our numerical experiments show the new algorithms to be superior in accuracy to, and often faster than, the original SchurPadÃ�Æ�Ã�Â© algorithm for computing matrix powers and more accurate than several alternative methods for computing the FrÃ�Æ�Ã�Â©chet derivative. They also show that reliable estimates of the condition number of $A^t$ are obtained by combining the algorithms with a matrix norm estimator.
Item Type:  Article 

Uncontrolled Keywords:  matrix power, fractional power, matrix root, FrÃ�Æ�Ã�Â©chet derivative, condition number, condition estimate, Schur decomposition, PadÃ�Æ�Ã�Â© approximation, PadÃ�Æ�Ã�Â© approximant, matrix logarithm, matrix exponential, 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:  Dr Lijing Lin 
Date Deposited:  18 Sep 2013 
Last Modified:  20 Oct 2017 14:13 
URI:  http://eprints.maths.manchester.ac.uk/id/eprint/2021 
Available Versions of this Item

An Improved SchurPadÃ© Algorithm for Fractional Powers of a Matrix and their FrÃ©chet Derivatives. (deposited 16 Jan 2013)

An Improved SchurPade Algorithm for Fractional Powers of a Matrix and their Frechet Derivatives. (deposited 03 May 2013)
 An Improved SchurPade Algorithm for Fractional Powers of a Matrix and their Frechet Derivatives. (deposited 18 Sep 2013) [Currently Displayed]

An Improved SchurPade Algorithm for Fractional Powers of a Matrix and their Frechet Derivatives. (deposited 03 May 2013)
Actions (login required)
View Item 