A fully parallel algorithm for the symmetric Eigenvalue problem
Journal Article
·
· SIAM J. Sci. Stat. Comput.; (United States)
In this paper the authors present a parallel algorithm for the symmetric algebraic eigenvalue problem. The algorithm is based upon a divide and conquer scheme suggested by Cuppen for computing the eigensystem of a symmetric tridiagonal matrix. They extend this idea to obtain a parallel algorithm that retains a number of active parallel processes that is greater than or equal to the initial number throughout the course of the computation. They give a new deflation technique which together with a robust root finding technique will assure computation of an eigensystem to full accuracy in the residuals and in the orthogonality of eigenvectors. A brief analysis of the numerical properties and sensitivity to round off error is presented to indicate where numerical difficulties may occur. The algorithm is able to exploit parallelism at all levels of the computation and is well suited to a variety of architectures. Computational results are presented for several machines. These results are very encouraging with respect to both accuracy and speedup. A surprising result is that the parallel algorithm, even when run in serial mode, can be significantly faster than the previously best sequential algorithm on large problems, and is effective on moderate size problems when run in serial mode.
- Research Organization:
- Mathematics and Computer Science Div., Argonne National Lab., 9700 South Cass Avenue, Argonne, IL
- DOE Contract Number:
- W-31109-ENG-38; AC05-84OR21400
- OSTI ID:
- 6849319
- Journal Information:
- SIAM J. Sci. Stat. Comput.; (United States), Journal Name: SIAM J. Sci. Stat. Comput.; (United States) Vol. 8:2; ISSN SIJCD
- Country of Publication:
- United States
- Language:
- English
Similar Records
Fully parallel algorithm for the symmetric eigenvalue problem
Implementation of a fully parallel algorithm for the symmetric eigenvalue problem
Parallel solution of the symmetric tridiagonal eigenproblem
Conference
·
Tue Dec 31 23:00:00 EST 1985
·
OSTI ID:5870703
Implementation of a fully parallel algorithm for the symmetric eigenvalue problem
Conference
·
Fri Aug 01 00:00:00 EDT 1986
·
OSTI ID:7018147
Parallel solution of the symmetric tridiagonal eigenproblem
Thesis/Dissertation
·
Sat Dec 31 23:00:00 EST 1988
·
OSTI ID:6155839