Abstract
In this paper, we present heuristic techniques for the reduction of the bandwidth of a sparse matrix as well as for the reduction of the cost of the associated Cholesky factorization. Our algorithms are inspired by the spectral method of Barnard, Pothen and Simon (1995), which derives a permutation for reducing the envelope-size of a sparse matrix by computing the second eigenvector of the associated Laplacian matrix. Two main modifications of that method are proposed and tested. The first is based on the experimental observation that it is often preferable to perform only few iterations of an iterative method converging to the second eigenvector; the second is the introduction of a weighted Laplacian. These simple ideas allow us to obtain a family of spectral methods that have been carefully tested on a set of matrices whose size ranges from few hundred to one million.
Similar content being viewed by others
References
G. Alaghband, Parallel pivoting combined with parallel reduction and fill-in control, Parallel Comput. 11 (1989) 210–221.
S.T. Barnard and H.D. Simon, A fast multilevel implementation of recursive spectral bisection for partitioning unstructured problems, Conc.: Pract. Exp. 6(2) (1994) 101–117.
S.T. Barnard, A. Pothen and H.D. Simon, A spectral algorithm for envelope reduction of sparse matrices, J. Numer. Lin. Alg. Appl. 2 (1995) 311–334.
M. Berry, B. Hendrickson and P. Raghavan, Sparse matrix reordering schemes for browsing hypertext, in: The Mathematics of Numerical Analysis, eds. S. Smale, J. Renegar, M. Shub, Lectures in Applied Mathematics, Vol. 32 (AMS, Providence, RI, 1996) pp. 99–123.
E.H. Cuthill and J. McKee, Reducing the bandwidth of sparse symmetric matrices, in: Proc. 24th Nat. Conf. Assoc. Comp. Mach. (ACM Publications, 1969) pp. 157–172.
T. Davis, Sparse matrix collection, http://www.cise.ufl.edu/∼davis/sparse/, ftp://ftp.cise.ufl.edu/pub/faculty/davis/matrices (1994). NA Digest, Vol. 92, No. 42 (October 16, 1994); NA Digest, Vol. 96, No. 28 (July 23, 1996); and NA Digest, Vol. 97, No. 23 (June 7, 1997).
E.F. D'Azevedo, P.A. Forsyth and W.P. Tang, Ordering methods for preconditioned conjugate gradient methods applied to unstructured grid problems, SIAM J. Matrix Anal. Appl. 13 (1992) 944–961.
I.S. Duff and G.A. Meurant, The effect of ordering on preconditioned conjugate gradients, BIT 29 (1989) 635–657.
A. Esposito, S. Fiorenzo Catalano, F. Malucelli and L. Tarricone, A new matrix bandwidth reduction algorithm, Operation Res. Lett. (1999), to appear.
A. Esposito, F. Malucelli and L. Tarricone, Bandwidth reduction of sparse symmetric and unsymmetric matrices: an experimental comparison of new heuristics, http://dvorak.istel.ing.unipg.it.
M. Fiedler, A property of non-negative symmetric matrices and its application to graph theory, Czech. Mat. J. 25 (1975) 619–633.
A. George and A. Pothen, An analysis of spectral envelope-reduction via quadratic assignment problems, SIAM J. Matrix Anal. Appl. 18(3) (1997) 706–732.
N. Gibbs, W. Poole and P. Stockmeyer, An algorithm for reducing the bandwidth and profile of a sparse matrix, SIAM J. Numer. Anal. 13 (1976) 235–251.
N. Gibbs, W. Poole and P. Stockmeyer, A comparison of several bandwidth and profile reduction algorithms, ACM Trans. Math. Software 2 (1976) 322–330.
M. Juvan and B. Mohar, Optimal linear labelings and eigenvalues of graphs, Discr. Appl. Math. 36 (1992) 153–168.
M. Juvan and B. Mohar, Laplace eigenvalues and bandwidth-type invariants of graphs, J. Graph Theory 17 (1993) 393–407.
J.W.H. Liu and H. Sherman, Comparative analysis of the Cuthill-Mckee and the reverse Cuthill-Mckee ordering algorithms for sparse matrices, SIAM J. Numer. Anal. 13 (1976) 198–213.
G. Kumfert and A. Pothen, Two improved algorithms for envelope and wavefront reduction, BIT 37(3) (1997) 559–590.
J.G. Lewis, Implementations of the Gibbs-Poole-Stockmeyer and Gibbs-King algorithms, ACM Trans. Math. Software 8 (1982) 180–189.
C. Papadimitriou, The NP-completeness of the bandwidth minimization problem, Computing 16 (1976) 263–270.
B.N. Parlett, The Symmetric Eigenvalue Problem (Prentice Hall, Englewood Cliffs, NJ, 1980).
Y. Saad, Numerical Methods for Large Eigenvalues Problems (Manchester Univ. Press, Manchester, 1992).
J.A. Scott, A new row ordering strategy for frontal solvers, Num. Lin. Alg. Appl. 6(3) (1999) 189–211.
S.W. Sloan, An algorithm for profile and wavefront reduction of sparse matrices, Int. J. Numer. Math. Engrg. 23 (1986) 239–251.
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Del Corso, G.M., Romani, F. Heuristic Spectral Techniques for the Reduction of Bandwidth and Work-Bound of Sparse Matrices. Numerical Algorithms 28, 117–136 (2001). https://doi.org/10.1023/A:1014082430392
Issue Date:
DOI: https://doi.org/10.1023/A:1014082430392