[Libmesh-devel] NewmarkSystem / LU - Solver

[Mathias N. <mne...@tu...>]

Hello,

I'm using the class NewmarkSystem. For solving the equation system I 
need to use a direct Solver.
For this reason I start my programm like this:

./programm -pc_type lu

The only drawback is that the Solver now calculates the LU-decomposition 
at every time step, as far as I know.
But this is just needed in the first time step and for further time 
steps only the forward and backward substitution is needed.

So, is there a easy way to handle this problem...

Thanks in advance,
Mathias N.