From 3692ab6d19671459071c7eb25ea423c7b45a8a25 Mon Sep 17 00:00:00 2001
From: Denis Davydov <davydden@gmail.com>
Date: Tue, 29 Mar 2016 23:02:58 +0200
Subject: [PATCH] use get_real_assert_zero_imag() instead of PetscRealPart()

---
 tests/mpi/ghost_01.cc                         | 12 ++++----
 tests/mpi/ghost_02.cc                         | 28 +++++++++----------
 tests/mpi/p4est_save_02.cc                    |  2 +-
 tests/mpi/p4est_save_03.cc                    |  4 +--
 tests/mpi/p4est_save_04.cc                    |  4 +--
 tests/mpi/periodicity_01.cc                   | 14 +++-------
 tests/mpi/petsc_bug_ghost_vector_01.cc        |  4 +--
 tests/mpi/petsc_distribute_01.cc              |  7 ++---
 tests/mpi/petsc_distribute_01_block.cc        | 14 ++++------
 tests/mpi/petsc_distribute_01_inhomogenous.cc | 19 ++++---------
 10 files changed, 45 insertions(+), 63 deletions(-)

diff --git a/tests/mpi/ghost_01.cc b/tests/mpi/ghost_01.cc
index 71578edf2b..5a917805c7 100644
--- a/tests/mpi/ghost_01.cc
+++ b/tests/mpi/ghost_01.cc
@@ -59,18 +59,18 @@ void test ()
   // check local values
   if (Utilities::MPI::this_mpi_process (MPI_COMM_WORLD) == 0)
     {
-      deallog << myid*2 << ":" << PetscRealPart(v(myid*2)) << std::endl;
-      deallog << myid*2+1 << ":" << PetscRealPart(v(myid*2+1)) << std::endl;
+      deallog << myid*2 << ":" << get_real_assert_zero_imag(v(myid*2)) << std::endl;
+      deallog << myid*2+1 << ":" << get_real_assert_zero_imag(v(myid*2+1)) << std::endl;
     }
 
-  Assert(PetscRealPart(v(myid*2)) == myid*4.0, ExcInternalError());
-  Assert(PetscRealPart(v(myid*2+1)) == myid*4.0+2.0, ExcInternalError());
+  Assert(get_real_assert_zero_imag(v(myid*2)) == myid*4.0, ExcInternalError());
+  Assert(get_real_assert_zero_imag(v(myid*2+1)) == myid*4.0+2.0, ExcInternalError());
 
 
   // check ghost values
   if (Utilities::MPI::this_mpi_process (MPI_COMM_WORLD) == 0)
-    deallog << "ghost: " << PetscRealPart(v(1)) << std::endl;
-  Assert(PetscRealPart(v(1)) == 2.0, ExcInternalError());
+    deallog << "ghost: " << get_real_assert_zero_imag(v(1)) << std::endl;
+  Assert(get_real_assert_zero_imag(v(1)) == 2.0, ExcInternalError());
 
   // done
   if (myid==0)
diff --git a/tests/mpi/ghost_02.cc b/tests/mpi/ghost_02.cc
index 63c02e9863..eaa877b5c6 100644
--- a/tests/mpi/ghost_02.cc
+++ b/tests/mpi/ghost_02.cc
@@ -52,10 +52,10 @@ void test ()
   vb = a;
   v2 = vb;
   if (Utilities::MPI::this_mpi_process (MPI_COMM_WORLD) == 0)
-    deallog << "ghost: " << PetscRealPart(v2(1)) << std::endl;
-  Assert(PetscRealPart(v2(1)) == 1.5, ExcInternalError());
-  Assert(PetscRealPart(v2(myid*2)) == 1.5, ExcInternalError());
-  Assert(PetscRealPart(v2(myid*2+1)) == 1.5, ExcInternalError());
+    deallog << "ghost: " << get_real_assert_zero_imag(v2(1)) << std::endl;
+  Assert(get_real_assert_zero_imag(v2(1)) == 1.5, ExcInternalError());
+  Assert(get_real_assert_zero_imag(v2(myid*2)) == 1.5, ExcInternalError());
+  Assert(get_real_assert_zero_imag(v2(myid*2+1)) == 1.5, ExcInternalError());
 
   // set local values
   vb(myid*2)=b;
