From: Benjamin K. <ben...@na...> - 2008-02-21 19:18:33
|
>> 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. -Ben |