<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Tue, Aug 30, 2016 at 7:01 AM, Florian Lindner <span dir="ltr"><<a href="mailto:mailinglists@xgm.de" target="_blank">mailinglists@xgm.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello,<br>
<br>
there is a FAQ and a Stackoverflow article about getting the condition number of a petsc matrix:<br>
<br>
<a href="http://www.mcs.anl.gov/petsc/documentation/faq.html#conditionnumber" rel="noreferrer" target="_blank">http://www.mcs.anl.gov/petsc/<wbr>documentation/faq.html#<wbr>conditionnumber</a><br>
<a href="http://scicomp.stackexchange.com/questions/34/how-can-i-estimate-the-condition-number-of-a-large-sparse-matrix-using-petsc" rel="noreferrer" target="_blank">http://scicomp.stackexchange.<wbr>com/questions/34/how-can-i-<wbr>estimate-the-condition-number-<wbr>of-a-large-sparse-matrix-<wbr>using-petsc</a><br>
<br>
Both tell me to add:<br>
<br>
-pc_type none -ksp_type gmres -ksp_monitor_singular_value -ksp_gmres_restart 1000<br>
<br>
to my options.<br>
<br>
I add the line to .petscrc but nothing happens, no additional output at all. I added -ksp_view, so my .petscrc looks<br>
like that:<br></blockquote><div><br></div><div>Each option must be on its own line .petscrc</div><div><br></div><div>  Thanks,</div><div><br></div><div>    Matt</div><div><br></div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
-ksp_view<br>
-pc_type none -ksp_type gmres -ksp_monitor_singular_value -ksp_gmres_restart 1000<br>
<br>
The complete output is below, but something I wonder about:<br>
<br>
GMRES: restart=30, shouldn't that be 1000<br>
<br>
And where can I read out the condition number approximation?<br>
<br>
Thanks,<br>
Florian<br>
<br>
<br>
KSP Object: 1 MPI processes<br>
  type: gmres<br>
    GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement<br>
    GMRES: happy breakdown tolerance 1e-30<br>
  maximum iterations=10000<br>
  tolerances:  relative=1e-09, absolute=1e-50, divergence=10000.<br>
  left preconditioning<br>
  using nonzero initial guess<br>
  using PRECONDITIONED norm type for convergence test<br>
PC Object: 1 MPI processes<br>
  type: none<br>
  linear system matrix = precond matrix:<br>
  Mat Object:  C   1 MPI processes<br>
    type: seqsbaij<br>
    rows=14403, cols=14403<br>
    total: nonzeros=1044787, allocated nonzeros=1123449<br>
    total number of mallocs used during MatSetValues calls =72016<br>
        block size is 1<br>
(0) 13:58:35 [precice::impl::<wbr>SolverInterfaceImpl]:395 in initialize: it 1 of 1 | dt# 1 | t 0 of 1 | dt 1 | max dt 1 |<br>
ongoing yes | dt complete no |<br>
(0) 13:58:35 [precice::impl::<wbr>SolverInterfaceImpl]:446 in advance: Iteration #1<br>
KSP Object: 1 MPI processes<br>
  type: gmres<br>
    GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement<br>
    GMRES: happy breakdown tolerance 1e-30<br>
  maximum iterations=10000<br>
  tolerances:  relative=1e-09, absolute=1e-50, divergence=10000.<br>
  left preconditioning<br>
  using nonzero initial guess<br>
  using PRECONDITIONED norm type for convergence test<br>
PC Object: 1 MPI processes<br>
  type: none<br>
  linear system matrix = precond matrix:<br>
  Mat Object:  C   1 MPI processes<br>
    type: seqsbaij<br>
    rows=14403, cols=14403<br>
    total: nonzeros=1044787, allocated nonzeros=1123449<br>
    total number of mallocs used during MatSetValues calls =72016<br>
        block size is 1<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature">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</div>
</div></div>