<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Fri, May 5, 2017 at 11:14 AM, Karin&NiKo <span dir="ltr"><<a href="mailto:niko.karin@gmail.com" target="_blank">niko.karin@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><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></blockquote><div><br></div><div>The FormJacobian() function that you pass in gets the current guess as an argument.</div><div><br></div><div>  Thanks,</div><div><br></div><div>    Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div> </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>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div>
</div></div>