论文信息 - A ScaLAPACK-Style Algorithm for Reducing a Regular Matrix Pair to Block Hessenberg-Triangular Form

A ScaLAPACK-Style Algorithm for Reducing a Regular Matrix Pair to Block Hessenberg-Triangular Form

A parallel algorithm for reduction of a regular matrix pair (A, B) to block Hessenberg-triangular form is presented. It is shown how a sequential elementwise algorithm can be reorganized in terms of blocked factorizations and matrix-matrix operations. Moreover, this LAPACK-style algorithm is straightforwardly extended to a parallel algorithm for a rectangular 2D processor grid using parallel kernels from ScaLAPACK. A hierarchical performance model is derived and used for algorithm analysis and selection of optimal blocking parameters and grid sizes.

Krister Dackland | Bo Kågström

[1] Krister Dackland,et al. Reduction of a Regular Matrix Pair (A, B) to Block Hessenberg Triangular Form , 1995, PARA.

[2] Gene H. Golub,et al. Matrix computations , 1983 .

[3] Jack J. Dongarra,et al. A Parallel Algorithm for the Reduction of a Nonsymmetric Matrix to Block Upper-Hessenberg Form , 1995, Parallel Comput..

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

[5] Wayne H. Enright,et al. A note on the efficient solution of matrix pencil systems , 1978 .

[6] Krister Dackland,et al. A Hierarchical Approach for Performance Analysis of ScaLAPACK-Based Routines Using the Distributed Linear Algebra Machine , 1996, PARA.

[7] G. Stewart,et al. An Algorithm for Generalized Matrix Eigenvalue Problems. , 1973 .

[8] Jack Dongarra,et al. ScaLAPACK Users' Guide , 1987 .