This should be good enough<br><br> <a href="http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/DA/DAGlobalToNaturalAllCreate.html">http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/DA/DAGlobalToNaturalAllCreate.html</a><br>
<br> Matt<br><br><div class="gmail_quote">On Thu, Aug 27, 2009 at 12:51 PM, nicolas aunai <span dir="ltr"><<a href="mailto:nicolas.aunai@gmail.com">nicolas.aunai@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hello,<br>
<br>
<br>
I have a problem using 'VecScatterCreateToZero'. I attach to this<br>
email a test code that illustrates the problem. I create a 2D DA then<br>
I ask for the coordinates (that have been set just before), which are<br>
returned in a global vector. I want to write the coordinates in a<br>
file, so I use 'VecScatterCreateToZero' to put all the values in<br>
sequential vector. It seems to work fine but sometimes (depending on<br>
the number of proc used) the coordinates in the sequential vectors are<br>
ordered differently.<br>
<br>
to run the code and see the problem :<br>
<br>
mpiexec -n 4 ./vecdacoord > problem.txt<br>
<br>
I have attached problem.txt too.<br>
<br>
Is someone can explain to me how I could obtain my coordinates in a<br>
natural order ?<br>
<br>
<br>
Thx<br>
</blockquote></div><br><br clear="all"><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>