In the DA we have some idea of topology, so we can make a sensible reordered local vector. The VecGhost is a lower level construct. All it does is maintain a vector of this form
[ local values | ghost values ] The update maps from local values to ghost values with a scatter. The VecGhostGetLocalForm() just pulls out the pointer. We could add code which maps into a different local vector using an ordering, but this can already be done with an AO as well. Matt On 6/20/06, Thomas Geenen <geenen at gmail.com> wrote: > > On Tuesday 20 June 2006 13:32, you wrote: > > > I now expect the vector to be in the local form. That means that the > > > positions > > > that I created as ghost points are filled with the solution > corresponding > > > to > > > the global solution. I expect that they are inserted in the position > in > > > my local vector corresponding to there corresponding global position. > > > > Not sure what you mean here. The local form is the same as the global > > vector > > in the owned spots, and then has the ghost points at the end. > oke so that means (4,5,3) ?? > So if I need the solution in the same positions as the ones that i filled > with > VecSetValuesLocal I will have to reorder the vector > What does VecGhostGetLocalForm do? inserts the global solution in the > ghost > points at the end? I thought VecGhostUpdateBegin VecGhostUpdateEnd did > that. > > thomas > > Matt > -- "Failure has a thousand explanations. Success doesn't need one" -- Sir Alec Guiness -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20060620/6c000784/attachment.htm>
