From bf2ed016a66d6daa7c70e5ade8deedd638ca84c7 Mon Sep 17 00:00:00 2001
From: "Toby D. Young" <tyoung@ippt.pan.pl>
Date: Sat, 9 Nov 2013 15:38:36 +0000
Subject: [PATCH] Explicitly compress petsc objects with a VectorOperation.

git-svn-id: https://svn.dealii.org/trunk@31602 0785d39b-7218-0410-832d-ea1e28bc413d
---
 tests/petsc/01.cc                                         | 2 +-
 tests/petsc/02.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/13.cc                                         | 2 +-
 tests/petsc/17.cc                                         | 2 +-
 tests/petsc/18.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/30.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/44.cc                                         | 8 ++++----
 tests/petsc/45.cc                                         | 4 ++--
 tests/petsc/46.cc                                         | 4 ++--
 tests/petsc/47.cc                                         | 6 +++---
 tests/petsc/48.cc                                         | 6 +++---
 tests/petsc/51.cc                                         | 4 ++--
 tests/petsc/52.cc                                         | 2 +-
 tests/petsc/55.cc                                         | 2 +-
 tests/petsc/56.cc                                         | 2 +-
 tests/petsc/57.cc                                         | 2 +-
 tests/petsc/58.cc                                         | 2 +-
 tests/petsc/59.cc                                         | 2 +-
 tests/petsc/60.cc                                         | 2 +-
 tests/petsc/61.cc                                         | 2 +-
 tests/petsc/66.cc                                         | 2 +-
 tests/petsc/67.cc                                         | 2 +-
 tests/petsc/68.cc                                         | 2 +-
 tests/petsc/69.cc                                         | 2 +-
 tests/petsc/deal_solver_01.cc                             | 6 +++---
 tests/petsc/deal_solver_02.cc                             | 6 +++---
 tests/petsc/deal_solver_03.cc                             | 6 +++---
 tests/petsc/deal_solver_04.cc                             | 6 +++---
 tests/petsc/deal_solver_05.cc                             | 6 +++---
 tests/petsc/full_matrix_01.cc                             | 2 +-
 tests/petsc/full_matrix_02.cc                             | 2 +-
 tests/petsc/full_matrix_05.cc                             | 2 +-
 tests/petsc/full_matrix_06.cc                             | 2 +-
 tests/petsc/full_matrix_07.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                      | 6 +++---
 tests/petsc/full_matrix_vector_02.cc                      | 6 +++---
 tests/petsc/full_matrix_vector_03.cc                      | 6 +++---
 tests/petsc/full_matrix_vector_04.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/parallel_sparse_matrix_01.cc                  | 2 +-
 tests/petsc/slowness_01.cc                                | 2 +-
 tests/petsc/slowness_02.cc                                | 2 +-
 tests/petsc/slowness_03.cc                                | 2 +-
 tests/petsc/slowness_04.cc                                | 2 +-
 tests/petsc/solver_01.cc                                  | 6 +++---
 tests/petsc/solver_02.cc                                  | 6 +++---
 tests/petsc/solver_03.cc                                  | 6 +++---
 tests/petsc/solver_03_mf.cc                               | 6 +++---
 tests/petsc/solver_03_precondition_boomeramg.cc           | 6 +++---
 tests/petsc/solver_03_precondition_boomeramg_symmetric.cc | 6 +++---
 tests/petsc/solver_03_precondition_eisenstat.cc           | 6 +++---
 tests/petsc/solver_03_precondition_icc.cc                 | 6 +++---
 tests/petsc/solver_03_precondition_ilu.cc                 | 6 +++---
 tests/petsc/solver_03_precondition_lu.cc                  | 6 +++---
 tests/petsc/solver_03_precondition_parasails.cc           | 6 +++---
 tests/petsc/solver_03_precondition_sor.cc                 | 6 +++---
 tests/petsc/solver_03_precondition_ssor.cc                | 6 +++---
 tests/petsc/solver_04.cc                                  | 6 +++---
 tests/petsc/solver_05.cc                                  | 6 +++---
 tests/petsc/solver_06.cc                                  | 6 +++---
 tests/petsc/solver_07.cc                                  | 6 +++---
 tests/petsc/solver_08.cc                                  | 6 +++---
 tests/petsc/solver_09.cc                                  | 6 +++---
 tests/petsc/solver_10.cc                                  | 6 +++---
 tests/petsc/solver_11.cc                                  | 6 +++---
 tests/petsc/solver_12.cc                                  | 6 +++---
 tests/petsc/solver_13.cc                                  | 6 +++---
 tests/petsc/sparse_direct_mumps.cc                        | 6 +++---
 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                    | 6 +++---
 tests/petsc/sparse_matrix_vector_02.cc                    | 6 +++---
 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 +++---
 115 files changed, 222 insertions(+), 222 deletions(-)

diff --git a/tests/petsc/01.cc b/tests/petsc/01.cc
index 5f3cc2ad78..bd4a692462 100644
--- a/tests/petsc/01.cc
+++ b/tests/petsc/01.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 ();
+  m.compress (VectorOperation::add);
 
   // then make sure we retrieve the same ones
   for (unsigned int i=0; i<m.m(); ++i)
diff --git a/tests/petsc/02.cc b/tests/petsc/02.cc
index 2dbcbcdb3b..e822870680 100644
--- a/tests/petsc/02.cc
+++ b/tests/petsc/02.cc
@@ -33,7 +33,7 @@ void test (PETScWrappers::SparseMatrix &m)
       if ((i+2*j+1) % 3 == 0)
         m.add (i,j, i*j*.5+.5);
 
-  m.compress ();
+  m.compress (VectorOperation::add);
 
   // 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 3991d6ecb0..da9d16bc4d 100644
--- a/tests/petsc/05.cc
+++ b/tests/petsc/05.cc
@@ -38,7 +38,7 @@ void test (PETScWrappers::SparseMatrix &m)
           ++counter;
         }
 
-  m.compress ();
+  m.compress VectorOperation::add();
 
   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 3cf3ad8b11..077103411b 100644
--- a/tests/petsc/06.cc
+++ b/tests/petsc/06.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 ();
+  m.compress (VectorOperation::add);
 
   // 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 c0bbc2a97b..e9f7649769 100644
--- a/tests/petsc/07.cc
+++ b/tests/petsc/07.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 ();
+  m.compress (VectorOperation::add);
 
   // 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 5a637d9bb3..540d479b7d 100644
--- a/tests/petsc/08.cc
+++ b/tests/petsc/08.cc
@@ -38,7 +38,7 @@ void test (PETScWrappers::SparseMatrix &m)
         }
   norm = std::sqrt(norm);
 
-  m.compress ();
+  m.compress (VectorOperation::add);
 
   // 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 15050cd795..e7e32c3418 100644
--- a/tests/petsc/09.cc
+++ b/tests/petsc/09.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 ();
+  m.compress (VectorOperation::add);
 
   // 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 ee2dab657f..3c1b7556b4 100644
--- a/tests/petsc/10.cc
+++ b/tests/petsc/10.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 ();
+  m.compress (VectorOperation::add);
 
   // 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 9f99bdba54..cc7213f923 100644
--- a/tests/petsc/11.cc
+++ b/tests/petsc/11.cc
@@ -30,7 +30,7 @@ void test (PETScWrappers::Vector &v)
   for (unsigned int i=0; i<v.size(); i+=1+i)
     v(i) = i;
 
