On Tue, Oct 23, 2012 at 6:17 PM, Jinquan Zhong <span dir="ltr"><<a href="mailto:jzhong@scsolutions.com" target="_blank">jzhong@scsolutions.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">
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">Any way to get the condition number.  I used
<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><code><span style="font-size:10.0pt">-ksp_type gmres -ksp_monitor_singular_value -ksp_gmres_restart 1000<u></u><u></u></span></code></p>
<p class="MsoNormal"><code><span style="font-size:10.0pt"><u></u> <u></u></span></code></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">It didn’t work.</span></p></blockquote><div><br></div><div>First, as a matter of course, remove "It didn't work" from your vocabulary. Please never utter those words on any mailing list ever again. It wastes everyone's time.</div>
</div><br><div>It is estimating the singular values of the <b>preconditioned</b> operator. Use -pc_type none to remove the preconditioner. Note that this is generally extremely inaccurate for the smallest eigenvalue, especially for an ill-conditioned problem and when you don't run to convergence. It can over-estimate the smallest singular value by many orders of magnitude. In general, it is difficult to estimate the smallest singular value without a <b>solver</b>. You can use SLEPc to accurately and efficiently compute the smallest singular values, provided you have a reasonable preconditioner (or direct solver) for the system.</div>