Fix call to compress in reinit of Trilinos sparse matrix

diff --git a/source/lac/trilinos_sparse_matrix.cc b/source/lac/trilinos_sparse_matrix.cc
index 01b51235e6b439e55394a4214322f0dbdfd529e5..2c464c1ead2ddbe98a9820dee818526828d8e9e4 100644 (file)
--- a/source/lac/trilinos_sparse_matrix.cc
+++ b/source/lac/trilinos_sparse_matrix.cc
@@ -788,8 +788,8 @@ namespace TrilinosWrappers
     else
       nonlocal_matrix.reset ();
 
-    compress(VectorOperation::insert);
     last_action = Zero;
+    compress(VectorOperation::insert);
   }
 
 
@@ -805,12 +805,14 @@ namespace TrilinosWrappers
     nonlocal_matrix_exporter.reset();
     matrix.reset (new Epetra_FECrsMatrix
                   (Copy, sparse_matrix.trilinos_sparsity_pattern(), false));
+
     if (sparse_matrix.nonlocal_matrix != 0)
       nonlocal_matrix.reset (new Epetra_CrsMatrix
                              (Copy, sparse_matrix.nonlocal_matrix->Graph()));
     else
       nonlocal_matrix.reset();
 
+    last_action = Zero;
     compress(VectorOperation::insert);
   }