<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><br><div><div>On Jul 10, 2012, at 10:51 AM, w_ang_temp wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div style="line-height:1.7;color:#000000;font-size:14px;font-family:arial"><div> </div>
<div> In my opinion, convergence in PETSc is decided by rtol, atol and dtol. The divergent hints just show that</div>
<div> </div>
<div>in the solving process it does not satisfy the rule. The "right" result may be different from the true result at</div>
<div> </div>
<div>the several back decimal places(I mean that they may be the same with four decimal places but may be not</div>
<div> </div>
<div>the same with more decimal places). </div></div></blockquote><div><br></div><div>No, If you have an rtol=1.e-30 and atol=1.e-300 then your solution will probably be correct to all decimal places even though PETSc will say you "diverged". You can't be too sure about how many digits of accuracy you have with simple linear algebra arguments unless your system is extremely well conditioned.</div><div><br></div><div>As Jed said if you are doing 10000 iterations are happy with the solution then your tolerances are probably too tight. Look at the residual history and see if you are stagnating.</div><br><blockquote type="cite"><div style="line-height:1.7;color:#000000;font-size:14px;font-family:arial">
<div> </div>
<div> Is it right?<br><br></div>
<div></div>
<div id="divNeteaseMailCard"></div>
<div><br></div>>At 2012-07-08 00:28:54,"Mark F. Adams" <<a href="mailto:mark.adams@columbia.edu">mark.adams@columbia.edu</a>> wrote:<br>
<blockquote style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" id="isReplyContent">>It sounds like your -ksp_rtol is too small. Experiment with looser tolerances until your solution is not "correct" to see >how much accuracy you want.
<div><br>
<div>
<div>>On Jul 7, 2012, at 12:15 PM, w_ang_temp wrote:</div><br class="Apple-interchange-newline">
<blockquote type="cite">
<div style="LINE-HEIGHT: 1.7; FONT-FAMILY: arial; COLOR: #000000; FONT-SIZE: 14px">
<div>> Maybe it is a problem of mathematical concept. I compare the result with the true result which is </div>
<div> </div>
<div>>computed and validated by other tools. I think it is right if I get the same result.</div>
<div></div>
<div></div>
<div><br></div>>>在 2012-07-08 00:03:21，"Matthew Knepley" <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> 写道：<br>
<blockquote style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" id="isReplyContent">>>On Sat, Jul 7, 2012 at 10:00 AM, w_ang_temp <span dir="ltr"><<a href="mailto:w_ang_temp@163.com" target="_blank">w_ang_temp@163.com</a>></span> wrote:<br>
<div class="gmail_quote">
<blockquote style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" class="gmail_quote">
<div style="LINE-HEIGHT: 1.7; FONT-FAMILY: arial; FONT-SIZE: 14px">
<div>>>Hello,</div>
<div> </div>
<div>>> I am a little puzzled that I get the right result while the converged reason says that 'Linear solve >>did not</div>
<div> </div>
<div>>>converge due to DIVERGED_ITS iterations 10000'. This infomation means that the iterations >reach >the maximum</div>
<div> </div>
<div>>>iterations. But the result is right now. So why says 'did not converge'? Can I think that the result is >>right and </div>
<div> </div>
<div>>>can be used?<br></div></div></blockquote>
<div>>>Obviously, your definition of "right" is not the same as the convergence tolerances you are using.</div>
<div><br></div>
<div>>> Matt</div>
<div> </div>
<blockquote style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" class="gmail_quote">
<div style="LINE-HEIGHT: 1.7; FONT-FAMILY: arial; FONT-SIZE: 14px">
<div><span style="LINE-HEIGHT: 1.7"> </span></div>
<div>>> Thanks.</div>
<div> </div>
<div>>> Jim</div></div><br><br><span title="neteasefooter"><span></span></span></blockquote></div><br><br clear="all">
<div><br></div>-- <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></blockquote></div><br><br><span title="neteasefooter"><span id="netease_mail_footer"></span></span></blockquote></div><br></div></blockquote></div><br><br><span title="neteasefooter"><span id="netease_mail_footer"></span></span></blockquote></div><br></body></html>