Re: [Libmesh-users] question about transient linear system

[John P. <jwp...@gm...>]

On Mon, Sep 4, 2017 at 12:46 PM, Michael Povolotskyi <mpo...@pu...>
wrote:

> Hello,
>
> I'm going to solve transient linear system.
>
> I need to set initial conditions.
>
> I'm going to code a special initialization function for this.
>
> In this function, do I need to set system.solution or
> system.current_local_solution?
>

system.solution and then call system.update() to set values in
current_local_solution. As far as I know, one never sets values in
current_local_solution directly...

Even better would be to call system.project_solution() with an appropriate
function pointer object. This will take care of both for you. See, for
example, the init_cd() function in transient/transient_ex1/transient_ex1.C.

-- 
John