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);
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++)
{
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);
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
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
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
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
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
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
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);
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);
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
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
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
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);
<< 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[])