On Wed, Jul 11, 2012 at 12:40 PM, Bao Kai <span dir="ltr"><<a href="mailto:paeanball@gmail.com" target="_blank">paeanball@gmail.com</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi, all,<br>
<br>
The following is the output from the solution of a Poisson equation<br>
from Darcy's law.<br>
<br>
To compute the condition number of matrix, I did not use PC and use<br>
GMRES KSP to do the test.<br>
<br>
It seems like that the condition number keep increasing during the<br>
iterative solution. Does this mean the matrix is ill-conditioned?<br></blockquote><div><br></div><div>Generally yes. Krylov methods take a long time to resolve the smallest</div><div>eigenvalues, so this approximation is not great.</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
For this test, it did not achieve convergence with 10000 iterations.<br>
<br>
When I use BJOCABI PC and BICGSTAB KSP, it generally takes about 600<br>
times iteration to get the iteration convergent.<br>
<br>
Any suggestion for improving the convergence rate will be much<br>
appreciated.  The solution of this equation has been the bottleneck of<br>
my code, it takes more than 90% of the total time.<br></blockquote><div><br></div><div>Try ML or GAMG.</div><div><br></div><div>    Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

Thank you very much.<br>
<br>
Best Regards,<br>
Kai<br>
</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>