@@ -69,26 +69,26 @@ void test ()
   // check local values
   if (Utilities::MPI::this_mpi_process (MPI_COMM_WORLD) == 0)
     {
-      deallog << myid*2 << ":" << PetscRealPart(v(myid*2)) << std::endl;
-      deallog << myid*2+1 << ":" << PetscRealPart(v(myid*2+1)) << std::endl;
+      deallog << myid*2 << ":" << get_real_assert_zero_imag(v(myid*2)) << std::endl;
+      deallog << myid*2+1 << ":" << get_real_assert_zero_imag(v(myid*2+1)) << std::endl;
     }
 
-  Assert(PetscRealPart(v(myid*2)) == myid*4.0, ExcInternalError());
-  Assert(PetscRealPart(v(myid*2+1)) == myid*4.0+2.0, ExcInternalError());
+  Assert(get_real_assert_zero_imag(v(myid*2)) == myid*4.0, ExcInternalError());
+  Assert(get_real_assert_zero_imag(v(myid*2+1)) == myid*4.0+2.0, ExcInternalError());
 
 
   // check ghost values
   if (Utilities::MPI::this_mpi_process (MPI_COMM_WORLD) == 0)
-    deallog << "ghost: " << PetscRealPart(v(1)) << std::endl;
-  Assert(PetscRealPart(v(1)) == 2.0, ExcInternalError());
+    deallog << "ghost: " << get_real_assert_zero_imag(v(1)) << std::endl;
+  Assert(get_real_assert_zero_imag(v(1)) == 2.0, ExcInternalError());
 
   //assignment from ghosted to ghosted
   v2 = v;
-  Assert(PetscRealPart(v2(1)) == 2.0, ExcInternalError());
-  Assert(PetscRealPart(v2(myid*2)) == myid*4.0, ExcInternalError());
-  Assert(PetscRealPart(v2(myid*2+1)) == myid*4.0+2.0, ExcInternalError());
+  Assert(get_real_assert_zero_imag(v2(1)) == 2.0, ExcInternalError());
+  Assert(get_real_assert_zero_imag(v2(myid*2)) == myid*4.0, ExcInternalError());
+  Assert(get_real_assert_zero_imag(v2(myid*2+1)) == myid*4.0+2.0, ExcInternalError());
   if (Utilities::MPI::this_mpi_process (MPI_COMM_WORLD) == 0)
-    deallog << "ghost: " << PetscRealPart(v2(1)) << std::endl;
+    deallog << "ghost: " << get_real_assert_zero_imag(v2(1)) << std::endl;
 
   // done
   if (myid==0)
diff --git a/tests/mpi/p4est_save_02.cc b/tests/mpi/p4est_save_02.cc
index 26592aab9a..05676d3e44 100644
--- a/tests/mpi/p4est_save_02.cc
+++ b/tests/mpi/p4est_save_02.cc
@@ -134,7 +134,7 @@ void test()
       {
         unsigned int idx = locally_owned_dofs.nth_index_in_set(i);
         //std::cout << '[' << idx << ']' << ' ' << solution(idx) << std::endl;
-        Assert(idx == PetscRealPart(solution(idx)),ExcInternalError());
+        Assert(idx == get_real_assert_zero_imag(solution(idx)),ExcInternalError());
 
       }
 
diff --git a/tests/mpi/p4est_save_03.cc b/tests/mpi/p4est_save_03.cc
index 723556e71f..bbadf28e9e 100644
--- a/tests/mpi/p4est_save_03.cc
+++ b/tests/mpi/p4est_save_03.cc
@@ -146,8 +146,8 @@ void test()
       {
         unsigned int idx = locally_owned_dofs.nth_index_in_set (i);
         //std::cout << '[' << idx << ']' << ' ' << solution(idx) << std::endl;
-        AssertThrow (idx == PetscRealPart(solution (idx)), ExcInternalError());
-        AssertThrow (2*idx == PetscRealPart(solution2 (idx)), ExcInternalError());
+        AssertThrow (idx == get_real_assert_zero_imag(solution (idx)), ExcInternalError());
+        AssertThrow (2*idx == get_real_assert_zero_imag(solution2 (idx)), ExcInternalError());
       }
 
     double norm = solution.l1_norm();
