<br><br>
<div class="gmail_quote">On Tue, Sep 15, 2009 at 8:29 PM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">
<div class="im">On Tue, Sep 15, 2009 at 6:04 PM, Ryan Yan <span dir="ltr"><<a href="mailto:vyan2000@gmail.com" target="_blank">vyan2000@gmail.com</a>></span> wrote:<br></div>
<div class="gmail_quote">
<div class="im">
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0pt 0pt 0pt 0.8ex; BORDER-LEFT: rgb(204,204,204) 1px solid">
<div>Thanks Jed. </div>
<div> </div>
<div>And in the term ||Ae||/||Ax||, PETSc impilictly assume that we use zero as initial guess and here the x really means exact solution, right?</div></blockquote></div>
<div><br>This has nothing to do with the initial guess. Ae is the residual, and x is the current guess. We do not have the exact solution,<br>but I guess A x^* = b, so that could have been used.<br><br></div></div></blockquote>
<div> </div>
<div> </div>
<div>Thanks Matt.</div>
<div> </div>
<div>I agree that in the print out ||Ae||/||Ax||, x here just represents the x*.</div>
<div> </div>
<div>Yan</div>
<div> </div>
<div> </div>
<div> </div>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">
<div class="gmail_quote">
<div> Matt<br> </div>
<div class="im">
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0pt 0pt 0pt 0.8ex; BORDER-LEFT: rgb(204,204,204) 1px solid">
<div><br></div><font color="#888888">
<div>Yan</div></font>
<div>
<div></div>
<div>
<div> </div>
<div><br><br> </div>
<div class="gmail_quote">On Tue, Sep 15, 2009 at 6:41 PM, Jed Brown <span dir="ltr"><<a href="mailto:jed@59a2.org" target="_blank">jed@59a2.org</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: rgb(204,204,204) 1px solid">
<div>Ryan Yan wrote:<br>> But can you say a little more about what happened in the first scenario, I<br>> did not quite fellow you. Let's say, If I set -ksp_rtol 1e-2, what does this<br>> mean in the left PC case. More specifically, which term will be considered<br>
> as a stop creteria for the KSP solve.<br><br></div>Look at the first column (preconditioned residual norm). Notice that<br>this is decreasing by 1e-2 and 1e-3 respectively in your examples.<br><font color="#888888"><br>
Jed<br><br></font></blockquote></div><br></div></div></blockquote></div></div><font color="#888888"><br><br clear="all"><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></font></blockquote></div><br>