On Mon, Sep 4, 2017 at 12:46 PM, Michael Povolotskyi <[email protected]> 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 ------------------------------------------------------------------------------ Check out the vibrant tech community on one of the world's most engaging tech sites, Slashdot.org! http://sdm.link/slashdot _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
