A Parallelizable Eigensolver for Real Diagonalizable Matrices with Real Eigenvalues

Abstract

In this paper, preliminary research results on a new algorithm for finding all the eigenvalues and eigenvectors of a real diagonalizable matrix with real eigenvalues are presented. The basic mathematical theory behind this approach is reviewed and is followed by a discussion of the numerical considerations of the actual implementation. The numerical algorithm has been tested on thousands of matrices on both a Cray-2 and an IBM RS/6000 Model 580 workstation. The results of these tests are presented. Finally, issues concerning the parallel implementation of the algorithm are discussed. The algorithm's heavy reliance on matrix--matrix multiplication, coupled with the divide and conquer nature of this algorithm, should yield a highly parallelizable algorithm.