From 2c42823fc4b8d19933b0c684ae055dbe106e634e Mon Sep 17 00:00:00 2001
From: Timo Heister <timo.heister@gmail.com>
Date: Tue, 24 Feb 2015 12:51:50 -0500
Subject: [PATCH] fix a few broken tests

Looks like somebody messed up the VectorOperation in compress() in many
places.
---
 tests/mpi/petsc_distribute_01_block.cc                    | 2 +-
 tests/petsc/01.cc                                         | 2 +-
 tests/petsc/05.cc                                         | 2 +-
 tests/petsc/06.cc                                         | 2 +-
 tests/petsc/07.cc                                         | 2 +-
 tests/petsc/08.cc                                         | 2 +-
 tests/petsc/09.cc                                         | 2 +-
 tests/petsc/10.cc                                         | 2 +-
 tests/petsc/11.cc                                         | 2 +-
 tests/petsc/12.cc                                         | 2 +-
 tests/petsc/17.cc                                         | 2 +-
 tests/petsc/19.cc                                         | 2 +-
 tests/petsc/20.cc                                         | 2 +-
 tests/petsc/21.cc                                         | 2 +-
 tests/petsc/22.cc                                         | 4 ++--
 tests/petsc/23.cc                                         | 4 ++--
 tests/petsc/24.cc                                         | 2 +-
 tests/petsc/25.cc                                         | 2 +-
 tests/petsc/26.cc                                         | 2 +-
 tests/petsc/27.cc                                         | 2 +-
 tests/petsc/28.cc                                         | 2 +-
 tests/petsc/31.cc                                         | 2 +-
 tests/petsc/32.cc                                         | 2 +-
 tests/petsc/33.cc                                         | 2 +-
 tests/petsc/34.cc                                         | 2 +-
 tests/petsc/35.cc                                         | 4 ++--
 tests/petsc/36.cc                                         | 4 ++--
 tests/petsc/37.cc                                         | 2 +-
 tests/petsc/38.cc                                         | 4 ++--
 tests/petsc/39.cc                                         | 4 ++--
 tests/petsc/40.cc                                         | 6 +++---
 tests/petsc/41.cc                                         | 4 ++--
 tests/petsc/42.cc                                         | 4 ++--
 tests/petsc/43.cc                                         | 6 +++---
 tests/petsc/45.cc                                         | 4 ++--
 tests/petsc/46.cc                                         | 4 ++--
 tests/petsc/47.cc                                         | 6 +++---
 tests/petsc/48.cc                                         | 6 +++---
 tests/petsc/51.cc                                         | 3 +--
 tests/petsc/52.cc                                         | 2 +-
 tests/petsc/55.cc                                         | 2 +-
 tests/petsc/56.cc                                         | 2 +-
 tests/petsc/64.cc                                         | 2 +-
 tests/petsc/deal_solver_01.cc                             | 4 +---
 tests/petsc/deal_solver_02.cc                             | 4 +---
 tests/petsc/deal_solver_03.cc                             | 4 +---
 tests/petsc/deal_solver_04.cc                             | 4 +---
 tests/petsc/deal_solver_05.cc                             | 4 +---
 tests/petsc/full_matrix_00.cc                             | 4 ++--
 tests/petsc/full_matrix_01.cc                             | 2 +-
 tests/petsc/full_matrix_05.cc                             | 2 +-
 tests/petsc/full_matrix_06.cc                             | 2 +-
 tests/petsc/full_matrix_08.cc                             | 2 +-
 tests/petsc/full_matrix_09.cc                             | 2 +-
 tests/petsc/full_matrix_10.cc                             | 2 +-
 tests/petsc/full_matrix_iterator_01.cc                    | 2 +-
 tests/petsc/full_matrix_vector_01.cc                      | 5 ++---
 tests/petsc/full_matrix_vector_02.cc                      | 5 ++---
 tests/petsc/full_matrix_vector_03.cc                      | 6 +++---
 tests/petsc/full_matrix_vector_05.cc                      | 6 +++---
 tests/petsc/full_matrix_vector_06.cc                      | 4 ++--
 tests/petsc/full_matrix_vector_07.cc                      | 6 +++---
 tests/petsc/solver_01.cc                                  | 4 +---
 tests/petsc/solver_02.cc                                  | 4 +---
 tests/petsc/solver_03.cc                                  | 4 +---
 tests/petsc/solver_03_mf.cc                               | 4 +---
 tests/petsc/solver_03_precondition_boomeramg.cc           | 4 +---
 tests/petsc/solver_03_precondition_boomeramg_symmetric.cc | 4 +---
 tests/petsc/solver_03_precondition_eisenstat.cc           | 4 +---
 tests/petsc/solver_03_precondition_icc.cc                 | 4 +---
 tests/petsc/solver_03_precondition_lu.cc                  | 4 +---
 tests/petsc/solver_03_precondition_parasails.cc           | 4 +---
 tests/petsc/solver_03_precondition_sor.cc                 | 4 +---
 tests/petsc/solver_03_precondition_ssor.cc                | 4 +---
 tests/petsc/solver_04.cc                                  | 4 +---
 tests/petsc/solver_05.cc                                  | 4 +---
 tests/petsc/solver_06.cc                                  | 4 +---
 tests/petsc/solver_07.cc                                  | 4 +---
 tests/petsc/solver_08.cc                                  | 4 +---
 tests/petsc/solver_09.cc                                  | 4 +---
 tests/petsc/solver_10.cc                                  | 4 +---
 tests/petsc/solver_11.cc                                  | 4 +---
 tests/petsc/solver_12.cc                                  | 4 +---
 tests/petsc/solver_13.cc                                  | 4 +---
 tests/petsc/sparse_direct_mumps.cc                        | 4 +---
 tests/petsc/sparse_matrix_01.cc                           | 4 ++--
 tests/petsc/sparse_matrix_02.cc                           | 2 +-
 tests/petsc/sparse_matrix_iterator_01.cc                  | 2 +-
 tests/petsc/sparse_matrix_vector_01.cc                    | 5 ++---
 tests/petsc/sparse_matrix_vector_02.cc                    | 5 ++---
 tests/petsc/sparse_matrix_vector_03.cc                    | 6 +++---
 tests/petsc/sparse_matrix_vector_04.cc                    | 6 +++---
 tests/petsc/sparse_matrix_vector_05.cc                    | 6 +++---
 tests/petsc/sparse_matrix_vector_06.cc                    | 4 ++--
 tests/petsc/sparse_matrix_vector_07.cc                    | 6 +++---
 95 files changed, 135 insertions(+), 196 deletions(-)

