<div dir="ltr">Thanks a lot for great support.<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jun 29, 2016 at 3:11 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
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.<br>
<div class="HOEnZb"><div class="h5"><br>
<br>
> On Jun 29, 2016, at 1:43 PM, ehsan sadrfaridpour <<a href="mailto:it.sadr@gmail.com">it.sadr@gmail.com</a>> wrote:<br>
><br>
> Hi,<br>
><br>
> I need to have access to most of elements of a parallel MPIAIJ matrix only from 1 process (rank 0).<br>
> I tried to copy or duplicate it to SEQAIJ, but I faced problems.<br>
><br>
> 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.<br>
><br>
> Best,<br>
> Ehsan<br>
><br>
><br>
<br>
</div></div></blockquote></div><br></div>