Minimum degree algorithm

In numerical analysis, the minimum degree algorithm is an algorithm used to permute the rows and columns of a symmetric sparse matrix before applying the Cholesky decomposition, to reduce the number of non-zeros in the Cholesky factor. This results in reduced storage requirements and means that the Cholesky factor can be applied with fewer arithmetic operations. (Sometimes it may also pertain to an incomplete Cholesky factor used as a preconditioner—for example, in the preconditioned conjugate gradient algorithm.)

Minimum degree algorithms are often used in the finite element method where the reordering of nodes can be carried out depending only on the topology of the mesh, rather than on the coefficients in the partial differential equation, resulting in efficiency savings when the same mesh is used for a variety of coefficient values.

Given a linear system

where A is an real symmetric sparse square matrix. The Cholesky factor L will typically suffer 'fill in', that is have more non-zeros than the upper triangle of A. We seek a permutation matrix P, so that the matrix , which is also symmetric, has the least possible fill in its Cholesky factor. We solve the reordered system

The problem of finding the best ordering is an NP-complete problem and is thus intractable, so heuristic methods are used instead. The minimum degree algorithm is derived from a method first proposed by Markowitz in 1959 for non-symmetric linear programming problems, which is loosely described as follows. At each step in Gaussian elimination row and column permutations are performed so as to minimize the number of off diagonal non-zeros in the pivot row and column. A symmetric version of Markowitz method was described by Tinney and Walker in 1967 and Rose later derived a graph theoretic version of the algorithm where the factorization is only simulated, and this was named the minimum degree algorithm. The graph referred to is the graph with n vertices, with vertices i and j connected by an edge when , and the degree is the degree of the vertices. A crucial aspect of such algorithms is a tie breaking strategy when there is a choice of renumbering resulting in the same degree.

A version of the minimum degree algorithm was implemented in the MATLAB function symmmd (where MMD stands for multiple minimum degree), but has now been superseded by a symmetric approximate multiple minimum degree function symamd, which is faster. This is confirmed by theoretical analysis, which shows that for graphs with n vertices and m edges, MMD has a tight upper bound of on its running time, whereas for AMD a tight bound of holds. Cummings, Fahrbach, and Fatehpuria designed an exact minimum degree algorithm with running time, and showed that no such algorithm can exist that runs in time , for any , assuming the strong exponential time hypothesis.

References

  • Cummings, Robert; Fahrbach, Matthew; Fatehpuria, Animesh (2021). "A fast minimum degree algorithm and matching lower bound". Proceedings of the 32nd Annual ACM-SIAM Symposium on Discrete Algorithms: 724–734. arXiv:1907.12119. doi:10.1137/1.9781611976465.45. ISBN 978-1-61197-646-5. S2CID 198968052.
  • George, Alan; Liu, Joseph (1989). "The evolution of the minimum degree ordering algorithm". SIAM Review. 31 (1): 1–19. doi:10.1137/1031001. JSTOR 2030845. OSTI 5686483.
  • Heggernes, P.; Eisenstat, S. C.; Kumfert, G.; Pothen, A. (2001), The Computational Complexity of the Minimum Degree Algorithm (PDF) (Technical report), Institute for Computer Applications in Science and Engineering
  • Markowitz, H. M. (1957). "The elimination form of the inverse and its application to linear programming". Management Science. 3 (3): 255–269. doi:10.1287/mnsc.3.3.255. JSTOR 2627454. Archived from the original on September 24, 2017.
  • Rose, D. J. (1972). "A graph-theoretic study of the numerical solution of sparse positive definite systems of linear equations". Graph Theory and Computing. Academic Press. pp. 183–217. ISBN 0-12-583850-6.
  • Tinney, W. F.; Walker, J. W. (1967). "Direct solution of sparse network equations by optimally ordered triangular factorization". Proc. IEEE. 55 (11): 1801–1809. doi:10.1109/PROC.1967.6011.
This article is issued from Wikipedia. The text is licensed under Creative Commons - Attribution - Sharealike. Additional terms may apply for the media files.