<div dir="ltr"><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Mar 27, 2017 at 10:42 AM, 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:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">1) sounds easier and matches your current code.</blockquote></div><br>Great! </div><div class="gmail_extra"><br></div><div class="gmail_extra">For matrix free TS, the form of F(t,u,u_t) =  <span style="color:rgb(80,0,80);font-size:12.800000190734863px">rho u*_t - div ( sigma) as written  "on paper" and should not </span></div><div class="gmail_extra"><span style="color:rgb(80,0,80);font-size:12.800000190734863px">contain time stepping factors like (1/2) that comes from CN?  Similarly for RHS = G(t,u). What about Jacobian of F (MATSHELL </span></div><div class="gmail_extra"><span style="color:rgb(80,0,80);font-size:12.800000190734863px">to solve momentum eqn) --- </span><span style="font-size:12.800000190734863px;color:rgb(80,0,80)">should that contain these factors?</span></div><div class="gmail_extra"><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div>--Amneet <br><br></div><div><br></div><div><br></div></div></div>
</div></div>