From: Denis Davydov <davydden@gmail.com>
Date: Sun, 6 Mar 2016 17:20:23 +0000 (+0100)
Subject: fix petsc_complex test
X-Git-Tag: v8.5.0-rc1~1217^2~7
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=991a7b57f6b83998f7fce2be1db821383242fc47;p=dealii.git

fix petsc_complex test
---

diff --git a/tests/petsc_complex/00.cc b/tests/petsc_complex/00.cc
index f00cdb481b..17781e5d09 100644
--- a/tests/petsc_complex/00.cc
+++ b/tests/petsc_complex/00.cc
@@ -1,7 +1,6 @@
 // ---------------------------------------------------------------------
-// $Id: 00.cc$
 //
-// Copyright (C) 2013 by the deal.II authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
@@ -14,21 +13,6 @@
 //
 // ---------------------------------------------------------------------
 
-
-// PETSc includes: We would like to do this cleanly
-//
-#include <petscksp.h>
-#include <petscsys.h>
-#include <petscmath.h>
-//
-// though a horrible warning message
-//
-// warning: imaginary constants are a GCC extension [enabled by default]
-//
-// appears. This seems to be related to the direct useage of
-// PetscScalar in accordance with gcc.gnu.org bug 7263.
-
-// deal.II includes
 #include "../tests.h"
 #include <deal.II/base/logstream.h>
 
@@ -159,16 +143,16 @@ void init_petsc_complex ()
   deallog << "OK" << std::endl;
 }
 
+
 int main (int argc, char **argv)
 {
   std::ofstream logfile ("output");
   dealii::deallog.attach (logfile);
-  dealii::deallog.depth_console (0);
   deallog.threshold_double(1.e-10);
 
   try
     {
-      PetscInitialize (&argc, &argv, (char *) 0, (char *) 0);
+      Utilities::MPI::MPI_InitFinalize mpi_initialization (argc, argv, 1);
       {
         // initialisation (zero and nonzero)
         init_petsc_complex ();
@@ -182,7 +166,6 @@ int main (int argc, char **argv)
         divide_petsc_complex_by_a_number ();
         multiply_petsc_complex_by_a_number ();
       }
-      PetscFinalize ();
     }
 
   catch (std::exception &exc)
@@ -210,8 +193,6 @@ int main (int argc, char **argv)
       return 1;
     }
 
-  logfile << std::endl;
-
   return 0;
 }
 
diff --git a/tests/petsc_complex/00.output b/tests/petsc_complex/00.output
index 8fc890b1ae..a7cb78f5fb 100644
--- a/tests/petsc_complex/00.output
+++ b/tests/petsc_complex/00.output
@@ -11,4 +11,3 @@ DEAL::Check PetscScalar division
 DEAL::OK
 DEAL::Check PetscScalar multiply
 DEAL::OK
-
diff --git a/tests/petsc_complex/01.cc b/tests/petsc_complex/01.cc
index 4dbcfec56c..61a78e9054 100644
--- a/tests/petsc_complex/01.cc
+++ b/tests/petsc_complex/01.cc
@@ -1,7 +1,6 @@
 // ---------------------------------------------------------------------
-// $Id: 01.cc $
 //
