Commit 48020626 authored by teo's avatar teo
Browse files

Fix a bug: when the matrix are recomputed also the stabilization terms must be recomputed.

parent 13993871
......@@ -1044,7 +1044,7 @@ updateSystem(const double alpha,
M_Displayer.leaderPrint(" F- Updating the stabilization terms ... ");
chrono.start();
if ( M_resetStab || !M_reuseStab || (M_matrStab.get() == 0) )
if ( M_resetStab || !M_reuseStab || (M_matrStab.get() == 0) || M_recomputeMatrix )
{
M_matrStab.reset ( new matrix_type(M_localMap) );
M_ipStab.apply( *M_matrStab, betaVec, false );
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment