<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Mon, Feb 22, 2016 at 6:02 PM, David Knezevic <span dir="ltr"><<a href="mailto:david.knezevic@akselos.com" target="_blank">david.knezevic@akselos.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Mon, Feb 22, 2016 at 6:57 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><br>
  Add -ksp_converged_reason<br></blockquote><div><br></div><div><br></div><div>That gives "DIVERGED_INDEFINITE_PC". The full output is below:</div></div></div></div></blockquote><div><br></div><div>Just for testing, replace CG with GMRES.</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"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><div><br></div><div><div>  105 KSP Residual norm 7.099534889984e-05 </div><div>  106 KSP Residual norm 6.615351528474e-05 </div><div>  107 KSP Residual norm 6.303646443688e-05 </div><div>  Linear solve converged due to CONVERGED_RTOL iterations 107</div><div>  NL step  5, |residual|_2 = 1.996487e+01</div><div>    0 KSP Residual norm 5.844456247091e+00 </div><div>    1 KSP Residual norm 1.308756943973e+00 </div><div>    2 KSP Residual norm 7.545133187463e-01 </div><div>    3 KSP Residual norm 1.346068942607e+00 </div><div>  Linear solve did not converge due to DIVERGED_INDEFINITE_PC iterations 4</div><div>SNES Object: 1 MPI processes</div><div>  type: newtonls</div><div>  maximum iterations=50, maximum function evaluations=10000</div><div>  tolerances: relative=1e-08, absolute=1e-50, solution=1e-08</div><div>  total number of linear solver iterations=535</div><div>  total number of function evaluations=11</div><div>  norm schedule ALWAYS</div><div>  SNESLineSearch Object:   1 MPI processes</div><div>    type: bt</div><div>      interpolation: cubic</div><div>      alpha=1.000000e-04</div><div>    maxstep=1.000000e+08, minlambda=1.000000e-12</div><div>    tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08</div><div>    maximum iterations=40</div><div>  KSP Object:   1 MPI processes</div><div>    type: cg</div><div>    maximum iterations=10000, initial guess is zero</div><div>    tolerances:  relative=1e-05, absolute=1e-50, divergence=10000</div><div>    left preconditioning</div><div>    using NATURAL norm type for convergence test</div><div>  PC Object:   1 MPI processes</div><div>    type: ilu</div><div>      ILU: out-of-place factorization</div><div>      0 levels of fill</div><div>      tolerance for zero pivot 2.22045e-14</div><div>      matrix ordering: natural</div><div>      factor fill ratio given 1, needed 1</div><div>        Factored matrix follows:</div><div>          Mat Object:           1 MPI processes</div><div>            type: seqaij</div><div>            rows=967608, cols=967608</div><div>            package used to perform factorization: petsc</div><div>            total: nonzeros=7.55539e+07, allocated nonzeros=7.55539e+07</div><div>            total number of mallocs used during MatSetValues calls =0</div><div>              using I-node routines: found 322536 nodes, limit used is 5</div><div>    linear system matrix = precond matrix:</div><div>    Mat Object:    ()     1 MPI processes</div><div>      type: seqaij</div><div>      rows=967608, cols=967608</div><div>      total: nonzeros=7.55539e+07, allocated nonzeros=7.55539e+07</div><div>      total number of mallocs used during MatSetValues calls =0</div><div>        using I-node routines: found 322536 nodes, limit used is 5</div><div>Number of nonlinear iterations: 5</div><div>Nonlinear solver convergence/divergence reason: DIVERGED_LINEAR_SOLVE</div></div><div><br></div><div><br></div><div><br></div></div></div></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="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>