diff --git a/tests/mpi/petsc_distribute_01_block.cc b/tests/mpi/petsc_distribute_01_block.cc
index 631a1d9dab..43fafab6b5 100644
--- a/tests/mpi/petsc_distribute_01_block.cc
+++ b/tests/mpi/petsc_distribute_01_block.cc
@@ -52,7 +52,7 @@ void test()
     vec.block(0)(i) = i;
   for (unsigned int i=vec.block(1).local_range().first; i<vec.block(1).local_range().second; ++i)
     vec.block(1)(i) = i;
-  vec.compress(VectorOperation::add);
+  vec.compress(VectorOperation::insert);
 
   // verify correctness so far
   {
diff --git a/tests/petsc/01.cc b/tests/petsc/01.cc
index a4de9a246a..3e29ed4b45 100644
--- a/tests/petsc/01.cc
+++ b/tests/petsc/01.cc
@@ -32,7 +32,7 @@ void test (PETScWrappers::SparseMatrix &m)
       if ((i+2*j+1) % 3 == 0)
         m.set (i,j, i*j*.5+.5);
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   // then make sure we retrieve the same ones
   for (unsigned int i=0; i<m.m(); ++i)
diff --git a/tests/petsc/05.cc b/tests/petsc/05.cc
index 27bb210372..07330cc82e 100644
--- a/tests/petsc/05.cc
+++ b/tests/petsc/05.cc
@@ -37,7 +37,7 @@ void test (PETScWrappers::SparseMatrix &m)
           ++counter;
         }
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   deallog << m.n_nonzero_elements() << std::endl;
   Assert (m.n_nonzero_elements() == counter,
diff --git a/tests/petsc/06.cc b/tests/petsc/06.cc
index aefef366b5..2c9d927e7c 100644
--- a/tests/petsc/06.cc
+++ b/tests/petsc/06.cc
@@ -32,7 +32,7 @@ void test (PETScWrappers::SparseMatrix &m)
       if ((i+2*j+1) % 3 == 0)
         m.set (i,j, i*j*.5+.5);
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   // compare against the exact value of the
   // l1-norm (max col-sum)
diff --git a/tests/petsc/07.cc b/tests/petsc/07.cc
index c61e90ecb5..90a9b745db 100644
--- a/tests/petsc/07.cc
+++ b/tests/petsc/07.cc
@@ -32,7 +32,7 @@ void test (PETScWrappers::SparseMatrix &m)
       if ((i+2*j+1) % 3 == 0)
         m.set (i,j, i*j*.5+.5);
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   // compare against the exact value of the
   // linfty-norm (max row-sum)
diff --git a/tests/petsc/08.cc b/tests/petsc/08.cc
index 2fb1dfa253..6fe6ec05ef 100644
--- a/tests/petsc/08.cc
+++ b/tests/petsc/08.cc
@@ -37,7 +37,7 @@ void test (PETScWrappers::SparseMatrix &m)
         }
   norm = std::sqrt(norm);
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   // compare against the exact value of the
   // l2-norm (max row-sum)
diff --git a/tests/petsc/09.cc b/tests/petsc/09.cc
index fb6b7a6e14..4914c88714 100644
--- a/tests/petsc/09.cc
+++ b/tests/petsc/09.cc
@@ -31,7 +31,7 @@ void test (PETScWrappers::SparseMatrix &m)
       if ((i+2*j+1) % 3 == 0)
         m.set (i,j, i*j*.5+.5);
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   // then multiply everything by 1.25 and
   // make sure we retrieve the values we
diff --git a/tests/petsc/10.cc b/tests/petsc/10.cc
index c0bc8b65ef..eeba560d9e 100644
--- a/tests/petsc/10.cc
+++ b/tests/petsc/10.cc
@@ -31,7 +31,7 @@ void test (PETScWrappers::SparseMatrix &m)
       if ((i+2*j+1) % 3 == 0)
         m.set (i,j, i*j*.5+.5);
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   // then divide everything by 4/3 and
   // make sure we retrieve the values we
diff --git a/tests/petsc/11.cc b/tests/petsc/11.cc
index 4cc621ed2d..0e8c572e93 100644
--- a/tests/petsc/11.cc
+++ b/tests/petsc/11.cc
@@ -29,7 +29,7 @@ void test (PETScWrappers::Vector &v)
   for (unsigned int i=0; i<v.size(); i+=1+i)
     v(i) = i;
 
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   Assert (v.size() == 100, ExcInternalError());
 
diff --git a/tests/petsc/12.cc b/tests/petsc/12.cc
index 74c51d4e7d..021e3e5e8a 100644
--- a/tests/petsc/12.cc
+++ b/tests/petsc/12.cc
@@ -36,7 +36,7 @@ void test (PETScWrappers::Vector &v)
       pattern[i] = true;
     }
 
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // check that they are ok, and this time
   // all of them
diff --git a/tests/petsc/17.cc b/tests/petsc/17.cc
index 3aa9d9951a..59f9ed3834 100644
--- a/tests/petsc/17.cc
+++ b/tests/petsc/17.cc
@@ -33,7 +33,7 @@ void test (PETScWrappers::Vector &v)
       v(i) = i;
       norm += std::fabs(1.*i);
     }
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // then check the norm
   Assert (v.l1_norm() == norm, ExcInternalError());
diff --git a/tests/petsc/19.cc b/tests/petsc/19.cc
index f04a700336..2f801f6488 100644
--- a/tests/petsc/19.cc
+++ b/tests/petsc/19.cc
@@ -33,7 +33,7 @@ void test (PETScWrappers::Vector &v)
       v(i) = i;
       norm = std::max(norm,fabs(i));
     }
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // then check the norm
   Assert (v.linfty_norm() == norm, ExcInternalError());
diff --git a/tests/petsc/20.cc b/tests/petsc/20.cc
index 1610b28891..6c845ce1ce 100644
--- a/tests/petsc/20.cc
+++ b/tests/petsc/20.cc
@@ -35,7 +35,7 @@ void test (PETScWrappers::Vector &v)
       v(i) = i;
       pattern[i] = true;
     }
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // multiply v with 5/4
   v *= 5./4.;
