[petsc-users] How to have a local copy (sequential) of a parallel matrix

[ehsan sadrfaridpour]

Thanks a lot for great support.

On Wed, Jun 29, 2016 at 3:11 PM, Barry Smith <bsmith at mcs.anl.gov> wrote:

>
>    MatGetSubmatrices() just have the first process request all the rows
> and columns and the others request none. You can use ISCreateStride() to
> create the ISs without having to make an array of all the indices.
>
>
> > On Jun 29, 2016, at 1:43 PM, ehsan sadrfaridpour <it.sadr at gmail.com>
> wrote:
> >
> > Hi,
> >
> > I need to have access to most of elements of a parallel MPIAIJ matrix
> only from 1 process (rank 0).
> > I tried to copy or duplicate it to SEQAIJ, but I faced problems.
> >
> > How can I have a local copy of a matrix which is distributed on multiple
> process? I don't want to update the matrix, and the read-only version of it
> would be enough.
> >
> > Best,
> > Ehsan
> >
> >
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160629/808b16f1/attachment.html>