<div dir="ltr"><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Jan 28, 2013 at 1:36 PM, Ling Zou <span dir="ltr"><<a href="mailto:lingzou80@gmail.com" target="_blank">lingzou80@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div id=":15m">I guess I got the impression that M is generally the same as A from reading the manual.(PETSc Users Manual, Reversion 3.3, page 71, under 4.1 Using KSP)<br>
<br>"<span style="background-color:rgb(255,255,153)"><span style="color:rgb(153,0,0)">Typically the preconditioning matrix (i.e., the matrix from which the preconditioner is to be constructed), Pmat, is the same as the matrix that defines the linear system, Amat</span></span>; however, occasionally these matrices differ (for instance, when a preconditioning matrix is obtained from a lower order method than that employed to form the linear system matrix)."</div>
</blockquote></div><br>I think you are getting confused by mixed notation. A PC in PETSc is an algorithm that takes a matrix (Pmat in the docs) and does some work to be able to apply an operation (named "M^{-1}" in your first email). This does *not* imply that M=Pmat, or that M is ever available or used, there is just a a linear operation named "M^{-1}".</div>
</div>