<p>Maybe a mistake implementing boundary conditions? </p>
<div class="gmail_quote">On Aug 31, 2011 7:15 PM, "Rongliang Chen" <<a href="mailto:rongliang.chan@gmail.com">rongliang.chan@gmail.com</a>> wrote:<br type="attribution">> Hello,<br>> <br>> I got the error message "Matrix is missing diagonal entry 23548!" when I<br>
> tried to run my code. I print out the 23548th row and find that the 23548th<br>> element of this row is one and the others are zero. Can anyone tell me<br>> 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>> ------------------------------------<br>
> [255]PETSC ERROR: Object is in wrong state!<br>> [255]PETSC ERROR: Matrix is missing diagonal entry 23548!<br>> [255]PETSC ERROR:<br>> ------------------------------------------------------------------------<br>
> [255]PETSC ERROR: Petsc Release Version 3.1.0, Patch 8, Thu Mar 17 13:37:48<br>> 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>> ------------------------------------------------------------------------<br>> [255]PETSC ERROR: /home/rchen/soft/fixedmesh/Cannula/./joab on a bgl-ibm-g<br>
> named R02M1N0 by Unknown Wed Aug 31 17:36:32 2011<br>> [255]PETSC ERROR: Libraries linked from<br>> /home/rchen/soft/petsc-3.1-p8-nodebug/bgl-ibm-goto-O3_440d/lib<br>> <br>> [255]PETSC ERROR:<br>> ------------------------------------------------------------------------<br>
> [255]PETSC ERROR: MatILUFactorSymbolic_SeqAIJ_ilu0() line 1627 in<br>> src/mat/impls/aij/seq/aijfact.c<br>> [255]PETSC ERROR: MatILUFactorSymbolic_SeqAIJ() line 1731 in<br>> src/mat/impls/aij/seq/aijfact.c<br>
> [255]PETSC ERROR: MatILUFactorSymbolic() line 5464 in<br>> src/mat/interface/matrix.c<br>> [255]PETSC ERROR: PCSetUp_ILU() line 204 in<br>> 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<br>> src/ksp/pc/impls/asm/asm.c<br>> [255]PETSC ERROR: PCSetUpOnBlocks() line 828 in<br>
> src/ksp/pc/interface/precon.c<br>> [255]PETSC ERROR: KSPSetUpOnBlocks() line 159 in<br>> 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<br>> /home/rlchen/rlchen/soft/fixedmesh/code_changed/joab.c<br>> [255]PETSC ERROR: SolveSteadyState() line 1362 in<br>
> /home/rlchen/rlchen/soft/fixedmesh/code_changed/joab.c<br>> [255]PETSC ERROR: ComputefixedBoundary() line 505 in<br>> /home/rlchen/rlchen/soft/fixedmesh/code_changed/joab.c<br>> [255]PETSC ERROR: JoabCtxSolve() line 165 in<br>
> /home/rlchen/rlchen/soft/fixedmesh/code_changed/joab.c<br>> [255]PETSC ERROR: main() line 107 in<br>> /home/rlchen/rlchen/soft/fixedmesh/code_changed/joab.c<br>> .....................<br></div>