On Thu, Oct 25, 2012 at 6:49 AM, Thomas Witkowski <span dir="ltr"><<a href="mailto:thomas.witkowski@tu-dresden.de" target="_blank">thomas.witkowski@tu-dresden.de</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">
Following the documentation it is allowed to set the rhs and solution vectors in KSPSolve to be the same. When my KSPSolve is based on a MatShell, must the MatShell operator function be implemented to handle also the case that both vectors are the same? Or can I assume, within MatShell, that both vectors, which are provided by KSPSolve, are really two different objects?</blockquote>
<div><br></div><div>This is taken care of by KSPSolve, so the vectors are always different in the implementations. Note that MatMult is typically used with different vectors internal to the KSP. The MatMult documentation requires the vectors to be different.</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="HOEnZb"><font color="#888888"><br>
<br>
Thomas<br>
</font></span></blockquote></div><br>