<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Tue, Sep 20, 2016 at 6:44 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>
to compute / approximate the condition number of a MATSBAIJ, I put<br>
<br>
-ksp_view<br>
<br>
# Conditon number estimate<br>
-pc_type none<br>
-ksp_type gmres<br>
-ksp_monitor_singular_value<br>
-ksp_gmres_restart 1000<br>
<br>
in my .petscrc<br>
<br>
Output is like:<br>
<br>
[...]<br>
566 KSP Residual norm 1.241765807317e-07 % max 7.020130499234e+02 min 6.429054752025e-04 max/min 1.091938203983e+06<br>
567 KSP Residual norm 1.219340847328e-07 % max 7.020130499278e+02 min 6.423976501709e-04 max/min 1.092801397609e+06<br>
568 KSP Residual norm 1.198886059519e-07 % max 7.020130499320e+02 min 6.419283878172e-04 max/min 1.093600257062e+06<br>
569 KSP Residual norm 1.178377018879e-07 % max 7.020130499362e+02 min 6.414517591235e-04 max/min 1.094412853268e+06<br>
KSP Object: 8 MPI processes<br>
  type: gmres<br>
    GMRES: restart=1000, 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: 8 MPI processes<br>
  type: none<br>
  linear system matrix = precond matrix:<br>
  Mat Object:  C   8 MPI processes<br>
    type: mpisbaij<br>
    rows=14404, cols=14404<br>
    total: nonzeros=1059188, allocated nonzeros=1244312<br>
    total number of mallocs used during MatSetValues calls =75272<br>
        block size is 1<br>
(0) 13:30:20 [precice::impl::<wbr>SolverInterfaceImpl]:380 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:30:20 [precice::mapping::<wbr>PetRadialBasisFctMapping]:519 in map: Mapping Data consistent from MeshA (ID 0) to MeshB<br>
(ID 1)<br>
  0 KSP Residual norm 1.178378120697e-07 % max 1.000000000000e+00 min 1.000000000000e+00 max/min 1.000000000000e+00<br>
KSP Object: 8 MPI processes<br>
  type: gmres<br>
    GMRES: restart=1000, 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: 8 MPI processes<br>
  type: none<br>
  linear system matrix = precond matrix:<br>
  Mat Object:  C   8 MPI processes<br>
    type: mpisbaij<br>
    rows=14404, cols=14404<br>
    total: nonzeros=1059188, allocated nonzeros=1244312<br>
    total number of mallocs used during MatSetValues calls =75272<br>
        block size is 1<br>
<br>
<br>
The approximate condition number is the max/min value, 1.094412853268e+06 here?<br></blockquote><div><br></div><div>Yes.</div><div><br></div><div>  Thanks,</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">
Just make sure my mathematical illerateracy does not spoil my report...<br>
<br>
<br>
Best thanks,<br>
Florian<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>