[petsc-dev] Vector operations on ghosted vectors

[Matthew Knepley]

On Wed, Feb 27, 2013 at 4:56 PM, John Mousel <john.mousel at gmail.com> wrote:

> Is there a possibility of adding a wrapper function around a few basic
> vector operations such as VecCopy, VecAXPY, VECAXPYPZ... to operate on
> ghosted vectors? I perform a lot of vector operations including the ghost
> region to avoid communication. It really clutters code to see stuff like
>
>     CALL VecGhostGetLocalForm(solver%u,uk,ierr)
>     CALL VecGhostGetLocalForm(solver%ukm1,ukm1,ierr)
>     CALL VecGhostGetLocalForm(solver%p,pk,ierr)
>     CALL VecGhostGetLocalForm(solver%pkm1,pkm1,ierr)
>     CALL VecCopy(uk,ukm1,ierr)
>     CALL VecCopy(pk,pkm1,ierr)
>     CALL VecGhostRestoreLocalForm(solver%u,uk,ierr)
>     CALL VecGhostRestoreLocalForm(solver%ukm1,ukm1,ierr)
>     CALL VecGhostRestoreLocalForm(solver%p,pk,ierr)
>     CALL VecGhostRestoreLocalForm(solver%pkm1,pkm1,ierr)
>
> Maybe there's a good reason for this not already existing, or maybe it
> does exist and I haven't found it, but I'd think something like
> VecCopyGhosted would be fairly useful.
>

Is there something we could do that would be better than you defining a
function VecCopyGhosted()?

  Thanks,

     Matt


> John
>
-- 
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-dev/attachments/20130227/9acf1462/attachment.html>