Hi all,<div><br></div><div>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?</div>
<div><br></div><div>Thank you!</div><div><br></div><div>Michel Cancelliere</div><div><br></div><div><br></div>