In this article the orthogonal decomposition of large sparse matrices on a hypercube multiprocessor is considered. The proposed algorithm offers a parallel implementation of the general row merging scheme for sparse Givens transformations recently developed by Joseph Liu. The proposed parallel algorithm is novel in several aspects. First, a new mapping strategy whose goal is to reduce the communication cost and balance the work load during the entire computing process is proposed. Second, a new sequential algorithm for merging two upper trapezoidal matrices (possibly of different dimensions) is described, wherein the order of computation is different from the standard Givens scheme, and is more suitable for parallel implementation. Third, it is shown that the hypercube network can be employed as a multi-loop multiprocessor. The performance of the parallel algorithm applied to a model problem is analyzed and computation/communication complexity results are presented. Finally it is shown that the parallel s...
[1]
O. Østerby,et al.
Direct Methods for Sparse Matrices
,
1983
.
[2]
Robert Schreiber,et al.
A New Implementation of Sparse Gaussian Elimination
,
1982,
TOMS.
[3]
George Ostrouchov.
Symbolic givens reduction and row-ordering in large sparse least squares problems
,
1987
.
[4]
Eleanor Chin-Hwa Lee Chu.
Orthogonal decomposition of dense and sparse matrices on multiprocessors
,
1988
.
[5]
A. George,et al.
Solution of sparse linear least squares problems using givens rotations
,
1980
.
[6]
M. Yannakakis.
Computing the Minimum Fill-in is NP^Complete
,
1981
.
[7]
Z. Zlatev.
Comparison of two pivotal strategies in sparse plane rotations
,
1980
.
[8]
E. Reingold,et al.
Combinatorial Algorithms: Theory and Practice
,
1977
.
[9]
G. A. Geist,et al.
A partitioning strategy for parallel sparse Cholesky factorization
,
1988
.