setting values in parallel vectors

[Khan, Irfan]

Excellent! Thank you

Irfan
----- Original Message -----
From: "Jed Brown" <jed at 59A2.org>
To: petsc-users at mcs.anl.gov
Sent: Thursday, April 2, 2009 3:18:06 PM GMT -05:00 US/Canada Eastern
Subject: Re: setting values in parallel vectors

On Thu 2009-04-02 15:09, Khan, Irfan wrote:
> Thank you, that was very helpful. Please do let me know if I understood this right. Generally VecSetValues()+VecGhostUpdateBegin/End() is faster than VecGetArray()+assign_array()+VecRestoreArray()+VecGhostUpdateBegin/End(). Also both these operation would be equivalent.

These operations would be equivalent.  The second option would be
faster, but only because it doesn't do a function call and
range-checking in an inner loop.  It would be nearly impossible to tell
the difference in a real code so don't worry about it, just use
whichever is more natural.

Jed