Dear Barry Sorry, I should complement that when I set DIFFERENT_NONZERO_PATTERN to the MatStructure, pastix works well. But it crash for SAME_NONZERO_PATTERN. However, I always use SAME_NONZERO_PATTERN before, which works for ksp solvers, mumps, superlu_dist, etc. Yours Gong Ding