Affiliation:
1. University of Zagreb, Croatia
Abstract
This article proposes an efficient algorithm for reducing matrices to generalized Hessenberg form by unitary similarity, and recommends using it as a preprocessor in a variety of applications. To illustrate its usefulness, two cases from control theory are analyzed in detail: a solution procedure for a sequence of shifted linear systems with multiple right hand sides (e.g. evaluating the transfer function of a MIMO LTI dynamical system at many points) and computation of the staircase form. The proposed algorithm for the generalized Hessenberg reduction uses two levels of aggregation of Householder reflectors, thus allowing efficient BLAS 3-based computation. Another level of aggregation is introduced when solving many shifted systems by processing the shifts in batches. Numerical experiments confirm that the proposed methods have superior efficiency.
Funder
Ministarstvo Obrazovanja, Znanosti i Sporta
Publisher
Association for Computing Machinery (ACM)
Subject
Applied Mathematics,Software
Reference26 articles.
1. Anderson E. Bai Z. Bischof C. Demmel J. Dongarra J. Croz J. D. Greenbaum A. Hammarling S. McKenny A. Ostrouchov S. and Sorensen D. 1992. LAPACK Users' Guide 2nd Ed. SIAM Philadelphia PA. Anderson E. Bai Z. Bischof C. Demmel J. Dongarra J. Croz J. D. Greenbaum A. Hammarling S. McKenny A. Ostrouchov S. and Sorensen D. 1992. LAPACK Users' Guide 2nd Ed. SIAM Philadelphia PA.
2. A note on shifted Hessenberg systems and frequency response computation
3. A framework for symmetric band reduction
4. Bischof C. and Tang P. T. P. 1991a. Generalized incremental condition estimation. Tech. rep. 32 LAPACK working note. http://www.netlib.org/lapack/lawnspdf/lawn32.pdf. Bischof C. and Tang P. T. P. 1991a. Generalized incremental condition estimation. Tech. rep. 32 LAPACK working note. http://www.netlib.org/lapack/lawnspdf/lawn32.pdf.
5. Bischof C. and Tang P. T. P. 1991b. Robust incremental condition estimation. Tech. rep. 33 LAPACK working note. http://www.netlib.org/lapack/lawnspdf/lawn33.pdf. Bischof C. and Tang P. T. P. 1991b. Robust incremental condition estimation. Tech. rep. 33 LAPACK working note. http://www.netlib.org/lapack/lawnspdf/lawn33.pdf.
Cited by
5 articles.
订阅此论文施引文献
订阅此论文施引文献,注册后可以免费订阅5篇论文的施引文献,订阅后可以查看论文全部施引文献