zero pivot

billy at dem.uminho.pt billy at dem.uminho.pt
Sun Nov 12 14:28:22 CST 2006



I am getting zero pivot errors. I know the matrix has a high condition number.


see
http://www.mcs.anl.gov/petsc/petsc-as/documentation/troubleshooting.html#ZeroPivot!
[12]PETSC ERROR: Zero pivot row 3556 value 6.56537e-13 tolerance 1e-12 * rs
0.727229!
[12]PETSC ERROR: MatLUFactorNumeric() line 2212 in src/mat/interface/matrix.c
[12]PETSC ERROR: PCSetUp_ILU() line 551 in src/ksp/pc/impls/factor/ilu/ilu.c
[12]PETSC ERROR: PCSetUp() line 798 in src/ksp/pc/interface/precon.c
[12]PETSC ERROR: KSPSetUp() line 234 in src/ksp/ksp/interface/itfunc.c
[12]PETSC ERROR: PCSetUpOnBlocks_ASM() line 224 in src/ksp/pc/impls/asm/asm.c
[12]PETSC ERROR: PCSetUpOnBlocks() line 831 in src/ksp/pc/interface/precon.c
[12]PETSC ERROR: KSPSetUpOnBlocks() line 158 in src/ksp/ksp/interface/itfunc.c
[12]PETSC ERROR: KSPSolve() line 335 in src/ksp/ksp/interface/itfunc.c


When I use -sub_pc_factor_shift_nonzero and -pc_factor_zeropivot 1E-18 I get:


[10]PETSC ERROR: Caught signal number 8 FPE: Floating Point Exception,probably
divide by zero
[10]PETSC ERROR: Try option -start_in_debugger or -on_error_attach_debugger
[10]PETSC ERROR: or try http://valgrind.org on linux to find memory corruption
errors
[10]PETSC ERROR: likely location of problem given in stack below
[10]PETSC ERROR: --------------- Stack Frames ---------------
[10]PETSC ERROR: Note: The EXACT line numbers in the stack are not available,
[10]PETSC ERROR:       INSTEAD the line number of the start of the function
[10]PETSC ERROR:       is given.
[10]PETSC ERROR: [10] VecDot_MPI line 25 src/vec/vec/impls/mpi/pbvec.c
[10]PETSC ERROR: [10] VecDot line 86 src/vec/vec/interface/rvector.c
[10]PETSC ERROR: [10] KSPSolve_BCGS line 29 src/ksp/ksp/impls/bcgs/bcgs.c
[10]PETSC ERROR: --------------------------------------------
[10]PETSC ERROR: User provided function() line 0 in unknown directory unknown file
[10]PETSC ERROR: Signal received!
[10]PETSC ERROR:  !

Is there anyway to solve this?




More information about the petsc-users mailing list