]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Setup tests to test periodic_face_map as well
authorDaniel Arndt <daniel.arndt@iwr.uni-heidelberg.de>
Tue, 22 Mar 2016 14:33:31 +0000 (15:33 +0100)
committerDaniel Arndt <daniel.arndt@iwr.uni-heidelberg.de>
Tue, 29 Mar 2016 11:24:09 +0000 (13:24 +0200)
tests/mpi/periodicity_01.cc
tests/mpi/periodicity_01.with_petsc=true.mpirun=3.output
tests/mpi/periodicity_01.with_petsc=true.mpirun=5.output
tests/mpi/periodicity_01.with_petsc=true.mpirun=7.output
tests/mpi/periodicity_02.cc
tests/mpi/periodicity_02.with_trilinos=true.mpirun=3.output
tests/mpi/periodicity_02.with_trilinos=true.mpirun=5.output
tests/mpi/periodicity_02.with_trilinos=true.mpirun=7.output
tests/mpi/periodicity_04.cc

index b64c039f85c6e1da49b2fbcc501509cca8fc0d2a..59924a2639dfb53e810c22a75a9a28d78f1eaebd 100644 (file)
@@ -328,16 +328,19 @@ namespace Step40
     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);
 
         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);
@@ -363,6 +366,9 @@ namespace Step40
     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++)
       for (unsigned int j=1; j< n_points; j++)
         {
@@ -372,21 +378,21 @@ namespace Step40
           Vector<double> 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);
index 4d8e21abc1aac1b3e2ee4edc74c31d2371fd1519..e758ab7dfd96dde25a72052d6b25ab3a8173fe58 100644 (file)
@@ -5,27 +5,27 @@ Cycle 0:
 DEAL::Starting value 0.0139546
 DEAL::Convergence step 8 value 0
    Solved in 8 iterations.
-0.500000 0.00000       -0.00549678
+0.531250 0.00000       -0.00531794
 
 Cycle 1:
 DEAL::Starting value 0.0223730
 DEAL::Convergence step 8 value 0
    Solved in 8 iterations.
-0.250000 0.00000       -0.00130251
-0.500000 0.00000       -0.00550617
-0.750000 0.00000       -0.00183466
+0.265625 0.00000       -0.00177028
+0.515625 0.00000       -0.00542131
+0.765625 0.00000       -0.00148613
 
 Cycle 2:
 DEAL::Starting value 0.0270972
 DEAL::Convergence step 8 value 0
    Solved in 8 iterations.
-0.125000 0.00000       0.00141483
-0.250000 0.00000       -0.00124254
-0.375000 0.00000       -0.00429173
-0.500000 0.00000       -0.00545458
-0.625000 0.00000       -0.00429054
-0.750000 0.00000       -0.00185255
-0.875000 0.00000       0.000356647
+0.132812 0.00000       0.00130870
+0.257812 0.00000       -0.00146441
+0.382812 0.00000       -0.00443592
+0.507812 0.00000       -0.00544040
+0.632812 0.00000       -0.00415112
+0.757812 0.00000       -0.00167411
+0.882812 0.00000       0.000448975
 
 Testing for dim=3
 
@@ -33,131 +33,131 @@ Cycle 0:
 DEAL::Starting value 0.00635044
 DEAL::Convergence step 7 value 0
    Solved in 7 iterations.
-0.500000 0.500000 0.00000      0.000839383
-0.500000 0.00000 0.500000      0.000542250
+0.531250 0.531250 0.00000      0.000730918
+0.531250 0.00000 0.531250      0.000460360
 
 Cycle 1:
 DEAL::Starting value 0.0136523
 DEAL::Convergence step 7 value 0
    Solved in 7 iterations.
-0.250000 0.250000 0.00000      -7.02574e-05
-0.250000 0.00000 0.250000      -8.84955e-05
-0.250000 0.500000 0.00000      2.30224e-05
-0.250000 0.00000 0.500000      -1.58988e-05
-0.250000 0.750000 0.00000      -5.50224e-05
-0.250000 0.00000 0.750000      -6.78044e-05
-0.500000 0.250000 0.00000      -0.000245815
-0.500000 0.00000 0.250000      -0.000367442
-0.500000 0.500000 0.00000      0.000854616
-0.500000 0.00000 0.500000      0.000545283
-0.500000 0.750000 0.00000      -7.60057e-05
-0.500000 0.00000 0.750000      -0.000116497
-0.750000 0.250000 0.00000      -0.000105394
-0.750000 0.00000 0.250000      -0.000112911
-0.750000 0.500000 0.00000      -2.69468e-06
-0.750000 0.00000 0.500000      -3.29142e-05
-0.750000 0.750000 0.00000      -7.97320e-05
-0.750000 0.00000 0.750000      -8.82296e-05
+0.265625 0.265625 0.00000      -7.01353e-05
+0.265625 0.00000 0.265625      -9.66620e-05
+0.265625 0.515625 0.00000      0.000111798
+0.265625 0.00000 0.515625      4.88063e-05
+0.265625 0.765625 0.00000      -6.75355e-05
+0.265625 0.00000 0.765625      -7.84425e-05
+0.515625 0.265625 0.00000      -0.000111934
+0.515625 0.00000 0.265625      -0.000242694
+0.515625 0.515625 0.00000      0.000800841
+0.515625 0.00000 0.515625      0.000505391
+0.515625 0.765625 0.00000      -0.000133271
+0.515625 0.00000 0.765625      -0.000148360
+0.765625 0.265625 0.00000      -8.91167e-05
+0.765625 0.00000 0.265625      -9.24145e-05
+0.765625 0.515625 0.00000      -7.45770e-05
+0.765625 0.00000 0.515625      -8.11949e-05
+0.765625 0.765625 0.00000      -8.35244e-05
+0.765625 0.00000 0.765625      -8.87256e-05
 
 Cycle 2:
 DEAL::Starting value 0.0250422
 DEAL::Convergence step 7 value 0
    Solved in 7 iterations.
