<div dir="ltr"><div>Dear PETSc team,</div><div><br></div><div>I am part of the development team of legacy fortran code with a tailored Newton's method. The software is already using PETSc's linear solvers and we enjoy it. Now I would like to evaluate the SNES solver.</div><div>I have already extracted a function in order to compute the Jacobian and another one to compute the residual. </div><div>But there is something I cannot figure out : at each Newton's iteration, our solver needs to know the unknowns value in order to compute the Jacobian. But the increment vector is computed within the SNES. </div><div>How can I synchronize PETSc's vector of unknowns  and mine? Is there some kind of SNESSetPostSolveShell ?</div><div><br></div><div>Thanks for developping PETSc,</div><div>Nicolas</div></div>