diff --git a/tests/mpi/p4est_save_04.cc b/tests/mpi/p4est_save_04.cc
index 57d9f278ba..28086449d6 100644
--- a/tests/mpi/p4est_save_04.cc
+++ b/tests/mpi/p4est_save_04.cc
@@ -88,7 +88,7 @@ void test()
         {
           unsigned int idx = locally_owned_dofs.nth_index_in_set (i);
           x (idx) = idx;
-          deallog << '[' << idx << ']' << ' ' << PetscRealPart(x(idx)) << std::endl;
+          deallog << '[' << idx << ']' << ' ' << get_real_assert_zero_imag(x(idx)) << std::endl;
         }
 
 
@@ -132,7 +132,7 @@ void test()
     for (unsigned int i = 0; i < locally_owned_dofs.n_elements(); ++i)
       {
         unsigned int idx = locally_owned_dofs.nth_index_in_set (i);
-        deallog << '[' << idx << ']' << ' ' << PetscRealPart(solution(idx)) << std::endl;
+        deallog << '[' << idx << ']' << ' ' << get_real_assert_zero_imag(solution(idx)) << std::endl;
       }
 
     deallog << "#cells = " << tr.n_global_active_cells() << std::endl;
diff --git a/tests/mpi/periodicity_01.cc b/tests/mpi/periodicity_01.cc
index fae79441fb..d645d559de 100644
--- a/tests/mpi/periodicity_01.cc
+++ b/tests/mpi/periodicity_01.cc
@@ -317,7 +317,7 @@ namespace Step40
 
     for (unsigned int i=0; i<value.size(); ++i)
       {
-        tmp[i]=PetscRealPart(value[i]);
+        tmp[i]=get_real_assert_zero_imag(value[i]);
         tmpI[i]=PetscImaginaryPart(value[i]);
       }
 #else
@@ -369,15 +369,13 @@ namespace Step40
 
         if (Utilities::MPI::this_mpi_process(mpi_communicator)==0)
           {
-            pcout << point1 << "\t" << PetscRealPart(value1[0]) << std::endl;
+            pcout << point1 << "\t" << get_real_assert_zero_imag(value1[0]) << std::endl;
             if (std::abs(value2[0]-value1[0])>1e-8)
               {
                 std::cout<<point1<< "\t" << value1[0] << std::endl;
                 std::cout<<point2<< "\t" << value2[0] << std::endl;
                 Assert(false, ExcInternalError());
               }
-            Assert(std::fabs(PetscImaginaryPart(value1[0]))<1e-10,
-                   ExcInternalError());
           }
       }
   }
@@ -421,24 +419,20 @@ namespace Step40
 
           if (Utilities::MPI::this_mpi_process(mpi_communicator)==0)
             {
-              pcout << point1 << "\t" << PetscRealPart(value1[0]) << std::endl;
+              pcout << point1 << "\t" << get_real_assert_zero_imag(value1[0]) << std::endl;
               if (std::abs(value2[0]-value1[0])>1e-8)
                 {
                   std::cout<<point1<< "\t" << value1[0] << std::endl;
                   std::cout<<point2<< "\t" << value2[0] << std::endl;
                   Assert(false, ExcInternalError());
                 }
-              Assert(std::fabs(PetscImaginaryPart(value1[0]))<1e-10,
-                     ExcInternalError());
-              pcout << point3 << "\t" << PetscRealPart(value3[0]) << std::endl;
+              pcout << point3 << "\t" << get_real_assert_zero_imag(value3[0]) << std::endl;
               if (std::abs(value4[0]-value3[0])>1e-8)
                 {
                   std::cout<<point3<< "\t" << value3[0] << std::endl;
                   std::cout<<point4<< "\t" << value4[0] << std::endl;
                   Assert(false, ExcInternalError());
                 }
-              Assert(std::fabs(PetscImaginaryPart(value3[0]))<1e-10,
-                     ExcInternalError());
             }
         }
   }