diff --git a/tests/petsc/21.cc b/tests/petsc/21.cc
index 036ed63cf4..eee63e666b 100644
--- a/tests/petsc/21.cc
+++ b/tests/petsc/21.cc
@@ -35,7 +35,7 @@ void test (PETScWrappers::Vector &v)
       v(i) = i;
       pattern[i] = true;
     }
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // multiply v with 3/4
   v /= 4./3.;
diff --git a/tests/petsc/22.cc b/tests/petsc/22.cc
index e6f2fc6982..821a1ff3bc 100644
--- a/tests/petsc/22.cc
+++ b/tests/petsc/22.cc
@@ -35,8 +35,8 @@ void test (PETScWrappers::Vector &v,
       v(i) = i;
     else
       w(i) = i;
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
 
   // make sure the scalar product is zero
   Assert (v*w == 0, ExcInternalError());
diff --git a/tests/petsc/23.cc b/tests/petsc/23.cc
index ef29358522..25f843160b 100644
--- a/tests/petsc/23.cc
+++ b/tests/petsc/23.cc
@@ -42,8 +42,8 @@ void test (PETScWrappers::Vector &v,
         }
     }
 
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
 
   // make sure the scalar product is zero
   Assert (v*w == product, ExcInternalError());
diff --git a/tests/petsc/24.cc b/tests/petsc/24.cc
index 3bbd795d0e..47ade91b2a 100644
--- a/tests/petsc/24.cc
+++ b/tests/petsc/24.cc
@@ -32,7 +32,7 @@ void test (PETScWrappers::Vector &v)
   for (unsigned int i=0; i<v.size(); ++i)
     if (i%3 == 0)
       v(i) = i+1.;
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // then clear it again and make sure the
   // vector is really empty
diff --git a/tests/petsc/25.cc b/tests/petsc/25.cc
index 993f4f4e0a..aa6715587b 100644
--- a/tests/petsc/25.cc
+++ b/tests/petsc/25.cc
@@ -30,7 +30,7 @@ void test (PETScWrappers::Vector &v)
   for (unsigned int i=0; i<v.size(); ++i)
     if (i%3 == 0)
       v(i) = i+1.;
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // then clear it again and make sure the
   // vector is really empty
diff --git a/tests/petsc/26.cc b/tests/petsc/26.cc
index 25309df9af..3ded098848 100644
--- a/tests/petsc/26.cc
+++ b/tests/petsc/26.cc
@@ -31,7 +31,7 @@ void test (PETScWrappers::Vector &v)
   for (unsigned int i=0; i<v.size(); ++i)
     if (i%3 == 0)
       v(i) = i+1.;
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   const unsigned int sz = v.size();
   v = 2;
diff --git a/tests/petsc/27.cc b/tests/petsc/27.cc
index 3c3632c06f..b773cb4933 100644
--- a/tests/petsc/27.cc
+++ b/tests/petsc/27.cc
@@ -30,7 +30,7 @@ void test (PETScWrappers::Vector &v)
   for (unsigned int i=0; i<v.size(); ++i)
     if (i%3 == 0)
       v(i) = i+1.;
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // then copy it
   PETScWrappers::Vector w (v.size());
diff --git a/tests/petsc/28.cc b/tests/petsc/28.cc
index 5e96a0a06d..130b4c73e5 100644
--- a/tests/petsc/28.cc
+++ b/tests/petsc/28.cc
@@ -31,7 +31,7 @@ void test (PETScWrappers::Vector &v)
   for (unsigned int i=0; i<v.size(); ++i)
     if (i%3 == 0)
       v(i) = i+1.;
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // then copy it to a vector of different
   // size
diff --git a/tests/petsc/31.cc b/tests/petsc/31.cc
index 875f12311c..0a1bcb31aa 100644
--- a/tests/petsc/31.cc
+++ b/tests/petsc/31.cc
@@ -33,7 +33,7 @@ void test (PETScWrappers::Vector &v)
       v(i) = i;
       norm += std::fabs(1.*i)*std::fabs(1.*i);
     }
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // then check the norm
   const double eps=typeid(PetscScalar)==typeid(double) ? 1e-14 : 1e-5;
diff --git a/tests/petsc/32.cc b/tests/petsc/32.cc
index d490ac0f3e..a372712006 100644
--- a/tests/petsc/32.cc
+++ b/tests/petsc/32.cc
@@ -33,7 +33,7 @@ void test (PETScWrappers::Vector &v)
       v(i) = i;
       sum += i;
     }
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // then check the norm
   const double eps=typeid(PetscScalar)==typeid(double) ? 1e-14 : 1e-5;
diff --git a/tests/petsc/33.cc b/tests/petsc/33.cc
index f367eb5594..62e686c748 100644
--- a/tests/petsc/33.cc
+++ b/tests/petsc/33.cc
@@ -33,7 +33,7 @@ void test (PETScWrappers::Vector &v)
       v(i) = i;
       sum += i*i*i;
     }
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // then check the norm
   const double eps=typeid(PetscScalar)==typeid(double) ? 1e-14 : 1e-5;
diff --git a/tests/petsc/34.cc b/tests/petsc/34.cc
index 707bbe2b0a..3d685b81dd 100644
--- a/tests/petsc/34.cc
+++ b/tests/petsc/34.cc
@@ -33,7 +33,7 @@ void test (PETScWrappers::Vector &v)
       v(i) = i;
       sum += i*i*i;
     }
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // set them to zero again
   v = 0;
diff --git a/tests/petsc/35.cc b/tests/petsc/35.cc
index 036130efdd..d3a1b0758a 100644
--- a/tests/petsc/35.cc
+++ b/tests/petsc/35.cc
@@ -36,8 +36,8 @@ void test (PETScWrappers::Vector &v,
         w(i) = i+1.;
     }
 
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
 
   v += w;
 
diff --git a/tests/petsc/36.cc b/tests/petsc/36.cc
index ff457604bb..718aaf7721 100644
--- a/tests/petsc/36.cc
+++ b/tests/petsc/36.cc
@@ -36,8 +36,8 @@ void test (PETScWrappers::Vector &v,
         w(i) = i+1.;
     }
 
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
 
   v -= w;
 
diff --git a/tests/petsc/37.cc b/tests/petsc/37.cc
index 91d11537e4..bf8522ae1e 100644
--- a/tests/petsc/37.cc
+++ b/tests/petsc/37.cc
@@ -29,7 +29,7 @@ void test (PETScWrappers::Vector &v)
   for (unsigned int i=0; i<v.size(); ++i)
     v(i) = i;
 
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   v.add (1.);
 
