Naslov
Efficient generalized Hessenberg form and applications

Autori
Bosner, Nela ; Bujanović, Zvonimir ; Drmač, Zlatko

Izvornik
Acm transactions on mathematical software (0098-3500) 39 (2013), 3; 19, 19

Vrsta, podvrsta i kategorija rada
Radovi u časopisima, članak, znanstveni

Ključne riječi
Hessenberg form ; shifted linear systems ; staircase algorithm ; transfer function

Sažetak
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 dynamic 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.

Izvorni jezik
Engleski

Znanstvena područja
Matematika

POVEZANOST RADA