[petsc-users] Matrix-free TS

Amneet Bhalla mail2amneet at gmail.com
Mon Mar 27 12:56:30 CDT 2017


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?


-- 
--Amneet
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20170327/c83ad4b4/attachment.html>


More information about the petsc-users mailing list