[petsc-dev] MatCreateSubMatrix question

[Jed Brown]

MatCreateSubMatrices

Mark Adams <mfadams at lbl.gov> writes:

> 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