<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Tue, Mar 1, 2016 at 8:56 AM, Jed Brown <span dir="ltr"><<a href="mailto:jed@jedbrown.org" target="_blank">jed@jedbrown.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><span class="">Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> writes:<br>
>> No.  You only need to call VecAssembly if you use VecSetValues.  But if<br>
>><br>
><br>
> I thought we had state increase on getarray() as well, but its just locking.<br>
<br>
</span>VecRestoreArray calls PetscObjectStateIncrease, but what does that has<br>
to do with VecAssembly?<br>
<span class=""><br>
>> you have a ghosted vector, you would normally use VecGetArray and never<br>
>> bother with VecSetValues.<br>
><br>
><br>
> Why would we assume people make the ghost values as well?<br>
<br>
</span>Huh?<br>
<div class=""><div class="h5"><br>
> However, if this is the only case, I don't think putting it in is a<br>
> great idea.<br>
</div></div></blockquote></div><br>I am trying to understand what you are saying here. I have a vector on 2 procs, divided as</div><div class="gmail_extra"><br></div><div class="gmail_extra">p0  |----interior---|--ghost--|<br clear="all"><div><br></div><div>p1  |----interior---|--ghost--|<br clear="all"></div><div><br></div><div>Now I calculate some values in the interior on both procs using VecGetArray(),</div><div class="gmail_extra"><br></div><div class="gmail_extra">p0  |xxxxinteriorxxx|--ghost--|<br clear="all"><div><br></div><div>p1  |xxxxinteriorxxx|--ghost--|<br clear="all"></div><div><br></div><div>What is wrong with calling VecGhostUpdate() here to fill in the ghosts? Is your contention</div><div>that I should have also calculated the ghosts?</div><div><br></div><div>  Matt</div><div><br></div></div>-- <br><div class="gmail_signature">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div>
</div></div>