-  v.compress ();
+  v.compress (VectorOperation::add);
 
   Assert (v.size() == 100, ExcInternalError());
 
diff --git a/tests/petsc/12.cc b/tests/petsc/12.cc
index a9954b0a8b..4e38d87fac 100644
--- a/tests/petsc/12.cc
+++ b/tests/petsc/12.cc
@@ -37,7 +37,7 @@ void test (PETScWrappers::Vector &v)
       pattern[i] = true;
     }
 
-  v.compress ();
+  v.compress (VectorOperation::add);
 
   // check that they are ok, and this time
   // all of them
diff --git a/tests/petsc/13.cc b/tests/petsc/13.cc
index 44be821ae3..ddc9632996 100644
--- a/tests/petsc/13.cc
+++ b/tests/petsc/13.cc
@@ -37,7 +37,7 @@ void test (PETScWrappers::Vector &v)
       pattern[i] = true;
     }
 
-  v.compress ();
+  v.compress (VectorOperation::add);
 
   // check that they are ok, and this time
   // all of them
diff --git a/tests/petsc/17.cc b/tests/petsc/17.cc
index 9501b1f042..cfb680d6cf 100644
--- a/tests/petsc/17.cc
+++ b/tests/petsc/17.cc
@@ -34,7 +34,7 @@ void test (PETScWrappers::Vector &v)
       v(i) = i;
       norm += std::fabs(1.*i);
     }
-  v.compress ();
+  v.compress (VectorOperation::add);
 
   // then check the norm
   Assert (v.l1_norm() == norm, ExcInternalError());
diff --git a/tests/petsc/18.cc b/tests/petsc/18.cc
index 3e3537cf30..6175fc23f1 100644
--- a/tests/petsc/18.cc
+++ b/tests/petsc/18.cc
@@ -34,7 +34,7 @@ void test (PETScWrappers::Vector &v)
       v(i) = i;
       norm += std::fabs(1.*i)*std::fabs(1.*i);
     }
-  v.compress ();
+  v.compress (VectorOperation::add);
 
   // then check the norm
   const double eps=typeid(PetscScalar)==typeid(double) ? 1e-14 : 1e-5;
diff --git a/tests/petsc/19.cc b/tests/petsc/19.cc
index 6de9a94180..7c5e074010 100644
--- a/tests/petsc/19.cc
+++ b/tests/petsc/19.cc
@@ -34,7 +34,7 @@ void test (PETScWrappers::Vector &v)
       v(i) = i;
       norm = std::max(norm,fabs(i));
     }
-  v.compress ();
+  v.compress (VectorOperation::add);
 
   // then check the norm
   Assert (v.linfty_norm() == norm, ExcInternalError());
diff --git a/tests/petsc/20.cc b/tests/petsc/20.cc
index aa303d23af..20bd818bad 100644
--- a/tests/petsc/20.cc
+++ b/tests/petsc/20.cc
@@ -36,7 +36,7 @@ void test (PETScWrappers::Vector &v)
       v(i) = i;
       pattern[i] = true;
     }
-  v.compress ();
+  v.compress (VectorOperation::add);
 
   // multiply v with 5/4
   v *= 5./4.;
diff --git a/tests/petsc/21.cc b/tests/petsc/21.cc
index 4b3602598b..9c0f7c134f 100644
--- a/tests/petsc/21.cc
+++ b/tests/petsc/21.cc
@@ -36,7 +36,7 @@ void test (PETScWrappers::Vector &v)
       v(i) = i;
       pattern[i] = true;
     }
-  v.compress ();
+  v.compress (VectorOperation::add);
 
   // multiply v with 3/4
   v /= 4./3.;
diff --git a/tests/petsc/22.cc b/tests/petsc/22.cc
index 593c88d366..4816c9fc8d 100644
--- a/tests/petsc/22.cc
+++ b/tests/petsc/22.cc
@@ -36,8 +36,8 @@ void test (PETScWrappers::Vector &v,
       v(i) = i;
     else
       w(i) = i;
-  v.compress ();
-  w.compress ();
+  v.compress (VectorOperation::add);
+  w.compress (VectorOperation::add);
 
   // 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 747ba74d3e..e277aae4f9 100644
--- a/tests/petsc/23.cc
+++ b/tests/petsc/23.cc
@@ -43,8 +43,8 @@ void test (PETScWrappers::Vector &v,
         }
     }
 
-  v.compress ();
-  w.compress ();
+  v.compress (VectorOperation::add);
+  w.compress (VectorOperation::add);
 
   // 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 dae86a02b1..837357831b 100644
--- a/tests/petsc/24.cc
+++ b/tests/petsc/24.cc
@@ -33,7 +33,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 ();
+  v.compress (VectorOperation::add);
 
   // 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 82a30e4adb..ee73c55135 100644
--- a/tests/petsc/25.cc
+++ b/tests/petsc/25.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 ();
+  v.compress (VectorOperation::add);
 
   // 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 2a26dfea28..e592f2b41f 100644
--- a/tests/petsc/26.cc
+++ b/tests/petsc/26.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 ();
+  v.compress (VectorOperation::add);
 
   const unsigned int sz = v.size();
   v = 2;
diff --git a/tests/petsc/27.cc b/tests/petsc/27.cc
index ea407ff938..a1bc82eaeb 100644
--- a/tests/petsc/27.cc
+++ b/tests/petsc/27.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 ();
+  v.compress (VectorOperation::add);
 
   // then copy it
   PETScWrappers::Vector w (v.size());
diff --git a/tests/petsc/28.cc b/tests/petsc/28.cc
index 278eb8fdfc..50ba45f847 100644
--- a/tests/petsc/28.cc
+++ b/tests/petsc/28.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 ();
+  v.compress (VectorOperation::add);
 
   // then copy it to a vector of different
   // size
diff --git a/tests/petsc/30.cc b/tests/petsc/30.cc
index 368a4d9194..1bbe5efea4 100644
--- a/tests/petsc/30.cc
+++ b/tests/petsc/30.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 ();
+  v.compress (VectorOperation::add);
 
   // then resize with setting to zero
   v.reinit (13, false);
diff --git a/tests/petsc/31.cc b/tests/petsc/31.cc
index 725303a2d0..276705a643 100644
--- a/tests/petsc/31.cc
+++ b/tests/petsc/31.cc
@@ -34,7 +34,7 @@ void test (PETScWrappers::Vector &v)
       v(i) = i;
       norm += std::fabs(1.*i)*std::fabs(1.*i);
     }
-  v.compress ();
+  v.compress (VectorOperation::add);
 
   // 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 1136971c1f..aada706965 100644
--- a/tests/petsc/32.cc
+++ b/tests/petsc/32.cc
@@ -34,7 +34,7 @@ void test (PETScWrappers::Vector &v)
       v(i) = i;
       sum += i;
     }
-  v.compress ();
+  v.compress (VectorOperation::add);
 
   // 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 7cdf00cdb2..17fae3e660 100644
--- a/tests/petsc/33.cc
+++ b/tests/petsc/33.cc
@@ -34,7 +34,7 @@ void test (PETScWrappers::Vector &v)
       v(i) = i;
       sum += i*i*i;
     }
