Changed the linear system solves back to GMRES; with ILU2 preconditioning, appears to still converge