> On Mar 27, 2017, at 12:56 PM, Amneet Bhalla <[email protected]> wrote:
> 
> 
> 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?

   None of the functions should contain factors associated with a particular 
time-stepping scheme.

   Barry

> 
> 
> -- 
> --Amneet 
> 
> 
> 

Reply via email to