-  v.compress ();
+  v.compress (VectorOperation::add);
 
   // 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 d2ce01c600..bc5ee0609c 100644
--- a/tests/petsc/34.cc
+++ b/tests/petsc/34.cc
@@ -34,7 +34,7 @@ void test (PETScWrappers::Vector &v)
       v(i) = i;
       sum += i*i*i;
     }
-  v.compress ();
+  v.compress (VectorOperation::add);
 
   // set them to zero again
   v = 0;
diff --git a/tests/petsc/35.cc b/tests/petsc/35.cc
index bf0b6bcd51..8742a8bd64 100644
--- a/tests/petsc/35.cc
+++ b/tests/petsc/35.cc
@@ -37,8 +37,8 @@ void test (PETScWrappers::Vector &v,
         w(i) = i+1.;
     }
 
-  v.compress ();
-  w.compress ();
+  v.compress (VectorOperation::add);
+  w.compress (VectorOperation::add);
 
   v += w;
 
diff --git a/tests/petsc/36.cc b/tests/petsc/36.cc
index a0eeea75ee..522f9cb116 100644
--- a/tests/petsc/36.cc
+++ b/tests/petsc/36.cc
@@ -37,8 +37,8 @@ void test (PETScWrappers::Vector &v,
         w(i) = i+1.;
     }
 
-  v.compress ();
-  w.compress ();
+  v.compress (VectorOperation::add);
+  w.compress (VectorOperation::add);
 
   v -= w;
 
diff --git a/tests/petsc/37.cc b/tests/petsc/37.cc
index ff11c6aabe..f513e665f6 100644
--- a/tests/petsc/37.cc
+++ b/tests/petsc/37.cc
@@ -30,7 +30,7 @@ void test (PETScWrappers::Vector &v)
   for (unsigned int i=0; i<v.size(); ++i)
     v(i) = i;
 
-  v.compress ();
+  v.compress (VectorOperation::add);
 
   v.add (1.);
 
diff --git a/tests/petsc/38.cc b/tests/petsc/38.cc
index aacb319b3d..ba288ac96f 100644
--- a/tests/petsc/38.cc
+++ b/tests/petsc/38.cc
@@ -34,8 +34,8 @@ void test (PETScWrappers::Vector &v,
       w(i) = i+1.;
     }
 
-  v.compress ();
-  w.compress ();
+  v.compress (VectorOperation::add);
+  w.compress (VectorOperation::add);
 
   v.add (w);
 
diff --git a/tests/petsc/39.cc b/tests/petsc/39.cc
index 662e3334f2..bed8660b2b 100644
--- a/tests/petsc/39.cc
+++ b/tests/petsc/39.cc
@@ -34,8 +34,8 @@ void test (PETScWrappers::Vector &v,
       w(i) = i+1.;
     }
 
-  v.compress ();
-  w.compress ();
+  v.compress (VectorOperation::add);
+  w.compress (VectorOperation::add);
 
   v.add (2, w);
 
diff --git a/tests/petsc/40.cc b/tests/petsc/40.cc
index 1ea51675ae..752a2765f8 100644
--- a/tests/petsc/40.cc
+++ b/tests/petsc/40.cc
@@ -36,9 +36,9 @@ void test (PETScWrappers::Vector &v,
       x(i) = i+2.;
     }
 
-  v.compress ();
-  w.compress ();
-  x.compress ();
+  v.compress (VectorOperation::add);
+  w.compress (VectorOperation::add);
+  x.compress (VectorOperation::add);
 
   v.add (2, w, 3, x);
 
diff --git a/tests/petsc/41.cc b/tests/petsc/41.cc
index feca17d765..69c98d621f 100644
--- a/tests/petsc/41.cc
+++ b/tests/petsc/41.cc
@@ -34,8 +34,8 @@ void test (PETScWrappers::Vector &v,
       w(i) = i+1.;
     }
 
-  v.compress ();
-  w.compress ();
+  v.compress (VectorOperation::add);
+  w.compress (VectorOperation::add);
 
   v.sadd (2, w);
 
diff --git a/tests/petsc/42.cc b/tests/petsc/42.cc
index caf0d11fd6..0591ff5fc7 100644
--- a/tests/petsc/42.cc
+++ b/tests/petsc/42.cc
@@ -34,8 +34,8 @@ void test (PETScWrappers::Vector &v,
       w(i) = i+1.;
     }
 
-  v.compress ();
-  w.compress ();
+  v.compress (VectorOperation::add);
+  w.compress (VectorOperation::add);
 
   v.sadd (3, 2, w);
 
diff --git a/tests/petsc/43.cc b/tests/petsc/43.cc
index 5675ffb352..216c4b0c4b 100644
--- a/tests/petsc/43.cc
+++ b/tests/petsc/43.cc
@@ -36,9 +36,9 @@ void test (PETScWrappers::Vector &v,
       x(i) = i+2.;
     }
 
-  v.compress ();
-  w.compress ();
-  x.compress ();
+  v.compress (VectorOperation::add);
+  w.compress (VectorOperation::add);
+  x.compress (VectorOperation::add);
 
   v.sadd (1.5, 2, w, 3, x);
 
diff --git a/tests/petsc/44.cc b/tests/petsc/44.cc
index f53690ca53..22ef05478c 100644
--- a/tests/petsc/44.cc
+++ b/tests/petsc/44.cc
@@ -38,10 +38,10 @@ void test (PETScWrappers::Vector &v,
       y(i) = i+3.;
     }
 
-  v.compress ();
-  w.compress ();
-  x.compress ();
-  y.compress ();
+  v.compress (VectorOperation::add);
+  w.compress (VectorOperation::add);
+  x.compress (VectorOperation::add);
+  y.compress (VectorOperation::add);
 
   v.sadd (1.5, 2, w, 3, x, 4, y);
 
diff --git a/tests/petsc/45.cc b/tests/petsc/45.cc
index c63f0f2787..4562d006e2 100644
--- a/tests/petsc/45.cc
+++ b/tests/petsc/45.cc
@@ -34,8 +34,8 @@ void test (PETScWrappers::Vector &v,
       w(i) = i+1.;
     }
 
-  v.compress ();
-  w.compress ();
+  v.compress (VectorOperation::add);
+  w.compress (VectorOperation::add);
 
   v.scale (w);
 
diff --git a/tests/petsc/46.cc b/tests/petsc/46.cc
index c063b75721..296c5e184d 100644
--- a/tests/petsc/46.cc
+++ b/tests/petsc/46.cc
@@ -34,8 +34,8 @@ void test (PETScWrappers::Vector &v,
       w(i) = i+1.;
     }
 
-  v.compress ();
-  w.compress ();
+  v.compress (VectorOperation::add);
+  w.compress (VectorOperation::add);
 
   v.equ (2, w);
 
diff --git a/tests/petsc/47.cc b/tests/petsc/47.cc
index d3065bf076..feb0c52086 100644
--- a/tests/petsc/47.cc
+++ b/tests/petsc/47.cc
@@ -36,9 +36,9 @@ void test (PETScWrappers::Vector &v,
       x(i) = i+2.;
     }
 
-  v.compress ();
-  w.compress ();
-  x.compress ();
+  v.compress (VectorOperation::add);
+  w.compress (VectorOperation::add);
+  x.compress (VectorOperation::add);
 
   v.equ (2, w, 3, x);
 
