[petsc-users] LU does not support zero diagonal entries?

[Fande Kong]

Hi all,

I am solving a saddle point problem using LU, and the problem has some
zeros at diagonal. I got error messages as follows:

[0]PETSC ERROR: --------------------- Error Message
--------------------------------------------------------------
[0]PETSC ERROR: Object is in wrong state
[0]PETSC ERROR: Matrix is missing diagonal entry 100
[0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for
trouble shooting.
[0]PETSC ERROR: Petsc Development GIT revision: v3.6-130-g2e90825  GIT
Date: 2015-07-13 18:57:36 -0500
[0]PETSC ERROR: ./ex9 on a arch-darwin-cxx-debug named
engr2-4-134-dhcp.int.colorado.edu by livia Mon Jul 27 14:28:23 2015
[0]PETSC ERROR: Configure options --with-clanguage=cxx
--with-shared-libraries=1 --download-fblaslapack=1
--with-mpi-dir=/Users/livia/math/mpich-3.1_install --download-parmetis=1
--download-metis=1 --with-64-bit-indices=1 --download-netcdf=1
--download-exodusii=1 --download-hdf5=1 --with-debugging=yes
--with-c2html=0 --download-hypre=1 --download-superlu_dist=1
PETSC_DIR=/Users/livia/math/petsc-mat PETSC_ARCH=arch-darwin-cxx-debug
[0]PETSC ERROR: #1 MatLUFactorSymbolic_SeqAIJ() line 298 in
/Users/livia/math/petsc-mat/src/mat/impls/aij/seq/aijfact.c
[0]PETSC ERROR: #2 MatLUFactorSymbolic() line 2894 in
/Users/livia/math/petsc-mat/src/mat/interface/matrix.c
[0]PETSC ERROR: #3 PCSetUp_LU() line 127 in
/Users/livia/math/petsc-mat/src/ksp/pc/impls/factor/lu/lu.c
[0]PETSC ERROR: #4 PCSetUp() line 982 in
/Users/livia/math/petsc-mat/src/ksp/pc/interface/precon.c
[0]PETSC ERROR: #5 KSPSetUp() line 332 in
/Users/livia/math/petsc-mat/src/ksp/ksp/interface/itfunc.c
[0]PETSC ERROR: #6 main() line 103 in
/Users/livia/math/petsc-mat/src/ksp/pc/examples/tests/ex9.c
[0]PETSC ERROR: PETSc Option Table entries:
[0]PETSC ERROR: -af ./test64/A_1
[0]PETSC ERROR: -bf ./test64/B_1
[0]PETSC ERROR: -bsf ./test64/IS_1
[0]PETSC ERROR: -ksp_type gmres
[0]PETSC ERROR: -matload_block_size 1
[0]PETSC ERROR: -pc_factor_mat_solver_package petsc
[0]PETSC ERROR: -pc_type lu
[0]PETSC ERROR: -rf ./test64/RHS_1
[0]PETSC ERROR: -vecload_block_size 1
[0]PETSC ERROR: ----------------End of Error Message -------send entire
error message to petsc-maint at mcs.anl.gov----------
application called MPI_Abort(MPI_COMM_WORLD, 73) - process 0
[cli_0]: aborting job:
application called MPI_Abort(MPI_COMM_WORLD, 73) - process 0

===================================================================================
=   BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES
=   PID 39186 RUNNING AT engr2-4-134-dhcp.int.colorado.edu
=   EXIT CODE: 73
=   CLEANING UP REMAINING PROCESSES
=   YOU CAN IGNORE THE BELOW CLEANUP MESSAGES
===================================================================================


I was wondering if petsc does permutation to avoid having zero pivot? There
are any options to turn on this function?

Thanks,

Fande Kong,
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20150727/d2ab1b6e/attachment-0001.html>