Given a sparse matrix and a row permutation matrix and a column permutation matrix , this subroutine performs the permutation . The nonzero entries of are stored by rows in a compact form and the user defines the permutation matrices and by index vectors of length .