-// Copyright (C) 2013 by the deal.II authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
@@ -31,7 +30,7 @@ void test (PETScWrappers::SparseMatrix &m)
       if ((i+2*j+1) % 3 == 0)
         m.set (i,j, std::complex<double> (0.,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_complex/02.cc b/tests/petsc_complex/02.cc
index dd70f49c83..ff4e8c87ab 100644
--- a/tests/petsc_complex/02.cc
+++ b/tests/petsc_complex/02.cc
@@ -1,7 +1,6 @@
 // ---------------------------------------------------------------------
-// $Id: 02.cc $
 //
-// Copyright (C) 2013 by the deal.II authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
diff --git a/tests/petsc_complex/04.cc b/tests/petsc_complex/04.cc
index b0a8d51e01..6312b33fbc 100644
--- a/tests/petsc_complex/04.cc
+++ b/tests/petsc_complex/04.cc
@@ -1,7 +1,6 @@
 // ---------------------------------------------------------------------
-// $Id: 04.cc $
 //
-// Copyright (C) 2013 by the deal.II authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
diff --git a/tests/petsc_complex/05.cc b/tests/petsc_complex/05.cc
index dd468a6e0b..f6e6422c10 100644
--- a/tests/petsc_complex/05.cc
+++ b/tests/petsc_complex/05.cc
@@ -1,7 +1,6 @@
 // ---------------------------------------------------------------------
-// $Id: 05.cc $
 //
-// Copyright (C) 2013 by the deal.II authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
@@ -36,7 +35,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_complex/06.cc b/tests/petsc_complex/06.cc
index e3cc027519..3790474535 100644
--- a/tests/petsc_complex/06.cc
+++ b/tests/petsc_complex/06.cc
@@ -1,7 +1,6 @@
 // ---------------------------------------------------------------------
-// $Id: 06.cc $
 //
-// Copyright (C) 2013 by the deal.II authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
@@ -31,7 +30,7 @@ void test (PETScWrappers::SparseMatrix &m)
       if ((i+2*j+1) % 3 == 0)
         m.set (i,j, std::complex<double> (0.,i*j*.5+.5));
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   // compare against the exact value of the l1-norm (max col-sum)
   deallog << m.l1_norm() << std::endl;
diff --git a/tests/petsc_complex/07.cc b/tests/petsc_complex/07.cc
index 3a1cafb1a9..75b106ebd4 100644
--- a/tests/petsc_complex/07.cc
+++ b/tests/petsc_complex/07.cc
@@ -1,7 +1,6 @@
 // ---------------------------------------------------------------------
-// $Id: 07.cc $
 //
-// Copyright (C) 2013 by the deal.II authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
@@ -32,7 +31,7 @@ void test (PETScWrappers::SparseMatrix &m)
       if ((i+2*j+1) % 3 == 0)
         m.set (i,j, std::complex<double> (0.,i*j*.5+.5));
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   // compare against the exact value of the linfty-norm (max row-sum).
   deallog << m.linfty_norm() << std::endl;
diff --git a/tests/petsc_complex/08.cc b/tests/petsc_complex/08.cc
index e34a00a744..3ab73e9b46 100644
--- a/tests/petsc_complex/08.cc
+++ b/tests/petsc_complex/08.cc
@@ -1,7 +1,6 @@
 // ---------------------------------------------------------------------
-// $Id: 08.cc $
 //
-// Copyright (C) 2013 by the deal.II authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
@@ -37,7 +36,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)
   deallog << m.frobenius_norm() << std::endl;
diff --git a/tests/petsc_complex/09.cc b/tests/petsc_complex/09.cc
index e04c644cc6..b4fb084a66 100644
--- a/tests/petsc_complex/09.cc
+++ b/tests/petsc_complex/09.cc
@@ -1,7 +1,6 @@
 // ---------------------------------------------------------------------
-// $Id: 09.cc $
 //
-// Copyright (C) 2013 by the deal.II authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
@@ -32,7 +31,7 @@ void test (PETScWrappers::SparseMatrix &m)
       if ((i+2*j+1) % 3 == 0)
         m.set (i,j, std::complex<double> (i*j*.5+.5,i*j*.5));
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   // then multiply everything by 1.25
   m *= 1.25;
diff --git a/tests/petsc_complex/10.cc b/tests/petsc_complex/10.cc
index 5d27dd617a..2e22c87c11 100644
--- a/tests/petsc_complex/10.cc
+++ b/tests/petsc_complex/10.cc
@@ -1,7 +1,6 @@
 // ---------------------------------------------------------------------
-// $id: 10.cc $
 //
-// Copyright (C) 2013 by the deal.II authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
@@ -31,7 +30,7 @@ void test (PETScWrappers::SparseMatrix &m)
       if ((i+2*j+1) % 3 == 0)
         m.set (i,j, std::complex<double> (0.,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 expect
diff --git a/tests/petsc_complex/11.cc b/tests/petsc_complex/11.cc
index 2b332d9d5b..07181f7bbe 100644
--- a/tests/petsc_complex/11.cc
+++ b/tests/petsc_complex/11.cc
@@ -1,7 +1,6 @@
 // ---------------------------------------------------------------------
-// $Id: 11.cc $
 //
-// Copyright (C) 2013 by the deal.II authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
@@ -30,7 +29,7 @@ void test (PETScWrappers::Vector &v)
   for (unsigned int k=0; k<v.size(); k+=1+k)
     v(k) = std::complex<double> (k,.5*k);
 
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   Assert (v.size() == 100, ExcInternalError());
 
diff --git a/tests/petsc_complex/12.cc b/tests/petsc_complex/12.cc
index 608fa3baee..d2ab3d5f4d 100644
--- a/tests/petsc_complex/12.cc
+++ b/tests/petsc_complex/12.cc
@@ -1,7 +1,6 @@
 // ---------------------------------------------------------------------
-// $Id: 12.cc $
 //
-// Copyright (C) 2013 by the deal.II authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
@@ -39,7 +38,7 @@ void test (PETScWrappers::Vector &v)
       pattern[k] = true;
     }
 
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // check that they are ok, and this time all of them
   for (unsigned int k=0; k<v.size(); ++k)
diff --git a/tests/petsc_complex/13.cc b/tests/petsc_complex/13.cc
index 03af432204..0d68cdf313 100644
--- a/tests/petsc_complex/13.cc
+++ b/tests/petsc_complex/13.cc
@@ -1,7 +1,6 @@
 // ---------------------------------------------------------------------
-// $Id: 13.cc $
 //
-// Copyright (C) 2013 by the deal.II authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
diff --git a/tests/petsc_complex/17.cc b/tests/petsc_complex/17.cc
index a4c85e11a7..637a8515fa 100644
--- a/tests/petsc_complex/17.cc
+++ b/tests/petsc_complex/17.cc
@@ -1,7 +1,6 @@
 // ---------------------------------------------------------------------
-// $Id: 17.cc $
 //
-// Copyright (C) 2013 by the deal.II authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
@@ -34,7 +33,7 @@ void test (PETScWrappers::Vector &v)
       v(k) = PetscScalar (k,2.*k);
       norm += std::fabs (1.*k+2.*k);
     }
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // then check the norm
   Assert (v.l1_norm() == norm, ExcInternalError());
diff --git a/tests/petsc_complex/18.cc b/tests/petsc_complex/18.cc
index c575c59fb5..37ef593f2d 100644
--- a/tests/petsc_complex/18.cc
+++ b/tests/petsc_complex/18.cc
@@ -1,7 +1,6 @@
 // ---------------------------------------------------------------------
-// $Id: 18.cc $
 //
-// Copyright (C) 2013 by the deal.II authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
@@ -38,7 +37,7 @@ void test (PETScWrappers::Vector &v)
       norm += std::fabs (1.*k*1.*k /*+ 1.*k*2*ki - 1.*k*2*ki*/ + 2*k*2.*k/*i*/);
     }
 
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // then check the l2-norm
   PetscReal l2_norm = v.l2_norm();
