Abstract
The paper describes a parallel algorithm for the LU factorization of sparse matrices on distributed memory machines by using SPMD as programming model and PVM as message passing interface. We address all the difficulties arising in sparse codes, as the fill-in or the dynamic movement of data inside the matrix. The cyclic distribution has been used to evenly distribute the elements onto a mesh of processors, whereas two local storage schemes are proposed: A semi-ordered and two-dimensional linked list, which fulfils better the requirements of the algorithm, and a compressed storage by rows, which behaves better in the use of memory. The properties of the code are extensively analyzed and execution times on the CRAY T3D are presented to illustrate the overall efficiency achieved by our methods.
This work was supported by the Ministry of Education and Science (CICYT) of Spain under project TIC92-0942-C03, by the Human Capital and Mobility programme of the European Union under proyect ERB4050P1921660, and by the Training and Research on Advanced Computing Systems (TRACS) at the Edinburgh Parallel Computing Centre (EPCC)
Preview
Unable to display preview. Download preview PDF.
References
R. Asenjo and E.L. Zapata. Sparse LU factorization on the Cray T3D. Technical report, Computer Architecture Dept. Univ. of Málaga, September 1994. Anonymous ftp: ftp.atc.ctima.uma.es (SpLUCray.ps).
T. A. Davis and P. C. Yew. A nondeterministic parallel algorithm for general unsymmetric sparse LU factorization. SIAM J. Matrix Anal. Appl., 11:383–402, 1990.
I. S. Duff and J. K. Reid. Some design features of a sparse matrix code. ACM Trans. Math. Software, 5:18–35, 1979.
I. S. Duff A. M. Erisman and J. K. Reid. Direct Methods for Sparse Matrices. Oxford University Press, Oxford, U.K., 1986.
G. H. Golub and C. F. Van Loan. Matrix Computation. The Johns Hopkins University Press, Baltimore, MD, 2 edition, 1989.
I. S. Duff R. G. Grimes and J. G. Lewis. Sparse matrix test problems. ACM Trans. Math. Software, 15:1–14, 1989.
A. Geist A. Beguelin J. Dongarra W. Jiang R. Manchek and V. Sunderam. PVM 3 User's Guide and Refernce Manual. Engineering Physic and Mathematics Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, May 1993.
H. M. Markowitz. The elimination form of the inverse and its application to linear programming. Management Sci., 3:255–269, 1957.
S. Booth J. Fisher N. MacDonald P. Maccallum E. Minty and A. Simpson. Parallel Programming on the Cray T3D. Edinburgh Parallel Computing Centre, University of Edinburgh, U.K., September 1994.
R. Asenjo L. F. Romero M. Ujaldón and E. L. Zapata. Sparse block and cyclic data distributions for matrix computations. In High Performance Computing: Technology and Application. Grandinetti et al. (Eds.) Elsevier Science, (to appear).
A. F. van der Stappen R. H. Bisseling and J. G. G. van de Vorst. Parallel sparse LU decomposition on a mesh network of transputers. SIAM J. Matrix Anal. Appl., 14(3):853–879, July 1993.
Z. Zlatev J. Wasniewski and K. Schaumburg. Y12M—Solution of Large and Sparse Systems of Linear Algebraic Equations. In Lecture Notes in Computer Science, number 121, pages 61–77, Berlin, 1981. Springer-Verlang.
Z. Zlatev. Computational Methods for General Sparse Matrices, Mathematics and Its Applications. 65. Kluwer Academin Publisher, Dordrecht, the Netherlands, 1991.
Author information
Authors and Affiliations
Editor information
Rights and permissions
Copyright information
© 1995 Springer-Verlag Berlin Heidelberg
About this paper
Cite this paper
Asenjo, R., Zapata, E.L. (1995). Sparse LU factorization of the Cray T3D. In: Hertzberger, B., Serazzi, G. (eds) High-Performance Computing and Networking. HPCN-Europe 1995. Lecture Notes in Computer Science, vol 919. Springer, Berlin, Heidelberg. https://doi.org/10.1007/BFb0046701
Download citation
DOI: https://doi.org/10.1007/BFb0046701
Published:
Publisher Name: Springer, Berlin, Heidelberg
Print ISBN: 978-3-540-59393-5
Online ISBN: 978-3-540-49242-9
eBook Packages: Springer Book Archive