Thank you  Jed.<br><br><div class="gmail_quote">On Sat, Dec 3, 2011 at 5:09 PM, Jed Brown <span dir="ltr">&lt;<a href="mailto:jedbrown@mcs.anl.gov">jedbrown@mcs.anl.gov</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div class="im"><br><br><div class="gmail_quote">On Sat, Dec 3, 2011 at 17:07, Hailong Xiao <span dir="ltr">&lt;<a href="mailto:xiaohl1986@gmail.com" target="_blank">xiaohl1986@gmail.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

So, I do not need to call the VecAssemblyBegin; and VecAssemblyEnd; routines in order to set the values of f?</blockquote></div><br></div><div>Only in order set values that you do not own.</div><div><br></div><div>Alternatively, you can set values in a local vector and DMDALocalToGlobalBegin/End.</div>

</blockquote></div><br><br clear="all"><br>-- <br>Hailong<br>