[petsc-dev] Ghost update in VecAssembly?

Jed Brown jed at jedbrown.org
Tue Mar 1 08:23:10 CST 2016


Matthew Knepley <knepley at gmail.com> writes:
>> ghost values.  Also, who uses VecGhost _and_ VecAssembly for the same
>> Vec?  That's quite wasteful; if you have local indexing, why not just
>> use VecGhostUpdate(...,SCATTER_REVERSE)?
>>
>
> I do not understand this argument. If you form the local part of each vector
> using VecGetArray(), then you must call both VecAssembly() 

No.  You only need to call VecAssembly if you use VecSetValues.  But if
you have a ghosted vector, you would normally use VecGetArray and never
bother with VecSetValues.

> and VecGhostUpdate(), which I think is the complaint.

-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 818 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20160301/e49afa36/attachment.sig>


More information about the petsc-dev mailing list