Abstract
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 axe block upper Hessenberg and upper triangular, respectively, and Q and Z orthogonal, may serve as a first 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 elements are compared. Specifically, the use of Givens rotations, Householder transformations, and combinations of the two. Performance results of the different variants are presented and compared to the LAPACK implementation DGGHRD, which indeed is unblocked.
Preview
Unable to display preview. Download preview PDF.
References
E. Anderson, Z. Bai, J. Demmel, J. Dongarra, J. Du Croz, A. Greenbaum, S. Hammarling, A. McKenney, S. Ostrouchov, and D. Sorensen. LAPACK Users' Guide, Second Edition. SIAM Press, Philadelphia, 1995.
M. Berry, J. Dongarra, and Y. Kim. A Highly Parallel Algorithm for the Reduction of a Nonsymmetric Matrix to Block Upper-Hessenberg Form. Technical Report, LAPACK Working Note 68, University of Tennessee, CS-94-221, 1994.
K. Dackland, E. Elmroth, and B. Kågström. Parallel Block Matrix Factorizations on the Shared Memory Multiprocessor IBM 3090 VF/600J. The International Journal of Supercomputer Applications, 6:69–97, 1992.
J. Dongarra, J. Du Croz, S. Hammarling, and Richard J. Hanson. An Extended Set of FORTRAN Basic Linear Algebra Subroutines. ACM Transactions on Mathematical Software, 14(1):1–17, March 1988.
Gene H. Golub and Charles F. Van Loan. Matrix Computations, 2nd ed. The John Hopkins University Press, Baltimore, Maryland, 1989.
S. Hammarling J. Dongarra, J. Du Croz and I. Duff. A Set of Level 3 Basic Linear Algebra Subprograms. ACM Trans. on Mathematical Software, 16:1–17, 1990.
B. Kågström and C. Van Loan. GEMM-Based Level 3 BLAS. Tech. rep. CTC91TR47, Department of Computer Science, Cornell University, 1989.
C. Lawson, R. Hanson, D. Kincaid, and F. Krogh. Basic Linear Algebra Subprograms for Fortran Usage. ACM Transactions on Mathematical Software, 5:308–323, 1979.
C. B. Moler and G. W. Stewart. An Algorithm for Generalized Matrix Eigenvalue Problems. SIAM J. Num. Anal., 10:241–256, 1973.
A. Pothen and P. Raghavan. Distributed Orthogonal Factorizations: Givens and Householder Algorithms. SIAM J. Sci. Statist. Comput, 10:1113–1134, 1989.
R. Schreiber and C. Van Loan. A Storage Efficient WY Representation for Products of Householder Transformations. SIAM J. Sci. Statist. Comput., 10:53–57, 1989.
Author information
Authors and Affiliations
Editor information
Rights and permissions
Copyright information
© 1996 Springer-Verlag Berlin Heidelberg
About this paper
Cite this paper
Dackland, K., Kågström, B. (1996). Reduction of a regular matrix pair (A, B) to block Hessenberg-triangular form. In: Dongarra, J., Madsen, K., Waśniewski, J. (eds) Applied Parallel Computing Computations in Physics, Chemistry and Engineering Science. PARA 1995. Lecture Notes in Computer Science, vol 1041. Springer, Berlin, Heidelberg. https://doi.org/10.1007/3-540-60902-4_15
Download citation
DOI: https://doi.org/10.1007/3-540-60902-4_15
Published:
Publisher Name: Springer, Berlin, Heidelberg
Print ISBN: 978-3-540-60902-5
Online ISBN: 978-3-540-49670-0
eBook Packages: Springer Book Archive