Hi, I'd like to try following strategy: For a nonlinear solver, at the first Newton step, do a complete LU factorization, and use the factorized matrix as preconditioner matrix in the following iteration. How to implement this in petsc? Gong Ding