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.