Reduction of a regular matrix pair (A, B) to block Hessenberg-triangular form [chapter]

Krister Dackland, Bo Kågström
1996 Lecture Notes in Computer Science  
An algorithm for reduction of a regular matrix pair (A; B) to block Hessenberg-triangular form is presented. This condensed form Q T (A; B)Z = (H; T), where H and T are block upper Hessenberg and upper triangular, respectively, and Q and Z orthogonal, may serve as a rst step in the solution of the generalized eigenvalue problem Ax = Bx. It is shown how an elementwise algorithm can be reorganized in terms of blocked factorizations and higher level BLAS operations. Several ways to annihilate
more » ... nts are compared. Speci cally, the use of Givens rotations, Householder transformations, and combinations of the two. Performance results of the di erent variants are presented and compared to the LAPACK implementation DGGHRD, which indeed is unblocked.
doi:10.1007/3-540-60902-4_15 fatcat:owos5z2w6radrlgrlyp75lw52q