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