[petsc-users] VecView doesn't work properly with DA global vectors

[Jed Brown]

You can set the "native" format to avoid doing this reordering.

http://www.mcs.anl.gov/petsc/petsc-2/snapshots/petsc-dev/docs/manualpages/Viewer/PetscViewerSetFormat.html

On Mar 22, 2011 2:28 AM, "Barry Smith" <bsmith at mcs.anl.gov> wrote:


On Mar 21, 2011, at 8:11 PM, Алексей Рязанов wrote:

> Hi
>
> Is it ok, that VecView doesn't work ...
  VecView() for DM global vectors present the vector entries in the global
natural ordering. Since PETSc divides the domain with DA into squareish
blocks when you print the values on a process (by calling DAVecGetArray) it
is not in the global natural ordering. Hence the two are different.

 Barry


> Thank you very much!
>
> Alexey Ryazanov.
>
>
>
> CODE:
>
> #include "petscksp.h"
> #incl...
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110322/7c0a7527/attachment.htm>