-0.125000 0.125000 0.00000      0.000714896
-0.125000 0.00000 0.125000      0.000732888
-0.125000 0.250000 0.00000      4.21241e-05
-0.125000 0.00000 0.250000      6.67244e-05
-0.125000 0.375000 0.00000      -0.000451000
-0.125000 0.00000 0.375000      -0.000348082
-0.125000 0.500000 0.00000      -0.000540492
-0.125000 0.00000 0.500000      -0.000389894
-0.125000 0.625000 0.00000      -0.000333853
-0.125000 0.00000 0.625000      -0.000227459
-0.125000 0.750000 0.00000      -4.30475e-05
-0.125000 0.00000 0.750000      -3.85607e-05
-0.125000 0.875000 0.00000      0.000218364
-0.125000 0.00000 0.875000      0.000150624
-0.250000 0.125000 0.00000      -0.000102974
-0.250000 0.00000 0.125000      -0.000107813
-0.250000 0.250000 0.00000      -7.06290e-05
-0.250000 0.00000 0.250000      -8.58914e-05
-0.250000 0.375000 0.00000      -2.80802e-05
-0.250000 0.00000 0.375000      -5.09368e-05
-0.250000 0.500000 0.00000      -2.22203e-06
-0.250000 0.00000 0.500000      -2.71771e-05
-0.250000 0.625000 0.00000      -1.09299e-05
-0.250000 0.00000 0.625000      -3.17431e-05
-0.250000 0.750000 0.00000      -5.07980e-05
-0.250000 0.00000 0.750000      -6.33648e-05
-0.250000 0.875000 0.00000      -9.87957e-05
-0.250000 0.00000 0.875000      -0.000103668
-0.375000 0.125000 0.00000      -0.00101184
-0.375000 0.00000 0.125000      -0.00104895
-0.375000 0.250000 0.00000      -0.000213964
-0.375000 0.00000 0.250000      -0.000284205
-0.375000 0.375000 0.00000      0.000442810
-0.375000 0.00000 0.375000      0.000274644
-0.375000 0.500000 0.00000      0.000616793
-0.375000 0.00000 0.500000      0.000394006
-0.375000 0.625000 0.00000      0.000372231
-0.375000 0.00000 0.625000      0.000207540
-0.375000 0.750000 0.00000      -5.19075e-05
-0.375000 0.00000 0.750000      -8.76469e-05
-0.375000 0.875000 0.00000      -0.000484075
-0.375000 0.00000 0.875000      -0.000427601
-0.500000 0.125000 0.00000      -0.00131821
-0.500000 0.00000 0.125000      -0.00136606
-0.500000 0.250000 0.00000      -0.000273020
-0.500000 0.00000 0.250000      -0.000364090
-0.500000 0.375000 0.00000      0.000588192
-0.500000 0.00000 0.375000      0.000370504
-0.500000 0.500000 0.00000      0.000816320
-0.500000 0.00000 0.500000      0.000527723
-0.500000 0.625000 0.00000      0.000496202
-0.500000 0.00000 0.625000      0.000283291
-0.500000 0.750000 0.00000      -5.86240e-05
-0.500000 0.00000 0.750000      -0.000103628
-0.500000 0.875000 0.00000      -0.000627027
-0.500000 0.00000 0.875000      -0.000552959
-0.625000 0.125000 0.00000      -0.000926278
-0.625000 0.00000 0.125000      -0.000960435
-0.625000 0.250000 0.00000      -0.000220905
-0.625000 0.00000 0.250000      -0.000286844
-0.625000 0.375000 0.00000      0.000368717
-0.625000 0.00000 0.375000      0.000216623
-0.625000 0.500000 0.00000      0.000530333
-0.625000 0.00000 0.500000      0.000330746
-0.625000 0.625000 0.00000      0.000314226
-0.625000 0.00000 0.625000      0.000166179
-0.625000 0.750000 0.00000      -6.93567e-05
-0.625000 0.00000 0.750000      -0.000103178
-0.625000 0.875000 0.00000      -0.000466912
-0.625000 0.00000 0.875000      -0.000419439
-0.750000 0.125000 0.00000      -0.000179286
-0.750000 0.00000 0.125000      -0.000185378
-0.750000 0.250000 0.00000      -0.000104709
-0.750000 0.00000 0.250000      -0.000119333
-0.750000 0.375000 0.00000      -3.05754e-05
-0.750000 0.00000 0.375000      -5.60127e-05
-0.750000 0.500000 0.00000      -1.89151e-06
-0.750000 0.00000 0.500000      -2.99975e-05
-0.750000 0.625000 0.00000      -2.44780e-05
-0.750000 0.00000 0.625000      -4.77630e-05
-0.750000 0.750000 0.00000      -7.84822e-05
-0.750000 0.00000 0.750000      -9.01734e-05
-0.750000 0.875000 0.00000      -0.000142966
-0.750000 0.00000 0.875000      -0.000142512
-0.875000 0.125000 0.00000      0.000456000
-0.875000 0.00000 0.125000      0.000475880
-0.875000 0.250000 0.00000      1.30647e-06
-0.875000 0.00000 0.250000      3.27430e-05
-0.875000 0.375000 0.00000      -0.000364914
-0.875000 0.00000 0.375000      -0.000280671
-0.875000 0.500000 0.00000      -0.000453633
-0.875000 0.00000 0.500000      -0.000336703
-0.875000 0.625000 0.00000      -0.000312918
-0.875000 0.00000 0.625000      -0.000230032
-0.875000 0.750000 0.00000      -7.87644e-05
-0.875000 0.00000 0.750000      -6.52054e-05
-0.875000 0.875000 0.00000      0.000142725
-0.875000 0.00000 0.875000      0.000106650
+0.132812 0.132812 0.00000      0.000646598
+0.132812 0.00000 0.132812      0.000661952
+0.132812 0.257812 0.00000      -1.72500e-06
+0.132812 0.00000 0.257812      2.56769e-05
+0.132812 0.382812 0.00000      -0.000451434
+0.132812 0.00000 0.382812      -0.000348091
+0.132812 0.507812 0.00000      -0.000513939
+0.132812 0.00000 0.507812      -0.000370142
+0.132812 0.632812 0.00000      -0.000305068
+0.132812 0.00000 0.632812      -0.000209183
+0.132812 0.757812 0.00000      -2.83548e-05
+0.132812 0.00000 0.757812      -3.09650e-05
+0.132812 0.882812 0.00000      0.000219467
+0.132812 0.00000 0.882812      0.000152889
+0.257812 0.132812 0.00000      -0.000165876
+0.257812 0.00000 0.132812      -0.000173604
+0.257812 0.257812 0.00000      -7.42803e-05
+0.257812 0.00000 0.257812      -9.44600e-05
+0.257812 0.382812 0.00000      1.16046e-05
+0.257812 0.00000 0.382812      -2.27690e-05
+0.257812 0.507812 0.00000      4.38914e-05
+0.257812 0.00000 0.507812      4.34583e-06
+0.257812 0.632812 0.00000      1.39282e-05
+0.257812 0.00000 0.632812      -1.65665e-05
+0.257812 0.757812 0.00000      -5.57917e-05
+0.257812 0.00000 0.757812      -6.88009e-05
+0.257812 0.882812 0.00000      -0.000130653
+0.257812 0.00000 0.882812      -0.000130155
+0.382812 0.132812 0.00000      -0.000995337
+0.382812 0.00000 0.132812      -0.00104505
+0.382812 0.257812 0.00000      -0.000164912
+0.382812 0.00000 0.257812      -0.000244631
+0.382812 0.382812 0.00000      0.000490431
+0.382812 0.00000 0.382812      0.000308958
+0.382812 0.507812 0.00000      0.000638078
+0.382812 0.00000 0.507812      0.000406714
+0.382812 0.632812 0.00000      0.000363883
+0.382812 0.00000 0.632812      0.000200367
+0.382812 0.757812 0.00000      -8.04884e-05
+0.382812 0.00000 0.757812      -0.000108485
+0.382812 0.882812 0.00000      -0.000527674
+0.382812 0.00000 0.882812      -0.000466275
+0.507812 0.132812 0.00000      -0.00125036
+0.507812 0.00000 0.132812      -0.00129986
+0.507812 0.257812 0.00000      -0.000203666
+0.507812 0.00000 0.257812      -0.000301637
+0.507812 0.382812 0.00000      0.000616785
+0.507812 0.00000 0.382812      0.000392221
+0.507812 0.507812 0.00000      0.000802115
+0.507812 0.00000 0.507812      0.000515750
+0.507812 0.632812 0.00000      0.000460227
+0.507812 0.00000 0.632812      0.000258189
+0.507812 0.757812 0.00000      -9.46937e-05
+0.507812 0.00000 0.757812      -0.000128448
+0.507812 0.882812 0.00000      -0.000656598
+0.507812 0.00000 0.882812      -0.000580311
+0.632812 0.132812 0.00000      -0.000846228
+0.632812 0.00000 0.132812      -0.000878904
+0.632812 0.257812 0.00000      -0.000169720
+0.632812 0.00000 0.257812      -0.000237495
+0.632812 0.382812 0.00000      0.000367537
+0.632812 0.00000 0.382812      0.000217482
+0.632812 0.507812 0.00000      0.000494364
+0.632812 0.00000 0.507812      0.000304766
+0.632812 0.632812 0.00000      0.000273516
+0.632812 0.00000 0.632812      0.000138834
+0.632812 0.757812 0.00000      -9.38392e-05
+0.632812 0.00000 0.757812      -0.000119091
+0.632812 0.882812 0.00000      -0.000470155
+0.632812 0.00000 0.882812      -0.000423658
+0.757812 0.132812 0.00000      -0.000128857
+0.757812 0.00000 0.132812      -0.000133514
+0.757812 0.257812 0.00000      -9.39660e-05
+0.757812 0.00000 0.257812      -0.000105571
+0.757812 0.382812 0.00000      -5.43154e-05
+0.757812 0.00000 0.382812      -7.19574e-05
+0.757812 0.507812 0.00000      -3.63902e-05
+0.757812 0.00000 0.507812      -5.47864e-05
+0.757812 0.632812 0.00000      -4.81446e-05
+0.757812 0.00000 0.632812      -6.35886e-05
+0.757812 0.757812 0.00000      -8.16497e-05
+0.757812 0.00000 0.757812      -9.07392e-05
+0.757812 0.882812 0.00000      -0.000121929
+0.757812 0.00000 0.882812      -0.000124171
+0.882812 0.132812 0.00000      0.000457188
+0.882812 0.00000 0.132812      0.000477415
+0.882812 0.257812 0.00000      -2.51848e-05
+0.882812 0.00000 0.257812      1.14568e-05
+0.882812 0.382812 0.00000      -0.000393924
+0.882812 0.00000 0.382812      -0.000301463
+0.882812 0.507812 0.00000      -0.000469132
+0.882812 0.00000 0.507812      -0.000344266
+0.882812 0.632812 0.00000      -0.000311403
+0.882812 0.00000 0.632812      -0.000228583
+0.882812 0.757812 0.00000      -6.69534e-05
+0.882812 0.00000 0.757812      -5.82098e-05
+0.882812 0.882812 0.00000      0.000172488
+0.882812 0.00000 0.882812      0.000134647
index 1cd6ca61ccf9c7c739e09c7f5186c6e06852da35..96e9a77053315fae79194c5ed37f5d0be07a0920 100644 (file)
@@ -5,27 +5,27 @@ Cycle 0:
 DEAL::Starting value 0.0140170
 DEAL::Convergence step 8 value 0
    Solved in 8 iterations.