diff --git a/tests/petsc_complex/19.cc b/tests/petsc_complex/19.cc
index 4bb0c1bc4e..169f908281 100644
--- a/tests/petsc_complex/19.cc
+++ b/tests/petsc_complex/19.cc
@@ -1,7 +1,6 @@
 // ---------------------------------------------------------------------
-// $Id: 19.cc $
 //
-// Copyright (C) 2013 by the deal.II authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
@@ -37,7 +36,7 @@ void test (PETScWrappers::Vector &v)
       norm = std::max (norm,std::fabs (std::sqrt ((k*k + 2.*k*2.*k))));
     }
 
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // then check the norm
   Assert (v.linfty_norm() == norm, ExcInternalError());
diff --git a/tests/petsc_complex/20.cc b/tests/petsc_complex/20.cc
index 0774b2cef2..c0bc9db93b 100644
--- a/tests/petsc_complex/20.cc
+++ b/tests/petsc_complex/20.cc
@@ -1,7 +1,6 @@
 // ---------------------------------------------------------------------
-// $Id: 23.cc 31349 2013-10-20 19:07:06Z maier $
 //
-// Copyright (C) 2004 - 2013 by the deal.II authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
@@ -45,8 +44,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_complex/element_access_00.cc b/tests/petsc_complex/element_access_00.cc
index a088479882..09fef7f132 100644
--- a/tests/petsc_complex/element_access_00.cc
+++ b/tests/petsc_complex/element_access_00.cc
@@ -1,8 +1,6 @@
-
 // ---------------------------------------------------------------------
