Personal tools
You are here: Home / openCMISS / Wiki / Libmesh Time-Dependent Problems
Log in

Forgot your password?

Libmesh Time-Dependent Problems

In version 0.5.0, you must reassemble your entire matrix and right-hand side at each time step. I do not need to speak of the inefficiency of this. The developer's speak of a the option of not doing this at each time step but I do not see how that is other than rewriting the code. Look on line 153 of 'implicit_system.C' for the assemble routine that zeros out the matrices and rhs's creating the need for reassembly.

The option 'System::assemble_before_solve' has been added in CVS to the "System class": 'include/solvers/system.h' revision 1.18 has:

date: 2006/02/16 22:17:57;  author: jwpeterson;  state: Exp;  lines: +21 -1
Added a boolean flag to the system class which allows you
to skip assembling during the solve step if you know that for some reason the
system has already been assembled e.g. during a previous call to assemble().
This can be useful if you are solving an implicit system with multiple right
hand sides where you will need to solve multiple times but only assemble once.

Use of the option is mentioned in the mailing lists: "Controlling when to assemble":