I have a parallel matrix and I want to extract parts of it locally to create small local matrices.
I use PETSC_COMM_SELF for the ISs, but that does not seem to be enough to tell MatCreateSubMatrix that I want a PETSC_COMM_SELF (sub)matrix. How should I do this? I am not doing any communication in this process. Could I just wrap the local part of the parallel matrix into a PETSC_COMM_SELF matrix? Thanks, Mark
