[EMAIL PROTECTED] writes: > >>> I want to solve a pde equation involving a parameter to be incremented > >>> sucessively. After obtaining a solution I want to use it as a guess > >>> for > >>> the system to be solved next. How can I implement this? > > > >> I'm not an expert on Libmesh, but I recall that PETSC allows one to > >> set the iinitial guess. If PETSc can be accessed through the libmesh > >> interface then it should be possible quite easily. > > > > Indeed it is -- this is in fact the default behavior in libMesh. Whatever > > is in the System.solution vector is handed off to the linear solver as the > > initial guess. You should be able to verify this by calling > > > > equation_systems.solve(); > > equation_systems.solve(); > > > > And seeing that the second solve converges in 0 iterations. > > > In fact I obtain: >
Is this being run in parallel by any chance? I wonder if you need a system.update() in between the solves so that during the second assembly you get the most recent values from current_local_solution? Any chance you could post a short version of your code that shows the problem? -J ------------------------------------------------------------------------- This SF.net email is sponsored by: Microsoft Defy all challenges. Microsoft(R) Visual Studio 2008. http://clk.atdmt.com/MRT/go/vse0120000070mrt/direct/01/ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
