Possibilities (not mutually exclusive)<div><br></div><div>1. the convergence tolerance is too tight to converge in double precision arithmetic, so the solver is stagnating (you could try --with-precision=__float128)</div><div>
<br></div><div>2. the tolerance is too tight to converge in 10000 iterations (either increase the number of iterations or use a better preconditioner</div><div><br></div><div>3. the tolerance is tighter than necessary to get a solution that is "right enough" for you (try loosening the tolerance)<br>
<br><div class="gmail_quote">On Tue, Jul 10, 2012 at 9:51 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><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="line-height:1.7;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>   </div>
<div>    Is it right?<br><br></div><div><div class="h5">
<div></div>
<div></div>
<div><br></div>>At 2012-07-08 00:28:54,"Mark F. Adams" <<a href="mailto:mark.adams@columbia.edu" target="_blank">mark.adams@columbia.edu</a>> wrote:<br>
<blockquote style="BORDER-LEFT:#ccc 1px solid;MARGIN:0px 0px 0px 0.8ex;PADDING-LEFT:1ex">>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>
<blockquote type="cite">
<div style="line-height:1.7;font-size:14px;font-family:arial">
<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" target="_blank">knepley@gmail.com</a>> 写道:<br>
<blockquote style="BORDER-LEFT:#ccc 1px solid;MARGIN:0px 0px 0px 0.8ex;PADDING-LEFT:1ex">>>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></span></span></blockquote></div><br></div></blockquote></div></div></div><br><br><span title="neteasefooter"><span></span></span></blockquote></div><br></div>