diff --git a/tests/petsc/38.cc b/tests/petsc/38.cc
index dc24186477..2817de384a 100644
--- a/tests/petsc/38.cc
+++ b/tests/petsc/38.cc
@@ -33,8 +33,8 @@ void test (PETScWrappers::Vector &v,
       w(i) = i+1.;
     }
 
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
 
   v.add (w);
 
diff --git a/tests/petsc/39.cc b/tests/petsc/39.cc
index 13ba1db9ed..3f85e85ab0 100644
--- a/tests/petsc/39.cc
+++ b/tests/petsc/39.cc
@@ -33,8 +33,8 @@ void test (PETScWrappers::Vector &v,
       w(i) = i+1.;
     }
 
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
 
   v.add (2, w);
 
diff --git a/tests/petsc/40.cc b/tests/petsc/40.cc
index 483edb8688..3918118b0e 100644
--- a/tests/petsc/40.cc
+++ b/tests/petsc/40.cc
@@ -35,9 +35,9 @@ void test (PETScWrappers::Vector &v,
       x(i) = i+2.;
     }
 
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
-  x.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
+  x.compress (VectorOperation::insert);
 
   v.add (2, w, 3, x);
 
diff --git a/tests/petsc/41.cc b/tests/petsc/41.cc
index 32d369d7ae..5c9eac41b1 100644
--- a/tests/petsc/41.cc
+++ b/tests/petsc/41.cc
@@ -33,8 +33,8 @@ void test (PETScWrappers::Vector &v,
       w(i) = i+1.;
     }
 
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
 
   v.sadd (2, w);
 
diff --git a/tests/petsc/42.cc b/tests/petsc/42.cc
index 319960ff5d..fe35044bb9 100644
--- a/tests/petsc/42.cc
+++ b/tests/petsc/42.cc
@@ -33,8 +33,8 @@ void test (PETScWrappers::Vector &v,
       w(i) = i+1.;
     }
 
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
 
   v.sadd (3, 2, w);
 
diff --git a/tests/petsc/43.cc b/tests/petsc/43.cc
index 7fa4fb03e6..1e137386bc 100644
--- a/tests/petsc/43.cc
+++ b/tests/petsc/43.cc
@@ -35,9 +35,9 @@ void test (PETScWrappers::Vector &v,
       x(i) = i+2.;
     }
 
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
-  x.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
+  x.compress (VectorOperation::insert);
 
   v.sadd (1.5, 2, w, 3, x);
 
diff --git a/tests/petsc/45.cc b/tests/petsc/45.cc
index 2feef69fc0..b5d70c399a 100644
--- a/tests/petsc/45.cc
+++ b/tests/petsc/45.cc
@@ -33,8 +33,8 @@ void test (PETScWrappers::Vector &v,
       w(i) = i+1.;
     }
 
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
 
   v.scale (w);
 
diff --git a/tests/petsc/46.cc b/tests/petsc/46.cc
index 5bf4056acb..a24193216b 100644
--- a/tests/petsc/46.cc
+++ b/tests/petsc/46.cc
@@ -33,8 +33,8 @@ void test (PETScWrappers::Vector &v,
       w(i) = i+1.;
     }
 
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
 
   v.equ (2, w);
 
diff --git a/tests/petsc/47.cc b/tests/petsc/47.cc
index 0763eaa0e8..c9ccac82b1 100644
--- a/tests/petsc/47.cc
+++ b/tests/petsc/47.cc
@@ -35,9 +35,9 @@ void test (PETScWrappers::Vector &v,
       x(i) = i+2.;
     }
 
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
-  x.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
+  x.compress (VectorOperation::insert);
 
   v.equ (2, w, 3, x);
 
diff --git a/tests/petsc/48.cc b/tests/petsc/48.cc
index eea0b002a9..97e658e7b6 100644
--- a/tests/petsc/48.cc
+++ b/tests/petsc/48.cc
@@ -35,9 +35,9 @@ void test (PETScWrappers::Vector &v,
       x(i) = i+2.;
     }
 
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
-  x.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
+  x.compress (VectorOperation::insert);
 
   v.ratio (w, x);
 
diff --git a/tests/petsc/51.cc b/tests/petsc/51.cc
index eb98b24895..2eb2ccea39 100644
--- a/tests/petsc/51.cc
+++ b/tests/petsc/51.cc
@@ -30,11 +30,10 @@ void test (PETScWrappers::Vector &v)
   for (unsigned int i=0; i<v.size(); ++i)
     if (i%3 == 0)
       v(i) = i+1.;
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // then copy it
   PETScWrappers::Vector w (v);
-  w.compress (VectorOperation::add);
 
   // make sure they're equal
   deallog << v *w << ' ' << v.l2_norm() * w.l2_norm()
diff --git a/tests/petsc/52.cc b/tests/petsc/52.cc
index 591944405a..1491e969b7 100644
--- a/tests/petsc/52.cc
+++ b/tests/petsc/52.cc
@@ -33,7 +33,7 @@ void test (PETScWrappers::SparseMatrix &m)
       if ((i+2*j+1) % 3 == 0)
         m.set (i,j, i*j*.5+.5);
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   // then make sure we retrieve the same ones
   for (unsigned int i=0; i<m.m(); ++i)
diff --git a/tests/petsc/55.cc b/tests/petsc/55.cc
index 0013f4da53..7949d49cf7 100644
--- a/tests/petsc/55.cc
+++ b/tests/petsc/55.cc
@@ -38,7 +38,7 @@ void test (PETScWrappers::Vector &v)
   for (unsigned int i=0; i<v.size(); i+=1+i)
     v(i) *= 2;
 
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // check that they are ok, and this time
   // all of them
diff --git a/tests/petsc/56.cc b/tests/petsc/56.cc
index 34b6eda159..7bd342ad93 100644
--- a/tests/petsc/56.cc
+++ b/tests/petsc/56.cc
@@ -38,7 +38,7 @@ void test (PETScWrappers::Vector &v)
   for (unsigned int i=0; i<v.size(); i+=1+i)
     v(i) /= 2;
 
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // check that they are ok, and this time
   // all of them