diff --git a/tests/petsc/48.cc b/tests/petsc/48.cc
index cb99358075..da45ea1145 100644
--- a/tests/petsc/48.cc
+++ b/tests/petsc/48.cc
@@ -36,9 +36,9 @@ void test (PETScWrappers::Vector &v,
       x(i) = i+2.;
     }
 
-  v.compress ();
-  w.compress ();
-  x.compress ();
+  v.compress (VectorOperation::add);
+  w.compress (VectorOperation::add);
+  x.compress (VectorOperation::add);
 
   v.ratio (w, x);
 
diff --git a/tests/petsc/51.cc b/tests/petsc/51.cc
index f7b826b64f..5cc1155cde 100644
--- a/tests/petsc/51.cc
+++ b/tests/petsc/51.cc
@@ -31,11 +31,11 @@ void test (PETScWrappers::Vector &v)
   for (unsigned int i=0; i<v.size(); ++i)
     if (i%3 == 0)
       v(i) = i+1.;
-  v.compress ();
+  v.compress (VectorOperation::add);
 
   // then copy it
   PETScWrappers::Vector w (v);
-  w.compress ();
+  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 3cf2e2e38d..754e61c7bd 100644
--- a/tests/petsc/52.cc
+++ b/tests/petsc/52.cc
@@ -34,7 +34,7 @@ void test (PETScWrappers::SparseMatrix &m)
       if ((i+2*j+1) % 3 == 0)
         m.set (i,j, i*j*.5+.5);
 
-  m.compress ();
+  m.compress (VectorOperation::add);
 
   // 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 a6dd450067..f50e927cf4 100644
--- a/tests/petsc/55.cc
+++ b/tests/petsc/55.cc
@@ -39,7 +39,7 @@ void test (PETScWrappers::Vector &v)
   for (unsigned int i=0; i<v.size(); i+=1+i)
     v(i) *= 2;
 
-  v.compress ();
+  v.compress (VectorOperation::add);
 
   // check that they are ok, and this time
   // all of them
diff --git a/tests/petsc/56.cc b/tests/petsc/56.cc
index 8183680678..67a49d5b77 100644
--- a/tests/petsc/56.cc
+++ b/tests/petsc/56.cc
@@ -39,7 +39,7 @@ void test (PETScWrappers::Vector &v)
   for (unsigned int i=0; i<v.size(); i+=1+i)
     v(i) /= 2;
 
-  v.compress ();
+  v.compress (VectorOperation::add);
 
   // check that they are ok, and this time
   // all of them
diff --git a/tests/petsc/57.cc b/tests/petsc/57.cc
index 4a20c3bfe8..b4bd1bfb03 100644
--- a/tests/petsc/57.cc
+++ b/tests/petsc/57.cc
@@ -36,7 +36,7 @@ void test (PETScWrappers::Vector &v)
       pattern[i] = true;
     }
 
-  v.compress ();
+  v.compress (VectorOperation::add);
 
   // check that the vector is really
   // non-negative
diff --git a/tests/petsc/58.cc b/tests/petsc/58.cc
index 5b195aeb5b..0dbf559611 100644
--- a/tests/petsc/58.cc
+++ b/tests/petsc/58.cc
@@ -38,7 +38,7 @@ void test (PETScWrappers::Vector &v)
       pattern[i] = true;
     }
 
-  v.compress ();
+  v.compress (VectorOperation::add);
 
   Vector<double> w (v);
   Vector<float>  x (v);
diff --git a/tests/petsc/59.cc b/tests/petsc/59.cc
index 3bb270d0fc..6e2b45c94a 100644
--- a/tests/petsc/59.cc
+++ b/tests/petsc/59.cc
@@ -38,7 +38,7 @@ void test (PETScWrappers::Vector &v)
       pattern[i] = true;
     }
 
-  v.compress ();
+  v.compress (VectorOperation::add);
 
   Vector<double> w (v);
   Vector<float>  x (v);
diff --git a/tests/petsc/60.cc b/tests/petsc/60.cc
index 77cba83030..2c532bc0a3 100644
--- a/tests/petsc/60.cc
+++ b/tests/petsc/60.cc
@@ -38,7 +38,7 @@ void test (PETScWrappers::Vector &v)
       pattern[i] = true;
     }
 
-  v.compress ();
+  v.compress (VectorOperation::add);
 
   Vector<double> w;
   w=v;
diff --git a/tests/petsc/61.cc b/tests/petsc/61.cc
index 75547b0ebd..d883b98132 100644
--- a/tests/petsc/61.cc
+++ b/tests/petsc/61.cc
@@ -38,7 +38,7 @@ void test (PETScWrappers::Vector &v)
       pattern[i] = true;
     }
 
-  v.compress ();
+  v.compress (VectorOperation::add);
 
   Vector<double> w;
   w=v;
diff --git a/tests/petsc/66.cc b/tests/petsc/66.cc
index a921906e92..16427b58b1 100644
--- a/tests/petsc/66.cc
+++ b/tests/petsc/66.cc
@@ -59,7 +59,7 @@ void test (PETScWrappers::MatrixBase &m)
       norm_sqr += s*s;
       ++nnz;
     }
-  m.compress ();
+  m.compress (VectorOperation::add);
 
   deallog << m.frobenius_norm() << ' ' << std::sqrt (norm_sqr)
           << std::endl;
diff --git a/tests/petsc/67.cc b/tests/petsc/67.cc
index 5653714578..6c18362cd9 100644
--- a/tests/petsc/67.cc
+++ b/tests/petsc/67.cc
@@ -61,7 +61,7 @@ void test (PETScWrappers::MatrixBase &m)
       norm_sqr += s*s;
       ++nnz;
     }
-  m.compress ();
+  m.compress (VectorOperation::add);
 
   deallog << m.frobenius_norm() << ' ' << std::sqrt (norm_sqr)
           << std::endl;
diff --git a/tests/petsc/68.cc b/tests/petsc/68.cc
index 676a852245..2719ce002a 100644
--- a/tests/petsc/68.cc
+++ b/tests/petsc/68.cc
@@ -59,7 +59,7 @@ void test (PETScWrappers::MatrixBase &m)
       norm_sqr += s*s;
       ++nnz;
     }
-  m.compress ();
+  m.compress (VectorOperation::add);
 
   deallog << m.frobenius_norm() << ' ' << std::sqrt (norm_sqr)
           << std::endl;
diff --git a/tests/petsc/69.cc b/tests/petsc/69.cc
index 55e76a9394..c01c3f0839 100644
--- a/tests/petsc/69.cc
+++ b/tests/petsc/69.cc
@@ -61,7 +61,7 @@ void test (PETScWrappers::MatrixBase &m)
       norm_sqr += s*s;
       ++nnz;
     }
-  m.compress ();
+  m.compress (VectorOperation::add);
 
   deallog << m.frobenius_norm() << ' ' << std::sqrt (norm_sqr)
           << std::endl;