-// $Id:
 //
-// Copyright (C) 2013 by the deal.II authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
@@ -39,7 +37,7 @@ void test_vector (PETScWrappers::Vector &v)
   for (unsigned int k=0; k<v.size(); ++k)
     v(k) = std::complex<double> (k,v.size()-k);
 
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // check that is what we get by casting PetscScalar to std::real()
   // and std::imag()
@@ -67,7 +65,7 @@ int main (int argc, char **argv)
 
   try
     {
-      PetscInitialize (&argc, &argv, (char *) 0, (char *) 0);
+      Utilities::MPI::MPI_InitFinalize mpi_initialization(argc, argv, 1);
       {
         PETScWrappers::Vector v (5);
         test_vector (v);
@@ -75,7 +73,6 @@ int main (int argc, char **argv)
         deallog << "vector:" << std::endl;
         v.print (logfile);
       }
-      PetscFinalize ();
     }
 
 
diff --git a/tests/petsc_complex/element_access_01.cc b/tests/petsc_complex/element_access_01.cc
index 31b0a4f940..69c28b990b 100644
--- a/tests/petsc_complex/element_access_01.cc
+++ b/tests/petsc_complex/element_access_01.cc
@@ -1,8 +1,6 @@
-
 // ---------------------------------------------------------------------
-// $Id:
 //
-// Copyright (C) 2013 by the deal.II authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
@@ -42,7 +40,7 @@ void test_matrix (PETScWrappers::SparseMatrix &m)
         m.set (k,l, std::complex<double> (el_r,el_i));
       }
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   for (unsigned int k=0; k<m.m(); ++k)
     for (unsigned int l=0; l<m.n(); ++l)
@@ -66,7 +64,7 @@ int main (int argc, char **argv)
 
   try
     {
-      PetscInitialize (&argc, &argv, (char *) 0, (char *) 0);
+      Utilities::MPI::MPI_InitFinalize mpi_initialization(argc, argv, 1);
       {
         PETScWrappers::SparseMatrix m (5,5,5);
         test_matrix (m);
@@ -74,7 +72,6 @@ int main (int argc, char **argv)
         deallog << "matrix:" << std::endl;
         m.print (logfile);
       }
-      PetscFinalize ();
     }
 
 
diff --git a/tests/petsc_complex/element_access_02.cc b/tests/petsc_complex/element_access_02.cc
index 0b66ad3801..02b4d04eda 100644
--- a/tests/petsc_complex/element_access_02.cc
+++ b/tests/petsc_complex/element_access_02.cc
@@ -1,8 +1,6 @@
-
 // ---------------------------------------------------------------------
-// $Id: element_access_02.cc
 //
-// Copyright (C) 2013 by the deal.II authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
@@ -35,7 +33,7 @@ void test (PETScWrappers::SparseMatrix &m)
       if (k>l)
         m.set (k,l, PetscScalar (k+l,-1.*(k+l)));
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   // check the matrix is correctly filled
   for (unsigned int k=0; k<m.m(); ++k)
@@ -55,7 +53,7 @@ int main (int argc, char **argv)
 
   try
     {
-      PetscInitialize (&argc, &argv, (char *) 0, (char *) 0);
+      Utilities::MPI::MPI_InitFinalize mpi_initialization(argc, argv, 1);
       {
         PETScWrappers::SparseMatrix m (10,10,10);
         test (m);
@@ -63,7 +61,6 @@ int main (int argc, char **argv)
         deallog << "matrix:" << std::endl;
         m.print (logfile);
       }
-      PetscFinalize ();
     }
 
   catch (std::exception &exc)
diff --git a/tests/petsc_complex/element_access_03.cc b/tests/petsc_complex/element_access_03.cc
index df1a0fc061..9d525755ae 100644
--- a/tests/petsc_complex/element_access_03.cc
+++ b/tests/petsc_complex/element_access_03.cc
@@ -1,8 +1,6 @@
-
 // ---------------------------------------------------------------------
-// $Id:
 //
