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.

John