  An algorithm for extracting identity submatrices of small rank and pivotal
units from large and sparse matrices is proposed. The procedure has already
been satisfactorily applied for solving the label switching problem in Bayesian
mixture models. Here we introduce it on its own and explore possible
applications in different contexts.
