On Fri, Jun 24, 2011 at 2:38 AM, Milan Mitrovic <span dir="ltr"><<a href="mailto:milanm@student.ethz.ch">milanm@student.ethz.ch</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
So, any ideas?<br>
<br>
when I call MatMult(A,b,x,info) the ghosts get updated. When I call<br>
KSPSolve(ksp,b,x,info) the ghosts dont get updated (the multiplication<br>
routine is called... I checked that)... It seems that KSP calles<br>
MatMult with different vectors... not the original b and x...<br>
</blockquote></div><br>Yes, that definitely happens. You cannot count on the same vector being<div>used internally. They are the same type (ghosted), but they may come</div><div>from VecDuplicate/Copy(), so anytime you need ghost values, you need</div>
<div>to update.</div><div><br></div><div> Matt<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>
</div>