diff --git a/tests/mpi/petsc_bug_ghost_vector_01.cc b/tests/mpi/petsc_bug_ghost_vector_01.cc
index 38d56d8585..fb07d83c35 100644
--- a/tests/mpi/petsc_bug_ghost_vector_01.cc
+++ b/tests/mpi/petsc_bug_ghost_vector_01.cc
@@ -229,14 +229,14 @@ void test ()
     // zero. check this first for the locally owned elements...
     for (unsigned int i=0; i<handler.n_dofs(); ++i)
       if (locally_owned_dofs.is_element(i))
-        AssertThrow (PetscRealPart(vector(i)) == 0,
+        AssertThrow (get_real_assert_zero_imag(vector(i)) == 0,
                      ExcInternalError());
     // ...end then also for the ghost elements
     for (unsigned int i=0; i<handler.n_dofs(); ++i)
       if (locally_relevant_dofs.is_element(i)
           &&
           !locally_owned_dofs.is_element(i))
-        AssertThrow (PetscRealPart(vector(i)) == 0,
+        AssertThrow (get_real_assert_zero_imag(vector(i)) == 0,
                      ExcInternalError());
   }
 }
diff --git a/tests/mpi/petsc_distribute_01.cc b/tests/mpi/petsc_distribute_01.cc
index 38c0120a16..87164feaa8 100644
--- a/tests/mpi/petsc_distribute_01.cc
+++ b/tests/mpi/petsc_distribute_01.cc
@@ -99,11 +99,11 @@ void test()
 
   // verify correctness
   if (myid != 0)
-    AssertThrow (PetscRealPart(vec(vec.local_range().first+10)) == vec.local_range().first-25,
+    AssertThrow (get_real_assert_zero_imag(vec(vec.local_range().first+10)) == vec.local_range().first-25,
                  ExcInternalError());
 
   if (myid != n_processes-1)
-    AssertThrow (PetscRealPart(vec(vec.local_range().first+90)) == vec.local_range().first+105,
+    AssertThrow (get_real_assert_zero_imag(vec(vec.local_range().first+90)) == vec.local_range().first+105,
                  ExcInternalError());
 
   for (unsigned int i=vec.local_range().first; i<vec.local_range().second; ++i)
@@ -113,8 +113,7 @@ void test()
           (i != vec.local_range().first+90))
         {
           PetscScalar val = vec(i);
-          AssertThrow (std::fabs(PetscRealPart(val) - i) <= 1e-6, ExcInternalError());
-          AssertThrow (std::fabs(PetscImaginaryPart(val)) <= 1e-6, ExcInternalError());
+          AssertThrow (std::fabs(get_real_assert_zero_imag(val) - i) <= 1e-6, ExcInternalError());
         }
     }
 
diff --git a/tests/mpi/petsc_distribute_01_block.cc b/tests/mpi/petsc_distribute_01_block.cc
index fc176c4afc..209517cfee 100644
--- a/tests/mpi/petsc_distribute_01_block.cc
+++ b/tests/mpi/petsc_distribute_01_block.cc
@@ -112,19 +112,19 @@ void test()
 
   // verify correctness
   if (myid != 0)
