Recursive formulation of the matrix Pade approximation in packed storage
The Extended Euclidean algorithm for matrix Pade approximants is applied to compute matrix Pade approximants when the coefficient matrices of the input matrix polynomial are triangular. The procedure given by Bjarne S. Anderson et al. for packing a triangular matrix in recursive packed storage is applied to pack a sequence of lower triangular matrices of a matrix polynomial in recursive packed storage. This recursive packed storage for a matrix polynomial is applied to compute matrix Pade approximants of the matrix polynomial using the Matrix NO Extended Euclidean algorithm in packed form. The CPU time and memory comparison, in computing the matrix Pade approximants of a matrix polynomial, between the packed case and the non-packed case are described in detail. (C) 2009 Elsevier Ltd. All rights reserved.
matrix pade approximants recursive packed storage