<div dir="ltr">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<div>
<br></div><div><div>    CALL VecGhostGetLocalForm(solver%u,uk,ierr)</div><div>    CALL VecGhostGetLocalForm(solver%ukm1,ukm1,ierr)</div><div>    CALL VecGhostGetLocalForm(solver%p,pk,ierr)</div><div>    CALL VecGhostGetLocalForm(solver%pkm1,pkm1,ierr)</div>
<div>    CALL VecCopy(uk,ukm1,ierr)</div><div>    CALL VecCopy(pk,pkm1,ierr)</div><div>    CALL VecGhostRestoreLocalForm(solver%u,uk,ierr)</div><div>    CALL VecGhostRestoreLocalForm(solver%ukm1,ukm1,ierr)</div><div>    CALL VecGhostRestoreLocalForm(solver%p,pk,ierr)</div>
<div>    CALL VecGhostRestoreLocalForm(solver%pkm1,pkm1,ierr)</div></div><div><br></div><div style>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.</div>
<div style><br></div><div style>John</div></div>