VecScatterCreateToZero

Matthew Knepley knepley at gmail.com
Thu Aug 27 13:20:37 CDT 2009


This should be good enough


http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/DA/DAGlobalToNaturalAllCreate.html

     Matt

On Thu, Aug 27, 2009 at 12:51 PM, nicolas aunai <nicolas.aunai at gmail.com>wrote:

> Hello,
>
>
> I have a problem using 'VecScatterCreateToZero'. I attach to this
> email a test code that illustrates the problem. I create a 2D DA then
> I ask for the coordinates (that have been set just before), which are
> returned in a global vector. I want to write the coordinates in a
> file, so I use 'VecScatterCreateToZero' to put all the values in
> sequential vector. It seems to work fine but sometimes (depending on
> the number of proc used) the coordinates in the sequential vectors are
> ordered differently.
>
> to run the code and see the problem :
>
> mpiexec -n 4 ./vecdacoord > problem.txt
>
> I have attached problem.txt too.
>
> Is someone can explain to me how I could obtain my coordinates in a
> natural order ?
>
>
> Thx
>



-- 
What most experimenters take for granted before they begin their experiments
is infinitely more interesting than any results to which their experiments
lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20090827/bae3aef9/attachment.htm>


More information about the petsc-users mailing list