<br><br><div class="gmail_quote">On Mon, Jan 28, 2013 at 12:43 PM, Jed Brown <span dir="ltr"><<a href="mailto:jedbrown@mcs.anl.gov" target="_blank">jedbrown@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr"><div class="gmail_extra"><div class="im"><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>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></div>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></blockquote><div><br>Hmmm...it's getting more complicated now. I guess I need read the manual more carefully and study the example codes more. By the way, could you redirect me to an example code or a tutorial to understand this better? <br>
<br>Ling<br></div></div><br>