<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">I never get to see the output. When it runs, I get:<div class=""><br class=""></div><div class=""><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class=""> 0 SNES Function norm 8.405398322914e-02 </span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[0]PETSC ERROR: Zero pivot in LU factorization: <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html#zeropivot" class="">http://www.mcs.anl.gov/petsc/documentation/faq.html#zeropivot</a></span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[0]PETSC ERROR: Zero pivot on row 0</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[0]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html" class="">http://www.mcs.anl.gov/petsc/documentation/faq.html</a> for trouble shooting.</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[0]PETSC ERROR: Petsc Release Version 3.6.1, unknown </span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[0]PETSC ERROR: ./blowup_refine on a arch-darwin-c-debug named gs_air by gideon Sat Sep 19 20:11:44 2015</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[0]PETSC ERROR: Configure options --download-mpich=yes --download-suitesparse=yes --download-superlu=yes --download-superlu_dist=yes --download-mumps=yes --download-sprng=yes --with-cxx=clang++ --with-cc=clang --with-fc=gfortran --download-metis=yes --download-parmetis=yes --download-scalapack=yes</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[0]PETSC ERROR: #1 MatSOR_SeqAIJ_Inode() line 2795 in /opt/petsc/src/mat/impls/aij/seq/inode.c</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[0]PETSC ERROR: #2 MatSOR_MPIAIJ() line 1496 in /opt/petsc/src/mat/impls/aij/mpi/mpiaij.c</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[0]PETSC ERROR: #3 MatSOR() line 3697 in /opt/petsc/src/mat/interface/matrix.c</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[0]PETSC ERROR: #4 PCApply_SOR() line 37 in /opt/petsc/src/ksp/pc/impls/sor/sor.c</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[0]PETSC ERROR: #5 PCApply() line 483 in /opt/petsc/src/ksp/pc/interface/precon.c</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[0]PETSC ERROR: #6 KSP_PCApply() line 242 in /opt/petsc/include/petsc/private/kspimpl.h</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[0]PETSC ERROR: #7 KSPInitialResidual() line 63 in /opt/petsc/src/ksp/ksp/interface/itres.c</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[0]PETSC ERROR: #8 KSPSolve_GMRES() line 235 in /opt/petsc/src/ksp/ksp/impls/gmres/gmres.c</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[0]PETSC ERROR: #9 KSPSolve() line 604 in /opt/petsc/src/ksp/ksp/interface/itfunc.c</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[0]PETSC ERROR: #10 KSPSolve_Chebyshev() line 381 in /opt/petsc/src/ksp/ksp/impls/cheby/cheby.c</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[0]PETSC ERROR: #11 KSPSolve() line 604 in /opt/petsc/src/ksp/ksp/interface/itfunc.c</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[0]PETSC ERROR: #12 PCMGMCycle_Private() line 19 in /opt/petsc/src/ksp/pc/impls/mg/mg.c</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[0]PETSC ERROR: #13 PCApply_MG() line 340 in /opt/petsc/src/ksp/pc/impls/mg/mg.c</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[0]PETSC ERROR: #14 PCApply() line 483 in /opt/petsc/src/ksp/pc/interface/precon.c</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[0]PETSC ERROR: #15 KSP_PCApply() line 242 in /opt/petsc/include/petsc/private/kspimpl.h</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[0]PETSC ERROR: #16 KSPInitialResidual() line 63 in /opt/petsc/src/ksp/ksp/interface/itres.c</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[0]PETSC ERROR: #17 KSPSolve_GMRES() line 235 in /opt/petsc/src/ksp/ksp/impls/gmres/gmres.c</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[0]PETSC ERROR: #18 KSPSolve() line 604 in /opt/petsc/src/ksp/ksp/interface/itfunc.c</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[0]PETSC ERROR: #19 SNESSolve_NEWTONLS() line 233 in /opt/petsc/src/snes/impls/ls/ls.c</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[0]PETSC ERROR: #20 SNESSolve() line 3906 in /opt/petsc/src/snes/interface/snes.c</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[1]PETSC ERROR: #1 VecNorm_MPI() line 63 in /opt/petsc/src/vec/vec/impls/mpi/pvec2.c</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[1]PETSC ERROR: #2 VecNorm() line 242 in /opt/petsc/src/vec/vec/interface/rvector.c</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[1]PETSC ERROR: #3 VecNormalize() line 337 in /opt/petsc/src/vec/vec/interface/rvector.c</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[1]PETSC ERROR: #4 KSPGMRESCycle() line 127 in /opt/petsc/src/ksp/ksp/impls/gmres/gmres.c</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[1]PETSC ERROR: #5 KSPSolve_GMRES() line 236 in /opt/petsc/src/ksp/ksp/impls/gmres/gmres.c</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[1]PETSC ERROR: #6 KSPSolve() line 604 in /opt/petsc/src/ksp/ksp/interface/itfunc.c</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[1]PETSC ERROR: #7 KSPSolve_Chebyshev() line 381 in /opt/petsc/src/ksp/ksp/impls/cheby/cheby.c</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[1]PETSC ERROR: #8 KSPSolve() line 604 in /opt/petsc/src/ksp/ksp/interface/itfunc.c</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[1]PETSC ERROR: #9 PCMGMCycle_Private() line 19 in /opt/petsc/src/ksp/pc/impls/mg/mg.c</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[1]PETSC ERROR: #10 PCApply_MG() line 340 in /opt/petsc/src/ksp/pc/impls/mg/mg.c</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[1]PETSC ERROR: #11 PCApply() line 483 in /opt/petsc/src/ksp/pc/interface/precon.c</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[1]PETSC ERROR: #12 KSP_PCApply() line 242 in /opt/petsc/include/petsc/private/kspimpl.h</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[1]PETSC ERROR: #13 KSPInitialResidual() line 63 in /opt/petsc/src/ksp/ksp/interface/itres.c</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[1]PETSC ERROR: #14 KSPSolve_GMRES() line 235 in /opt/petsc/src/ksp/ksp/impls/gmres/gmres.c</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[1]PETSC ERROR: #15 KSPSolve() line 604 in /opt/petsc/src/ksp/ksp/interface/itfunc.c</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[1]PETSC ERROR: #16 SNESSolve_NEWTONLS() line 233 in /opt/petsc/src/snes/impls/ls/ls.c</span></font></div><div style="margin: 0px;" class=""><font face="Menlo" class=""><span style="font-size: 12px;" class="">[1]PETSC ERROR: #17 SNESSolve() line 3906 in /opt/petsc/src/snes/interface/snes.c</span></font></div><div style="margin: 0px; font-family: Menlo; min-height: 16px;" class=""><br class=""></div><div class="">and then it just sits there.</div><div class=""><br class=""></div><div class="">
<span class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px;">-gideon</span>
</div>
<br class=""><div><blockquote type="cite" class=""><div class="">On Sep 19, 2015, at 7:18 PM, Matthew Knepley <<a href="mailto:knepley@gmail.com" class="">knepley@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" style="font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><div class="gmail_extra"><div class="gmail_quote">On Sat, Sep 19, 2015 at 6:10 PM, Gideon Simpson<span class="Apple-converted-space"> </span><span dir="ltr" class=""><<a href="mailto:gideon.simpson@gmail.com" target="_blank" class="">gideon.simpson@gmail.com</a>></span><span class="Apple-converted-space"> </span>wrote:<br class=""><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;">It gives zero pivots on:<br class=""></blockquote><div class=""><br class=""></div><div class="">Send -ksp_view</div><div class=""><br class=""></div><div class=""> Matt</div><div class=""> </div><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;">-ksp_type gmres -pc_type gamg -mg_coarse_sub_pc_factor_shift_type NONZERO<br class=""><br class="">It has no problem with:<br class=""><br class="">-ksp_type gmres -pc_type bjacobi -sub_pc_factor_shift_type NONZERO<br class=""><br class="">-gideon<br class=""><br class="">> On Sep 19, 2015, at 12:37 AM, Barry Smith <<a href="mailto:bsmith@mcs.anl.gov" class="">bsmith@mcs.anl.gov</a>> wrote:<br class="">><br class="">><br class="">>> On Sep 18, 2015, at 10:42 PM, Gideon Simpson <<a href="mailto:gideon.simpson@gmail.com" class="">gideon.simpson@gmail.com</a>> wrote:<br class="">>><br class="">>> I’m getting a zero pivot in LU factorization error when I use the gamg preconditioner. I know there is some issue with my matrix, as I get the same error if I use pc_type lu in serial, or if I use pc_type bjacobi in parallel, both of which are corrected by adding -pc_factor_shift_type NONZERO -sub_pc_factor_shift_type NONZERO, respectively. However, I can’t quite seem to figure out the analogous command to use for the gamg preconditioner.<br class="">><br class="">> You need to figure out the prefix to the option that is added by gamg to set it. One way is to run with the gamg options and -help and grep for factor_shift.<br class="">><br class="">> The option is like something like -mg_coarse_pc_factor_shift_type nonzero or -mg_coarse_sub_pc_factor_shift_type<br class="">><br class="">> Barry<br class="">><br class="">>><br class="">>> -gideon<br class="">>><br class="">><br class=""><br class=""></blockquote></div><br class=""><br clear="all" class=""><div class=""><br class=""></div>--<span class="Apple-converted-space"> </span><br class=""><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 class="">-- Norbert Wiener</div></div></div></div></blockquote></div><br class=""></div></body></html>