<div dir="ltr"><div dir="ltr">On Wed, Jun 10, 2020 at 9:08 AM Eda Oktay <<a href="mailto:eda.oktay@metu.edu.tr">eda.oktay@metu.edu.tr</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear Matt,<br>
<br>
Matthew Knepley <<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>>, 10 Haz 2020 Çar, 16:03 tarihinde<br>
şunu yazdı:<br>
><br>
> On Wed, Jun 10, 2020 at 8:56 AM Eda Oktay <<a href="mailto:eda.oktay@metu.edu.tr" target="_blank">eda.oktay@metu.edu.tr</a>> wrote:<br>
>><br>
>> Hi all,<br>
>><br>
>> I am trying to get all the rows of a parallel matrix as individual<br>
>> vectors. For instance, if I have 72*4 matrix, I want to get 72<br>
>> different vectors having size 4.<br>
>><br>
>> As far as I understood, MatGetRow is only for local rows, so<br>
>> MatGetOwnershipRange is used, however, when I tried this one, I<br>
>> couldn't get the whole and desired row vectors.<br>
>><br>
>> In MatGetRow explanation, it is written that I should use<br>
>> MatCreateSubMatrices first, then use MatGetRow. But I couldn't<br>
>> understand to which extent I should create submatrices. I just need to<br>
>> have all 72 rows as 72 different vectors each having 4 elements.<br>
><br>
><br>
> 1) For sparse matrices, the storage is always divided by row, so that values can only be retrieved for local rows with MatGetRow()<br>
><br>
> 2) Is this matrix sparse? It sounds like it is dense.<br>
<br>
Matrix is dense.<br>
<br>
><br>
> 3) Are you asking to get all matrix values on all processes? If so, I think the easiest thing to do is first wrap a Vec around the<br>
>     values, then use VecScatterToAll(), then wrap each one in a MatDense again.<br>
<br>
Yes, I want all row vectors on all processes. In a dense matrix,<br>
should I still wrap a Vec around the values? I know I should use<br>
scatter but I couldn't even wrap a Vec around them.<br></blockquote><div><br></div><div>I would do</div><div><br></div><div>  MatGetSize(&N);</div><div>  MatGetLocalSize(&m);</div><div>  <a href="https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatDenseGetArray.html">https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatDenseGetArray.html</a></div><div>  <create vector of local size m*N></div><div>  <a href="https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Vec/VecPlaceArray.html">https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Vec/VecPlaceArray.html</a><br></div><div>  <a href="https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Vec/VecScatterCreateToAll.html">https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Vec/VecScatterCreateToAll.html</a></div><div>  <do scatter></div><div>  <a href="https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Vec/VecResetArray.html#VecResetArray">https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Vec/VecResetArray.html#VecResetArray</a></div><div>  <a href="https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatCreateDense.html">https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatCreateDense.html</a></div><div>  <use it></div><div>  <destroy matrix></div><div>  <destroy vector from CreateToAll></div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
Thanks so much!<br>
<br>
Eda<br>
<br>
><br>
>   Thanks,<br>
><br>
>      Matt<br>
><br>
>><br>
>> Thanks!<br>
>><br>
>> Eda<br>
><br>
><br>
><br>
> --<br>
> What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
> -- Norbert Wiener<br>
><br>
> <a href="https://www.cse.buffalo.edu/~knepley/" rel="noreferrer" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>