On Thu, Apr 30, 2009 at 7:41 AM, Jed Brown <span dir="ltr"><<a href="mailto:jed@59a2.org">jed@59a2.org</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<div class="im">Matthew Knepley wrote:<br>
<br>
> We could make a copy and do the same thing for GetArray. We don't because<br>
> it is more expensive. The same argument applies here. Some people will<br>
> complain about any unnecessary expense.<br>
<br>
</div>It *is* expensive with GetArray, so that would be pointless.  In this<br>
case, it does *not* cost extra, you just set the state to the maximum.<br><div class="im"></div></blockquote><div><br>I misunderstood. I thought you were checking for element maxima, which would<br>tell you if something changed. You are just checking state numbers? I did not think<br>
these had any guarantees of synchronization.<br> </div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div class="im">
> The local form is meaningless after restore so its state is unimportant.<br>
<br>
</div>VecGhostGetLocalForm(g,&lf);<br>
VecNorm(lf,NORM_1,&nrm);<br>
VecGhostRestoreLocalForm(g,&lf);<br>
VecSet(g,0);<br>
VecGhostGetLocalForm(g,&lf);<br>
VecNorm(lf,NORM_1,&nrm); /* will erroneously use cached norm */<br>
<br>
You might say that the user should not do this (the ghost values are<br>
stale, and the norm of a local vector has dubious value), but it doesn't<br>
cost anything to at least do it in a consistent way.  This is why I<br>
prefer to syncronize state on both GetLocalForm and RestoreLocalForm.</blockquote><div><br>You are not guaranteed to get the same local vector.<br><br>   Matt<br> </div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<font color="#888888"><br>
Jed<br></font></blockquote></div><br>-- <br>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<br>