diff --git a/tests/petsc/64.cc b/tests/petsc/64.cc
index 0f17ad5f6d..59e804ea68 100644
--- a/tests/petsc/64.cc
+++ b/tests/petsc/64.cc
@@ -34,8 +34,8 @@ template<typename MatrixType>
 void test (MatrixType &m)
 {
   m.add(0,0,1);
+  m.compress(VectorOperation::add);
   m = 0;
-  m.compress();
 
   Assert(fabs(m.frobenius_norm())<1e-15, ExcInternalError());
 
diff --git a/tests/petsc/deal_solver_01.cc b/tests/petsc/deal_solver_01.cc
index 51fe084280..c9ad7e4c87 100644
--- a/tests/petsc/deal_solver_01.cc
+++ b/tests/petsc/deal_solver_01.cc
@@ -83,13 +83,11 @@ int main(int argc, char **argv)
     FDMatrix testproblem(size, size);
     PETScWrappers::SparseMatrix  A(dim, dim, 5);
     testproblem.five_point(A);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
 
     GrowingVectorMemory<PETScWrappers::Vector> mem;
     SolverCG<PETScWrappers::Vector> solver(control,mem);
diff --git a/tests/petsc/deal_solver_02.cc b/tests/petsc/deal_solver_02.cc
index 0ef008f3e7..07ca880f34 100644
--- a/tests/petsc/deal_solver_02.cc
+++ b/tests/petsc/deal_solver_02.cc
@@ -84,13 +84,11 @@ int main(int argc, char **argv)
     FDMatrix testproblem(size, size);
     PETScWrappers::SparseMatrix  A(dim, dim, 5);
     testproblem.five_point(A);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
 
     GrowingVectorMemory<PETScWrappers::Vector> mem;
     SolverBicgstab<PETScWrappers::Vector> solver(control,mem);
diff --git a/tests/petsc/deal_solver_03.cc b/tests/petsc/deal_solver_03.cc
index 3edb73d612..27daaebc5f 100644
--- a/tests/petsc/deal_solver_03.cc
+++ b/tests/petsc/deal_solver_03.cc
@@ -89,9 +89,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     GrowingVectorMemory<PETScWrappers::Vector> mem;
     SolverGMRES<PETScWrappers::Vector> solver(control,mem);
diff --git a/tests/petsc/deal_solver_04.cc b/tests/petsc/deal_solver_04.cc
index d3f932664b..99c19b0b6e 100644
--- a/tests/petsc/deal_solver_04.cc
+++ b/tests/petsc/deal_solver_04.cc
@@ -88,9 +88,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     GrowingVectorMemory<PETScWrappers::Vector> mem;
     SolverMinRes<PETScWrappers::Vector> solver(control,mem);
diff --git a/tests/petsc/deal_solver_05.cc b/tests/petsc/deal_solver_05.cc
index ce225df08d..72bd9bdade 100644
--- a/tests/petsc/deal_solver_05.cc
+++ b/tests/petsc/deal_solver_05.cc
@@ -87,9 +87,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     GrowingVectorMemory<PETScWrappers::Vector> mem;
     SolverQMRS<PETScWrappers::Vector> solver(control,mem);
diff --git a/tests/petsc/full_matrix_00.cc b/tests/petsc/full_matrix_00.cc
index 824d1fbe6a..c9914e2127 100644
--- a/tests/petsc/full_matrix_00.cc
+++ b/tests/petsc/full_matrix_00.cc
@@ -34,7 +34,7 @@ void test (PETScWrappers::FullMatrix &m)
     for (unsigned int j=0; j<m.n(); ++j)
       m.set (i, j, i+2*j);
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   // things we know
   Assert (m.m() == 3, ExcInternalError());
@@ -57,7 +57,7 @@ void test (PETScWrappers::FullMatrix &m)
   for (unsigned int i=0; i<m.m(); ++i)
     for (unsigned int j=0; j<m.n(); ++j)
       m.set (i, j, j+2*i);
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
   
   // things we know
   Assert (m.m() == 5, ExcInternalError());
diff --git a/tests/petsc/full_matrix_01.cc b/tests/petsc/full_matrix_01.cc
index a19fb62e18..f3975ca740 100644
--- a/tests/petsc/full_matrix_01.cc
+++ b/tests/petsc/full_matrix_01.cc
@@ -32,7 +32,7 @@ void test (PETScWrappers::FullMatrix &m)
       if ((i+2*j+1) % 3 == 0)
         m.set (i,j, i*j*.5+.5);
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   // then make sure we retrieve the same ones
   for (unsigned int i=0; i<m.m(); ++i)
diff --git a/tests/petsc/full_matrix_05.cc b/tests/petsc/full_matrix_05.cc
index 776152ba56..d3c45b6947 100644
--- a/tests/petsc/full_matrix_05.cc
+++ b/tests/petsc/full_matrix_05.cc
@@ -37,7 +37,7 @@ void test (PETScWrappers::FullMatrix &m)
           ++counter;
         }
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   // check n_nonzero_elements then
   // output a dummy number.
diff --git a/tests/petsc/full_matrix_06.cc b/tests/petsc/full_matrix_06.cc
index a7b4a2b3fc..76bc6c60d3 100644
--- a/tests/petsc/full_matrix_06.cc
+++ b/tests/petsc/full_matrix_06.cc
@@ -32,7 +32,7 @@ void test (PETScWrappers::FullMatrix &m)
       if ((i+2*j+1) % 3 == 0)
         m.set (i,j, i*j*.5+.5);
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   // compare against the exact value of the
   // l1-norm (max col-sum)
diff --git a/tests/petsc/full_matrix_08.cc b/tests/petsc/full_matrix_08.cc
index d6f0fdcbe5..3ef0233d99 100644
--- a/tests/petsc/full_matrix_08.cc
+++ b/tests/petsc/full_matrix_08.cc
@@ -37,7 +37,7 @@ void test (PETScWrappers::FullMatrix &m)
         }
   norm = std::sqrt(norm);
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   // compare against the exact value of the
   // l2-norm (max row-sum)
diff --git a/tests/petsc/full_matrix_09.cc b/tests/petsc/full_matrix_09.cc
index 744fd1cc6b..513ad58f09 100644
--- a/tests/petsc/full_matrix_09.cc
+++ b/tests/petsc/full_matrix_09.cc
@@ -31,7 +31,7 @@ void test (PETScWrappers::FullMatrix &m)
       if ((i+2*j+1) % 3 == 0)
         m.set (i,j, i*j*.5+.5);
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   // then multiply everything by 1.25 and
   // make sure we retrieve the values we
