A Givens-Weight Representation for Rank Structured Matrices

In this paper we introduce a Givens-weight representation for rank structured matrices, where the rank structure is defined by certain submatrices starting from the bottom left or upper right matrix corner being of low rank. We proceed in two steps. First, we introduce a unitary-weight representation. This representation will be compared to the (block) quasiseparable representations introduced by P. Dewilde and A.-J. van der Veen [Time-varying Systems and Computations, Kluwer Academic Publishers, Boston, 1998]. More specifically, we show that our unitary-weight representations are theoretically equivalent to the so-called block quasiseparable representations in input or output normal form introduced by Dewilde and van der Veen [Time varying Systems and Computations, Kluwer Academic Publishers, Boston, 1998]. Next, we move from the unitary-weight to the Givens-weight representation. We then provide some basic algorithms for the unitary/Givens-weight representation, showing how to obtain such a representation for a dense matrix by means of numerical approximation. We also show how to “swap” the representation and how to reduce the number of parameters of the representation, whenever appropriate. As such, these results will become the basis for algorithms on unitary/Givens-weight representations to be described in subsequent papers.

[1]  Gene H. Golub,et al.  Matrix computations (3rd ed.) , 1996 .

[2]  Gene H. Golub,et al.  Matrix computations , 1983 .

[3]  Marc Van Barel,et al.  A Hessenberg Reduction Algorithm for Rank Structured Matrices , 2007, SIAM J. Matrix Anal. Appl..

[4]  Marc Van Barel,et al.  A QR-Based Solver for Rank Structured Matrices , 2008, SIAM J. Matrix Anal. Appl..

[5]  Eugene E. Tyrtyshnikov Mosaic Ranks for Weakly Semiseparable Matrices , 1999, LSSC.

[6]  I. Gohberg,et al.  A modification of the Dewilde-van der Veen method for inversion of finite structured matrices , 2002 .

[7]  P. Dewilde,et al.  Time-Varying Systems and Computations , 1998 .

[8]  A. V. D. Veen,et al.  Inner-outer factorization and the inversion of locally finite systems of equations , 2000 .

[9]  Marc Van Barel,et al.  Orthogonal similarity transformation into block-semiseparable matrices of semiseparability rank k , 2005, Numer. Linear Algebra Appl..

[10]  Raf Vandebril,et al.  A note on the representation and definition of semiseparable matrices , 2005, Numer. Linear Algebra Appl..

[11]  Israel Koltracht,et al.  Linear complexity algorithm for semiseparable matrices , 1985 .

[12]  I. Gohberg,et al.  On generators of quasiseparable finite block matrices , 2005 .

[13]  Alle-Jan van der Veen,et al.  Fast Stable Solver for Sequentially Semi-separable Linear Systems of Equations , 2002, HiPC.

[14]  Marc Van Barel,et al.  Rank structures preserved by the QR-algorithm: The singular case , 2006 .

[15]  Raf Vandebril,et al.  An implicit QR algorithm for symmetric semiseparable matrices , 2005, Numer. Linear Algebra Appl..

[16]  I. Gohberg,et al.  On a new class of structured matrices , 1999 .

[17]  S. Delvaux,et al.  Structures preserved by the QR-algorithm , 2006 .

[18]  Thomas Kailath,et al.  Linear complexity algorithms for semiseparable matrices , 1985 .