Thanks for the suggestions, I'm going to try both the possibilities. Lorenzo
Il giorno 12/ott/2009, alle ore 18.01, John Peterson ha scritto: > On Mon, Oct 12, 2009 at 10:57 AM, Roy Stogner <[email protected] > > wrote: >> >> On Mon, 12 Oct 2009, Lorenzo Botti wrote: >> >>> I'm solving the incompressible NS equations with a fully implicit >>> dg code >>> and I would like >>> to know if there is an easy way to solve for the increment instead >>> of the >>> solution itself. >> >>> I'm using a transient_system but maybe there is a best framework >>> to do that. >> >> The simplest change might be for you to use a >> TransientSystem<NonlinearImplicitSystem> and solve it with the >> PetscNonlinearSolver; that expects you to build just the global >> residual and Jacobian for your time step. A more radical change >> might >> be for you to create a FEMSystem subclass; that framework doesn't >> even >> require you to worry about the time stepping, just to build the >> differential and the algebraic parts of element and side residuals in >> separate methods. > > FYI: FEMSystem example in ex18. > > -- > John ------------------------------------------------------------------------------ Come build with us! The BlackBerry(R) Developer Conference in SF, CA is the only developer event you need to attend this year. Jumpstart your developing skills, take BlackBerry mobile applications to market and stay ahead of the curve. Join us from November 9 - 12, 2009. Register now! http://p.sf.net/sfu/devconference _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
