Dear all, How can I reorder a sparse matrix in Petsc? Sth like in a dirty MATLAB like pseudo-code: A = diag([1 3 4 5]); neworder = [2 4 1 3]; B = A(neworder, neworder) I know I will destroy the sparsity pattern, but for my application there is no other way, the matrices should be reordered for some partitioning. Best regards, Umut