<div dir="ltr"><div dir="ltr">On Wed, Oct 18, 2023 at 5:55 AM Enrico <<a href="mailto:degregori@dkrz.de">degregori@dkrz.de</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">Hello,<br>
<br>
I'm trying to use Petsc to solve a linear system in an application. I'm <br>
using the coordinate format to define the matrix and the vector (it <br>
should work better on GPU but at the moment every test is on CPU). After <br>
the call to VecSetValuesCOO, I've noticed that the vector is storing the <br>
data in a different way from my application. For example with two <br>
processes in the application<br>
<br>
process 0 owns cells 2, 3, 4<br>
<br>
process 1 owns cells 0, 1, 5<br>
<br>
But in the vector data structure of Petsc<br>
<br>
process 0 owns cells 0, 1, 2<br>
<br>
process 1 owns cells 3, 4, 5<br>
<br>
This is in principle not a big issue, but after solving the linear <br>
system I get the solution vector x and I want to get the values in the <br>
correct processes. Is there a way to get vector values from other <br>
processes or to get a mapping so that I can do it myself?<br></blockquote><div><br></div><div>By definition, PETSc vectors and matrices own contiguous row blocks. If you want to have another,</div><div>global ordering, we support that with <a href="https://petsc.org/main/manualpages/AO/">https://petsc.org/main/manualpages/AO/</a></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">
Cheers,<br>
Enrico Degregori<br>
</blockquote></div><br clear="all"><div><br></div><span class="gmail_signature_prefix">-- </span><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>