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>