-// Copyright (C) 2013 by the deal.II authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
@@ -43,7 +41,7 @@ void test_matrix (PETScWrappers::SparseMatrix &m)
     for (unsigned int l=0; l<m.n(); ++l)
       m.set (k,l, PetscScalar ((k+l)*r,(k+l)*i));
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   // Check elements have the correct value
   for (unsigned int k=0; k<m.m(); ++k)
@@ -75,7 +73,7 @@ int main (int argc, char **argv)
 
   try
     {
-      PetscInitialize (&argc, &argv, (char *) 0, (char *) 0);
+      Utilities::MPI::MPI_InitFinalize mpi_initialization(argc, argv, 1);
       {
         PETScWrappers::SparseMatrix m (5,5,5);
         test_matrix (m);
@@ -83,7 +81,6 @@ int main (int argc, char **argv)
         deallog << "matrix:" << std::endl;
         m.print (logfile);
       }
-      PetscFinalize ();
     }
 
   catch (std::exception &exc)
diff --git a/tests/petsc_complex/fe_get_function_values.cc b/tests/petsc_complex/fe_get_function_values.cc
index 999e4f0d91..9aca85dff6 100644
--- a/tests/petsc_complex/fe_get_function_values.cc
+++ b/tests/petsc_complex/fe_get_function_values.cc
@@ -1,7 +1,6 @@
 // ---------------------------------------------------------------------
-// $Id: testsuite.html 32446 2014-05-21 14:31:22Z davydov $
 //
-// Copyright (C) 2013 by the deal.II Authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
@@ -170,8 +169,8 @@ int main (int argc, char *argv[])
   deallog.depth_console(0);
   deallog.threshold_double(1.e-10);
 
-  Utilities::MPI::MPI_InitFinalize mpi_initialization (argc, argv);
-
-  test ();
-
-}
\ No newline at end of file
+  Utilities::MPI::MPI_InitFinalize mpi_initialization(argc, argv, 1);
+  {
+    test ();
+  }
+}
diff --git a/tests/petsc_complex/parallel_sparse_matrix_01.cc b/tests/petsc_complex/parallel_sparse_matrix_01.cc
index 503bb3f2ff..5eb530915e 100644
--- a/tests/petsc_complex/parallel_sparse_matrix_01.cc
+++ b/tests/petsc_complex/parallel_sparse_matrix_01.cc
@@ -1,7 +1,6 @@
 // ---------------------------------------------------------------------
-// $Id: testsuite.html 32446 2014-02-10 14:31:22Z bangerth $
 //
-// Copyright (C) 2013 by the deal.II Authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
@@ -113,7 +112,7 @@ void test (const unsigned int poly_degree = 1)
       if (cell->is_locally_owned())
         {
           fe_values.reinit (cell);
-          cell_mass_matrix    = PetscScalar(0);
+          cell_mass_matrix    = PetscScalar();
 
           for (unsigned int q_point=0; q_point<n_q_points; ++q_point)
             for (unsigned int i=0; i<dofs_per_cell; ++i)
@@ -160,8 +159,8 @@ void test (const unsigned int poly_degree = 1)
   vector.compress(VectorOperation::insert);
   constraints.distribute(vector);
 
-  deallog<<"symmetric: "<<mass_matrix.is_symmetric()<<std::endl;
-  deallog<<"hermitian: "<<mass_matrix.is_hermitian()<<std::endl;
+  deallog<<"symmetric: "<<(mass_matrix.is_symmetric()==PETSC_TRUE)<<std::endl;
+  deallog<<"hermitian: "<<(mass_matrix.is_hermitian()==PETSC_TRUE)<<std::endl;
 
   PETScWrappers::MPI::Vector tmp(vector);
   mass_matrix.vmult (tmp, vector);
@@ -192,8 +191,8 @@ int main (int argc, char *argv[])
   deallog.depth_console(0);
   deallog.threshold_double(1.e-10);
 
-  Utilities::MPI::MPI_InitFinalize mpi_initialization (argc, argv);
-
-  test<2>();
-
-}
\ No newline at end of file
+  Utilities::MPI::MPI_InitFinalize mpi_initialization(argc, argv, 1);
+  {
+    test<2>();
+  }
+}
diff --git a/tests/petsc_complex/solver_real_01.cc b/tests/petsc_complex/solver_real_01.cc
index 9c045ca958..c5126c999e 100644
--- a/tests/petsc_complex/solver_real_01.cc
+++ b/tests/petsc_complex/solver_real_01.cc
@@ -1,7 +1,6 @@
 // ---------------------------------------------------------------------
