<div dir="ltr">Sorry, I forgot to say that my computational domain is decomposed with a parallel library using MPI, and the particles are defined with a distributed vector. So, the entries of the matrix are basically from that distributed vector.<div><br></div><div>Thanks,</div><div>Massoud</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Dec 12, 2016 at 1:21 AM, Massoud Rezavand <span dir="ltr"><<a href="mailto:msdrezavand@gmail.com" target="_blank">msdrezavand@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class=""><div dir="ltr"><span style="font-size:12.8px">Thanks,</span><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">as I already discussed with you, the matrix is coming from SPH discretization, which is not fixed on a grid and is changing over time.</div></div></span><div class="gmail_extra"><br><div class="gmail_quote"><span class="">On Mon, Dec 12, 2016 at 1:10 AM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br></span><div><div class="h5"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span><br>
> On Dec 11, 2016, at 6:04 PM, Massoud Rezavand <<a href="mailto:msdrezavand@gmail.com" target="_blank">msdrezavand@gmail.com</a>> wrote:<br>
><br>
> Thank you very much,<br>
><br>
> So, if I am using PetscSplitOwnership() and then MatGetOwnershipRange() to be prepared for preallocation, then MatSetSizes(A, local_size, local_size, N, N) should be called with the calculated  local_size from PetscSplitOwnership() ?<br>
<br>
</span>   Confusion from the two responses. You cannot use MatGetOwnershipRange() for preallocation.<br>
<br>
Without preallocation:<br>
<span><br>
> > PetscInt local_size = PETSC_DECIDE;<br>
> ><br>
> > MatSetSizes(A, local_size, local_size, N, N);<br>
<br>
</span>     MatGetOwnershipRanges(...)<br>
<br>
With preallocation:<br>
<span>> ><br>
> ><br>
> > 2)<br>
> ><br>
> > PetscInt local_size = PETSC_DECIDE;<br>
> ><br>
> > PetscSplitOwnership(PETSC_COMM<wbr>_WORLD, &local_size, &N);<br>
> ><br>
> > MPI_Scan(&local_size, &end_row, 1, MPIU_INT, MPI_SUM, PETSC_COMM_WORLD);<br>
> > begin_row = end_row - local_size;<br>
<br>
</span>    MatMPIAIJSetPreallocation(....<wbr>.).<br>
<br>
<br>
But note that normally if the matrix comes from a discretization on a grid you would not use either approach above. The parallel layout of the grid would determine the local sizes and you won't not obtain them with PetscSplitOwnership() or local_size = PETSC_DECIDE;<br>
<br>
Where is your matrix coming from?<br>
<span class="m_-420401091346126432HOEnZb"><font color="#888888"><br>
   Barry<br>
</font></span><div class="m_-420401091346126432HOEnZb"><div class="m_-420401091346126432h5"><br>
<br>
<br>
> ><br>
> ><br>
<br>
<br>
><br>
> Thanks,<br>
> Massoud<br>
><br>
><br>
> On Mon, Dec 12, 2016 at 12:35 AM, Jed Brown <<a href="mailto:jed@jedbrown.org" target="_blank">jed@jedbrown.org</a>> wrote:<br>
> Massoud Rezavand <<a href="mailto:msdrezavand@gmail.com" target="_blank">msdrezavand@gmail.com</a>> writes:<br>
><br>
> > Dear PETSc team,<br>
> ><br>
> > What is the difference between the following two methods to get the local<br>
> > dimensions of a square matrix A? If they do the same, which one is<br>
> > recommended? Should I use MPI_Scan after both?<br>
><br>
> I would typically use 1 because it's fewer calls and automatically uses<br>
> the correct communicator.  You can use MatGetOwnershipRange() instead of<br>
> manually using MPI_Scan.<br>
><br>
> > 1)<br>
> ><br>
> > PetscInt local_size = PETSC_DECIDE;<br>
> ><br>
> > MatSetSizes(A, local_size, local_size, N, N);<br>
> ><br>
> ><br>
> > 2)<br>
> ><br>
> > PetscInt local_size = PETSC_DECIDE;<br>
> ><br>
> > PetscSplitOwnership(PETSC_COMM<wbr>_WORLD, &local_size, &N);<br>
> ><br>
> > MPI_Scan(&local_size, &end_row, 1, MPIU_INT, MPI_SUM, PETSC_COMM_WORLD);<br>
> > begin_row = end_row - local_size;<br>
> ><br>
> ><br>
> > Thanks in advance,<br>
> > Massoud<br>
><br>
<br>
</div></div></blockquote></div></div></div><br></div>
</blockquote></div><br></div>