Reduction of a general matrix to tridiagonal form using a hypercube multiprocessor

An algorithm for reducing a nonsymmetric matrix to tridiagonal form as a first step toward finding its eigenvalues is described. The algorithm uses a variation of threshold pivoting, where at each step, the pivot is chosen to minimize the maximum entry in the transformation matrix that reduces the next column and row of the matrix. Situations are given where the tridiagonalization process breaks down, and two recovery methods are presented for these situations. Although no existing tridiagonalization algorithm is guaranteed to succeed, this algorithm is found to be very robust and fast in practice. A gradual loss of similarity is also observed as the order of the matrix increases.

[1]  G. A. Geist,et al.  Stabilized Gaussian reduction of an arbitrary matrix to tridiagonal form , 1989 .

[2]  Christopher Strachey,et al.  The Reduction of a Matrix to Codiagonal Form by Eliminations , 1961, Comput. J..

[3]  Jack Dongarra,et al.  Computing the eigenvalues and eigenvectors of a general matrix by reduction to general tridiagonal form , 1990 .

[4]  Daniel L. Boley Practical iterative methods for large scale computations : proceedings of the Minnesota Supercomputer Institute Workshop on Practical Iterative Methods for Large Scale Computations, Minneapolis, 23-25 October 1988 , 1989 .

[5]  C. Donald LaBudde,et al.  The reduction of an arbitrary real square matrix to tridiagonal form using similarity transformations , 1963 .

[6]  J. Cullum,et al.  Lanczos algorithms for large symmetric eigenvalue computations , 1985 .

[7]  B. S. Garbow,et al.  Matrix Eigensystem Routines — EISPACK Guide , 1974, Lecture Notes in Computer Science.

[8]  Robert Todd Gregory,et al.  A collection of matrices for testing computational algorithms , 1969 .

[9]  J. G. F. Francis,et al.  The QR Transformation - Part 2 , 1962, Comput. J..

[10]  Self-equivalent flows associated with the generalized eigenvalue problem , 1989 .

[11]  J. H. Wilkinson,et al.  IMPROVING THE ACCURACY OF COMPUTED EIGENVALUES AND EIGENVECTORS , 1983 .

[12]  P. Lax Mathematics and computing , 1986 .

[13]  J. H. Wilkinson The algebraic eigenvalue problem , 1966 .

[14]  S. Kaniel,et al.  The ELR Method for Computing the Eigenvalues of a General Matrix , 1981 .