-// $Id: solver_real_01.cc $
 //
-// Copyright (C) 2013 by the deal.II authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
@@ -85,9 +84,9 @@ 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);
+    f.compress (VectorOperation::insert);
+    u.compress (VectorOperation::insert);
 
     PETScWrappers::SolverRichardson solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
diff --git a/tests/petsc_complex/solver_real_02.cc b/tests/petsc_complex/solver_real_02.cc
index 28f4e0484b..8d952c9285 100644
--- a/tests/petsc_complex/solver_real_02.cc
+++ b/tests/petsc_complex/solver_real_02.cc
@@ -1,7 +1,6 @@
 // ---------------------------------------------------------------------
-// $Id: solver_real_02.cc $
 //
-// Copyright (C) 2013 by the deal.II authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
@@ -83,9 +82,9 @@ 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);
+    f.compress (VectorOperation::insert);
+    u.compress (VectorOperation::insert);
 
     PETScWrappers::SolverChebychev solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
diff --git a/tests/petsc_complex/solver_real_03.cc b/tests/petsc_complex/solver_real_03.cc
index 0b69601e07..cc5406ac71 100644
--- a/tests/petsc_complex/solver_real_03.cc
+++ b/tests/petsc_complex/solver_real_03.cc
@@ -1,7 +1,6 @@
 // ---------------------------------------------------------------------
-// $Id: solver_real_03.cc $
 //
-// Copyright (C) 2013 by the deal.II authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
@@ -82,9 +81,9 @@ 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);
+    f.compress (VectorOperation::insert);
+    u.compress (VectorOperation::insert);
 
     PETScWrappers::SolverCG solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
diff --git a/tests/petsc_complex/solver_real_03_mf.cc b/tests/petsc_complex/solver_real_03_mf.cc
index b3197e70ea..1bea00f78f 100644
--- a/tests/petsc_complex/solver_real_03_mf.cc
+++ b/tests/petsc_complex/solver_real_03_mf.cc
@@ -1,7 +1,6 @@
 // ---------------------------------------------------------------------
-// $Id: solver_03_mf.cc $
 //
-// Copyright (C) 2013 by the deal.II authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
@@ -92,9 +91,9 @@ 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);
+    f.compress (VectorOperation::insert);
+    u.compress (VectorOperation::insert);
 
     PETScWrappers::SolverCG solver(control);
     PETScWrappers::PreconditionNone preconditioner(A);
diff --git a/tests/petsc_complex/solver_real_04.cc b/tests/petsc_complex/solver_real_04.cc
index 9ace49088c..eed403f0b8 100644
--- a/tests/petsc_complex/solver_real_04.cc
+++ b/tests/petsc_complex/solver_real_04.cc
@@ -1,7 +1,6 @@
 // ---------------------------------------------------------------------
-// $Id: solver_real_04.cc $
 //
-// Copyright (C) 2013 by the deal.II authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
@@ -82,9 +81,9 @@ 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);
+    f.compress (VectorOperation::insert);
+    u.compress (VectorOperation::insert);
 
     PETScWrappers::SolverBiCG solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
diff --git a/tests/petsc_complex/sparse_matrix_01.cc b/tests/petsc_complex/sparse_matrix_01.cc
index 7d7a749b1c..46f783f61d 100644
--- a/tests/petsc_complex/sparse_matrix_01.cc
+++ b/tests/petsc_complex/sparse_matrix_01.cc
@@ -1,7 +1,6 @@
 // ---------------------------------------------------------------------
-// $Id sparse_matrix_01.cc $
 //
-// Copyright (C) 2013 by the deal.II authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
@@ -33,13 +32,13 @@ void test ()
     for (unsigned int l=0; l<=k; ++l)
       m.set (k,l, k+2*l);
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   PETScWrappers::SparseMatrix m2(s,s,s);
   m2.set(0,1,5.0);
   for (unsigned int k=0; k<m2.m(); ++k)
     m2.set(k,k,PetscScalar(1.0+k,-1.0-k));
