Dear Barry, Trouble you again. I found that my nonlinear solver failed to convergence with superlu_dist (even broken some times), but works well with mumps and pastix. When I change SAME_NONZERO_PATTERN to DIFFERENT_NONZERO_PATTERN, everything goes well. Will you please try to fix this? Gong Ding