-0.500000 0.00000       -0.00549678
+0.531250 0.00000       -0.00531794
 
 Cycle 1:
 DEAL::Starting value 0.0219353
 DEAL::Convergence step 8 value 0
    Solved in 8 iterations.
-0.250000 0.00000       -0.00130251
-0.500000 0.00000       -0.00550617
-0.750000 0.00000       -0.00183466
+0.265625 0.00000       -0.00177028
+0.515625 0.00000       -0.00542131
+0.765625 0.00000       -0.00148613
 
 Cycle 2:
 DEAL::Starting value 0.0268844
 DEAL::Convergence step 8 value 0
    Solved in 8 iterations.
-0.125000 0.00000       0.00141483
-0.250000 0.00000       -0.00124254
-0.375000 0.00000       -0.00429173
-0.500000 0.00000       -0.00545458
-0.625000 0.00000       -0.00429054
-0.750000 0.00000       -0.00185255
-0.875000 0.00000       0.000356647
+0.132812 0.00000       0.00130870
+0.257812 0.00000       -0.00146441
+0.382812 0.00000       -0.00443592
+0.507812 0.00000       -0.00544040
+0.632812 0.00000       -0.00415112
+0.757812 0.00000       -0.00167411
+0.882812 0.00000       0.000448975
 
 Testing for dim=3
 
@@ -33,131 +33,131 @@ Cycle 0:
 DEAL::Starting value 0.00632368
 DEAL::Convergence step 7 value 0
    Solved in 7 iterations.
-0.500000 0.500000 0.00000      0.000839383
-0.500000 0.00000 0.500000      0.000542250
+0.531250 0.531250 0.00000      0.000730918
+0.531250 0.00000 0.531250      0.000460360
 
 Cycle 1:
 DEAL::Starting value 0.0137083
 DEAL::Convergence step 7 value 0
    Solved in 7 iterations.
-0.250000 0.250000 0.00000      -7.02574e-05
-0.250000 0.00000 0.250000      -8.84955e-05
-0.250000 0.500000 0.00000      2.30224e-05
-0.250000 0.00000 0.500000      -1.58988e-05
-0.250000 0.750000 0.00000      -5.50224e-05
-0.250000 0.00000 0.750000      -6.78044e-05
-0.500000 0.250000 0.00000      -0.000245815
-0.500000 0.00000 0.250000      -0.000367442
-0.500000 0.500000 0.00000      0.000854616
-0.500000 0.00000 0.500000      0.000545283
-0.500000 0.750000 0.00000      -7.60057e-05
-0.500000 0.00000 0.750000      -0.000116497
-0.750000 0.250000 0.00000      -0.000105394
-0.750000 0.00000 0.250000      -0.000112911
-0.750000 0.500000 0.00000      -2.69468e-06
-0.750000 0.00000 0.500000      -3.29142e-05
-0.750000 0.750000 0.00000      -7.97320e-05
-0.750000 0.00000 0.750000      -8.82296e-05
+0.265625 0.265625 0.00000      -7.01353e-05
+0.265625 0.00000 0.265625      -9.66620e-05
+0.265625 0.515625 0.00000      0.000111798
+0.265625 0.00000 0.515625      4.88063e-05
+0.265625 0.765625 0.00000      -6.75355e-05
+0.265625 0.00000 0.765625      -7.84425e-05
+0.515625 0.265625 0.00000      -0.000111934
+0.515625 0.00000 0.265625      -0.000242694
+0.515625 0.515625 0.00000      0.000800841
+0.515625 0.00000 0.515625      0.000505391
+0.515625 0.765625 0.00000      -0.000133271
+0.515625 0.00000 0.765625      -0.000148360
+0.765625 0.265625 0.00000      -8.91167e-05
+0.765625 0.00000 0.265625      -9.24145e-05
+0.765625 0.515625 0.00000      -7.45770e-05
+0.765625 0.00000 0.515625      -8.11949e-05
+0.765625 0.765625 0.00000      -8.35244e-05
+0.765625 0.00000 0.765625      -8.87256e-05
 
 Cycle 2:
 DEAL::Starting value 0.0248658
 DEAL::Convergence step 8 value 0
    Solved in 8 iterations.