-    AssertThrow (PetscRealPart(vec(vec.block(0).local_range().first+10)) == vec.block(0).local_range().first-25,
+    AssertThrow (get_real_assert_zero_imag(vec(vec.block(0).local_range().first+10)) == vec.block(0).local_range().first-25,
                  ExcInternalError());
 
   if (myid != n_processes-1)
-    AssertThrow (PetscRealPart(vec(vec.block(0).local_range().first+90)) == vec.block(0).local_range().first+105,
+    AssertThrow (get_real_assert_zero_imag(vec(vec.block(0).local_range().first+90)) == vec.block(0).local_range().first+105,
                  ExcInternalError());
 
   if (myid != 0)
-    AssertThrow (PetscRealPart(vec(vec.block(0).size()+vec.block(1).local_range().first+10)) == vec.block(1).local_range().first-25,
+    AssertThrow (get_real_assert_zero_imag(vec(vec.block(0).size()+vec.block(1).local_range().first+10)) == vec.block(1).local_range().first-25,
                  ExcInternalError());
 
   if (myid != n_processes-1)
-    AssertThrow (PetscRealPart(vec(vec.block(0).size()+vec.block(1).local_range().first+90)) == vec.block(1).local_range().first+105,
+    AssertThrow (get_real_assert_zero_imag(vec(vec.block(0).size()+vec.block(1).local_range().first+90)) == vec.block(1).local_range().first+105,
                  ExcInternalError());
 
 
@@ -135,8 +135,7 @@ void test()
           (i != vec.block(0).local_range().first+90))
         {
           PetscScalar val = vec.block(0)(i);
-          AssertThrow (std::fabs(PetscRealPart(val) - i) <= 1e-6, ExcInternalError());
-          AssertThrow (PetscImaginaryPart(val) == 0.0, ExcInternalError());
+          AssertThrow (std::fabs(get_real_assert_zero_imag(val) - i) <= 1e-6, ExcInternalError());
         }
     }
   for (unsigned int i=vec.block(1).local_range().first; i<vec.block(1).local_range().second; ++i)
@@ -146,8 +145,7 @@ void test()
           (i != vec.block(1).local_range().first+90))
         {
           PetscScalar val = vec.block(1)(i);
-          AssertThrow (std::fabs(PetscRealPart(val) - i) <= 1e-6, ExcInternalError());
-          AssertThrow ( PetscImaginaryPart(val) == 0.0, ExcInternalError());
+          AssertThrow (std::fabs(get_real_assert_zero_imag(val) - i) <= 1e-6, ExcInternalError());
         }
     }
 
diff --git a/tests/mpi/petsc_distribute_01_inhomogenous.cc b/tests/mpi/petsc_distribute_01_inhomogenous.cc
index eaa56a3904..bb92a31878 100644
--- a/tests/mpi/petsc_distribute_01_inhomogenous.cc
+++ b/tests/mpi/petsc_distribute_01_inhomogenous.cc
@@ -93,20 +93,12 @@ void test()
 
   // verify correctness
   if (myid != 0)
-    {
-      AssertThrow (PetscRealPart(vec(vec.local_range().first+10)) == vec.local_range().first-25+1,
-                   ExcInternalError());
-      AssertThrow (PetscImaginaryPart(vec(vec.local_range().first+10)) == 0.,
-                   ExcInternalError());
-    }
+    AssertThrow (get_real_assert_zero_imag(vec(vec.local_range().first+10)) == vec.local_range().first-25+1,
+                 ExcInternalError());
 
   if (myid != n_processes-1)
-    {
-      AssertThrow (PetscRealPart(vec(vec.local_range().first+90)) == vec.local_range().first+105-1,
-                   ExcInternalError());
-      AssertThrow (PetscImaginaryPart(vec(vec.local_range().first+90)) == 0.,
-                   ExcInternalError());
-    }
+    AssertThrow (get_real_assert_zero_imag(vec(vec.local_range().first+90)) == vec.local_range().first+105-1,
+                 ExcInternalError());
 
   for (unsigned int i=vec.local_range().first; i<vec.local_range().second; ++i)
     {
@@ -115,8 +107,7 @@ void test()
           (i != vec.local_range().first+90))
         {
           PetscScalar val = vec(i);
-          AssertThrow (std::fabs(PetscRealPart(val) - i) <= 1e-6, ExcInternalError());
-          AssertThrow (std::fabs(PetscImaginaryPart(val)) <= 1e-6, ExcInternalError());
+          AssertThrow (std::fabs(get_real_assert_zero_imag(val) - i) <= 1e-6, ExcInternalError());
         }
     }
 
-- 
2.39.5