A multiprocessor scheme for the singular value decomposition
Conference
·
OSTI ID:5576354
We present a multiprocessor scheme for determining the singular value decomposition of rectangular matrices in which the number of rows is substantially large (or smaller) than the number of columns. In this scheme, we perform an initial QR factorization on the tall matrix (either A or A/sup T/) using a multiprocessor block Householder algorithm. We then use a parallel one-sided Jacobi scheme to orthogonalize the columns of the upper triangular matrix R to yield the factorization RV = U..sigma.., from which the desired singular value decomposition is obtained. Preliminary experiments on an Alliant FX/8 computer system with 8 processors indicate speedups near 5 for our scheme over an optimized implementation of the Linpack/Eispack routines which perform the classical bi-diagonalization technique. Our scheme performs exceptionally well for rank deficient matrices as well as for those rectangular matrices having clustered or multiple singular values, and may be well suited for applications such as real-time signal processing. We present performance results on the Alliant FX/8 and Cray X-MP/48 computer systems with particular emphasis on speedups obtained for our scheme over classical SVD algorithms.
- Research Organization:
- Illinois Univ., Urbana (USA). Center for Supercomputing Research and Development
- DOE Contract Number:
- FG02-85ER25001
- OSTI ID:
- 5576354
- Report Number(s):
- DOE/ER/25001-52; CONF-871251-2; ON: DE88003591
- Country of Publication:
- United States
- Language:
- English
Similar Records
Multiprocessor Jacobi algorithms for dense symmetric eigenvalue and singular value decompositions
Multiprocessor algorithm for the symmetric tridiagonal eigenvalue problem
Optimizing matrix operations on a parallel multiprocessor with a memory hierarchy
Conference
·
Tue Feb 18 23:00:00 EST 1986
·
OSTI ID:7162999
Multiprocessor algorithm for the symmetric tridiagonal eigenvalue problem
Technical Report
·
Tue Jan 07 23:00:00 EST 1986
·
OSTI ID:7008081
Optimizing matrix operations on a parallel multiprocessor with a memory hierarchy
Technical Report
·
Thu Feb 27 23:00:00 EST 1986
·
OSTI ID:7204163