Fast Parallel Algorithms for Solving Triangular Systems of Linear Equations on the Hypercube
โ Scribed by O.H. Ibarra; M.H. Kim
- Publisher
- Elsevier Science
- Year
- 1994
- Tongue
- English
- Weight
- 864 KB
- Volume
- 20
- Category
- Article
- ISSN
- 0743-7315
No coin nor oath required. For personal study only.
โฆ Synopsis
This paper presents efficient hypercube algorithms for solving triangular systems of linear equations by using various matrix partitioning and mapping schemes. Recently, several parallel algorithms have been developed for this problem. In these algorithms, the triangular solver is treated as the second stage of Gauss elimination. Thus, the triangular matrix is distributed by columns (or rows) in a wrap fashion since it is likely that the matrix is distributed this way after an LU decomposition has been done on the matrix. However, the efficiency of the algorithms is low. Our motivation is to develop various data partitioning and mapping schemes for hypercube algorithms by treating the triangular solver as an independent problem. Theoretically, the computation time of our best algorithm is (\left((12 p+1) n^{2}+36 p^{3}-\right.) (\left.28 p^{2}\right) /\left(24 p^{2}\right)), and an upper bound on the communication time is (2 \alpha p \log p(\log n-\log p)+2 \alpha(\log n-\log p-1) \log p+(c n / p-) 2c) ((2 \log p-1)+\log p(c n-c-\alpha)), where (\alpha) is the (communication startup time)/(one entry scanning time), (c) is a constant, (n) is the order of the triangular system and (p) is the number of nodes in the hypercube. Experimental results show that the algorithm is efficient. The efficiency of the algorithm is 0.945 when (p=2, n=) 513, and 0.93 when (p=8, n=1025). 1994 Academic Press. Inc.
๐ SIMILAR VOLUMES
A parallel algorithm for solving the Poisson equation with either Dirichlet or Neumann conditions is presented. The solver follows some of the principles introduced in a previous fast algorithm for evaluating singular integral transforms by Daripa et al. Here we present recursive relations in Fourie
## Abstract An algorithm based on a small matrix approach to the solution of a system of inhomogeneous linear algebraic equations is developed and tested in this short communication. The solution is assumed to lie in an initial subspace and the dimension of the subspace is augmented iteratively by
The order of the matrices involved in several algebraic problems decreases during the solution process. In these cases, parallel algorithms which use adaptive solving block sizes offer better performance results than the ones obtained on parallel algorithms using traditional constant block sizes. Re