No.<br><br><div class="gmail_quote">On Mon, Dec 14, 2009 at 9:40 AM, Jarunan Panyasantisuk <span dir="ltr">&lt;<a href="mailto:jarunan@ascomp.ch">jarunan@ascomp.ch</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Thank you. I see it.<br>
I also see in the output file that sometimes the initial residual<br>
||r_0|| = ||b|| is simply rnorm ||r|| �at iteration 0. Does it change<br>
the reference during the computation?<br>
<br>
Jarunan<br>
<div class="im"><br>
<br>
Matthew Knepley wrote:<br>
&gt; You can see in the info output what is happening. The test is<br>
&gt;<br>
&gt; � rnorm &gt;= ksp-&gt;divtol*ksp-&gt;rnorm0<br>
&gt;<br>
&gt; Here ||r_0|| = ||b|| = 0.00331536 and ||r|| = 513.979, and thus it results<br>
&gt; in divergence. It appears that ILU(0) is very bad for this problem.<br>
&gt;<br>
&gt; � �Matt<br>
&gt;<br>
&gt; On Mon, Dec 14, 2009 at 9:07 AM, &lt;<a href="mailto:jarunan@ascomp.ch">jarunan@ascomp.ch</a><br>
</div><div class="im">&gt; &lt;mailto:<a href="mailto:jarunan@ascomp.ch">jarunan@ascomp.ch</a>&gt;&gt; wrote:<br>
&gt;<br>
&gt; � � Hi,<br>
&gt;<br>
&gt; � � Thank you for your advice. With this mail I attach the output file<br>
&gt; � � run.log.<br>
&gt; � � My grid is 2-level grid. The coarse grid is computed first then<br>
&gt; � � the fine grid. The problem occur in the fine grid. You can see at<br>
&gt; � � lines 3098, 4552 and 6127 where it is diverged.<br>
&gt;<br>
&gt; � � Regards,<br>
&gt; � � Jarunan<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; � � Quoting Matthew Knepley &lt;<a href="mailto:knepley@gmail.com">knepley@gmail.com</a><br>
</div>&gt; � � &lt;mailto:<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>&gt;&gt;:<br>
<div class="im">&gt;<br>
&gt; � � � � Run with -ksp_view -ksp_monitor -info and send us all the output.<br>
&gt;<br>
&gt; � � � � �Matt<br>
&gt;<br>
&gt; � � � � On Mon, Dec 14, 2009 at 8:36 AM, &lt;<a href="mailto:jarunan@ascomp.ch">jarunan@ascomp.ch</a><br>
</div><div><div></div><div class="h5">&gt; � � � � &lt;mailto:<a href="mailto:jarunan@ascomp.ch">jarunan@ascomp.ch</a>&gt;&gt; wrote:<br>
&gt;<br>
&gt;<br>
&gt; � � � � � � Dear Petsc Team,<br>
&gt;<br>
&gt; � � � � � � There is something not clear for me about divergence<br>
&gt; � � � � � � tolerance. Firstly, my<br>
&gt; � � � � � � computation diverges with the reason of dtol (-4), which<br>
&gt; � � � � � � was set by default<br>
&gt; � � � � � � (10000). When I increased dtol to 1.0e+8, it became<br>
&gt; � � � � � � converged. However, the<br>
&gt; � � � � � � relative residual where it was diverged is less than the<br>
&gt; � � � � � � default dtol. I<br>
&gt; � � � � � � would like to know, how it is exactly calculated.<br>
&gt;<br>
&gt; � � � � � � ----------I describe my problem more in detail<br>
&gt; � � � � � � here.-------------------<br>
&gt; � � � � � � As I set the tolerances for the solver as below, my<br>
&gt; � � � � � � problem diverges with<br>
&gt; � � � � � � converged reason -4 (from divergence tolerance)<br>
&gt;<br>
&gt; � � � � � � call KSPSetTolerances(ksp,1.0e-06,1.0e-20,<br>
&gt; � � � � � � PETSC_DEFAULT_DOUBLE_PRECISION,10000,ierr)<br>
&gt;<br>
&gt; � � � � � � I monitor the residual. it printed the information below<br>
&gt; � � � � � � and quit the<br>
&gt; � � � � � � computation.(I delete the true residual norm and Ax/Ae.<br>
&gt; � � � � � � Please tell me if<br>
&gt; � � � � � � you need them.)<br>
&gt;<br>
&gt; � � � � � � �0 KSP preconditioned resid norm 5.139790011051e+02<br>
&gt;<br>
&gt; � � � � � � With -ksp_view, it shows that default value of dtol is<br>
&gt; � � � � � � 10000. I increased<br>
&gt; � � � � � � dtol to 1.0e+8, then the solver solved the problem. With<br>
&gt; � � � � � � -ksp_monitor_true_residual, the first 5 iterations read:<br>
&gt;<br>
&gt; � � � � � � �0 KSP preconditioned resid norm 5.139790011051e+02<br>
&gt; � � � � � � �1 KSP preconditioned resid norm 1.840751662590e+01<br>
&gt; � � � � � � �2 KSP preconditioned resid norm 3.653339750400e+00<br>
&gt; � � � � � � �3 KSP preconditioned resid norm 2.538794826933e+00<br>
&gt; � � � � � � �4 KSP preconditioned resid norm 2.299844715275e-01<br>
&gt;<br>
&gt; � � � � � � Here, from iteration 0 to iteration 1,the relative<br>
&gt; � � � � � � residual is 0.03580,<br>
&gt; � � � � � � much less than 1.0e+4 the default dtol value. So, why was<br>
&gt; � � � � � � it diverged at the<br>
&gt; � � � � � � first time?<br>
&gt;<br>
&gt; � � � � � � ----------------------------------------------------------------------------<br>
&gt;<br>
&gt; � � � � � � Regards,<br>
&gt; � � � � � � Jarunan<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; � � � � � � --<br>
&gt; � � � � � � Jarunan Panyasantisuk<br>
&gt; � � � � � � Development Engineer<br>
&gt; � � � � � � ASCOMP GmbH, Technoparkstr. 1<br>
&gt; � � � � � � CH-8005 Zurich, Switzerland<br>
&gt; � � � � � � Phone : +41 44 445 4072<br>
&gt; � � � � � � Fax � : +41 44 445 4075<br>
</div></div>&gt; � � � � � � E-mail: <a href="mailto:jarunan@ascomp.ch">jarunan@ascomp.ch</a> &lt;mailto:<a href="mailto:jarunan@ascomp.ch">jarunan@ascomp.ch</a>&gt;<br>
&gt; � � � � � � <a href="http://www.ascomp.ch" target="_blank">www.ascomp.ch</a> &lt;<a href="http://www.ascomp.ch" target="_blank">http://www.ascomp.ch</a>&gt;<br>
<div class="im">&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; � � � � --<br>
&gt; � � � � What most experimenters take for granted before they begin<br>
&gt; � � � � their experiments<br>
&gt; � � � � is infinitely more interesting than any results to which their<br>
&gt; � � � � experiments<br>
&gt; � � � � lead.<br>
&gt; � � � � -- Norbert Wiener<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; � � --<br>
&gt; � � Jarunan Panyasantisuk<br>
&gt; � � Development Engineer<br>
&gt; � � ASCOMP GmbH, Technoparkstr. 1<br>
&gt; � � CH-8005 Zurich, Switzerland<br>
&gt; � � Phone : +41 44 445 4072<br>
&gt; � � Fax � : +41 44 445 4075<br>
</div>&gt; � � E-mail: <a href="mailto:jarunan@ascomp.ch">jarunan@ascomp.ch</a> &lt;mailto:<a href="mailto:jarunan@ascomp.ch">jarunan@ascomp.ch</a>&gt;<br>
&gt; � � <a href="http://www.ascomp.ch" target="_blank">www.ascomp.ch</a> &lt;<a href="http://www.ascomp.ch" target="_blank">http://www.ascomp.ch</a>&gt;<br>
<div><div></div><div class="h5">&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; --<br>
&gt; What most experimenters take for granted before they begin their<br>
&gt; experiments is infinitely more interesting than any results to which<br>
&gt; their experiments lead.<br>
&gt; -- Norbert Wiener<br>
<br>
<br>
--<br>
Jarunan Panyasantisuk<br>
Development Engineer<br>
ASCOMP GmbH, Technoparkstr. 1<br>
CH-8005 Zurich, Switzerland<br>
Phone : +41 44 445 4072<br>
Fax � : +41 44 445 4075<br>
E-mail: <a href="mailto:jarunan@ascomp.ch">jarunan@ascomp.ch</a><br>
<a href="http://www.ascomp.ch" target="_blank">www.ascomp.ch</a><br>
<br>
<br>
</div></div></blockquote></div><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>