<div style="line-height:1.7;color:#000000;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>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>    Thanks.</DIV>
<DIV>                                                       Jim</DIV>
<DIV><BR> </DIV>
<DIV></DIV>
<DIV id="divNeteaseMailCard"></DIV>
<DIV><BR></DIV>>>ÔÚ 2012-07-11 22:19:51£¬"Jed Brown" <jedbrown@mcs.anl.gov> Ð´µÀ£º<BR>
<BLOCKQUOTE style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" id="isReplyContent">
<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"><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 id="netease_mail_footer"></span></span>