<div dir="ltr">Hi all,<div><br></div><div>I am solving a saddle point problem using LU, and the problem has some zeros at diagonal. I got error messages as follows:</div><div><br></div><div><div>[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------</div><div>[0]PETSC ERROR: Object is in wrong state</div><div>[0]PETSC ERROR: Matrix is missing diagonal entry 100</div><div>[0]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html">http://www.mcs.anl.gov/petsc/documentation/faq.html</a> for trouble shooting.</div><div>[0]PETSC ERROR: Petsc Development GIT revision: v3.6-130-g2e90825  GIT Date: 2015-07-13 18:57:36 -0500</div><div>[0]PETSC ERROR: ./ex9 on a arch-darwin-cxx-debug named <a href="http://engr2-4-134-dhcp.int.colorado.edu">engr2-4-134-dhcp.int.colorado.edu</a> by livia Mon Jul 27 14:28:23 2015</div><div>[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</div><div>[0]PETSC ERROR: #1 MatLUFactorSymbolic_SeqAIJ() line 298 in /Users/livia/math/petsc-mat/src/mat/impls/aij/seq/aijfact.c</div><div>[0]PETSC ERROR: #2 MatLUFactorSymbolic() line 2894 in /Users/livia/math/petsc-mat/src/mat/interface/matrix.c</div><div>[0]PETSC ERROR: #3 PCSetUp_LU() line 127 in /Users/livia/math/petsc-mat/src/ksp/pc/impls/factor/lu/lu.c</div><div>[0]PETSC ERROR: #4 PCSetUp() line 982 in /Users/livia/math/petsc-mat/src/ksp/pc/interface/precon.c</div><div>[0]PETSC ERROR: #5 KSPSetUp() line 332 in /Users/livia/math/petsc-mat/src/ksp/ksp/interface/itfunc.c</div><div>[0]PETSC ERROR: #6 main() line 103 in /Users/livia/math/petsc-mat/src/ksp/pc/examples/tests/ex9.c</div><div>[0]PETSC ERROR: PETSc Option Table entries:</div><div>[0]PETSC ERROR: -af ./test64/A_1</div><div>[0]PETSC ERROR: -bf ./test64/B_1</div><div>[0]PETSC ERROR: -bsf ./test64/IS_1</div><div>[0]PETSC ERROR: -ksp_type gmres</div><div>[0]PETSC ERROR: -matload_block_size 1</div><div>[0]PETSC ERROR: -pc_factor_mat_solver_package petsc</div><div>[0]PETSC ERROR: -pc_type lu</div><div>[0]PETSC ERROR: -rf ./test64/RHS_1</div><div>[0]PETSC ERROR: -vecload_block_size 1</div><div>[0]PETSC ERROR: ----------------End of Error Message -------send entire error message to petsc-maint@mcs.anl.gov----------</div><div>application called MPI_Abort(MPI_COMM_WORLD, 73) - process 0</div><div>[cli_0]: aborting job:</div><div>application called MPI_Abort(MPI_COMM_WORLD, 73) - process 0</div><div><br></div><div>===================================================================================</div><div>=   BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES</div><div>=   PID 39186 RUNNING AT <a href="http://engr2-4-134-dhcp.int.colorado.edu">engr2-4-134-dhcp.int.colorado.edu</a></div><div>=   EXIT CODE: 73</div><div>=   CLEANING UP REMAINING PROCESSES</div><div>=   YOU CAN IGNORE THE BELOW CLEANUP MESSAGES</div><div>===================================================================================</div></div><div><br></div><div><br></div><div>I was wondering if petsc does permutation to avoid having zero pivot? There are any options to turn on this function?</div><div><br></div><div>Thanks,</div><div><br></div><div>Fande Kong,</div></div>