Hello,<br><br>I got the error message "Matrix is missing diagonal entry 23548!" when I tried to run my code. I print out the 23548th row and find that the 23548th element of this row is one and the others are zero. Can anyone tell me what's maybe the problem and how to fix this? Thanks.<br>
<br>Best Regards,<br>Rongliang<br><br>/bin/bash: SHELL: readonly variable<br>/bin/bash: PATH: readonly variable<br>[255]PETSC ERROR: --------------------- Error Message ------------------------------------<br>[255]PETSC ERROR: Object is in wrong state!<br>
[255]PETSC ERROR: Matrix is missing diagonal entry 23548!<br>[255]PETSC ERROR: ------------------------------------------------------------------------<br>[255]PETSC ERROR: Petsc Release Version 3.1.0, Patch 8, Thu Mar 17 13:37:48 CDT 2011<br>
[255]PETSC ERROR: See docs/changes/index.html for recent updates.<br>[255]PETSC ERROR: See docs/faq.html for hints about trouble shooting.<br>[255]PETSC ERROR: See docs/index.html for manual pages.<br>[255]PETSC ERROR: ------------------------------------------------------------------------<br>
[255]PETSC ERROR: /home/rchen/soft/fixedmesh/Cannula/./joab on a bgl-ibm-g named R02M1N0 by Unknown Wed Aug 31 17:36:32 2011<br>[255]PETSC ERROR: Libraries linked from /home/rchen/soft/petsc-3.1-p8-nodebug/bgl-ibm-goto-O3_440d/lib<br>
<br>[255]PETSC ERROR: ------------------------------------------------------------------------<br>[255]PETSC ERROR: MatILUFactorSymbolic_SeqAIJ_ilu0() line 1627 in src/mat/impls/aij/seq/aijfact.c<br>[255]PETSC ERROR: MatILUFactorSymbolic_SeqAIJ() line 1731 in src/mat/impls/aij/seq/aijfact.c<br>
[255]PETSC ERROR: MatILUFactorSymbolic() line 5464 in src/mat/interface/matrix.c<br>[255]PETSC ERROR: PCSetUp_ILU() line 204 in src/ksp/pc/impls/factor/ilu/ilu.c<br>[255]PETSC ERROR: PCSetUp() line 795 in src/ksp/pc/interface/precon.c<br>
[255]PETSC ERROR: KSPSetUp() line 237 in src/ksp/ksp/interface/itfunc.c<br>[255]PETSC ERROR: PCSetUpOnBlocks_ASM() line 335 in src/ksp/pc/impls/asm/asm.c<br>[255]PETSC ERROR: PCSetUpOnBlocks() line 828 in src/ksp/pc/interface/precon.c<br>
[255]PETSC ERROR: KSPSetUpOnBlocks() line 159 in src/ksp/ksp/interface/itfunc.c<br>[255]PETSC ERROR: KSPSolve() line 354 in src/ksp/ksp/interface/itfunc.c<br>[255]PETSC ERROR: SNES_KSPSolve() line 2944 in src/snes/interface/snes.c<br>
[255]PETSC ERROR: SNESSolve() line 2255 in src/snes/interface/snes.c<br>[255]PETSC ERROR: TimeStep() line 1462 in /home/rlchen/rlchen/soft/fixedmesh/code_changed/joab.c<br>[255]PETSC ERROR: SolveSteadyState() line 1362 in /home/rlchen/rlchen/soft/fixedmesh/code_changed/joab.c<br>
[255]PETSC ERROR: ComputefixedBoundary() line 505 in /home/rlchen/rlchen/soft/fixedmesh/code_changed/joab.c<br>[255]PETSC ERROR: JoabCtxSolve() line 165 in /home/rlchen/rlchen/soft/fixedmesh/code_changed/joab.c<br>[255]PETSC ERROR: main() line 107 in /home/rlchen/rlchen/soft/fixedmesh/code_changed/joab.c<br>
.....................<br>