Hi all, 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. Sometimes with non-linear fluxes building the residual on the rhs is easier than summing it to the product J*U^(n-1) as you do in ex13. The problem is that petsc_linear_solver use to initialize the solution vector to the previous solution and clearly the previous increment is not a good guess. I'm using a transient_system but maybe there is a best framework to do that.
Thanks Lorenzo ------------------------------------------------------------------------------ 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
