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

Reply via email to