diff --git a/tests/petsc/full_matrix_10.cc b/tests/petsc/full_matrix_10.cc
index 2f1088d211..d0b87e0d3d 100644
--- a/tests/petsc/full_matrix_10.cc
+++ b/tests/petsc/full_matrix_10.cc
@@ -31,7 +31,7 @@ void test (PETScWrappers::FullMatrix &m)
       if ((i+2*j+1) % 3 == 0)
         m.set (i,j, i*j*.5+.5);
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   // then divide everything by 4/3 and
   // make sure we retrieve the values we
diff --git a/tests/petsc/full_matrix_iterator_01.cc b/tests/petsc/full_matrix_iterator_01.cc
index 21fdf5414e..108d4c13f4 100644
--- a/tests/petsc/full_matrix_iterator_01.cc
+++ b/tests/petsc/full_matrix_iterator_01.cc
@@ -29,7 +29,7 @@ void test ()
   m.set (0,0,1);
   m.set (1,1,2);
   m.set (1,2,3);
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
   PETScWrappers::FullMatrix::const_iterator i = m.begin();
   deallog << i->row() << ' ' << i->column() << ' ' << i->value() << std::endl;
   ++i;
diff --git a/tests/petsc/full_matrix_vector_01.cc b/tests/petsc/full_matrix_vector_01.cc
index 9685ee4382..ef9a9a8304 100644
--- a/tests/petsc/full_matrix_vector_01.cc
+++ b/tests/petsc/full_matrix_vector_01.cc
@@ -36,9 +36,8 @@ void test (PETScWrappers::Vector &v,
   for (unsigned int i=0; i<v.size(); ++i)
     v(i) = i;
 
-  m.compress (VectorOperation::add);
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
+  v.compress (VectorOperation::insert);
 
   // w:=Mv
   m.vmult (w,v);
diff --git a/tests/petsc/full_matrix_vector_02.cc b/tests/petsc/full_matrix_vector_02.cc
index 12d2454b6d..b63166e148 100644
--- a/tests/petsc/full_matrix_vector_02.cc
+++ b/tests/petsc/full_matrix_vector_02.cc
@@ -36,9 +36,8 @@ void test (PETScWrappers::Vector &v,
   for (unsigned int i=0; i<v.size(); ++i)
     v(i) = i;
 
-  m.compress (VectorOperation::add);
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
+  v.compress (VectorOperation::insert);
 
   // w:=Mv
   m.Tvmult (w,v);
diff --git a/tests/petsc/full_matrix_vector_03.cc b/tests/petsc/full_matrix_vector_03.cc
index a9c1519431..5abc75fa6d 100644
--- a/tests/petsc/full_matrix_vector_03.cc
+++ b/tests/petsc/full_matrix_vector_03.cc
@@ -39,9 +39,9 @@ void test (PETScWrappers::Vector &v,
       w(i) = i;
     }
 
-  m.compress (VectorOperation::add);
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
 
   // w:=Mv
   m.vmult_add (w,v);
diff --git a/tests/petsc/full_matrix_vector_05.cc b/tests/petsc/full_matrix_vector_05.cc
index d5dd8bc86f..27982b5dad 100644
--- a/tests/petsc/full_matrix_vector_05.cc
+++ b/tests/petsc/full_matrix_vector_05.cc
@@ -39,9 +39,9 @@ void test (PETScWrappers::Vector &v,
       w(i) = i+1;
     }
 
-  m.compress (VectorOperation::add);
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
 
   // <w,Mv>
   const PetscScalar s = m.matrix_scalar_product (w,v);
diff --git a/tests/petsc/full_matrix_vector_06.cc b/tests/petsc/full_matrix_vector_06.cc
index a93f05b9d0..5a70894fcc 100644
--- a/tests/petsc/full_matrix_vector_06.cc
+++ b/tests/petsc/full_matrix_vector_06.cc
@@ -35,8 +35,8 @@ void test (PETScWrappers::Vector &v)
   for (unsigned int i=0; i<v.size(); ++i)
     v(i) = i;
 
-  m.compress (VectorOperation::add);
-  v.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
+  v.compress (VectorOperation::insert);
 
   // <w,Mv>
   const PetscScalar s = m.matrix_norm_square (v);
diff --git a/tests/petsc/full_matrix_vector_07.cc b/tests/petsc/full_matrix_vector_07.cc
index a8226aa8b0..001829635e 100644
--- a/tests/petsc/full_matrix_vector_07.cc
+++ b/tests/petsc/full_matrix_vector_07.cc
@@ -40,9 +40,9 @@ void test (PETScWrappers::Vector &v,
       w(i) = i+1;
     }
 
-  m.compress (VectorOperation::add);
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
 
   // x=w-Mv
   const double s = m.residual (x, v, w);
diff --git a/tests/petsc/solver_01.cc b/tests/petsc/solver_01.cc
index 0294903e55..0a48809690 100644
--- a/tests/petsc/solver_01.cc
+++ b/tests/petsc/solver_01.cc
@@ -81,9 +81,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverRichardson solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
diff --git a/tests/petsc/solver_02.cc b/tests/petsc/solver_02.cc
index 372fe4acc7..58de09df2e 100644
--- a/tests/petsc/solver_02.cc
+++ b/tests/petsc/solver_02.cc
@@ -84,9 +84,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverChebychev solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
diff --git a/tests/petsc/solver_03.cc b/tests/petsc/solver_03.cc
index c8ba190a51..4c839e3c80 100644
--- a/tests/petsc/solver_03.cc
+++ b/tests/petsc/solver_03.cc
@@ -82,9 +82,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverCG solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
diff --git a/tests/petsc/solver_03_mf.cc b/tests/petsc/solver_03_mf.cc
index ad3b8eae4c..f651d2bc45 100644
--- a/tests/petsc/solver_03_mf.cc
+++ b/tests/petsc/solver_03_mf.cc
@@ -80,9 +80,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverCG solver(control);
     PETScWrappers::PreconditionNone preconditioner(A);
diff --git a/tests/petsc/solver_03_precondition_boomeramg.cc b/tests/petsc/solver_03_precondition_boomeramg.cc
index 3952dd18c1..404b096793 100644
--- a/tests/petsc/solver_03_precondition_boomeramg.cc
+++ b/tests/petsc/solver_03_precondition_boomeramg.cc
@@ -83,9 +83,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverCG solver(control);
     PETScWrappers::PreconditionBoomerAMG preconditioner(A);
diff --git a/tests/petsc/solver_03_precondition_boomeramg_symmetric.cc b/tests/petsc/solver_03_precondition_boomeramg_symmetric.cc
index c301083b4a..ce04b203bd 100644
--- a/tests/petsc/solver_03_precondition_boomeramg_symmetric.cc
+++ b/tests/petsc/solver_03_precondition_boomeramg_symmetric.cc
@@ -84,9 +84,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverCG solver(control);
     PETScWrappers::PreconditionBoomerAMG preconditioner(A, true);
diff --git a/tests/petsc/solver_03_precondition_eisenstat.cc b/tests/petsc/solver_03_precondition_eisenstat.cc
index 8528015206..5747cf0e35 100644
--- a/tests/petsc/solver_03_precondition_eisenstat.cc
+++ b/tests/petsc/solver_03_precondition_eisenstat.cc
@@ -82,9 +82,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverCG solver(control);
     PETScWrappers::PreconditionEisenstat preconditioner(A);
diff --git a/tests/petsc/solver_03_precondition_icc.cc b/tests/petsc/solver_03_precondition_icc.cc
index 6f3fa7708a..7a04bc7770 100644
--- a/tests/petsc/solver_03_precondition_icc.cc
+++ b/tests/petsc/solver_03_precondition_icc.cc
@@ -82,9 +82,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverCG solver(control);
     PETScWrappers::PreconditionICC preconditioner(A);
diff --git a/tests/petsc/solver_03_precondition_lu.cc b/tests/petsc/solver_03_precondition_lu.cc
index 3bb0a590b4..b7b88a0f37 100644
--- a/tests/petsc/solver_03_precondition_lu.cc
+++ b/tests/petsc/solver_03_precondition_lu.cc
@@ -89,9 +89,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverCG solver(control);
     PETScWrappers::PreconditionLU preconditioner(A);
diff --git a/tests/petsc/solver_03_precondition_parasails.cc b/tests/petsc/solver_03_precondition_parasails.cc
index def46b1520..269fdaae8f 100644
--- a/tests/petsc/solver_03_precondition_parasails.cc
+++ b/tests/petsc/solver_03_precondition_parasails.cc
@@ -82,9 +82,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverCG solver(control);
     PETScWrappers::PreconditionParaSails preconditioner(A);
diff --git a/tests/petsc/solver_03_precondition_sor.cc b/tests/petsc/solver_03_precondition_sor.cc
index ceab998ba0..91edc6de6d 100644
--- a/tests/petsc/solver_03_precondition_sor.cc
+++ b/tests/petsc/solver_03_precondition_sor.cc
@@ -82,9 +82,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverCG solver(control);
     PETScWrappers::PreconditionSOR preconditioner(A);
diff --git a/tests/petsc/solver_03_precondition_ssor.cc b/tests/petsc/solver_03_precondition_ssor.cc
index 5f47ab3d90..9856f14bb7 100644
--- a/tests/petsc/solver_03_precondition_ssor.cc
+++ b/tests/petsc/solver_03_precondition_ssor.cc
@@ -82,9 +82,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverCG solver(control);
     PETScWrappers::PreconditionSSOR preconditioner(A);
diff --git a/tests/petsc/solver_04.cc b/tests/petsc/solver_04.cc
index 895845c74e..0f57a698b6 100644
--- a/tests/petsc/solver_04.cc
+++ b/tests/petsc/solver_04.cc
@@ -82,9 +82,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverBiCG solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
diff --git a/tests/petsc/solver_05.cc b/tests/petsc/solver_05.cc
index c88c9e4786..d3bcd3d20a 100644
--- a/tests/petsc/solver_05.cc
+++ b/tests/petsc/solver_05.cc
@@ -82,9 +82,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverGMRES solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
diff --git a/tests/petsc/solver_06.cc b/tests/petsc/solver_06.cc
index 1a81d8e151..95c65ac3af 100644
--- a/tests/petsc/solver_06.cc
+++ b/tests/petsc/solver_06.cc
@@ -85,9 +85,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverBicgstab solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
diff --git a/tests/petsc/solver_07.cc b/tests/petsc/solver_07.cc
index 6d677ade05..0a568fc940 100644
--- a/tests/petsc/solver_07.cc
+++ b/tests/petsc/solver_07.cc
@@ -82,9 +82,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverCGS solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
diff --git a/tests/petsc/solver_08.cc b/tests/petsc/solver_08.cc
index 4b749fd338..87286242f8 100644
--- a/tests/petsc/solver_08.cc
+++ b/tests/petsc/solver_08.cc
@@ -82,9 +82,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverTFQMR solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
diff --git a/tests/petsc/solver_09.cc b/tests/petsc/solver_09.cc
index 08520c11ba..01e44da878 100644
--- a/tests/petsc/solver_09.cc
+++ b/tests/petsc/solver_09.cc
@@ -86,9 +86,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverTCQMR solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
diff --git a/tests/petsc/solver_10.cc b/tests/petsc/solver_10.cc
index 5fc7cd7efc..3197a839fe 100644
--- a/tests/petsc/solver_10.cc
+++ b/tests/petsc/solver_10.cc
@@ -82,9 +82,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverCR solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
diff --git a/tests/petsc/solver_11.cc b/tests/petsc/solver_11.cc
index c22613b528..5f8a6c81d9 100644
--- a/tests/petsc/solver_11.cc
+++ b/tests/petsc/solver_11.cc
@@ -86,9 +86,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverLSQR solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
diff --git a/tests/petsc/solver_12.cc b/tests/petsc/solver_12.cc
index b0c4032d39..79a6e9d43e 100644
--- a/tests/petsc/solver_12.cc
+++ b/tests/petsc/solver_12.cc
@@ -75,9 +75,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverPreOnly solver(control);
 
diff --git a/tests/petsc/solver_13.cc b/tests/petsc/solver_13.cc
index 50909db7cd..5deb2efdd8 100644
--- a/tests/petsc/solver_13.cc
+++ b/tests/petsc/solver_13.cc
@@ -73,9 +73,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  u(dim);
     u = 0.;
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
     PETScWrappers::Vector t=u;
 
     PETScWrappers::SolverPreOnly solver(control);
diff --git a/tests/petsc/sparse_direct_mumps.cc b/tests/petsc/sparse_direct_mumps.cc
index 4670fb127e..3432a94fc6 100644
--- a/tests/petsc/sparse_direct_mumps.cc
+++ b/tests/petsc/sparse_direct_mumps.cc
@@ -57,9 +57,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  u(dim);
     u = 0.;
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     SolverControl cn;
     PETScWrappers::SparseDirectMUMPS solver(cn);
diff --git a/tests/petsc/sparse_matrix_01.cc b/tests/petsc/sparse_matrix_01.cc
index 3d39fb34b2..41b4d0d97a 100644
--- a/tests/petsc/sparse_matrix_01.cc
+++ b/tests/petsc/sparse_matrix_01.cc
@@ -33,13 +33,13 @@ void test ()
     for (unsigned int j=0; j<=i; ++j)
       m.set (i,j, i+2*j);
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   PETScWrappers::SparseMatrix m2(s,s,s);
   m2.set(0,1,5.0);
   for (unsigned int i=0; i<m2.m(); ++i)
     m2.set(i,i,1.0+ i);
-  m2.compress (VectorOperation::add);
+  m2.compress (VectorOperation::insert);
 
   // we now print the matrix m,
   // print after adding m2, and then subtract m2 again
diff --git a/tests/petsc/sparse_matrix_02.cc b/tests/petsc/sparse_matrix_02.cc
index 92329cfc3e..acab595b6f 100644
--- a/tests/petsc/sparse_matrix_02.cc
+++ b/tests/petsc/sparse_matrix_02.cc
@@ -34,7 +34,7 @@ void test ()
       m.set (i,j, i+2*j);
 
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   {
 
diff --git a/tests/petsc/sparse_matrix_iterator_01.cc b/tests/petsc/sparse_matrix_iterator_01.cc
index e9091f7e8c..740bc0d846 100644
--- a/tests/petsc/sparse_matrix_iterator_01.cc
+++ b/tests/petsc/sparse_matrix_iterator_01.cc
@@ -29,7 +29,7 @@ void test ()
   m.set (0,0,1);
   m.set (1,1,2);
   m.set (1,2,3);
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
   PETScWrappers::SparseMatrix::const_iterator i = m.begin();
   deallog << i->row() << ' ' << i->column() << ' ' << i->value() << std::endl;
   ++i;
diff --git a/tests/petsc/sparse_matrix_vector_01.cc b/tests/petsc/sparse_matrix_vector_01.cc
index e510513ac2..d54aa90c13 100644
--- a/tests/petsc/sparse_matrix_vector_01.cc
+++ b/tests/petsc/sparse_matrix_vector_01.cc
@@ -36,9 +36,8 @@ void test (PETScWrappers::Vector &v,
   for (unsigned int i=0; i<v.size(); ++i)
     v(i) = i;
 
-  m.compress (VectorOperation::add);
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
+  v.compress (VectorOperation::insert);
 
   // w:=Mv
   m.vmult (w,v);
diff --git a/tests/petsc/sparse_matrix_vector_02.cc b/tests/petsc/sparse_matrix_vector_02.cc
index bc98ad15b8..39d9cc5744 100644
--- a/tests/petsc/sparse_matrix_vector_02.cc
+++ b/tests/petsc/sparse_matrix_vector_02.cc
@@ -36,9 +36,8 @@ void test (PETScWrappers::Vector &v,
   for (unsigned int i=0; i<v.size(); ++i)
     v(i) = i;
 
-  m.compress (VectorOperation::add);
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
+  v.compress (VectorOperation::insert);
 
   // w:=Mv
   m.Tvmult (w,v);
diff --git a/tests/petsc/sparse_matrix_vector_03.cc b/tests/petsc/sparse_matrix_vector_03.cc
index 9f2df38acd..ec264536ba 100644
--- a/tests/petsc/sparse_matrix_vector_03.cc
+++ b/tests/petsc/sparse_matrix_vector_03.cc
@@ -39,9 +39,9 @@ void test (PETScWrappers::Vector &v,
       w(i) = i;
     }
 
-  m.compress (VectorOperation::add);
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
 
   // w:=Mv
   m.vmult_add (w,v);
diff --git a/tests/petsc/sparse_matrix_vector_04.cc b/tests/petsc/sparse_matrix_vector_04.cc
index 287cdd6a1e..ce746f37eb 100644
--- a/tests/petsc/sparse_matrix_vector_04.cc
+++ b/tests/petsc/sparse_matrix_vector_04.cc
@@ -39,9 +39,9 @@ void test (PETScWrappers::Vector &v,
       w(i) = i;
     }
 
-  m.compress (VectorOperation::add);
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
 
   // w:=Mv
   m.Tvmult_add (w,v);
diff --git a/tests/petsc/sparse_matrix_vector_05.cc b/tests/petsc/sparse_matrix_vector_05.cc
index 49fa71e3ac..d08b2642f1 100644
--- a/tests/petsc/sparse_matrix_vector_05.cc
+++ b/tests/petsc/sparse_matrix_vector_05.cc
@@ -39,9 +39,9 @@ void test (PETScWrappers::Vector &v,
       w(i) = i+1;
     }
 
-  m.compress (VectorOperation::add);
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
 
   // <w,Mv>
   const PetscScalar s = m.matrix_scalar_product (w,v);
diff --git a/tests/petsc/sparse_matrix_vector_06.cc b/tests/petsc/sparse_matrix_vector_06.cc
index d589fe11fd..a84f749f09 100644
--- a/tests/petsc/sparse_matrix_vector_06.cc
+++ b/tests/petsc/sparse_matrix_vector_06.cc
@@ -35,8 +35,8 @@ void test (PETScWrappers::Vector &v)
   for (unsigned int i=0; i<v.size(); ++i)
     v(i) = i;
 
-  m.compress (VectorOperation::add);
-  v.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
+  v.compress (VectorOperation::insert);
 
   // <w,Mv>
   const PetscScalar s = m.matrix_norm_square (v);
diff --git a/tests/petsc/sparse_matrix_vector_07.cc b/tests/petsc/sparse_matrix_vector_07.cc
index 26041ff6bc..35fc56c870 100644
--- a/tests/petsc/sparse_matrix_vector_07.cc
+++ b/tests/petsc/sparse_matrix_vector_07.cc
@@ -40,9 +40,9 @@ void test (PETScWrappers::Vector &v,
       w(i) = i+1;
     }
 
-  m.compress (VectorOperation::add);
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
 
   // x=w-Mv
   const double s = m.residual (x, v, w);
-- 
2.39.5