MC58: Estimate rank and find independent rows/columns of a sparse unsymmetric or rectangular matrix

To estimate the rank and find a nonsingular submatrix of an unsymmetric or rectangular sparse matrix \(\mathbf{A}\) using Gaussian elimination. The main entry performs a sparse \(\mathbf{LU}\) factorization of the matrix optionally using rook pivoting. The factors are not returned.