Hybridizing Nested Dissection and Halo Approximate Minimum Degree for Efficient Sparse Matrix Ordering

Minimum Degree and Nested Dissection are the two most popular reordering schemes used to reduce ll-in and operation count when factoring and solving sparse matrices. Most of the state-of-the-art ordering packages hybridize these methods by performing incomplete Nested Dissection and ordering by Minimum Degree the subgraphs associated with the leaves of the separation tree, but to date only loose couplings have been achieved, resulting in poorer performance than could have been expected. This paper presents a tight coupling of the Nested Dissection and Halo Approximate Minimum Degree algorithms, which allows the MinimumDegree algorithm to use exact degrees on the boundaries of the subgraphs passed to it, and to yield back not only the ordering of the nodes of the subgraph, but also the amalgamated assembly subtrees, for e cient block computations. Experimental results show the performance improvement, both in terms of ll-in reduction and concurrency during numerical factorization.

[1]  Joseph W. H. Liu The Minimum Degree Ordering with Constraints , 1989 .

[2]  Patrick R. Amestoy,et al.  An Approximate Minimum Degree Ordering Algorithm , 1996, SIAM J. Matrix Anal. Appl..

[3]  Joseph W. H. Liu,et al.  Robust Ordering of Sparse Matrices using Multisection , 1998 .

[4]  Joseph W. H. Liu,et al.  Modification of the minimum-degree algorithm by multiple elimination , 1985, TOMS.

[5]  Vipin Kumar,et al.  A Fast and High Quality Multilevel Scheme for Partitioning Irregular Graphs , 1998, SIAM J. Sci. Comput..

[6]  Michael T. Heath,et al.  Sparse Cholesky factorization on a local-memory multiprocessor , 1988 .

[7]  John B. Shoven,et al.  I , Edinburgh Medical and Surgical Journal.

[8]  H. Wilf,et al.  Direct Solutions of Sparse Network Equations by Optimally Ordered Triangular Factorization , 1967 .

[9]  Jean Roman,et al.  Sparse Matrix Ordering with SCOTCH , 1997, HPCN Europe.

[10]  J. Pasciak,et al.  Computer solution of large sparse positive definite systems , 1982 .

[11]  Cleve Ashcraft,et al.  Compressed Graphs and the Minimum Degree Algorithm , 1995, SIAM J. Sci. Comput..

[12]  Robert Schreiber,et al.  Scalability of Sparse Direct Solvers , 1993 .

[13]  Pascal Hénon,et al.  A Mapping and Scheduling Algorithm for Parallel Sparse Fan-In Numerical Factorization , 1999, Euro-Par.

[14]  Alan George,et al.  The Evolution of the Minimum Degree Ordering Algorithm , 1989, SIAM Rev..

[15]  Kevin Barraclough,et al.  I and i , 2001, BMJ : British Medical Journal.

[16]  Norman E. Gibbs,et al.  A Comparison of Several Bandwidth and Profile Reduction Algorithms , 1976, TOMS.

[17]  P. Charrier,et al.  Algorithmique et calculs de complexité pour un solveur de type dissections emboîtées , 1989 .

[18]  Patrick Amestoy,et al.  Hybridizing Nested Dissection and Halo Approximate Minimum Degree for Efficient Sparse Matrix Ordering , 1999, IPPS/SPDP Workshops.

[19]  D. Rose,et al.  Generalized nested dissection , 1977 .

[20]  Bruce Hendrickson,et al.  Improving the Run Time and Quality of Nested Dissection Ordering , 1998, SIAM J. Sci. Comput..

[21]  Vipin Kumar,et al.  Scalable Parallel Algorithms for Sparse Linear Systems , 1997 .

[22]  Joseph W. H. Liu,et al.  A Comparison of Three Column-Based Distributed Sparse Factorization Schemes. , 1990 .

[23]  Bruce Hendrickson,et al.  Effective Sparse Matrix Ordering: Just Around the BEND , 1997, PPSC.

[24]  Patrick R. Amestoy,et al.  Multifrontal parallel distributed symmetric and unsymmetric solvers , 2000 .