<div dir="ltr"><div class="gmail_quote"><div dir="ltr">On Fri, Dec 7, 2018 at 10:05 AM Yingjie Wu via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div dir="ltr"><div>Dear Petsc developer:</div><div>Hi,</div><div> I recently begin a new program, which is a two-dimensional nonlinear PDEs
 problem, and solved by SNES.  However, there is an error. The log is :</div><div><br></div><div><div style="margin-left:40px">Linear solve did not converge due to DIVERGED_PCSETUP_FAILED iterations 0<br>                 PCSETUP_FAILED due to FACTOR_NUMERIC_ZEROPIVOT <br>Nonlinear solve did not converge due to DIVERGED_LINEAR_SOLVE iterations 0<br>SNES Object: 1 MPI processes<br>  type: newtonls<br>  maximum iterations=50, maximum function evaluations=10000<br>  tolerances: relative=1e-08, absolute=1e-50, solution=1e-08<br>  total number of linear solver iterations=0<br>  total number of function evaluations=18903<br>  norm schedule ALWAYS<br>  Jacobian is built using finite differences one column at a time<br>  SNESLineSearch Object: 1 MPI processes<br>    type: bt<br>      interpolation: cubic<br>      alpha=1.000000e-04<br>    maxstep=1.000000e+08, minlambda=1.000000e-12<br>    tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08<br>    maximum iterations=40<br>  KSP Object: 1 MPI processes<br>    type: gmres<br>      restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement<br>      happy breakdown tolerance 1e-30<br>    maximum iterations=10000, initial guess is zero<br>    tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.<br>    left preconditioning<br>    using PRECONDITIONED norm type for convergence test<br>  PC Object: 1 MPI processes<br>    type: lu<br>      out-of-place factorization<br>      tolerance for zero pivot 2.22045e-14<br>      matrix ordering: nd<br>      factor fill ratio given 5., needed 20.7926<br>        Factored matrix follows:<br>          Mat Object: 1 MPI processes<br>            type: seqaij<br>            rows=18901, cols=18901<br>            package used to perform factorization: petsc<br>            total: nonzeros=3029644, allocated nonzeros=3029644<br>            total number of mallocs used during MatSetValues calls =0<br>              not using I-node routines<br>    linear system matrix = precond matrix:<br>    Mat Object: 1 MPI processes<br>      type: seqaij<br>      rows=18901, cols=18901<br>      total: nonzeros=145708, allocated nonzeros=341390<br>      total number of mallocs used during MatSetValues calls =16459<br>        not using I-node routines</div><br></div><div>I used "-snes_fd -pc_type lu" to test the program.<br></div><div>The error shows that there are some zero elements in the matrix. What might cause this problem?</div></div></div></blockquote><div><br></div><div>Error in the residual? It looks like changes in some variable do not effect it.</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 dir="ltr"><div>Thanks,</div><div>Yingjie<br></div></div></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>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><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>