On Tue, Aug 7, 2012 at 8:56 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 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><br>Hello, Jed</div>
<div>    I am a little puzzled about the "How can I determine the condition number of a matrix?" in the</div>
<div>FAQ. The details are as follows:<br> <br>    "...For larger matrices, you can run with <em>-pc_type none -ksp_type gmres -ksp_monitor_singular_value </em></div>
<div><em>-ksp_gmres_restart 1000</em> to get approximations to the condition number of the operator. This will </div>
<div>generally be accurate for the largest singular values, but may overestimate the smallest singular </div>
<div>value unless the method has converged. Make sure to avoid restarts...".<br>    <br>    Since the command "<em>-pc_type none -ksp_type gmres -ksp_monitor_singular_value</em> </div>
<div><em>-ksp_gmres_restart 1000</em>" using restart parameter, so, what does "Make sure to avoid restarts" mean?</div></div></blockquote><div><br></div><div>Do not believe the answer after the restart.</div>
<div> </div><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>Besides, is the parameter "<em>1000</em>" a reference value which will be changed in specific project in order</div>
<div>to be converged? Or must use <em>1000 </em>to get the condition number?</div></div></blockquote><div><br></div><div>The former.</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">
<div style="line-height:1.7;font-size:14px;font-family:arial"><div> </div>
<div>    Thanks.</div>
<div>                                                       Jim</div>
<div><br> </div>
<div></div>
<div></div>
<div><br></div>>>在 2012-07-11 22:19:51，"Jed Brown" <<a href="mailto:jedbrown@mcs.anl.gov" target="_blank">jedbrown@mcs.anl.gov</a>> 写道：<br>
<blockquote style="BORDER-LEFT:#ccc 1px solid;MARGIN:0px 0px 0px 0.8ex;PADDING-LEFT:1ex">
<div class="gmail_quote">>>On Wed, Jul 11, 2012 at 8: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 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> </div>
<div>>>Hello,</div>
<div>    </div>
<div>>>    Is there function in PETSc for judging the properties of the matrix A ? Such as positive definitiveness and </div>
<div> </div>
<div>>>conditional number. I know that there are several functions begin with MatIs+, but it has limited amounts. </div></div></blockquote>
<div><br></div>
<div>>The Ritz values are estimates of the eigenvalues that you can obtain using KSPSolve() with ->ksp_compute_eigenvalues or -ksp_plot_eigenvalues. Note that these are estimates for the _preconditioned_ >operator, so you should use -pc_type none if you want estimates for the original operator.</div>

<div><br></div>
<div><a href="http://www.mcs.anl.gov/petsc/documentation/faq.html#conditionnumber" target="_blank"><font color="#000000">></font>http://www.mcs.anl.gov/petsc/documentation/faq.html#conditionnumber</a></div></div></blockquote>
</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>