On Tue, Nov 15, 2011 at 5:43 PM, Jed Brown <span dir="ltr"><<a href="mailto:jedbrown@mcs.anl.gov">jedbrown@mcs.anl.gov</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div class="im"><div class="gmail_quote">On Tue, Nov 15, 2011 at 17:37, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

  b) FormFunction and FormJacobian in a SNES could be used by an enclosing TS to do pseudo-timestepping</blockquote></div><br></div><div>This can't work. What equation would the FormFunction and FormJacobian correspond to?</div>
</blockquote><div><br></div><div>What are you on? It is F(x) = 0. This embedding will not work for all problems, but it will work for many, and be easy. If you</div><div>have a more complicated system, by all means us TS, however I am pointing out that there are useful regimes where this</div>
<div>kind of embedding would work.</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>The TS *defines* the equation that needs to be solved. The user provides information according to an interface. The TS interface is more general, and really is required for transient things. The SNESFormFunction isn't enough because you haven't identified positivity and constraints. Add that transient term to all variables equally and you are going to see some confusing stuff like huge conservation error and incorrect boundary conditions, such that you could easily end up on the wrong branch or mis-diagnose a limit cycle versus steady state.</div>

</blockquote></div><br><br clear="all"><div><br></div>-- <br>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<br>