Nested-dissection orderings for sparse LU with partial pivoting (2000) [9 citations — 2 self]
Abstract:
Abstract. We describe the implementation and performance of a novel fill-minimization ordering technique for sparse LU factorization with partial pivoting. The technique was proposed by Gilbert and Schreiber in 1980 but never implemented and tested. Like other techniques for ordering sparse matrices for LU with partial pivoting, our new method preorders the columns of the matrix (the row permutation is chosen by the pivoting sequence during the numerical factorization). Also like other methods, the column permutation Q that we select is a permutation that minimizes the fill in the Cholesky factor of Q T A T AQ. Unlike existing column-ordering techniques, which all rely on minimum-degree heuristics, our new method is based on a nested-dissection ordering of A T A. Our algorithm, however, never computes a representation of A T A, which can be expensive. We only work with a representation of A itself. Our experiments demonstrate that the method is efficient and that it can reduce fill significanly relative to the best existing methods. The method reduces the LU running time on some very large matrices (tens of millions of nonzeros in the factors) by more than a factor of 2. 1
Citations
| 141 | A Supernodal Approach to Sparse Partial Pivoting – Demmel, Eisenstat, et al. - 1999 |
| 108 | The evolution of the minimum degree ordering algorithm – George, Liu - 1989 |
| 89 | An unsymmetric-pattern multifrontal method for sparse LU factorization – Davis, Duff - 1997 |
| 49 | Computing the block triangular form of a sparse matrix – Pothen, Fan - 1990 |
| 48 | A combined unifrontal/multifrontal method for unsymmetric sparse matrices – Davis, Duff - 1999 |
| 36 | Predicting structure in sparse matrix computations – GILBERT - 1986 |
| 31 | Sparse Gaussian Elimination on High Performance Computers – Li - 1996 |
| 30 | SPOOLES: An Object-oriented Sparse Matrix Library, in – Ashcraft, Grimes - 1999 |
| 14 | Graph Separator Theorems and Sparse Gaussian Elimination – Gilbert - 1980 |
| 9 | On the complexity of sparse QR and LU factorization of finite-element matrices – George, Ng - 1988 |
| 6 | An Approximate Minimum Degree Column Ordering Algorithm – Larimore - 1998 |
| 2 | Nested dissection with partial pivoting – Gilbert, Schreiber - 1982 |
