From: Martin Kronbichler <kronbichler@lnm.mw.tum.de>
Date: Fri, 1 Apr 2016 06:18:41 +0000 (+0200)
Subject: Merge pull request #2395 from masterleinad/periodic_neighbor
X-Git-Tag: v8.5.0-rc1~1151
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=f29876fed0862ab1ccb565e140996506968c9c46;p=dealii.git

Merge pull request #2395 from masterleinad/periodic_neighbor

Periodic neighbors
---

f29876fed0862ab1ccb565e140996506968c9c46
diff --cc tests/mpi/periodicity_01.cc
index bfa9e4eb11,59924a2639..aef92fd4c8
--- a/tests/mpi/periodicity_01.cc
+++ b/tests/mpi/periodicity_01.cc
@@@ -333,16 -328,19 +333,19 @@@ namespace Step4
      for (unsigned int i = 0; i<cycle; i++)
        n_points*=2;
  
+     //don't test exactly at the support points, since point_value is not stable there
+     const double eps = 1./(16.*n_points);
+ 
      for (unsigned int i=1; i< n_points; i++)
        {
 -        Vector<double> value1(1);
 -        Vector<double> value2(1);
 +        Vector<PetscScalar> value1(1);
 +        Vector<PetscScalar> value2(1);
  
          Point <2> point1;
-         point1(0)=1.*i/n_points;
+         point1(0)=1.*i/n_points+eps;
          point1(1)=0.;
          Point <2> point2;
-         point2(0)=1.*i/n_points;
+         point2(0)=1.*i/n_points+eps;
          point2(1)=1.;
  
          get_point_value (point1, 0, value1);
@@@ -371,27 -372,27 +377,27 @@@
      for (unsigned int i=1; i< n_points; i++)
        for (unsigned int j=1; j< n_points; j++)
          {
 -          Vector<double> value1(1);
 -          Vector<double> value2(1);
 -          Vector<double> value3(1);
 -          Vector<double> value4(1);
 +          Vector<PetscScalar> value1(1);
 +          Vector<PetscScalar> value2(1);
 +          Vector<PetscScalar> value3(1);
 +          Vector<PetscScalar> value4(1);
  
            Point <3> point1;
-           point1(0)=1.*i/n_points;
-           point1(1)=1.*j/n_points;
-           point1(2)=0.;
+           point1(0)=1.*i/n_points+eps;
+           point1(1)=1.*j/n_points+eps;
+           point1(2)=0;
            Point <3> point2;
-           point2(0)=1.*i/n_points;
-           point2(1)=1.*j/n_points;
+           point2(0)=1.*i/n_points+eps;
+           point2(1)=1.*j/n_points+eps;
            point2(2)=1.;
            Point <3> point3;
-           point3(0)=1.*i/n_points;
+           point3(0)=1.*i/n_points+eps;
            point3(1)=0.;
-           point3(2)=1.*j/n_points;;
+           point3(2)=1.*j/n_points+eps;
            Point <3> point4;
-           point4(0)=1.*i/n_points;
+           point4(0)=1.*i/n_points+eps;
            point4(1)=1.;
-           point4(2)=1.*j/n_points;;
+           point4(2)=1.*j/n_points+eps;
  
            get_point_value (point1, 0, value1);
            get_point_value (point2, 0, value2);