-0.125000 0.125000 0.00000      0.000714896
-0.125000 0.00000 0.125000      0.000732888
-0.125000 0.250000 0.00000      4.21241e-05
-0.125000 0.00000 0.250000      6.67244e-05
-0.125000 0.375000 0.00000      -0.000451000
-0.125000 0.00000 0.375000      -0.000348082
-0.125000 0.500000 0.00000      -0.000540492
-0.125000 0.00000 0.500000      -0.000389894
-0.125000 0.625000 0.00000      -0.000333853
-0.125000 0.00000 0.625000      -0.000227459
-0.125000 0.750000 0.00000      -4.30475e-05
-0.125000 0.00000 0.750000      -3.85607e-05
-0.125000 0.875000 0.00000      0.000218364
-0.125000 0.00000 0.875000      0.000150624
-0.250000 0.125000 0.00000      -0.000102974
-0.250000 0.00000 0.125000      -0.000107813
-0.250000 0.250000 0.00000      -7.06290e-05
-0.250000 0.00000 0.250000      -8.58914e-05
-0.250000 0.375000 0.00000      -2.80802e-05
-0.250000 0.00000 0.375000      -5.09368e-05
-0.250000 0.500000 0.00000      -2.22203e-06
-0.250000 0.00000 0.500000      -2.71771e-05
-0.250000 0.625000 0.00000      -1.09299e-05
-0.250000 0.00000 0.625000      -3.17431e-05
-0.250000 0.750000 0.00000      -5.07980e-05
-0.250000 0.00000 0.750000      -6.33648e-05
-0.250000 0.875000 0.00000      -9.87957e-05
-0.250000 0.00000 0.875000      -0.000103668
-0.375000 0.125000 0.00000      -0.00101184
-0.375000 0.00000 0.125000      -0.00104895
-0.375000 0.250000 0.00000      -0.000213964
-0.375000 0.00000 0.250000      -0.000284205
-0.375000 0.375000 0.00000      0.000442810
-0.375000 0.00000 0.375000      0.000274644
-0.375000 0.500000 0.00000      0.000616793
-0.375000 0.00000 0.500000      0.000394006
-0.375000 0.625000 0.00000      0.000372231
-0.375000 0.00000 0.625000      0.000207540
-0.375000 0.750000 0.00000      -5.19075e-05
-0.375000 0.00000 0.750000      -8.76469e-05
-0.375000 0.875000 0.00000      -0.000484075
-0.375000 0.00000 0.875000      -0.000427601
-0.500000 0.125000 0.00000      -0.00131821
-0.500000 0.00000 0.125000      -0.00136606
-0.500000 0.250000 0.00000      -0.000273020
-0.500000 0.00000 0.250000      -0.000364090
-0.500000 0.375000 0.00000      0.000588192
-0.500000 0.00000 0.375000      0.000370504
-0.500000 0.500000 0.00000      0.000816320
-0.500000 0.00000 0.500000      0.000527723
-0.500000 0.625000 0.00000      0.000496202
-0.500000 0.00000 0.625000      0.000283291
-0.500000 0.750000 0.00000      -5.86240e-05
-0.500000 0.00000 0.750000      -0.000103628
-0.500000 0.875000 0.00000      -0.000627027
-0.500000 0.00000 0.875000      -0.000552959
-0.625000 0.125000 0.00000      -0.000926278
-0.625000 0.00000 0.125000      -0.000960435
-0.625000 0.250000 0.00000      -0.000220905
-0.625000 0.00000 0.250000      -0.000286844
-0.625000 0.375000 0.00000      0.000368717
-0.625000 0.00000 0.375000      0.000216623
-0.625000 0.500000 0.00000      0.000530333
-0.625000 0.00000 0.500000      0.000330746
-0.625000 0.625000 0.00000      0.000314226
-0.625000 0.00000 0.625000      0.000166179
-0.625000 0.750000 0.00000      -6.93567e-05
-0.625000 0.00000 0.750000      -0.000103178
-0.625000 0.875000 0.00000      -0.000466912
-0.625000 0.00000 0.875000      -0.000419439
-0.750000 0.125000 0.00000      -0.000179286
-0.750000 0.00000 0.125000      -0.000185378
-0.750000 0.250000 0.00000      -0.000104709
-0.750000 0.00000 0.250000      -0.000119333
-0.750000 0.375000 0.00000      -3.05754e-05
-0.750000 0.00000 0.375000      -5.60127e-05
-0.750000 0.500000 0.00000      -1.89151e-06
-0.750000 0.00000 0.500000      -2.99975e-05
-0.750000 0.625000 0.00000      -2.44780e-05
-0.750000 0.00000 0.625000      -4.77630e-05
-0.750000 0.750000 0.00000      -7.84822e-05
-0.750000 0.00000 0.750000      -9.01734e-05
-0.750000 0.875000 0.00000      -0.000142966
-0.750000 0.00000 0.875000      -0.000142512
-0.875000 0.125000 0.00000      0.000456000
-0.875000 0.00000 0.125000      0.000475880
-0.875000 0.250000 0.00000      1.30647e-06
-0.875000 0.00000 0.250000      3.27430e-05
-0.875000 0.375000 0.00000      -0.000364914
-0.875000 0.00000 0.375000      -0.000280671
-0.875000 0.500000 0.00000      -0.000453633
-0.875000 0.00000 0.500000      -0.000336703
-0.875000 0.625000 0.00000      -0.000312918
-0.875000 0.00000 0.625000      -0.000230032
-0.875000 0.750000 0.00000      -7.87644e-05
-0.875000 0.00000 0.750000      -6.52054e-05
-0.875000 0.875000 0.00000      0.000142725
-0.875000 0.00000 0.875000      0.000106650
+0.132812 0.132812 0.00000      0.000646598
+0.132812 0.00000 0.132812      0.000661952
+0.132812 0.257812 0.00000      -1.72500e-06
+0.132812 0.00000 0.257812      2.56769e-05
+0.132812 0.382812 0.00000      -0.000451434
+0.132812 0.00000 0.382812      -0.000348091
+0.132812 0.507812 0.00000      -0.000513939
+0.132812 0.00000 0.507812      -0.000370142
+0.132812 0.632812 0.00000      -0.000305068
+0.132812 0.00000 0.632812      -0.000209183
+0.132812 0.757812 0.00000      -2.83548e-05
+0.132812 0.00000 0.757812      -3.09650e-05
+0.132812 0.882812 0.00000      0.000219467
+0.132812 0.00000 0.882812      0.000152889
+0.257812 0.132812 0.00000      -0.000165876
+0.257812 0.00000 0.132812      -0.000173604
+0.257812 0.257812 0.00000      -7.42803e-05
+0.257812 0.00000 0.257812      -9.44600e-05
+0.257812 0.382812 0.00000      1.16046e-05
+0.257812 0.00000 0.382812      -2.27690e-05
+0.257812 0.507812 0.00000      4.38914e-05
+0.257812 0.00000 0.507812      4.34583e-06
+0.257812 0.632812 0.00000      1.39282e-05
+0.257812 0.00000 0.632812      -1.65665e-05
+0.257812 0.757812 0.00000      -5.57917e-05
+0.257812 0.00000 0.757812      -6.88009e-05
+0.257812 0.882812 0.00000      -0.000130653
+0.257812 0.00000 0.882812      -0.000130155
+0.382812 0.132812 0.00000      -0.000995337
+0.382812 0.00000 0.132812      -0.00104505
+0.382812 0.257812 0.00000      -0.000164912
+0.382812 0.00000 0.257812      -0.000244631
+0.382812 0.382812 0.00000      0.000490431
+0.382812 0.00000 0.382812      0.000308958
+0.382812 0.507812 0.00000      0.000638078
+0.382812 0.00000 0.507812      0.000406714
+0.382812 0.632812 0.00000      0.000363883
+0.382812 0.00000 0.632812      0.000200367
+0.382812 0.757812 0.00000      -8.04884e-05
+0.382812 0.00000 0.757812      -0.000108485
+0.382812 0.882812 0.00000      -0.000527674
+0.382812 0.00000 0.882812      -0.000466275
+0.507812 0.132812 0.00000      -0.00125036
+0.507812 0.00000 0.132812      -0.00129986
+0.507812 0.257812 0.00000      -0.000203666
+0.507812 0.00000 0.257812      -0.000301637
+0.507812 0.382812 0.00000      0.000616785
+0.507812 0.00000 0.382812      0.000392221
+0.507812 0.507812 0.00000      0.000802115
+0.507812 0.00000 0.507812      0.000515750
+0.507812 0.632812 0.00000      0.000460227
+0.507812 0.00000 0.632812      0.000258189
+0.507812 0.757812 0.00000      -9.46937e-05
+0.507812 0.00000 0.757812      -0.000128448
+0.507812 0.882812 0.00000      -0.000656598
+0.507812 0.00000 0.882812      -0.000580311
+0.632812 0.132812 0.00000      -0.000846228
+0.632812 0.00000 0.132812      -0.000878904
+0.632812 0.257812 0.00000      -0.000169720
+0.632812 0.00000 0.257812      -0.000237495
+0.632812 0.382812 0.00000      0.000367537
+0.632812 0.00000 0.382812      0.000217482
+0.632812 0.507812 0.00000      0.000494364
+0.632812 0.00000 0.507812      0.000304766
+0.632812 0.632812 0.00000      0.000273516
+0.632812 0.00000 0.632812      0.000138834
+0.632812 0.757812 0.00000      -9.38392e-05
+0.632812 0.00000 0.757812      -0.000119091
+0.632812 0.882812 0.00000      -0.000470155
+0.632812 0.00000 0.882812      -0.000423658
+0.757812 0.132812 0.00000      -0.000128857
+0.757812 0.00000 0.132812      -0.000133514
+0.757812 0.257812 0.00000      -9.39660e-05
+0.757812 0.00000 0.257812      -0.000105571
+0.757812 0.382812 0.00000      -5.43154e-05
+0.757812 0.00000 0.382812      -7.19574e-05
+0.757812 0.507812 0.00000      -3.63902e-05
+0.757812 0.00000 0.507812      -5.47864e-05
+0.757812 0.632812 0.00000      -4.81446e-05
+0.757812 0.00000 0.632812      -6.35886e-05
+0.757812 0.757812 0.00000      -8.16497e-05
+0.757812 0.00000 0.757812      -9.07392e-05
+0.757812 0.882812 0.00000      -0.000121929
+0.757812 0.00000 0.882812      -0.000124171
+0.882812 0.132812 0.00000      0.000457188
+0.882812 0.00000 0.132812      0.000477415
+0.882812 0.257812 0.00000      -2.51848e-05
+0.882812 0.00000 0.257812      1.14568e-05
+0.882812 0.382812 0.00000      -0.000393924
+0.882812 0.00000 0.382812      -0.000301463
+0.882812 0.507812 0.00000      -0.000469132
+0.882812 0.00000 0.507812      -0.000344266
+0.882812 0.632812 0.00000      -0.000311403
+0.882812 0.00000 0.632812      -0.000228583
+0.882812 0.757812 0.00000      -6.69534e-05
+0.882812 0.00000 0.757812      -5.82098e-05
+0.882812 0.882812 0.00000      0.000172488
+0.882812 0.00000 0.882812      0.000134647
index ea81594e516d6421231a1b3c388533b4fc036951..39e50ca71c8504563df190c6ba057c82e23620aa 100644 (file)
@@ -5,27 +5,27 @@ Cycle 0:
 DEAL::Starting value 0.0139631
 DEAL::Convergence step 8 value 0
    Solved in 8 iterations.
-0.500000 0.00000       -0.00549678
+0.531250 0.00000       -0.00531794
 
 Cycle 1:
 DEAL::Starting value 0.0222274
 DEAL::Convergence step 8 value 0
    Solved in 8 iterations.
-0.250000 0.00000       -0.00130251
-0.500000 0.00000       -0.00550617
-0.750000 0.00000       -0.00183466
+0.265625 0.00000       -0.00177028
+0.515625 0.00000       -0.00542131
+0.765625 0.00000       -0.00148613
 
 Cycle 2:
 DEAL::Starting value 0.0268966
 DEAL::Convergence step 8 value 0
    Solved in 8 iterations.
-0.125000 0.00000       0.00141483
-0.250000 0.00000       -0.00124254
-0.375000 0.00000       -0.00429173
-0.500000 0.00000       -0.00545458
-0.625000 0.00000       -0.00429054
-0.750000 0.00000       -0.00185255
-0.875000 0.00000       0.000356647
+0.132812 0.00000       0.00130870
+0.257812 0.00000       -0.00146441
+0.382812 0.00000       -0.00443592
+0.507812 0.00000       -0.00544040
+0.632812 0.00000       -0.00415112
+0.757812 0.00000       -0.00167411
+0.882812 0.00000       0.000448975
 
 Testing for dim=3
 
@@ -33,131 +33,131 @@ Cycle 0:
 DEAL::Starting value 0.00630053
 DEAL::Convergence step 7 value 0
    Solved in 7 iterations.
-0.500000 0.500000 0.00000      0.000839383
-0.500000 0.00000 0.500000      0.000542250
+0.531250 0.531250 0.00000      0.000730918
+0.531250 0.00000 0.531250      0.000460360
 
 Cycle 1:
 DEAL::Starting value 0.0135106
 DEAL::Convergence step 7 value 0
    Solved in 7 iterations.
-0.250000 0.250000 0.00000      -7.02574e-05
-0.250000 0.00000 0.250000      -8.84955e-05
-0.250000 0.500000 0.00000      2.30224e-05
-0.250000 0.00000 0.500000      -1.58988e-05
-0.250000 0.750000 0.00000      -5.50224e-05
-0.250000 0.00000 0.750000      -6.78044e-05
-0.500000 0.250000 0.00000      -0.000245815
-0.500000 0.00000 0.250000      -0.000367442
-0.500000 0.500000 0.00000      0.000854616
-0.500000 0.00000 0.500000      0.000545283
-0.500000 0.750000 0.00000      -7.60057e-05
-0.500000 0.00000 0.750000      -0.000116497
-0.750000 0.250000 0.00000      -0.000105394
-0.750000 0.00000 0.250000      -0.000112911
-0.750000 0.500000 0.00000      -2.69468e-06
-0.750000 0.00000 0.500000      -3.29142e-05
-0.750000 0.750000 0.00000      -7.97320e-05
-0.750000 0.00000 0.750000      -8.82296e-05
+0.265625 0.265625 0.00000      -7.01353e-05
+0.265625 0.00000 0.265625      -9.66620e-05
+0.265625 0.515625 0.00000      0.000111798
+0.265625 0.00000 0.515625      4.88063e-05
+0.265625 0.765625 0.00000      -6.75355e-05
+0.265625 0.00000 0.765625      -7.84425e-05
+0.515625 0.265625 0.00000      -0.000111934
+0.515625 0.00000 0.265625      -0.000242694
+0.515625 0.515625 0.00000      0.000800841
+0.515625 0.00000 0.515625      0.000505391
+0.515625 0.765625 0.00000      -0.000133271
+0.515625 0.00000 0.765625      -0.000148360
+0.765625 0.265625 0.00000      -8.91167e-05
+0.765625 0.00000 0.265625      -9.24145e-05
+0.765625 0.515625 0.00000      -7.45770e-05
+0.765625 0.00000 0.515625      -8.11949e-05
+0.765625 0.765625 0.00000      -8.35244e-05
+0.765625 0.00000 0.765625      -8.87256e-05
 
 Cycle 2:
 DEAL::Starting value 0.0248824
 DEAL::Convergence step 8 value 0
    Solved in 8 iterations.
