Parallel frontal solvers for large sparse linear systems

Many applications in science and engineering give rise to large sparse linear systems of equations that need to be solved as efficiently as possible. As the size of the problems of interest increases, it can become necessary to consider exploiting multiprocessors to solve these systems. We report on the design and development of parallel frontal solvers for the numerical solution of large sparse linear systems. Three codes have been developed for the mathematical software library HSL (www.cse.clrc.ac.uk/Activity/HSL). The first is for unsymmetric finite-element problems; the second is for symmetric positive definite finite-element problems; and the third is for highly unsymmetric linear systems such as those that arise in chemical process engineering. In each case, the problem is subdivided into a small number of loosely connected subproblems and a frontal method is then applied to each of the subproblems in parallel. We discuss how our software is designed to achieve the goals of portability, ease of use, efficiency, and flexibility, and illustrate the performance using problems arising from real applications.

[1]  Jennifer A. Scott,et al.  Two-stage ordering for unsymmetric parallel row-by-row frontal solvers , 2001 .

[2]  I. Duff,et al.  Direct Methods for Sparse Matrices , 1987 .

[3]  I. Duff,et al.  The Use of multiple fronts in Gaussian elimination , 1994 .

[4]  Ivo Marek,et al.  Czech–US workshop on iterative methods and parallel computing , 1999 .

[5]  I. Duff,et al.  The Design of MA48, a code for the direct solution of sparse unsymmetric linear systems of equations , 1995 .

[6]  John K. Reid,et al.  The design of MA48: a code for the direct solution of sparse unsymmetric linear systems of equations , 1996, TOMS.

[7]  Jennifer A. Scott,et al.  MA42 - A new frontal code for solving sparse unsymmetric systems , 1993 .

[8]  E. Cuthill,et al.  Reducing the bandwidth of sparse symmetric matrices , 1969, ACM '69.

[9]  Jennifer A. Scott,et al.  The design of a new frontal code for solving sparse, unsymmetric systems , 1996, TOMS.

[10]  Jennifer A. Scott,et al.  Row ordering for frontal solvers in chemical process engineering , 2000 .

[12]  Jennifer A. Scott,et al.  A new row ordering strategy for frontal solvers , 1999, Numer. Linear Algebra Appl..

[13]  Jennifer A. Scott,et al.  A parallel frontal solver for finite element applications , 2001 .

[14]  R. J. Blake,et al.  A multilevel unsymmetric matrix ordering algorithm for parallel process simulation , 2000 .

[15]  Mark A. Stadtherr,et al.  A parallel block frontal solver for large scale process simulation : Reordering effects , 1997 .

[16]  H. V. D. Vorst,et al.  A parallel linear system solver for circuit simulation problems , 2000 .

[17]  Iain S. Duff,et al.  Design features of a frontal code for solving sparse unsymmetric linear systems out-of-core , 1984 .

[18]  Message P Forum,et al.  MPI: A Message-Passing Interface Standard , 1994 .

[19]  Jennifer A. Scott.  The design of a portable parallel frontal solver for chemical process engineering problems , 2001 .

[20]  Bruce Hendrickson,et al.  The Chaco user`s guide. Version 1.0 , 1993 .

[21]  John Reid,et al.  Ordering symmetric sparse matrices for small profile and wavefront , 1999 .

[22]  E. Lui,et al.  A parallel frontal solver on the alliant FX/80 , 1991 .

[23]  Gene H. Golub,et al.  A parallel balance scheme for banded linear systems , 2001, Numer. Linear Algebra Appl..

[24]  Jennifer A. Scott,et al.  ELEMENT RESEQUENCING FOR USE WITH A MULTIPLE FRONT ALGORITHM , 1996 .

[25]  Jennifer A. Scott.  On ordering elements for a frontal solver , 1999 .

[26]  P. Hood,et al.  Frontal solution program for unsymmetric matrices , 1976 .

[27]  Bruce M. Irons,et al.  A frontal solution program for finite element analysis , 1970 .

[28]  Scott W. Sloan,et al.  A FORTRAN program for profile and wavefront reduction , 1989 .

[29]  Iain S. Duff,et al.  Concurrent Multifrontal Methods: Shared Memory, Cache, and Frontwidth Issues , 1987 .

[30]  Jack J. Dongarra,et al.  A set of level 3 basic linear algebra subprograms , 1990, TOMS.

[31]  Henk A. van der Vorst,et al.  A parallel linear system solver for circuit simulation problems , 2000, Numer. Linear Algebra Appl..

[32]  Jennifer A. Scott,et al.  A parallel direct solver for large sparse highly unsymmetric linear systems , 2004, TOMS.