-  m2.compress(VectorOperation::add);
+  m2.compress(VectorOperation::insert);
 
   // we now print the matrix m, it is all real. Then print after
   // adding m2 which is complex, and then subtract m2 again to get the
diff --git a/tests/petsc_complex/vector_02.cc b/tests/petsc_complex/vector_02.cc
index 33586ea2e9..9cacfcbe95 100644
--- a/tests/petsc_complex/vector_02.cc
+++ b/tests/petsc_complex/vector_02.cc
@@ -1,7 +1,6 @@
 // ---------------------------------------------------------------------
-// $Id vector_02.cc $
 //
-// Copyright (C) 2013 by the deal.II authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
@@ -33,7 +32,7 @@ void test ()
   for (unsigned int k=0; k<v.size(); ++k)
     v(k) = k;
 
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   PETScWrappers::Vector v2(s);
   for (unsigned int k=0; k<v2.size(); ++k)
@@ -44,7 +43,7 @@ void test ()
       v2(k) = PetscScalar (k,-1.0*k);
     }
 
-  v2.compress(VectorOperation::add);
+  v2.compress(VectorOperation::insert);
 
   // we now print the vector v, it is all real. Then print after
   // adding v2 which is complex, and then subtract v2 again to get the
diff --git a/tests/petsc_complex/vector_assign_01.cc b/tests/petsc_complex/vector_assign_01.cc
index d8f85360c0..c952cf3569 100644
--- a/tests/petsc_complex/vector_assign_01.cc
+++ b/tests/petsc_complex/vector_assign_01.cc
@@ -1,7 +1,6 @@
 // ---------------------------------------------------------------------
-// $Id: vector_assign_01.cc $
 //
-// Copyright (C) 2013 by the deal.II authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
diff --git a/tests/petsc_complex/vector_assign_02.cc b/tests/petsc_complex/vector_assign_02.cc
index 64c20333b0..ffc0923934 100644
--- a/tests/petsc_complex/vector_assign_02.cc
+++ b/tests/petsc_complex/vector_assign_02.cc
@@ -1,7 +1,6 @@
 // ---------------------------------------------------------------------
-// $Id: vector_assign_02.cc $
 //
-// Copyright (C) 2013 by the deal.II authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
diff --git a/tests/petsc_complex/vector_equality_1.cc b/tests/petsc_complex/vector_equality_1.cc
index d9d14e9a7c..ff9953e85e 100644
--- a/tests/petsc_complex/vector_equality_1.cc
+++ b/tests/petsc_complex/vector_equality_1.cc
@@ -1,7 +1,6 @@
 // ---------------------------------------------------------------------
-// $Id: vector_equality_1.cc $
 //
-// Copyright (C) 2013 by the deal.II authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
diff --git a/tests/petsc_complex/vector_equality_2.cc b/tests/petsc_complex/vector_equality_2.cc
index 53aa85cce2..814f8043bb 100644
--- a/tests/petsc_complex/vector_equality_2.cc
+++ b/tests/petsc_complex/vector_equality_2.cc
@@ -1,7 +1,6 @@
 // ---------------------------------------------------------------------
-// $Id: vector_equality_2.cc $
 //
-// Copyright (C) 2013 by the deal.II authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
diff --git a/tests/petsc_complex/vector_print.cc b/tests/petsc_complex/vector_print.cc
index 7d7b64cf1b..071fa8b79a 100644
--- a/tests/petsc_complex/vector_print.cc
+++ b/tests/petsc_complex/vector_print.cc
@@ -1,7 +1,6 @@
 // ---------------------------------------------------------------------
-// $Id: vector_print.cc $
 //
-// Copyright (C) 2013 by the deal.II authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
diff --git a/tests/petsc_complex/vector_wrap_01.cc b/tests/petsc_complex/vector_wrap_01.cc
index 2f520c6f47..d9c1f4c5e1 100644
--- a/tests/petsc_complex/vector_wrap_01.cc
+++ b/tests/petsc_complex/vector_wrap_01.cc
@@ -1,7 +1,6 @@
 // ---------------------------------------------------------------------
-// $Id: vector_wrap_01.cc $
 //
-// Copyright (C) 2013 by the deal.II authors
+// Copyright (C) 2004 - 2016 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //