On Mon, Mar 27, 2017 at 10:42 AM, Matthew Knepley <[email protected]> wrote:
> 1) sounds easier and matches your current code. Great! For matrix free TS, the form of F(t,u,u_t) = rho u*_t - div ( sigma) as written "on paper" and should not contain time stepping factors like (1/2) that comes from CN? Similarly for RHS = G(t,u). What about Jacobian of F (MATSHELL to solve momentum eqn) --- should that contain these factors? -- --Amneet