diff --git a/tests/petsc/deal_solver_01.cc b/tests/petsc/deal_solver_01.cc
index 1e73402ba0..88df05b761 100644
--- a/tests/petsc/deal_solver_01.cc
+++ b/tests/petsc/deal_solver_01.cc
@@ -86,9 +86,9 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress ();
-    f.compress ();
-    u.compress ();
+    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 46f04c6fa5..1a96781ed8 100644
--- a/tests/petsc/deal_solver_02.cc
+++ b/tests/petsc/deal_solver_02.cc
@@ -88,9 +88,9 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress ();
-    f.compress ();
-    u.compress ();
+    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 831ebb38ec..795c13704a 100644
--- a/tests/petsc/deal_solver_03.cc
+++ b/tests/petsc/deal_solver_03.cc
@@ -88,9 +88,9 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress ();
-    f.compress ();
-    u.compress ();
+    A.compress (VectorOperation::add);
+    f.compress (VectorOperation::add);
+    u.compress (VectorOperation::add);
 
     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 f4a88b8bc4..378775dd98 100644
--- a/tests/petsc/deal_solver_04.cc
+++ b/tests/petsc/deal_solver_04.cc
@@ -87,9 +87,9 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress ();
-    f.compress ();
-    u.compress ();
+    A.compress (VectorOperation::add);
+    f.compress (VectorOperation::add);
+    u.compress (VectorOperation::add);
 
     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 033b498db6..f6bf2ebd84 100644
--- a/tests/petsc/deal_solver_05.cc
+++ b/tests/petsc/deal_solver_05.cc
@@ -86,9 +86,9 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress ();
-    f.compress ();
-    u.compress ();
+    A.compress (VectorOperation::add);
+    f.compress (VectorOperation::add);
+    u.compress (VectorOperation::add);
 
     GrowingVectorMemory<PETScWrappers::Vector> mem;
     SolverQMRS<PETScWrappers::Vector> solver(control,mem);
diff --git a/tests/petsc/full_matrix_01.cc b/tests/petsc/full_matrix_01.cc
index 6125db2d92..697edcecab 100644
--- a/tests/petsc/full_matrix_01.cc
+++ b/tests/petsc/full_matrix_01.cc
@@ -33,7 +33,7 @@ void test (PETScWrappers::FullMatrix &m)
       if ((i+2*j+1) % 3 == 0)
         m.set (i,j, i*j*.5+.5);
 
-  m.compress ();
+  m.compress (VectorOperation::add);
 
   // then make sure we retrieve the same ones
   for (unsigned int i=0; i<m.m(); ++i)
diff --git a/tests/petsc/full_matrix_02.cc b/tests/petsc/full_matrix_02.cc
index 629ba76c68..1e9c279b1c 100644
--- a/tests/petsc/full_matrix_02.cc
+++ b/tests/petsc/full_matrix_02.cc
@@ -33,7 +33,7 @@ void test (PETScWrappers::FullMatrix &m)
       if ((i+2*j+1) % 3 == 0)
         m.add (i,j, i*j*.5+.5);
 
-  m.compress ();
+  m.compress (VectorOperation::add);
 
   // 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 36851aa8fc..a59679105f 100644
--- a/tests/petsc/full_matrix_05.cc
+++ b/tests/petsc/full_matrix_05.cc
@@ -38,7 +38,7 @@ void test (PETScWrappers::FullMatrix &m)
           ++counter;
         }
 
-  m.compress ();
+  m.compress (VectorOperation::add);
 
   // 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 941c32f8fe..7c90c423e8 100644
--- a/tests/petsc/full_matrix_06.cc
+++ b/tests/petsc/full_matrix_06.cc
@@ -33,7 +33,7 @@ void test (PETScWrappers::FullMatrix &m)
       if ((i+2*j+1) % 3 == 0)
         m.set (i,j, i*j*.5+.5);
 
-  m.compress ();
+  m.compress (VectorOperation::add);
 
   // compare against the exact value of the
   // l1-norm (max col-sum)
diff --git a/tests/petsc/full_matrix_07.cc b/tests/petsc/full_matrix_07.cc
index 29d8e07e0c..98277237c5 100644
--- a/tests/petsc/full_matrix_07.cc
+++ b/tests/petsc/full_matrix_07.cc
@@ -33,7 +33,7 @@ void test (PETScWrappers::FullMatrix &m)
       if ((i+2*j+1) % 3 == 0)
         m.set (i,j, i*j*.5+.5);
 
-  m.compress ();
+  m.compress (VectorOperation::add);
 
   // compare against the exact value of the
   // linfty-norm (max row-sum)
diff --git a/tests/petsc/full_matrix_08.cc b/tests/petsc/full_matrix_08.cc
index f3c6201ccb..5e741375e1 100644
--- a/tests/petsc/full_matrix_08.cc
+++ b/tests/petsc/full_matrix_08.cc
@@ -38,7 +38,7 @@ void test (PETScWrappers::FullMatrix &m)
         }
   norm = std::sqrt(norm);
 
-  m.compress ();
+  m.compress (VectorOperation::add);
 
   // 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 a7b0bb8a08..22781e7e10 100644
--- a/tests/petsc/full_matrix_09.cc
+++ b/tests/petsc/full_matrix_09.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 ();
+  m.compress (VectorOperation::add);
 
   // 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 203f39f0bc..f66fc56133 100644
--- a/tests/petsc/full_matrix_10.cc
+++ b/tests/petsc/full_matrix_10.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 ();
+  m.compress (VectorOperation::add);
 
   // 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 94d29eb3cf..a5c1ad4c26 100644
--- a/tests/petsc/full_matrix_iterator_01.cc
+++ b/tests/petsc/full_matrix_iterator_01.cc
@@ -30,7 +30,7 @@ void test ()
   m.set (0,0,1);
   m.set (1,1,2);
   m.set (1,2,3);
-  m.compress ();
+  m.compress (VectorOperation::add);
   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 daabd47328..392252810e 100644
