<div dir="ltr">Apologies, in the end it seems this was more of a Firedrake question: with the help of Lawrence Mitchell, I now believe I should simply intercept SNESFormFunction().<br><div><div class="gmail_extra"><br><div class="gmail_quote">On 24 March 2016 at 17:39, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span><br>
> On Mar 24, 2016, at 10:18 AM, Andrew McRae <<a href="mailto:A.T.T.McRae@bath.ac.uk" target="_blank">A.T.T.McRae@bath.ac.uk</a>> wrote:<br>
><br>
> I have a finite element discretisation of the following nonlinear equation:<br>
><br>
> m*(phi_xx * phi_yy - phi_xy^2) = const,<br>
><br>
> solving for phi.  Unfortunately, the function m depends on phi in a complicated way -- let's assume I need to call my own function to handle this.<br>
<br>
</span>  Andrew<br>
<br>
   So you are actually solving<br>
<br>
  m(phi)*(phi_xx * phi_yy - phi_xy^2) - const = 0<br>
<br>
  with finite elements for phi?<br>
<br>
<br>
    What are you providing for a Jacobian?<br></blockquote><div><br></div><div>The Jacobian I give treats m as being independent of phi, so just whatever you get from linearising det(Hessian(phi)).<br> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<span><br>
<br>
><br>
> I'm using PETSc's SNES  in Python via petsc4py, within the wider environment of the software Firedrake.<br>
><br>
> Currently I'm hacking in the m update (and various output diagnostics) by writing a Python function "fakemonitor" and calling snes.setMonitor(fakemonitor).  This allows me to update m each nonlinear iteration.<br>
<br>
</span>    Hmm, I don't understand this. It sounds like you are passing (phi_xx * phi_yy - phi_xy^2) or something to SNES as the SNESFormFunction()? Why is this? Why not pass the entire function to SNES?<br></blockquote><div><br></div><div>I was passing in m(phi^n)(phi_xx * phi_yy - phi_xy^2) - const, i.e., m was effectively frozen from the last nonlinear iteration.  As stated above, I think it's as simple as arranging for m to be updated whenever SNESFormFunction() is called, which involves hacking Firedrake code but not PETSc code.<br><br></div><div>Thanks,<br></div><div>Andrew<br></div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<span><font color="#888888"><br>
  Barry<br>
</font></span><div><div><br>
><br>
> While this is better than nothing, there's still some problems: if I use e.g. snes_linesearch_type: "l2", the fnorms for lambda = 1.0, 0.5 and 0.0 are calculated without updating m, and so the step length taken is (seemingly) far from optimal.  I tried adding a damping parameter, but all this does is change the lambdas used to generate the quadratic fit; it doesn't actually make the step length smaller.<br>
><br>
> Is there some cleaner way to do what I want, perhaps by intercepting the fnorm calculation to update m, rather than abusing a custom monitor routine?<br>
><br>
> Thanks,<br>
> Andrew<br>
<br>
</div></div></blockquote></div><br></div></div></div>