Hi Matt, Thank you for your suggestion. It solves the problem.
Khai On Mon, Mar 23, 2015 at 8:47 AM, Matthew Knepley <[email protected]> wrote: > On Mon, Mar 23, 2015 at 7:41 AM, Khai Pham <[email protected]> wrote: > >> Dear Petsc-Users, >> >> I received the below error message. There could be some error in the >> matrix assembly and I try to find it out. Based on the error message, could >> you give me any hint where should I focus on for finding the bug ? Thanks >> > > If you want to use ILU, we require that a diagonal element exists in the > matrix. If it is 0, you must explicitly > set the 0. > > Thanks, > > Matt > > >> Khai >> >> >> [2]PETSC ERROR: --------------------- Error Message >> -------------------------------------------------------------- >> [2]PETSC ERROR: Object is in wrong state >> [2]PETSC ERROR: Matrix is missing diagonal entry 1 >> [2]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html >> for trouble shooting. >> [2]PETSC ERROR: Petsc Release Version 3.5.1, Jul, 24, 2014 >> [2]PETSC ERROR: ./3DMMF on a arch-linux2-c-debug named khp by khp Mon Mar >> 23 08:24:40 2015 >> [2]PETSC ERROR: Configure options --with-cc=gcc --with-cxx=g++ >> --with-fc=gfortran --download-fblaslapack --download-mpich --download-hdf5 >> --download-metis --download-parmetis >> [2]PETSC ERROR: #1 MatILUFactorSymbolic_SeqAIJ() line 1738 in >> /home/khp/bin/petsc-3.5.1/src/mat/impls/aij/seq/aijfact.c >> [2]PETSC ERROR: #2 MatILUFactorSymbolic() line 6203 in >> /home/khp/bin/petsc-3.5.1/src/mat/interface/matrix.c >> [2]PETSC ERROR: #3 PCSetUp_ILU() line 206 in >> /home/khp/bin/petsc-3.5.1/src/ksp/pc/impls/factor/ilu/ilu.c >> [2]PETSC ERROR: #4 PCSetUp() line 902 in >> /home/khp/bin/petsc-3.5.1/src/ksp/pc/interface/precon.c >> [2]PETSC ERROR: #5 KSPSetUp() line 305 in >> /home/khp/bin/petsc-3.5.1/src/ksp/ksp/interface/itfunc.c >> [2]PETSC ERROR: #6 PCSetUpOnBlocks_BJacobi_Singleblock() line 651 in >> /home/khp/bin/petsc-3.5.1/src/ksp/pc/impls/bjacobi/bjacobi.c >> [2]PETSC ERROR: #7 PCSetUpOnBlocks() line 935 in >> /home/khp/bin/petsc-3.5.1/src/ksp/pc/interface/precon.c >> [2]PETSC ERROR: #8 KSPSetUpOnBlocks() line 166 in >> /home/khp/bin/petsc-3.5.1/src/ksp/ksp/interface/itfunc.c >> [2]PETSC ERROR: #9 KSPSolve() line 418 in >> /home/khp/bin/petsc-3.5.1/src/ksp/ksp/interface/itfunc.c >> [2]PETSC ERROR: #10 SNESSolve_NEWTONLS() line 232 in >> /home/khp/bin/petsc-3.5.1/src/snes/impls/ls/ls.c >> [2]PETSC ERROR: #11 SNESSolve() line 3743 in >> /home/khp/bin/petsc-3.5.1/src/snes/interface/snes.c >> [2]PETSC ERROR: #12 PF_solver() line 63 in src/PF_solver.cpp >> [1]PETSC ERROR: --------------------- Error Message >> -------------------------------------------------------------- >> [1]PETSC ERROR: Object is in wrong state >> [1]PETSC ERROR: Matrix is missing diagonal entry 0 >> [1]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html >> for trouble shooting. >> [1]PETSC ERROR: Petsc Release Version 3.5.1, Jul, 24, 2014 >> [1]PETSC ERROR: ./3DMMF on a arch-linux2-c-debug named khp by khp Mon Mar >> 23 08:24:40 2015 >> [1]PETSC ERROR: Configure options --with-cc=gcc --with-cxx=g++ >> --with-fc=gfortran --download-fblaslapack --download-mpich --download-hdf5 >> --download-metis --download-parmetis >> [1]PETSC ERROR: #1 MatILUFactorSymbolic_SeqAIJ() line 1738 in >> /home/khp/bin/petsc-3.5.1/src/mat/impls/aij/seq/aijfact.c >> [1]PETSC ERROR: #2 MatILUFactorSymbolic() line 6203 in >> /home/khp/bin/petsc-3.5.1/src/mat/interface/matrix.c >> [1]PETSC ERROR: #3 PCSetUp_ILU() line 206 in >> /home/khp/bin/petsc-3.5.1/src/ksp/pc/impls/factor/ilu/ilu.c >> [1]PETSC ERROR: #4 PCSetUp() line 902 in >> /home/khp/bin/petsc-3.5.1/src/ksp/pc/interface/precon.c >> [1]PETSC ERROR: #5 KSPSetUp() line 305 in >> /home/khp/bin/petsc-3.5.1/src/ksp/ksp/interface/itfunc.c >> [1]PETSC ERROR: #6 PCSetUpOnBlocks_BJacobi_Singleblock() line 651 in >> /home/khp/bin/petsc-3.5.1/src/ksp/pc/impls/bjacobi/bjacobi.c >> [1]PETSC ERROR: #7 PCSetUpOnBlocks() line 935 in >> /home/khp/bin/petsc-3.5.1/src/ksp/pc/interface/precon.c >> [1]PETSC ERROR: #8 KSPSetUpOnBlocks() line 166 in >> /home/khp/bin/petsc-3.5.1/src/ksp/ksp/interface/itfunc.c >> [1]PETSC ERROR: #9 KSPSolve() line 418 in >> /home/khp/bin/petsc-3.5.1/src/ksp/ksp/interface/itfunc.c >> [1]PETSC ERROR: #10 SNESSolve_NEWTONLS() line 232 in >> /home/khp/bin/petsc-3.5.1/src/snes/impls/ls/ls.c >> [1]PETSC ERROR: #11 SNESSolve() line 3743 in >> /home/khp/bin/petsc-3.5.1/src/snes/interface/snes.c >> > > > > -- > What most experimenters take for granted before they begin their > experiments is infinitely more interesting than any results to which their > experiments lead. > -- Norbert Wiener >