--- a/tests/petsc/full_matrix_vector_01.cc
+++ b/tests/petsc/full_matrix_vector_01.cc
@@ -37,9 +37,9 @@ void test (PETScWrappers::Vector &v,
   for (unsigned int i=0; i<v.size(); ++i)
     v(i) = i;
 
-  m.compress ();
-  v.compress ();
-  w.compress ();
+  m.compress (VectorOperation::add);
+  v.compress (VectorOperation::add);
+  w.compress (VectorOperation::add);
 
   // 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 6bf332a6f1..b84340e73c 100644
--- a/tests/petsc/full_matrix_vector_02.cc
+++ b/tests/petsc/full_matrix_vector_02.cc
@@ -37,9 +37,9 @@ void test (PETScWrappers::Vector &v,
   for (unsigned int i=0; i<v.size(); ++i)
     v(i) = i;
 
-  m.compress ();
-  v.compress ();
-  w.compress ();
+  m.compress (VectorOperation::add);
+  v.compress (VectorOperation::add);
+  w.compress (VectorOperation::add);
 
   // 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 1309f3cb23..f1315a707b 100644
--- a/tests/petsc/full_matrix_vector_03.cc
+++ b/tests/petsc/full_matrix_vector_03.cc
@@ -40,9 +40,9 @@ void test (PETScWrappers::Vector &v,
       w(i) = i;
     }
 
-  m.compress ();
-  v.compress ();
-  w.compress ();
+  m.compress (VectorOperation::add);
+  v.compress (VectorOperation::add);
+  w.compress (VectorOperation::add);
 
   // w:=Mv
   m.vmult_add (w,v);
diff --git a/tests/petsc/full_matrix_vector_04.cc b/tests/petsc/full_matrix_vector_04.cc
index 4ed5c1b261..eaad49e52e 100644
--- a/tests/petsc/full_matrix_vector_04.cc
+++ b/tests/petsc/full_matrix_vector_04.cc
@@ -40,9 +40,9 @@ void test (PETScWrappers::Vector &v,
       w(i) = i;
     }
 
-  m.compress ();
-  v.compress ();
-  w.compress ();
+  m.compress (VectorOperation::add);
+  v.compress (VectorOperation::add);
+  w.compress (VectorOperation::add);
 
   // w:=Mv
   m.Tvmult_add (w,v);
diff --git a/tests/petsc/full_matrix_vector_05.cc b/tests/petsc/full_matrix_vector_05.cc
index 1d67d92f3d..565405adfa 100644
--- a/tests/petsc/full_matrix_vector_05.cc
+++ b/tests/petsc/full_matrix_vector_05.cc
@@ -40,9 +40,9 @@ void test (PETScWrappers::Vector &v,
       w(i) = i+1;
     }
 
-  m.compress ();
-  v.compress ();
-  w.compress ();
+  m.compress (VectorOperation::add);
+  v.compress (VectorOperation::add);
+  w.compress (VectorOperation::add);
 
   // <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 5514768e60..07b61100e7 100644
--- a/tests/petsc/full_matrix_vector_06.cc
+++ b/tests/petsc/full_matrix_vector_06.cc
@@ -36,8 +36,8 @@ void test (PETScWrappers::Vector &v)
   for (unsigned int i=0; i<v.size(); ++i)
     v(i) = i;
 
-  m.compress ();
-  v.compress ();
+  m.compress (VectorOperation::add);
+  v.compress (VectorOperation::add);
 
   // <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 4eead31ae6..93f4876507 100644
--- a/tests/petsc/full_matrix_vector_07.cc
+++ b/tests/petsc/full_matrix_vector_07.cc
@@ -41,9 +41,9 @@ void test (PETScWrappers::Vector &v,
       w(i) = i+1;
     }
 
-  m.compress ();
-  v.compress ();
-  w.compress ();
+  m.compress (VectorOperation::add);
+  v.compress (VectorOperation::add);
+  w.compress (VectorOperation::add);
 
   // x=w-Mv
   const double s = m.residual (x, v, w);
diff --git a/tests/petsc/parallel_sparse_matrix_01.cc b/tests/petsc/parallel_sparse_matrix_01.cc
index fc1cfb4e94..24f5ba8316 100644
--- a/tests/petsc/parallel_sparse_matrix_01.cc
+++ b/tests/petsc/parallel_sparse_matrix_01.cc
@@ -112,7 +112,7 @@ void test ()
     for (unsigned int j=0; j<csp.row_length(i); ++j)
       m.add (i, csp.column_number(i,j), 1.);
 
-  m.compress ();
+  m.compress (VectorOperation::add);
 
   deallog << "OK" << std::endl;
 }
diff --git a/tests/petsc/slowness_01.cc b/tests/petsc/slowness_01.cc
index 7290ec6caa..09176095b3 100644
--- a/tests/petsc/slowness_01.cc
+++ b/tests/petsc/slowness_01.cc
@@ -63,7 +63,7 @@ void test ()
             sparsity.add(global, global+N);
           }
       }
-  sparsity.compress();
+  sparsity.compress (VectorOperation::add);
 
 
   // next build the sparse matrix itself
diff --git a/tests/petsc/slowness_02.cc b/tests/petsc/slowness_02.cc
index 9611d27990..a5b4c63d8a 100644
--- a/tests/petsc/slowness_02.cc
+++ b/tests/petsc/slowness_02.cc
@@ -63,7 +63,7 @@ void test ()
             matrix.add(global, global+N, -1);
           }
       }
-  matrix.compress ();
+  matrix.compress (VectorOperation::add);
 
 
   // then do a single matrix-vector
diff --git a/tests/petsc/slowness_03.cc b/tests/petsc/slowness_03.cc
index d004153eb0..8b3687f811 100644
--- a/tests/petsc/slowness_03.cc
+++ b/tests/petsc/slowness_03.cc
@@ -66,7 +66,7 @@ void test ()
             matrix.add(global, global+N, -1);
           }
       }
-  matrix.compress ();
+  matrix.compress (VectorOperation::add);
 
   // then do a single matrix-vector
   // multiplication with subsequent formation
diff --git a/tests/petsc/slowness_04.cc b/tests/petsc/slowness_04.cc
index 221ff38714..295bc36966 100644
--- a/tests/petsc/slowness_04.cc
+++ b/tests/petsc/slowness_04.cc
@@ -97,7 +97,7 @@ void test ()
             matrix.add(global, global+N, rand());
           }
       }
-  matrix.compress ();
+  matrix.compress (VectorOperation::add);
 
   // then do a single matrix-vector
   // multiplication with subsequent formation
diff --git a/tests/petsc/solver_01.cc b/tests/petsc/solver_01.cc
index 5c68a79f41..9a7c99f66a 100644
--- a/tests/petsc/solver_01.cc
+++ b/tests/petsc/solver_01.cc
@@ -80,9 +80,9 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress ();
-    f.compress ();
-    u.compress ();
+    A.compress (VectorOperation::add);
+    f.compress (VectorOperation::add);
+    u.compress (VectorOperation::add);
 
     PETScWrappers::SolverRichardson solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
diff --git a/tests/petsc/solver_02.cc b/tests/petsc/solver_02.cc
index baf08b9793..e534c96a45 100644
--- a/tests/petsc/solver_02.cc
+++ b/tests/petsc/solver_02.cc
@@ -80,9 +80,9 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress ();
-    f.compress ();
-    u.compress ();
+    A.compress (VectorOperation::add);
+    f.compress (VectorOperation::add);
+    u.compress (VectorOperation::add);
 
     PETScWrappers::SolverChebychev solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
diff --git a/tests/petsc/solver_03.cc b/tests/petsc/solver_03.cc
index 57b908cc00..74524db363 100644
--- a/tests/petsc/solver_03.cc
+++ b/tests/petsc/solver_03.cc
@@ -81,9 +81,9 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress ();
-    f.compress ();
-    u.compress ();
+    A.compress (VectorOperation::add);
+    f.compress (VectorOperation::add);
+    u.compress (VectorOperation::add);
 
     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 84428afb54..4debf5c5c5 100644
--- a/tests/petsc/solver_03_mf.cc
+++ b/tests/petsc/solver_03_mf.cc
@@ -79,9 +79,9 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress ();
-    f.compress ();
-    u.compress ();
+    A.compress (VectorOperation::add);
+    f.compress (VectorOperation::add);
+    u.compress (VectorOperation::add);
 
     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 8c452a78b1..e775a8e41e 100644
--- a/tests/petsc/solver_03_precondition_boomeramg.cc
+++ b/tests/petsc/solver_03_precondition_boomeramg.cc
@@ -82,9 +82,9 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress ();
-    f.compress ();
-    u.compress ();
+    A.compress (VectorOperation::add);
+    f.compress (VectorOperation::add);
+    u.compress (VectorOperation::add);
 
     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 73a442f1a5..782e461eea 100644
