Thank you.
I found that before calling system.update() I have to call
system.solution->close()
Is this technically correct?
Michael.
On 09/05/2017 11:17 AM, John Peterson wrote:
On Mon, Sep 4, 2017 at 12:46 PM, Michael Povolotskyi
<mpovo...@purdue.edu <mailto:mpovo...@purdue.edu>> 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
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users
