[petsc-users] Matrix-free TS

[Barry Smith]

> On Mar 27, 2017, at 12:56 PM, Amneet Bhalla <mail2amneet at gmail.com> wrote:
> 
> 
> On Mon, Mar 27, 2017 at 10:42 AM, Matthew Knepley <knepley at gmail.com> 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 
> 
> 
>