--- a/tests/petsc/solver_03_precondition_boomeramg_symmetric.cc
+++ b/tests/petsc/solver_03_precondition_boomeramg_symmetric.cc
@@ -83,9 +83,9 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress ();
-    f.compress ();
-    u.compress ();
+    A.compress (VectorOperation::add);
+    f.compress (VectorOperation::add);
+    u.compress (VectorOperation::add);
 
     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 6551a25c05..70a6b99b12 100644
--- a/tests/petsc/solver_03_precondition_eisenstat.cc
+++ b/tests/petsc/solver_03_precondition_eisenstat.cc
@@ -81,9 +81,9 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress ();
-    f.compress ();
-    u.compress ();
+    A.compress (VectorOperation::add);
+    f.compress (VectorOperation::add);
+    u.compress (VectorOperation::add);
 
     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 6849d74979..628bee413f 100644
--- a/tests/petsc/solver_03_precondition_icc.cc
+++ b/tests/petsc/solver_03_precondition_icc.cc
@@ -81,9 +81,9 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress ();
-    f.compress ();
-    u.compress ();
+    A.compress (VectorOperation::add);
+    f.compress (VectorOperation::add);
+    u.compress (VectorOperation::add);
 
     PETScWrappers::SolverCG solver(control);
     PETScWrappers::PreconditionICC preconditioner(A);
diff --git a/tests/petsc/solver_03_precondition_ilu.cc b/tests/petsc/solver_03_precondition_ilu.cc
index faad508f0a..477852b01d 100644
--- a/tests/petsc/solver_03_precondition_ilu.cc
+++ b/tests/petsc/solver_03_precondition_ilu.cc
@@ -81,9 +81,9 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress ();
-    f.compress ();
-    u.compress ();
+    A.compress (VectorOperation::add);
+    f.compress (VectorOperation::add);
+    u.compress (VectorOperation::add);
 
     PETScWrappers::SolverCG solver(control);
     PETScWrappers::PreconditionILU preconditioner(A);
diff --git a/tests/petsc/solver_03_precondition_lu.cc b/tests/petsc/solver_03_precondition_lu.cc
index beabcf434a..ca0af140cb 100644
--- a/tests/petsc/solver_03_precondition_lu.cc
+++ b/tests/petsc/solver_03_precondition_lu.cc
@@ -88,9 +88,9 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress ();
-    f.compress ();
-    u.compress ();
+    A.compress (VectorOperation::add);
+    f.compress (VectorOperation::add);
+    u.compress (VectorOperation::add);
 
     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 05ffee26f7..3e8698282a 100644
--- a/tests/petsc/solver_03_precondition_parasails.cc
+++ b/tests/petsc/solver_03_precondition_parasails.cc
@@ -81,9 +81,9 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress ();
-    f.compress ();
-    u.compress ();
+    A.compress (VectorOperation::add);
+    f.compress (VectorOperation::add);
+    u.compress (VectorOperation::add);
 
     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 f0700eebc5..4a7893f6f4 100644
--- a/tests/petsc/solver_03_precondition_sor.cc
+++ b/tests/petsc/solver_03_precondition_sor.cc
@@ -81,9 +81,9 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress ();
-    f.compress ();
-    u.compress ();
+    A.compress (VectorOperation::add);
+    f.compress (VectorOperation::add);
+    u.compress (VectorOperation::add);
 
     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 67fd339d47..64ffcf8a55 100644
--- a/tests/petsc/solver_03_precondition_ssor.cc
+++ b/tests/petsc/solver_03_precondition_ssor.cc
@@ -81,9 +81,9 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress ();
-    f.compress ();
-    u.compress ();
+    A.compress (VectorOperation::add);
+    f.compress (VectorOperation::add);
+    u.compress (VectorOperation::add);
 
     PETScWrappers::SolverCG solver(control);
     PETScWrappers::PreconditionSSOR preconditioner(A);
diff --git a/tests/petsc/solver_04.cc b/tests/petsc/solver_04.cc
index 2d9819fba4..57fc9afb53 100644
--- a/tests/petsc/solver_04.cc
+++ b/tests/petsc/solver_04.cc
@@ -81,9 +81,9 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress ();
-    f.compress ();
-    u.compress ();
+    A.compress (VectorOperation::add);
+    f.compress (VectorOperation::add);
+    u.compress (VectorOperation::add);
 
     PETScWrappers::SolverBiCG solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
diff --git a/tests/petsc/solver_05.cc b/tests/petsc/solver_05.cc
index 15f6466e8c..be72a50a57 100644
--- a/tests/petsc/solver_05.cc
+++ b/tests/petsc/solver_05.cc
@@ -81,9 +81,9 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress ();
-    f.compress ();
-    u.compress ();
+    A.compress (VectorOperation::add);
+    f.compress (VectorOperation::add);
+    u.compress (VectorOperation::add);
 
     PETScWrappers::SolverGMRES solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
diff --git a/tests/petsc/solver_06.cc b/tests/petsc/solver_06.cc
index e19bddaf0c..bccd1a5c76 100644
--- a/tests/petsc/solver_06.cc
+++ b/tests/petsc/solver_06.cc
@@ -81,9 +81,9 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress ();
-    f.compress ();
-    u.compress ();
+    A.compress (VectorOperation::add);
+    f.compress (VectorOperation::add);
+    u.compress (VectorOperation::add);
 
     PETScWrappers::SolverBicgstab solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
diff --git a/tests/petsc/solver_07.cc b/tests/petsc/solver_07.cc
index 12d44a6ae2..30b25fa356 100644
--- a/tests/petsc/solver_07.cc
+++ b/tests/petsc/solver_07.cc
@@ -81,9 +81,9 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress ();
-    f.compress ();
-    u.compress ();
+    A.compress (VectorOperation::add);
+    f.compress (VectorOperation::add);
+    u.compress (VectorOperation::add);
 
     PETScWrappers::SolverCGS solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
diff --git a/tests/petsc/solver_08.cc b/tests/petsc/solver_08.cc
index 56b38ef60d..6d680f1448 100644
--- a/tests/petsc/solver_08.cc
+++ b/tests/petsc/solver_08.cc
@@ -81,9 +81,9 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress ();
-    f.compress ();
-    u.compress ();
+    A.compress (VectorOperation::add);
+    f.compress (VectorOperation::add);
+    u.compress (VectorOperation::add);
 
     PETScWrappers::SolverTFQMR solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
diff --git a/tests/petsc/solver_09.cc b/tests/petsc/solver_09.cc
index 3297bd8fe2..884195c278 100644
--- a/tests/petsc/solver_09.cc
+++ b/tests/petsc/solver_09.cc
@@ -85,9 +85,9 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress ();
-    f.compress ();
-    u.compress ();
+    A.compress (VectorOperation::add);
+    f.compress (VectorOperation::add);
+    u.compress (VectorOperation::add);
 
     PETScWrappers::SolverTCQMR solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
diff --git a/tests/petsc/solver_10.cc b/tests/petsc/solver_10.cc
index df81e7dbd2..9e4c41e08e 100644
--- a/tests/petsc/solver_10.cc
+++ b/tests/petsc/solver_10.cc
@@ -81,9 +81,9 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress ();
-    f.compress ();
-    u.compress ();
+    A.compress (VectorOperation::add);
+    f.compress (VectorOperation::add);
+    u.compress (VectorOperation::add);
 
     PETScWrappers::SolverCR solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