-0.125000 0.125000 0.00000      0.000714896
-0.125000 0.00000 0.125000      0.000732888
-0.125000 0.250000 0.00000      4.21241e-05
-0.125000 0.00000 0.250000      6.67244e-05
-0.125000 0.375000 0.00000      -0.000451000
-0.125000 0.00000 0.375000      -0.000348082
-0.125000 0.500000 0.00000      -0.000540492
-0.125000 0.00000 0.500000      -0.000389894
-0.125000 0.625000 0.00000      -0.000333853
-0.125000 0.00000 0.625000      -0.000227459
-0.125000 0.750000 0.00000      -4.30475e-05
-0.125000 0.00000 0.750000      -3.85607e-05
-0.125000 0.875000 0.00000      0.000218364
-0.125000 0.00000 0.875000      0.000150624
-0.250000 0.125000 0.00000      -0.000102974
-0.250000 0.00000 0.125000      -0.000107813
-0.250000 0.250000 0.00000      -7.06290e-05
-0.250000 0.00000 0.250000      -8.58914e-05
-0.250000 0.375000 0.00000      -2.80802e-05
-0.250000 0.00000 0.375000      -5.09368e-05
-0.250000 0.500000 0.00000      -2.22203e-06
-0.250000 0.00000 0.500000      -2.71771e-05
-0.250000 0.625000 0.00000      -1.09299e-05
-0.250000 0.00000 0.625000      -3.17431e-05
-0.250000 0.750000 0.00000      -5.07980e-05
-0.250000 0.00000 0.750000      -6.33648e-05
-0.250000 0.875000 0.00000      -9.87957e-05
-0.250000 0.00000 0.875000      -0.000103668
-0.375000 0.125000 0.00000      -0.00101184
-0.375000 0.00000 0.125000      -0.00104895
-0.375000 0.250000 0.00000      -0.000213964
-0.375000 0.00000 0.250000      -0.000284205
-0.375000 0.375000 0.00000      0.000442810
-0.375000 0.00000 0.375000      0.000274644
-0.375000 0.500000 0.00000      0.000616793
-0.375000 0.00000 0.500000      0.000394006
-0.375000 0.625000 0.00000      0.000372231
-0.375000 0.00000 0.625000      0.000207540
-0.375000 0.750000 0.00000      -5.19075e-05
-0.375000 0.00000 0.750000      -8.76469e-05
-0.375000 0.875000 0.00000      -0.000484075
-0.375000 0.00000 0.875000      -0.000427601
-0.500000 0.125000 0.00000      -0.00131821
-0.500000 0.00000 0.125000      -0.00136606
-0.500000 0.250000 0.00000      -0.000273020
-0.500000 0.00000 0.250000      -0.000364090
-0.500000 0.375000 0.00000      0.000588192
-0.500000 0.00000 0.375000      0.000370504
-0.500000 0.500000 0.00000      0.000816320
-0.500000 0.00000 0.500000      0.000527723
-0.500000 0.625000 0.00000      0.000496202
-0.500000 0.00000 0.625000      0.000283291
-0.500000 0.750000 0.00000      -5.86240e-05
-0.500000 0.00000 0.750000      -0.000103628
-0.500000 0.875000 0.00000      -0.000627027
-0.500000 0.00000 0.875000      -0.000552959
-0.625000 0.125000 0.00000      -0.000926278
-0.625000 0.00000 0.125000      -0.000960435
-0.625000 0.250000 0.00000      -0.000220905
-0.625000 0.00000 0.250000      -0.000286844
-0.625000 0.375000 0.00000      0.000368717
-0.625000 0.00000 0.375000      0.000216623
-0.625000 0.500000 0.00000      0.000530333
-0.625000 0.00000 0.500000      0.000330746
-0.625000 0.625000 0.00000      0.000314226
-0.625000 0.00000 0.625000      0.000166179
-0.625000 0.750000 0.00000      -6.93567e-05
-0.625000 0.00000 0.750000      -0.000103178
-0.625000 0.875000 0.00000      -0.000466912
-0.625000 0.00000 0.875000      -0.000419439
-0.750000 0.125000 0.00000      -0.000179286
-0.750000 0.00000 0.125000      -0.000185378
-0.750000 0.250000 0.00000      -0.000104709
-0.750000 0.00000 0.250000      -0.000119333
-0.750000 0.375000 0.00000      -3.05754e-05
-0.750000 0.00000 0.375000      -5.60127e-05
-0.750000 0.500000 0.00000      -1.89151e-06
-0.750000 0.00000 0.500000      -2.99975e-05
-0.750000 0.625000 0.00000      -2.44780e-05
-0.750000 0.00000 0.625000      -4.77630e-05
-0.750000 0.750000 0.00000      -7.84822e-05
-0.750000 0.00000 0.750000      -9.01734e-05
-0.750000 0.875000 0.00000      -0.000142966
-0.750000 0.00000 0.875000      -0.000142512
-0.875000 0.125000 0.00000      0.000456000
-0.875000 0.00000 0.125000      0.000475880
-0.875000 0.250000 0.00000      1.30647e-06
-0.875000 0.00000 0.250000      3.27430e-05
-0.875000 0.375000 0.00000      -0.000364914
-0.875000 0.00000 0.375000      -0.000280671
-0.875000 0.500000 0.00000      -0.000453633
-0.875000 0.00000 0.500000      -0.000336703
-0.875000 0.625000 0.00000      -0.000312918
-0.875000 0.00000 0.625000      -0.000230032
-0.875000 0.750000 0.00000      -7.87644e-05
-0.875000 0.00000 0.750000      -6.52054e-05
-0.875000 0.875000 0.00000      0.000142725
-0.875000 0.00000 0.875000      0.000106650
+0.132812 0.132812 0.00000      0.000646598
+0.132812 0.00000 0.132812      0.000661952
+0.132812 0.257812 0.00000      -1.72500e-06
+0.132812 0.00000 0.257812      2.56769e-05
+0.132812 0.382812 0.00000      -0.000451434
+0.132812 0.00000 0.382812      -0.000348091
+0.132812 0.507812 0.00000      -0.000513939
+0.132812 0.00000 0.507812      -0.000370142
+0.132812 0.632812 0.00000      -0.000305068
+0.132812 0.00000 0.632812      -0.000209183
+0.132812 0.757812 0.00000      -2.83548e-05
+0.132812 0.00000 0.757812      -3.09650e-05
+0.132812 0.882812 0.00000      0.000219467
+0.132812 0.00000 0.882812      0.000152889
+0.257812 0.132812 0.00000      -0.000165876
+0.257812 0.00000 0.132812      -0.000173604
+0.257812 0.257812 0.00000      -7.42803e-05
+0.257812 0.00000 0.257812      -9.44600e-05
+0.257812 0.382812 0.00000      1.16046e-05
+0.257812 0.00000 0.382812      -2.27690e-05
+0.257812 0.507812 0.00000      4.38914e-05
+0.257812 0.00000 0.507812      4.34583e-06
+0.257812 0.632812 0.00000      1.39282e-05
+0.257812 0.00000 0.632812      -1.65665e-05
+0.257812 0.757812 0.00000      -5.57917e-05
+0.257812 0.00000 0.757812      -6.88009e-05
+0.257812 0.882812 0.00000      -0.000130653
+0.257812 0.00000 0.882812      -0.000130155
+0.382812 0.132812 0.00000      -0.000995337
+0.382812 0.00000 0.132812      -0.00104505
+0.382812 0.257812 0.00000      -0.000164912
+0.382812 0.00000 0.257812      -0.000244631
+0.382812 0.382812 0.00000      0.000490431
+0.382812 0.00000 0.382812      0.000308958
+0.382812 0.507812 0.00000      0.000638078
+0.382812 0.00000 0.507812      0.000406714
+0.382812 0.632812 0.00000      0.000363883
+0.382812 0.00000 0.632812      0.000200367
+0.382812 0.757812 0.00000      -8.04884e-05
+0.382812 0.00000 0.757812      -0.000108485
+0.382812 0.882812 0.00000      -0.000527674
+0.382812 0.00000 0.882812      -0.000466275
+0.507812 0.132812 0.00000      -0.00125036
+0.507812 0.00000 0.132812      -0.00129986
+0.507812 0.257812 0.00000      -0.000203666
+0.507812 0.00000 0.257812      -0.000301637
+0.507812 0.382812 0.00000      0.000616785
+0.507812 0.00000 0.382812      0.000392221
+0.507812 0.507812 0.00000      0.000802115
+0.507812 0.00000 0.507812      0.000515750
+0.507812 0.632812 0.00000      0.000460227
+0.507812 0.00000 0.632812      0.000258189
+0.507812 0.757812 0.00000      -9.46937e-05
+0.507812 0.00000 0.757812      -0.000128448
+0.507812 0.882812 0.00000      -0.000656598
+0.507812 0.00000 0.882812      -0.000580311
+0.632812 0.132812 0.00000      -0.000846228
+0.632812 0.00000 0.132812      -0.000878904
+0.632812 0.257812 0.00000      -0.000169720
+0.632812 0.00000 0.257812      -0.000237495
+0.632812 0.382812 0.00000      0.000367537
+0.632812 0.00000 0.382812      0.000217482
+0.632812 0.507812 0.00000      0.000494364
+0.632812 0.00000 0.507812      0.000304766
+0.632812 0.632812 0.00000      0.000273516
+0.632812 0.00000 0.632812      0.000138834
+0.632812 0.757812 0.00000      -9.38392e-05
+0.632812 0.00000 0.757812      -0.000119091
+0.632812 0.882812 0.00000      -0.000470155
+0.632812 0.00000 0.882812      -0.000423658
+0.757812 0.132812 0.00000      -0.000128857
+0.757812 0.00000 0.132812      -0.000133514
+0.757812 0.257812 0.00000      -9.39660e-05
+0.757812 0.00000 0.257812      -0.000105571
+0.757812 0.382812 0.00000      -5.43154e-05
+0.757812 0.00000 0.382812      -7.19574e-05
+0.757812 0.507812 0.00000      -3.63902e-05
+0.757812 0.00000 0.507812      -5.47864e-05
+0.757812 0.632812 0.00000      -4.81446e-05
+0.757812 0.00000 0.632812      -6.35886e-05
+0.757812 0.757812 0.00000      -8.16497e-05
+0.757812 0.00000 0.757812      -9.07392e-05
+0.757812 0.882812 0.00000      -0.000121929
+0.757812 0.00000 0.882812      -0.000124171
+0.882812 0.132812 0.00000      0.000457188
+0.882812 0.00000 0.132812      0.000477415
+0.882812 0.257812 0.00000      -2.51848e-05
+0.882812 0.00000 0.257812      1.14568e-05
+0.882812 0.382812 0.00000      -0.000393924
+0.882812 0.00000 0.382812      -0.000301463
+0.882812 0.507812 0.00000      -0.000469132
+0.882812 0.00000 0.507812      -0.000344266
+0.882812 0.632812 0.00000      -0.000311403
+0.882812 0.00000 0.632812      -0.000228583
+0.882812 0.757812 0.00000      -6.69534e-05
+0.882812 0.00000 0.757812      -5.82098e-05
+0.882812 0.882812 0.00000      0.000172488
+0.882812 0.00000 0.882812      0.000134647
index 04e513b7756fd42f145e54c2c7b54f85dd42e5b1..31387eec908601a94eb23765d74334a3c2db3c08 100644 (file)
@@ -589,6 +589,9 @@ namespace Step22
     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(3);
@@ -596,9 +599,9 @@ namespace Step22
 
         Point<2> point1;
         point1(0)=0;
-        point1(1)=.5*(1.+1.*i/n_points);
+        point1(1)=.5*(1.+1.*i/n_points+eps);
         Point<2> point2;
-        point2(0)=.5*(1.+1.*i/n_points);
+        point2(0)=.5*(1.+1.*i/n_points+eps);
         point2(1)=0.;
 
         get_point_value (point1, 0, value1);
index aa60359f7aa061484456e712bf8e90707a08e058..53d4955a34a3ea44ca046f0082d200f609a4ffc6 100644 (file)
@@ -2,39 +2,39 @@ Refinement cycle 0
    Assembling...
    Computing preconditioner...
    Solving...
-0.00000 0.625000       -1.60000        0.00000
-0.00000 0.750000       -1.33333        0.00000
-0.00000 0.875000       -1.14286        0.00000
+0.00000 0.632812       -1.58049        0.00000
+0.00000 0.757812       -1.31971        0.00000
+0.00000 0.882812       -1.13281        0.00000
 
 Refinement cycle 1
    Assembling...
    Computing preconditioner...
    Solving...
-0.00000 0.562500       -1.77778        0.00000
-0.00000 0.625000       -1.60000        0.00000
-0.00000 0.687500       -1.45455        0.00000
-0.00000 0.750000       -1.33333        5.49010e-13
-0.00000 0.812500       -1.23077        2.80537e-13
-0.00000 0.875000       -1.14286        0.00000
-0.00000 0.937500       -1.06667        0.00000
+0.00000 0.566406       -1.76557        0.00000
+0.00000 0.628906       -1.59010        0.00000
+0.00000 0.691406       -1.44635        0.00000
+0.00000 0.753906       -1.32644        0.00000
+0.00000 0.816406       -1.22489        0.00000
+0.00000 0.878906       -1.13782        0.00000
+0.00000 0.941406       -1.06222        0.00000
 
 Refinement cycle 2
    Assembling...
    Computing preconditioner...
    Solving...
-0.00000 0.531250       -1.88235        0.00000
-0.00000 0.562500       -1.77778        0.00000
-0.00000 0.593750       -1.68421        0.00000
-0.00000 0.625000       -1.60000        0.00000
-0.00000 0.656250       -1.52381        0.00000
-0.00000 0.687500       -1.45455        0.00000
-0.00000 0.718750       -1.39130        1.95566e-17
-0.00000 0.750000       -1.33333        5.49063e-13
-0.00000 0.781250       -1.28000        0.00000
-0.00000 0.812500       -1.23077        2.80586e-13
-0.00000 0.843750       -1.18519        1.69980e-17
-0.00000 0.875000       -1.14286        0.00000
-0.00000 0.906250       -1.10357        0.00000
-0.00000 0.937500       -1.06667        0.00000
-0.00000 0.968750       -1.03214        0.00000
+0.00000 0.533203       -1.87547        0.00000
+0.00000 0.564453       -1.77163        0.00000
+0.00000 0.595703       -1.67869        0.00000
+0.00000 0.626953       -1.59504        0.00000
+0.00000 0.658203       -1.51928        0.00000
+0.00000 0.689453       -1.45044        0.00000
+0.00000 0.720703       -1.38753        0.00000
+0.00000 0.751953       -1.32988        0.00000
+0.00000 0.783203       -1.27680        0.00000
+0.00000 0.814453       -1.22782        0.00000
+0.00000 0.845703       -1.18244        0.00000
+0.00000 0.876953       -1.14033        0.00000
+0.00000 0.908203       -1.10120        0.00000
+0.00000 0.939453       -1.06444        0.00000
+0.00000 0.970703       -1.03006        0.00000
 
index aa60359f7aa061484456e712bf8e90707a08e058..53d4955a34a3ea44ca046f0082d200f609a4ffc6 100644 (file)
@@ -2,39 +2,39 @@ Refinement cycle 0
    Assembling...
    Computing preconditioner...
    Solving...
-0.00000 0.625000       -1.60000        0.00000
-0.00000 0.750000       -1.33333        0.00000
-0.00000 0.875000       -1.14286        0.00000
+0.00000 0.632812       -1.58049        0.00000
+0.00000 0.757812       -1.31971        0.00000
+0.00000 0.882812       -1.13281        0.00000
 
 Refinement cycle 1
    Assembling...
    Computing preconditioner...
    Solving...
-0.00000 0.562500       -1.77778        0.00000
-0.00000 0.625000       -1.60000        0.00000
-0.00000 0.687500       -1.45455        0.00000
-0.00000 0.750000       -1.33333        5.49010e-13
-0.00000 0.812500       -1.23077        2.80537e-13
-0.00000 0.875000       -1.14286        0.00000
-0.00000 0.937500       -1.06667        0.00000
+0.00000 0.566406       -1.76557        0.00000
+0.00000 0.628906       -1.59010        0.00000
+0.00000 0.691406       -1.44635        0.00000
+0.00000 0.753906       -1.32644        0.00000
+0.00000 0.816406       -1.22489        0.00000
+0.00000 0.878906       -1.13782        0.00000
+0.00000 0.941406       -1.06222        0.00000
 
 Refinement cycle 2
    Assembling...
    Computing preconditioner...
    Solving...
