[petsc-users] Question about matrix permutation

[Michel Cancelliere]

Hi all,

I'm implementing a preconditioner that involves to extract some submatrices,
the system matrix is of type seqbaij, and i need to permute some rows and
columns before extract the submatrices. I'm trying to use MatPermute() but i
should need create and additional matrix that i don't want (and it seems not
to work with baij matrices), I took a look to MatGetOrdering but not sure if
it fit to my case. there is any alternative to use something like
MatPermute() in a seqbaij matrix but without create a new matrix?

Thank you!

Michel Cancelliere
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20100129/4bfae1b7/attachment.htm>