Evan : > PCFactorSetMatSolverPackage(pc_fetd_dt, MATSOLVERMUMPS); You use mumps' Cholesky.
> PCFactorSetShiftType(pc_fetd_dt,MAT_SHIFT_POSITIVE_DEFINITE); This option only works for PETSc Cholesky, not mumps. I'm not aware of any option in mumps to add a shift to zero pivot. I'll send a request to mumps developer about it. Hong