-0.00000 0.531250       -1.88235        0.00000
-0.00000 0.562500       -1.77778        0.00000
-0.00000 0.593750       -1.68421        0.00000
-0.00000 0.625000       -1.60000        0.00000
-0.00000 0.656250       -1.52381        0.00000
-0.00000 0.687500       -1.45455        0.00000
-0.00000 0.718750       -1.39130        1.95566e-17
-0.00000 0.750000       -1.33333        5.49063e-13
-0.00000 0.781250       -1.28000        0.00000
-0.00000 0.812500       -1.23077        2.80586e-13
-0.00000 0.843750       -1.18519        1.69980e-17
-0.00000 0.875000       -1.14286        0.00000
-0.00000 0.906250       -1.10357        0.00000
-0.00000 0.937500       -1.06667        0.00000
-0.00000 0.968750       -1.03214        0.00000
+0.00000 0.533203       -1.87547        0.00000
+0.00000 0.564453       -1.77163        0.00000
+0.00000 0.595703       -1.67869        0.00000
+0.00000 0.626953       -1.59504        0.00000
+0.00000 0.658203       -1.51928        0.00000
+0.00000 0.689453       -1.45044        0.00000
+0.00000 0.720703       -1.38753        0.00000
+0.00000 0.751953       -1.32988        0.00000
+0.00000 0.783203       -1.27680        0.00000
+0.00000 0.814453       -1.22782        0.00000
+0.00000 0.845703       -1.18244        0.00000
+0.00000 0.876953       -1.14033        0.00000
+0.00000 0.908203       -1.10120        0.00000
+0.00000 0.939453       -1.06444        0.00000
+0.00000 0.970703       -1.03006        0.00000
 
index aa60359f7aa061484456e712bf8e90707a08e058..53d4955a34a3ea44ca046f0082d200f609a4ffc6 100644 (file)
@@ -2,39 +2,39 @@ Refinement cycle 0
    Assembling...
    Computing preconditioner...
    Solving...
-0.00000 0.625000       -1.60000        0.00000
-0.00000 0.750000       -1.33333        0.00000
-0.00000 0.875000       -1.14286        0.00000
+0.00000 0.632812       -1.58049        0.00000
+0.00000 0.757812       -1.31971        0.00000
+0.00000 0.882812       -1.13281        0.00000
 
 Refinement cycle 1
    Assembling...
    Computing preconditioner...
    Solving...
-0.00000 0.562500       -1.77778        0.00000
-0.00000 0.625000       -1.60000        0.00000
-0.00000 0.687500       -1.45455        0.00000
-0.00000 0.750000       -1.33333        5.49010e-13
-0.00000 0.812500       -1.23077        2.80537e-13
-0.00000 0.875000       -1.14286        0.00000
-0.00000 0.937500       -1.06667        0.00000
+0.00000 0.566406       -1.76557        0.00000
+0.00000 0.628906       -1.59010        0.00000
+0.00000 0.691406       -1.44635        0.00000
+0.00000 0.753906       -1.32644        0.00000
+0.00000 0.816406       -1.22489        0.00000
+0.00000 0.878906       -1.13782        0.00000
+0.00000 0.941406       -1.06222        0.00000
 
 Refinement cycle 2
    Assembling...
    Computing preconditioner...
    Solving...
-0.00000 0.531250       -1.88235        0.00000
-0.00000 0.562500       -1.77778        0.00000
-0.00000 0.593750       -1.68421        0.00000
-0.00000 0.625000       -1.60000        0.00000
-0.00000 0.656250       -1.52381        0.00000
-0.00000 0.687500       -1.45455        0.00000
-0.00000 0.718750       -1.39130        1.95566e-17
-0.00000 0.750000       -1.33333        5.49063e-13
-0.00000 0.781250       -1.28000        0.00000
-0.00000 0.812500       -1.23077        2.80586e-13
-0.00000 0.843750       -1.18519        1.69980e-17
-0.00000 0.875000       -1.14286        0.00000
-0.00000 0.906250       -1.10357        0.00000
-0.00000 0.937500       -1.06667        0.00000
-0.00000 0.968750       -1.03214        0.00000
+0.00000 0.533203       -1.87547        0.00000
+0.00000 0.564453       -1.77163        0.00000
+0.00000 0.595703       -1.67869        0.00000
+0.00000 0.626953       -1.59504        0.00000
+0.00000 0.658203       -1.51928        0.00000
+0.00000 0.689453       -1.45044        0.00000
+0.00000 0.720703       -1.38753        0.00000
+0.00000 0.751953       -1.32988        0.00000
+0.00000 0.783203       -1.27680        0.00000
+0.00000 0.814453       -1.22782        0.00000
+0.00000 0.845703       -1.18244        0.00000
+0.00000 0.876953       -1.14033        0.00000
+0.00000 0.908203       -1.10120        0.00000
+0.00000 0.939453       -1.06444        0.00000
+0.00000 0.970703       -1.03006        0.00000
 
index d43d231f55b1c2aff85d38086a3f394250563bb2..7211af713562b2ae444fb4de6b68a5a9799012b4 100644 (file)
@@ -173,6 +173,7 @@ void check
   generate_grid(triangulation, orientation);
   set_periodicity(triangulation, reverse);
 
+  //first without refinement
   FE_Q<dim> fe(1);
   DoFHandler<dim> dof_handler(triangulation);
   dof_handler.distribute_dofs(fe);
@@ -230,6 +231,50 @@ void check
             << std::endl;
   AssertThrow(n_constraints == n_expected_constraints,
               ExcDimensionMismatch(n_constraints, n_expected_constraints));
+
+  //now refine and check if the neighboring faces are correctly found
+  typename Triangulation<dim>::active_cell_iterator cell;
+  for (cell = triangulation.begin_active(); cell != triangulation.end(); ++cell)
+    if (cell->is_locally_owned() && cell->center()(dim-1)>0)
+      cell->set_refine_flag();
+
+  triangulation.execute_coarsening_and_refinement();
+
+  typedef std::pair<typename Triangulation<dim>::cell_iterator, unsigned int> CellFace;
+  const typename std::map<CellFace, std::pair<CellFace, std::bitset<3> > > &face_map = triangulation.periodic_face_map;
+  typename std::map<CellFace, std::pair<CellFace, std::bitset<3> > >::const_iterator it;
+  int sum_of_pairs_local = face_map.size();
+  int sum_of_pairs_global;
+  MPI_Allreduce(&sum_of_pairs_local, &sum_of_pairs_global, 1, MPI_INT, MPI_SUM, triangulation.get_communicator());
+  Assert(sum_of_pairs_global>0, ExcInternalError());
+  for (it = face_map.begin(); it != face_map.end(); ++it)
+    {
+      const typename Triangulation<dim>::cell_iterator cell_1 = it->first.first;
+      const unsigned int face_no_1 = it->first.second;
+      const typename Triangulation<dim>::cell_iterator cell_2 = it->second.first.first;
+      const unsigned int face_no_2 = it->second.first.second;
+      const Point<dim> face_center_1 = cell_1->face(face_no_1)->center();
+      const Point<dim> face_center_2 = cell_2->face(face_no_2)->center();
+      Assert(std::min(std::abs(face_center_1(dim-1)-3.), std::abs(face_center_1(dim-1)+3.))<1.e-8, ExcInternalError());
+      Assert(std::min(std::abs(face_center_2(dim-1)-3.), std::abs(face_center_2(dim-1)+3.))<1.e-8, ExcInternalError());
+      if (cell_1->level()== cell_2->level())
+        for (unsigned int c = 0; c<dim-1; ++c)
+          if (std::abs(face_center_1(c)-face_center_2(c))>1.e-8)
+            {
+              std::cout << "face_center_1: " << face_center_1 << std::endl;
+              std::cout << "face_center_2: " << face_center_2 << std::endl;
+              for (auto it = triangulation.periodic_face_map.begin();
+                   it !=  triangulation.periodic_face_map.end(); ++it)
+                {
+                  std::cout << "The cell with center " << it->first.first->center()
+                            << " has on face " << it->first.second
+                            << " as neighbor the cell with center " << it->second.first.first->center()
+                            << " on face " << it->second.first.second
+                            << std::endl;
+                }
+              Assert(false, ExcInternalError());
+            }
+    }
 }
 
 int main (int argc, char *argv[])

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.