diff --git a/tests/petsc/solver_11.cc b/tests/petsc/solver_11.cc
index a7e585eb33..02c1e52981 100644
--- a/tests/petsc/solver_11.cc
+++ b/tests/petsc/solver_11.cc
@@ -82,9 +82,9 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress ();
-    f.compress ();
-    u.compress ();
+    A.compress (VectorOperation::add);
+    f.compress (VectorOperation::add);
+    u.compress (VectorOperation::add);
 
     PETScWrappers::SolverLSQR solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
diff --git a/tests/petsc/solver_12.cc b/tests/petsc/solver_12.cc
index 79e44183a9..39dc5468af 100644
--- a/tests/petsc/solver_12.cc
+++ b/tests/petsc/solver_12.cc
@@ -74,9 +74,9 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress ();
-    f.compress ();
-    u.compress ();
+    A.compress (VectorOperation::add);
+    f.compress (VectorOperation::add);
+    u.compress (VectorOperation::add);
 
     PETScWrappers::SolverPreOnly solver(control);
 
diff --git a/tests/petsc/solver_13.cc b/tests/petsc/solver_13.cc
index b28e88cafe..a567e11a63 100644
--- a/tests/petsc/solver_13.cc
+++ b/tests/petsc/solver_13.cc
@@ -72,9 +72,9 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  u(dim);
     u = 0.;
     f = 1.;
-    A.compress ();
-    f.compress ();
-    u.compress ();
+    A.compress VectorOperation::add();
+    f.compress VectorOperation::add();
+    u.compress VectorOperation::add();
     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 01956e4e02..4ca1387710 100644
--- a/tests/petsc/sparse_direct_mumps.cc
+++ b/tests/petsc/sparse_direct_mumps.cc
@@ -80,9 +80,9 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  u(dim);
     u = 0.;
     f = 1.;
-    A.compress ();
-    f.compress ();
-    u.compress ();
+    A.compress (VectorOperation::add);
+    f.compress (VectorOperation::add);
+    u.compress (VectorOperation::add);
 
     SolverControl cn;
     PETScWrappers::SparseDirectMUMPS solver(cn);
diff --git a/tests/petsc/sparse_matrix_01.cc b/tests/petsc/sparse_matrix_01.cc
index d1482625ba..767f140e60 100644
--- a/tests/petsc/sparse_matrix_01.cc
+++ b/tests/petsc/sparse_matrix_01.cc
@@ -34,13 +34,13 @@ void test ()
     for (unsigned int j=0; j<=i; ++j)
       m.set (i,j, i+2*j);
 
-  m.compress ();
+  m.compress (VectorOperation::add);
 
   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();
+  m2.compress (VectorOperation::add);
 
   // 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 f4073e5e32..493bfda620 100644
--- a/tests/petsc/sparse_matrix_02.cc
+++ b/tests/petsc/sparse_matrix_02.cc
@@ -35,7 +35,7 @@ void test ()
       m.set (i,j, i+2*j);
 
 
-  m.compress();
+  m.compress (VectorOperation::add);
 
   {
 
diff --git a/tests/petsc/sparse_matrix_iterator_01.cc b/tests/petsc/sparse_matrix_iterator_01.cc
index b070cc721c..e9a1444dfd 100644
--- a/tests/petsc/sparse_matrix_iterator_01.cc
+++ b/tests/petsc/sparse_matrix_iterator_01.cc
@@ -30,7 +30,7 @@ void test ()
   m.set (0,0,1);
   m.set (1,1,2);
   m.set (1,2,3);
-  m.compress ();
+  m.compress (VectorOperation::add);
   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 6562499418..96094736a0 100644
--- a/tests/petsc/sparse_matrix_vector_01.cc
+++ b/tests/petsc/sparse_matrix_vector_01.cc
@@ -37,9 +37,9 @@ void test (PETScWrappers::Vector &v,
   for (unsigned int i=0; i<v.size(); ++i)
     v(i) = i;
 
-  m.compress ();
-  v.compress ();
-  w.compress ();
+  m.compress (VectorOperation::add);
+  v.compress (VectorOperation::add);
+  w.compress (VectorOperation::add);
 
   // 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 b127bb49ee..a0d345cc78 100644
--- a/tests/petsc/sparse_matrix_vector_02.cc
+++ b/tests/petsc/sparse_matrix_vector_02.cc
@@ -37,9 +37,9 @@ void test (PETScWrappers::Vector &v,
   for (unsigned int i=0; i<v.size(); ++i)
     v(i) = i;
 
-  m.compress ();
-  v.compress ();
-  w.compress ();
+  m.compress (VectorOperation::add);
+  v.compress (VectorOperation::add);
+  w.compress (VectorOperation::add);
 
   // 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 667d48d3c3..09728f409c 100644
--- a/tests/petsc/sparse_matrix_vector_03.cc
+++ b/tests/petsc/sparse_matrix_vector_03.cc
@@ -40,9 +40,9 @@ void test (PETScWrappers::Vector &v,
       w(i) = i;
     }
 
-  m.compress ();
-  v.compress ();
-  w.compress ();
+  m.compress (VectorOperation::add);
+  v.compress (VectorOperation::add);
+  w.compress (VectorOperation::add);
 
   // 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 357e719c5e..6b311421a9 100644
--- a/tests/petsc/sparse_matrix_vector_04.cc
+++ b/tests/petsc/sparse_matrix_vector_04.cc
@@ -40,9 +40,9 @@ void test (PETScWrappers::Vector &v,
       w(i) = i;
     }
 
-  m.compress ();
-  v.compress ();
-  w.compress ();
+  m.compress (VectorOperation::add);
+  v.compress (VectorOperation::add);
+  w.compress (VectorOperation::add);
 
   // 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 d17bad38c1..6f83c1c0bb 100644
--- a/tests/petsc/sparse_matrix_vector_05.cc
+++ b/tests/petsc/sparse_matrix_vector_05.cc
@@ -40,9 +40,9 @@ void test (PETScWrappers::Vector &v,
       w(i) = i+1;
     }
 
-  m.compress ();
-  v.compress ();
-  w.compress ();
+  m.compress (VectorOperation::add);
+  v.compress (VectorOperation::add);
+  w.compress (VectorOperation::add);
 
   // <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 5901150030..144d246bf8 100644
--- a/tests/petsc/sparse_matrix_vector_06.cc
+++ b/tests/petsc/sparse_matrix_vector_06.cc
@@ -36,8 +36,8 @@ void test (PETScWrappers::Vector &v)
   for (unsigned int i=0; i<v.size(); ++i)
     v(i) = i;
 
-  m.compress ();
-  v.compress ();
+  m.compress (VectorOperation::add);
+  v.compress (VectorOperation::add);
 
   // <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 45bd1b80c8..b4bfef4701 100644
--- a/tests/petsc/sparse_matrix_vector_07.cc
+++ b/tests/petsc/sparse_matrix_vector_07.cc
@@ -41,9 +41,9 @@ void test (PETScWrappers::Vector &v,
       w(i) = i+1;
     }
 
-  m.compress ();
-  v.compress ();
-  w.compress ();
+  m.compress (VectorOperation::add);
+  v.compress (VectorOperation::add);
+  w.compress (VectorOperation::add);
 
   // x=w-Mv
   const double s = m.residual (x, v, w);
-- 
2.39.5