Cuthill–McKee algorithm

In the mathematical subfield of matrix theory the Cuthill–McKee algorithm is an algorithm to reduce the bandwidth of sparse symmetric matrices. The reverse Cuthill–McKee algorithm (RCM) is the same algorithm but with the resulting index numbers reversed. In practice this is generally a better solution.

Algorithm

Given a symmetric n×n matrix we visualize the matrix as the adjacency matrix of a graph. The Cuthill-McKee algorithm is then a relabeling of the vertices of the graph to reduce the bandwidth of the adjacency matrix.

The algorithm produces an ordered n-tuple R of vertices which is the new order of the vertices.

First we choose a peripheral vertex x and set R:= ({x}).

Then for i=1,2,... we iterate the following steps while |R| < n

• Construct the adjacency set ${\displaystyle A_{i}}$  of ${\displaystyle R_{i}}$  (with ${\displaystyle R_{i}}$  the i-th component of R) and exclude the vertices we already have in R

${\displaystyle A_{i}:=\operatorname {Adj} (R_{i})\setminus R}$ 

• Sort ${\displaystyle A_{i}}$  with ascending vertex order.
• Append ${\displaystyle A_{i}}$  to the Result set R.

In other words, number the vertices according to a particular breadth-first traversal where neighboring vertices are visited in order from lowest to highest vertex order.

References

E. Cuthill and J. McKee. Reducing the bandwidth of sparse symmetric matrices In Proc. 24th Nat. Conf. ACM, pages 157-172, 1969.

External links

• Cuthill–McKee documentation for the Boost C++ Libraries.
• A detailed description of the Cuthill–McKee algorithm.