Aggregations of Elementary Transformations

The techniques of aggregating transformations have been used in the development of blocked algorithms, in particular, the compact representations of aggregated Householder transformation have been wellknown. We show in this paper that the aggregation techniques can be integrated in a framework for aggregating general elementary transformations. The framework consists of three basic compact forms, as extensions of the three compact forms for aggregated Householder transformations, and two componentwise representations of the compact forms. The componentwise representations, introduced in this paper, are a central link to various computational procedures (known and new) for the compact forms, and provide a stepping stone for numerical analysis of blocked algorithms. We show also the connection between the computation of the compact forms and the inverse of general nonsingular triangular matrices.

[1]  Åke Björck,et al.  Numerical methods for least square problems , 1996 .

[2]  H. Walker Implementation of the GMRES method using householder transformations , 1988 .

[3]  Charles L. Lawson,et al.  Basic Linear Algebra Subprograms for Fortran Usage , 1979, TOMS.

[4]  Nicholas J. Higham,et al.  INVERSE PROBLEMS NEWSLETTER , 1991 .

[5]  Alston S. Householder,et al.  The Theory of Matrices in Numerical Analysis , 1964 .

[6]  C. H. Bischof,et al.  A summary of block schemes for reducing a general matrix to Hessenberg form , 1993 .

[7]  C. Loan,et al.  A Storage-Efficient $WY$ Representation for Products of Householder Transformations , 1989 .

[8]  Robert Schreiber,et al.  Block Algorithms for Parallel Machines , 1988 .

[9]  G. Stewart Introduction to matrix computations , 1973 .

[10]  Christopher C. Paige,et al.  Loss and Recapture of Orthogonality in the Modified Gram-Schmidt Algorithm , 1992, SIAM J. Matrix Anal. Appl..

[11]  Christian H. Bischof,et al.  A Basis-Kernel Representation of Orthogonal Matrices , 1995, SIAM J. Matrix Anal. Appl..

[12]  Christian H. Bischof,et al.  The WY representation for products of householder matrices , 1985, PPSC.

[13]  C. Puglisi Modification of the householder method based on the compact WY representation , 1992 .

[14]  Jack J. Dongarra,et al.  A set of level 3 basic linear algebra subprograms , 1990, TOMS.

[15]  J. Navarro-Pedreño Numerical Methods for Least Squares Problems , 1996 .

[16]  John B. Shoven,et al.  I , Edinburgh Medical and Surgical Journal.