]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Set the correct boundary description on the first cell. This avoids NaN entries for...
authorkronbichler <kronbichler@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 18 Jan 2013 13:58:19 +0000 (13:58 +0000)
committerkronbichler <kronbichler@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 18 Jan 2013 13:58:19 +0000 (13:58 +0000)
git-svn-id: https://svn.dealii.org/trunk@28126 0785d39b-7218-0410-832d-ea1e28bc413d

42 files changed:
tests/fe/cell_similarity_01.cc
tests/fe/cell_similarity_02.cc
tests/fe/cell_similarity_03.cc
tests/fe/cell_similarity_04.cc
tests/fe/cell_similarity_05.cc
tests/fe/cell_similarity_05/cmp/generic
tests/fe/cell_similarity_06.cc
tests/fe/cell_similarity_06/cmp/generic
tests/fe/cell_similarity_07.cc
tests/fe/cell_similarity_07/cmp/generic
tests/fe/cell_similarity_08.cc
tests/fe/cell_similarity_08/cmp/generic
tests/fe/cell_similarity_09.cc
tests/fe/cell_similarity_10.cc
tests/fe/cell_similarity_dgp_monomial_01.cc
tests/fe/cell_similarity_dgp_monomial_02.cc
tests/fe/cell_similarity_dgp_monomial_03.cc
tests/fe/cell_similarity_dgp_monomial_04.cc
tests/fe/cell_similarity_dgp_monomial_05.cc
tests/fe/cell_similarity_dgp_monomial_05/cmp/generic
tests/fe/cell_similarity_dgp_monomial_06.cc
tests/fe/cell_similarity_dgp_monomial_06/cmp/generic
tests/fe/cell_similarity_dgp_monomial_07.cc
tests/fe/cell_similarity_dgp_monomial_07/cmp/generic
tests/fe/cell_similarity_dgp_monomial_08.cc
tests/fe/cell_similarity_dgp_monomial_08/cmp/generic
tests/fe/cell_similarity_dgp_monomial_09.cc
tests/fe/cell_similarity_dgp_monomial_10.cc
tests/fe/cell_similarity_dgp_nonparametric_01.cc
tests/fe/cell_similarity_dgp_nonparametric_02.cc
tests/fe/cell_similarity_dgp_nonparametric_03.cc
tests/fe/cell_similarity_dgp_nonparametric_04.cc
tests/fe/cell_similarity_dgp_nonparametric_05.cc
tests/fe/cell_similarity_dgp_nonparametric_05/cmp/generic
tests/fe/cell_similarity_dgp_nonparametric_06.cc
tests/fe/cell_similarity_dgp_nonparametric_06/cmp/generic
tests/fe/cell_similarity_dgp_nonparametric_07.cc
tests/fe/cell_similarity_dgp_nonparametric_07/cmp/generic
tests/fe/cell_similarity_dgp_nonparametric_08.cc
tests/fe/cell_similarity_dgp_nonparametric_08/cmp/generic
tests/fe/cell_similarity_dgp_nonparametric_09.cc
tests/fe/cell_similarity_dgp_nonparametric_10.cc

index 1c5341256a492f34959afee7168cb84e56c08464..a8d03db821f206ad8318842f9f06a90ae0e928c3 100644 (file)
@@ -148,7 +148,7 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary;
+  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
   tr.set_boundary (1, boundary);
 
                                   // set boundary id on cell 1
index 1aa167efb56c4928d3e033b7b10b6bc1fc548b9a..dea47632c76381f34d5b523c263e5ff8313185ff 100644 (file)
@@ -148,7 +148,7 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary;
+  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
   tr.set_boundary (1, boundary);
 
                                   // set boundary id on cell 1
index 7f3d6bbd0ad1872e4d53530ed65db597bc1ffbb8..cd502338027de3e606b42e1b421db9afb78e633a 100644 (file)
@@ -148,7 +148,7 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary;
+  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
   tr.set_boundary (1, boundary);
 
                                   // set boundary id on cell 1
index 270c4cd68e294853df55a887ee0c3733089da0cd..4942892a85e1db480a7934f87b064039ab51a41f 100644 (file)
@@ -148,7 +148,7 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary;
+  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
   tr.set_boundary (1, boundary);
 
                                   // set boundary id on cell 1
index 702edb4bb98d34cb7c9ce7b568602f5d3820d7ad..bdc67860345e23d2340cb85ffa7fa4aa8aee402e 100644 (file)
@@ -148,7 +148,7 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary;
+  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
   tr.set_boundary (1, boundary);
 
                                   // set boundary id on cell 1
index 441b8428134dae60dbe36818c916bd5f2a070973..ba39a73840137444e9c19ff55f36d5d39cd1c8df 100644 (file)
@@ -5,15 +5,15 @@ DEAL::Mass matrices are not equal.
 DEAL::Laplace matrices are not equal.
 DEAL::cell=0
 DEAL::mass_matrix:
-DEAL::0.2371 0.07167 0.06719 0.02582 
-DEAL::0.07167 0.04962 0.02582 0.03610 
-DEAL::0.06719 0.02582 0.03170 0.03162 
-DEAL::0.02582 0.03610 0.03162 0.09477 
+DEAL::0.2392 0.1111 0.1196 0.05556 
+DEAL::0.1111 0.2052 0.05556 0.1026 
+DEAL::0.1196 0.05556 0.2392 0.1111 
+DEAL::0.05556 0.1026 0.1111 0.2052 
 DEAL::laplace_matrix:
-DEAL::4.337 -0.2381 -2.153 -1.946 
-DEAL::-0.2381 0.9387 -0.5940 -0.1066 
-DEAL::-2.153 -0.5940 1.899 0.8483 
-DEAL::-1.946 -0.1066 0.8483 1.205 
+DEAL::1.128 -0.4478 -0.3766 -0.3037 
+DEAL::-0.4478 0.8880 -0.3037 -0.1365 
+DEAL::-0.3766 -0.3037 1.128 -0.4478 
+DEAL::-0.3037 -0.1365 -0.4478 0.8880 
 DEAL::cell=1
 DEAL::mass_matrix:
 DEAL::0.1111 0.05556 0.05556 0.02778 
@@ -31,23 +31,23 @@ DEAL::Mass matrices are not equal.
 DEAL::Laplace matrices are not equal.
 DEAL::cell=0
 DEAL::mass_matrix:
-DEAL::0.07572 0.02718 0.02414 0.009471 0.02414 0.009471 0.008538 0.003898 
-DEAL::0.02718 0.03298 0.009471 0.01375 0.009471 0.01375 0.003898 0.007053 
-DEAL::0.02414 0.009471 0.02083 0.01071 0.008538 0.003898 0.01001 0.006119 
-DEAL::0.009471 0.01375 0.01071 0.02201 0.003898 0.007053 0.006119 0.0144
-DEAL::0.02414 0.009471 0.008538 0.003898 0.02083 0.01071 0.01001 0.006119 
-DEAL::0.009471 0.01375 0.003898 0.007053 0.01071 0.02201 0.006119 0.0144
-DEAL::0.008538 0.003898 0.01001 0.006119 0.01001 0.006119 0.01922 0.01377 
-DEAL::0.003898 0.007053 0.006119 0.01447 0.006119 0.01447 0.01377 0.03585 
+DEAL::0.07543 0.03584 0.03771 0.01792 0.03771 0.01792 0.01886 0.008961 
+DEAL::0.03584 0.06794 0.01792 0.03397 0.01792 0.03397 0.008961 0.01699 
+DEAL::0.03771 0.01792 0.07543 0.03584 0.01886 0.008961 0.03771 0.01792 
+DEAL::0.01792 0.03397 0.03584 0.06794 0.008961 0.01699 0.01792 0.0339
+DEAL::0.03771 0.01792 0.01886 0.008961 0.07543 0.03584 0.03771 0.01792 
+DEAL::0.01792 0.03397 0.008961 0.01699 0.03584 0.06794 0.01792 0.0339
+DEAL::0.01886 0.008961 0.03771 0.01792 0.03771 0.01792 0.07543 0.03584 
+DEAL::0.008961 0.01699 0.01792 0.03397 0.01792 0.03397 0.03584 0.06794 
 DEAL::laplace_matrix:
-DEAL::0.7840 0.04931 0.04652 -0.2138 0.04652 -0.2138 -0.2672 -0.2315 
-DEAL::0.04931 0.3584 -0.04035 -0.04210 -0.04035 -0.04210 -0.09114 -0.1517 
-DEAL::0.04652 -0.04035 0.2547 -0.006678 -0.03606 -0.07879 -0.03423 -0.1051 
-DEAL::-0.2138 -0.04210 -0.006678 0.3292 -0.07879 -0.05266 0.002197 0.06269 
-DEAL::0.04652 -0.04035 -0.03606 -0.07879 0.2547 -0.006678 -0.03423 -0.1051 
-DEAL::-0.2138 -0.04210 -0.07879 -0.05266 -0.006678 0.3292 0.002197 0.06269 
-DEAL::-0.2672 -0.09114 -0.03423 0.002197 -0.03423 0.002197 0.3052 0.1173 
-DEAL::-0.2315 -0.1517 -0.1051 0.06269 -0.1051 0.06269 0.1173 0.3509 
+DEAL::1.090 -0.1204 -0.1704 -0.1779 -0.1704 -0.1779 -0.2178 -0.05474 
+DEAL::-0.1204 0.7538 -0.1779 -0.06053 -0.1779 -0.06053 -0.05474 -0.1017 
+DEAL::-0.1704 -0.1779 1.090 -0.1204 -0.2178 -0.05474 -0.1704 -0.1779 
+DEAL::-0.1779 -0.06053 -0.1204 0.7538 -0.05474 -0.1017 -0.1779 -0.06053 
+DEAL::-0.1704 -0.1779 -0.2178 -0.05474 1.090 -0.1204 -0.1704 -0.1779 
+DEAL::-0.1779 -0.06053 -0.05474 -0.1017 -0.1204 0.7538 -0.1779 -0.06053 
+DEAL::-0.2178 -0.05474 -0.1704 -0.1779 -0.1704 -0.1779 1.090 -0.1204 
+DEAL::-0.05474 -0.1017 -0.1779 -0.06053 -0.1779 -0.06053 -0.1204 0.7538 
 DEAL::cell=1
 DEAL::mass_matrix:
 DEAL::0.03704 0.01852 0.01852 0.009259 0.01852 0.009259 0.009259 0.004630 
index 7ab965f39e318233431f446fa8ed24ba0c5238aa..9d8856d27f56d83e7c76ae7e1e37a9518cd783ce 100644 (file)
@@ -148,7 +148,7 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary;
+  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
   tr.set_boundary (1, boundary);
 
                                   // set boundary id on cell 1
index 4f68dc5fabcd55feb36745bbb0397c31516cc96e..a0f11d3b7aa90a2d54d81f3c2be52e3ef151bdd2 100644 (file)
@@ -5,25 +5,25 @@ DEAL::Mass matrices are not equal.
 DEAL::Laplace matrices are not equal.
 DEAL::cell=0
 DEAL::mass_matrix:
-DEAL::0.02634 -0.007964 -0.007466 0.002869 0.03775 -0.01019 0.03675 -0.009193 0.05512 
-DEAL::-0.007964 0.005514 0.002869 -0.004011 -0.01019 0.003005 -0.004900 -0.002283 -0.01437 
-DEAL::-0.007466 0.002869 0.003522 -0.003513 -0.007887 -0.001288 -0.009193 1.858e-05 -0.01835 
-DEAL::0.002869 -0.004011 -0.003513 0.01053 -0.001288 0.01304 -0.002283 0.01403 0.02350 
-DEAL::0.03775 -0.01019 -0.007887 -0.001288 0.05973 -0.02295 0.05512 -0.01835 0.07355 
-DEAL::-0.01019 0.003005 -0.001288 0.01304 -0.02295 0.03209 -0.01437 0.02350 0.01827 
-DEAL::0.03675 -0.004900 -0.009193 -0.002283 0.05512 -0.01437 0.06371 -0.02295 0.08151 
-DEAL::-0.009193 -0.002283 1.858e-05 0.01403 -0.01835 0.02350 -0.02295 0.02810 0.01030 
-DEAL::0.05512 -0.01437 -0.01835 0.02350 0.07355 0.01827 0.08151 0.01030 0.1836 
+DEAL::0.02658 -0.01235 -0.01329 0.006173 0.02658 -0.01235 0.02846 -0.01423 0.02846 
+DEAL::-0.01235 0.02281 0.006173 -0.01140 -0.01235 0.02281 0.02092 -0.01046 0.02092 
+DEAL::-0.01329 0.006173 0.02658 -0.01235 0.02658 -0.01235 -0.01423 0.02846 0.02846 
+DEAL::0.006173 -0.01140 -0.01235 0.02281 -0.01235 0.02281 -0.01046 0.02092 0.02092 
+DEAL::0.02658 -0.01235 0.02658 -0.01235 0.1063 -0.04938 0.02846 0.02846 0.1139 
+DEAL::-0.01235 0.02281 -0.01235 0.02281 -0.04938 0.09122 0.02092 0.02092 0.08368 
+DEAL::0.02846 0.02092 -0.01423 -0.01046 0.02846 0.02092 0.09877 -0.04938 0.09877 
+DEAL::-0.01423 -0.01046 0.02846 0.02092 0.02846 0.02092 -0.04938 0.09877 0.09877 
+DEAL::0.02846 0.02092 0.02846 0.02092 0.1139 0.08368 0.09877 0.09877 0.3951 
 DEAL::laplace_matrix:
-DEAL::1.086 -0.0007442 0.2950 -0.2162 -0.5114 0.5990 -0.007372 1.084 -2.328 
-DEAL::-0.0007442 0.7950 -0.06600 -0.2749 0.1645 1.362 -1.277 0.2042 -0.9067 
-DEAL::0.2950 -0.06600 2.485 -0.9357 -4.333 0.8656 0.1882 5.057 -3.557 
-DEAL::-0.2162 -0.2749 -0.9357 0.9632 1.633 -0.8646 0.5105 -1.743 0.9273 
-DEAL::-0.5114 0.1645 -4.333 1.633 7.947 -1.400 -0.5993 -8.994 6.09
-DEAL::0.5990 1.362 0.8656 -0.8646 -1.400 4.384 -2.653 1.631 -3.923 
-DEAL::-0.007372 -1.277 0.1882 0.5105 -0.5993 -2.653 2.792 -0.1247 1.17
-DEAL::1.084 0.2042 5.057 -1.743 -8.994 1.631 -0.1247 12.09 -9.205 
-DEAL::-2.328 -0.9067 -3.557 0.9273 6.093 -3.923 1.171 -9.205 11.73 
+DEAL::1.189 -0.1638 -0.2365 -0.03375 0.1812 0.2293 0.09840 0.2461 -1.510 
+DEAL::-0.1638 0.9131 -0.03375 -0.2186 0.1547 0.2309 -0.1113 0.2229 -0.9941 
+DEAL::-0.2365 -0.03375 1.189 -0.1638 0.1812 0.2293 0.2461 0.09840 -1.510 
+DEAL::-0.03375 -0.2186 -0.1638 0.9131 0.1547 0.2309 0.2229 -0.1113 -0.9941 
+DEAL::0.1812 0.1547 0.1812 0.1547 0.8301 0.1281 -0.1432 -0.1432 -1.34
+DEAL::0.2293 0.2309 0.2293 0.2309 0.1281 1.241 -0.2245 -0.2245 -1.841 
+DEAL::0.09840 -0.1113 0.2461 0.2229 -0.1432 -0.2245 0.8657 -0.04194 -0.912
+DEAL::0.2461 0.2229 0.09840 -0.1113 -0.1432 -0.2245 -0.04194 0.8657 -0.9121 
+DEAL::-1.510 -0.9941 -1.510 -0.9941 -1.343 -1.841 -0.9121 -0.9121 10.02 
 DEAL::cell=1
 DEAL::mass_matrix:
 DEAL::0.01235 -0.006173 -0.006173 0.003086 0.01235 -0.006173 0.01235 -0.006173 0.01235 
@@ -51,61 +51,61 @@ DEAL::Mass matrices are not equal.
 DEAL::Laplace matrices are not equal.
 DEAL::cell=0
 DEAL::mass_matrix:
-DEAL::0.002805 -0.001006 -0.0008940 0.0003508 -0.0008940 0.0003508 0.0003162 -0.0001444 0.003821 -0.001311 0.003596 -0.001086 -0.001156 0.0004129 -0.001086 0.0003437 0.003821 -0.001311 -0.001156 0.0004129 0.005331 -0.001797 0.005020 -0.001485 0.005020 -0.001485 0.007068 
-DEAL::-0.001006 0.001221 0.0003508 -0.0005092 0.0003508 -0.0005092 -0.0001444 0.0002612 -0.001311 0.001424 0.0004298 -0.0003167 0.0004129 -0.0004958 -0.0003167 0.0002337 -0.001311 0.001424 0.0004129 -0.0004958 -0.001797 0.001857 0.0002261 -0.0001660 0.0002261 -0.0001660 0.0001202 
-DEAL::-0.0008940 0.0003508 0.0007714 -0.0003967 0.0003162 -0.0001444 -0.0003708 0.0002266 -0.0002452 -9.174e-05 -0.001086 0.0007495 -0.0001092 0.0001646 0.0003437 -0.0002883 -0.001156 0.0004129 0.0008012 -0.0003400 -0.0007088 0.0001457 -0.001485 0.0009223 -0.0006739 0.0001107 -0.001126 
-DEAL::0.0003508 -0.0005092 -0.0003967 0.0008152 -0.0001444 0.0002612 0.0002266 -0.0005358 -9.174e-05 0.0006122 -0.0003167 0.0008371 0.0001646 -0.0005491 0.0002337 -0.0006182 0.0004129 -0.0004958 -0.0003400 0.0005589 0.0001457 0.0001262 -0.0001660 0.0004378 0.001041 -0.0007690 0.0005436 
-DEAL::-0.0008940 0.0003508 0.0003162 -0.0001444 0.0007714 -0.0003967 -0.0003708 0.0002266 -0.001156 0.0004129 -0.001086 0.0003437 0.0008012 -0.0003400 0.0007495 -0.0002883 -0.0002452 -9.174e-05 -0.0001092 0.0001646 -0.0007088 0.0001457 -0.0006739 0.0001107 -0.001485 0.0009223 -0.001126 
-DEAL::0.0003508 -0.0005092 -0.0001444 0.0002612 -0.0003967 0.0008152 0.0002266 -0.0005358 0.0004129 -0.0004958 -0.0003167 0.0002337 -0.0003400 0.0005589 0.0008371 -0.0006182 -9.174e-05 0.0006122 0.0001646 -0.0005491 0.0001457 0.0001262 0.001041 -0.0007690 -0.0001660 0.0004378 0.0005436 
-DEAL::0.0003162 -0.0001444 -0.0003708 0.0002266 -0.0003708 0.0002266 0.0007118 -0.0005099 -0.0001092 0.0001646 0.0003437 -0.0002883 0.0006821 -0.0005666 -0.0002883 0.0004038 -0.0001092 0.0001646 0.0006821 -0.0005666 0.001146 -0.0008040 0.0001107 0.0002310 0.0001107 0.0002310 0.0006834 
-DEAL::-0.0001444 0.0002612 0.0002266 -0.0005358 0.0002266 -0.0005358 -0.0005099 0.001328 0.0001646 -0.0005491 0.0002337 -0.0006182 -0.0005666 0.001584 -0.0006182 0.001636 0.0001646 -0.0005491 -0.0005666 0.001584 -0.0008040 0.002070 -0.0007690 0.002035 -0.0007690 0.002035 0.002532 
-DEAL::0.003821 -0.001311 -0.0002452 -9.174e-05 -0.001156 0.0004129 -0.0001092 0.0001646 0.007152 -0.002806 0.005020 -0.0006739 -0.002530 0.001155 -0.001485 0.0001107 0.005331 -0.001797 -0.0007088 0.0001457 0.009245 -0.003303 0.007068 -0.001126 0.008691 -0.002750 0.01188 
-DEAL::-0.001311 0.001424 -9.174e-05 0.0006122 0.0004129 -0.0004958 0.0001646 -0.0005491 -0.002806 0.004073 0.0002261 0.001041 0.001155 -0.002090 -0.0001660 -0.0007690 -0.001797 0.001857 0.0001457 0.0001262 -0.003303 0.003967 0.0001202 0.0005436 0.002534 -0.001870 0.001328 
-DEAL::0.003596 0.0004298 -0.001086 -0.0003167 -0.001086 -0.0003167 0.0003437 0.0002337 0.005020 0.0002261 0.008052 -0.002806 -0.001485 -0.0001660 -0.002806 0.001155 0.005020 0.0002261 -0.001485 -0.0001660 0.007068 0.0001202 0.01049 -0.003303 0.01049 -0.003303 0.01438 
-DEAL::-0.001086 -0.0003167 0.0007495 0.0008371 0.0003437 0.0002337 -0.0002883 -0.0006182 -0.0006739 0.001041 -0.002806 0.003173 0.0001107 -0.0007690 0.001155 -0.001813 -0.001485 -0.0001660 0.0009223 0.0004378 -0.001126 0.0005436 -0.003303 0.002720 0.0007339 -0.001317 -0.001165 
-DEAL::-0.001156 0.0004129 -0.0001092 0.0001646 0.0008012 -0.0003400 0.0006821 -0.0005666 -0.002530 0.001155 -0.001485 0.0001107 0.002967 -0.001813 0.0009223 0.0002310 -0.0007088 0.0001457 0.001146 -0.0008040 0.0008737 -0.001317 -0.001126 0.0006834 -0.002750 0.002307 -0.0008857 
-DEAL::0.0004129 -0.0004958 0.0001646 -0.0005491 -0.0003400 0.0005589 -0.0005666 0.001584 0.001155 -0.002090 -0.0001660 -0.0007690 -0.001813 0.004286 0.0004378 0.002035 0.0001457 0.0001262 -0.0008040 0.002070 -0.001317 0.004393 0.0005436 0.002532 -0.001870 0.004946 0.006152 
-DEAL::-0.001086 -0.0003167 0.0003437 0.0002337 0.0007495 0.0008371 -0.0002883 -0.0006182 -0.001485 -0.0001660 -0.002806 0.001155 0.0009223 0.0004378 0.003173 -0.001813 -0.0006739 0.001041 0.0001107 -0.0007690 -0.001126 0.0005436 0.0007339 -0.001317 -0.003303 0.002720 -0.001165 
-DEAL::0.0003437 0.0002337 -0.0002883 -0.0006182 -0.0002883 -0.0006182 0.0004038 0.001636 0.0001107 -0.0007690 0.001155 -0.001813 0.0002310 0.002035 -0.001813 0.004079 0.0001107 -0.0007690 0.0002310 0.002035 0.0006834 0.002532 -0.001317 0.004532 -0.001317 0.004532 0.00643
-DEAL::0.003821 -0.001311 -0.001156 0.0004129 -0.0002452 -9.174e-05 -0.0001092 0.0001646 0.005331 -0.001797 0.005020 -0.001485 -0.0007088 0.0001457 -0.0006739 0.0001107 0.007152 -0.002806 -0.002530 0.001155 0.009245 -0.003303 0.008691 -0.002750 0.007068 -0.001126 0.01188 
-DEAL::-0.001311 0.001424 0.0004129 -0.0004958 -9.174e-05 0.0006122 0.0001646 -0.0005491 -0.001797 0.001857 0.0002261 -0.0001660 0.0001457 0.0001262 0.001041 -0.0007690 -0.002806 0.004073 0.001155 -0.002090 -0.003303 0.003967 0.002534 -0.001870 0.0001202 0.0005436 0.001328 
-DEAL::-0.001156 0.0004129 0.0008012 -0.0003400 -0.0001092 0.0001646 0.0006821 -0.0005666 -0.0007088 0.0001457 -0.001485 0.0009223 0.001146 -0.0008040 0.0001107 0.0002310 -0.002530 0.001155 0.002967 -0.001813 0.0008737 -0.001317 -0.002750 0.002307 -0.001126 0.0006834 -0.0008857 
-DEAL::0.0004129 -0.0004958 -0.0003400 0.0005589 0.0001646 -0.0005491 -0.0005666 0.001584 0.0001457 0.0001262 -0.0001660 0.0004378 -0.0008040 0.002070 -0.0007690 0.002035 0.001155 -0.002090 -0.001813 0.004286 -0.001317 0.004393 -0.001870 0.004946 0.0005436 0.002532 0.006152 
-DEAL::0.005331 -0.001797 -0.0007088 0.0001457 -0.0007088 0.0001457 0.001146 -0.0008040 0.009245 -0.003303 0.007068 -0.001126 0.0008737 -0.001317 -0.001126 0.0006834 0.009245 -0.003303 0.0008737 -0.001317 0.02024 -0.009239 0.01188 -0.0008857 0.01188 -0.0008857 0.02200 
-DEAL::-0.001797 0.001857 0.0001457 0.0001262 0.0001457 0.0001262 -0.0008040 0.002070 -0.003303 0.003967 0.0001202 0.0005436 -0.001317 0.004393 0.0005436 0.002532 -0.003303 0.003967 -0.001317 0.004393 -0.009239 0.01672 0.001328 0.006152 0.001328 0.006152 0.01496 
-DEAL::0.005020 0.0002261 -0.001485 -0.0001660 -0.0006739 0.001041 0.0001107 -0.0007690 0.007068 0.0001202 0.01049 -0.003303 -0.001126 0.0005436 0.0007339 -0.001317 0.008691 0.002534 -0.002750 -0.001870 0.01188 0.001328 0.02245 -0.009239 0.01438 -0.001165 0.02642 
-DEAL::-0.001485 -0.0001660 0.0009223 0.0004378 0.0001107 -0.0007690 0.0002310 0.002035 -0.001126 0.0005436 -0.003303 0.002720 0.0006834 0.002532 -0.001317 0.004532 -0.002750 -0.001870 0.002307 0.004946 -0.0008857 0.006152 -0.009239 0.01451 -0.001165 0.006432 0.01053 
-DEAL::0.005020 0.0002261 -0.0006739 0.001041 -0.001485 -0.0001660 0.0001107 -0.0007690 0.008691 0.002534 0.01049 0.0007339 -0.002750 -0.001870 -0.003303 -0.001317 0.007068 0.0001202 -0.001126 0.0005436 0.01188 0.001328 0.01438 -0.001165 0.02245 -0.009239 0.02642 
-DEAL::-0.001485 -0.0001660 0.0001107 -0.0007690 0.0009223 0.0004378 0.0002310 0.002035 -0.002750 -0.001870 -0.003303 -0.001317 0.002307 0.004946 0.002720 0.004532 -0.001126 0.0005436 0.0006834 0.002532 -0.0008857 0.006152 -0.001165 0.006432 -0.009239 0.01451 0.01053 
-DEAL::0.007068 0.0001202 -0.001126 0.0005436 -0.001126 0.0005436 0.0006834 0.002532 0.01188 0.001328 0.01438 -0.001165 -0.0008857 0.006152 -0.001165 0.006432 0.01188 0.001328 -0.0008857 0.006152 0.02200 0.01496 0.02642 0.01053 0.02642 0.01053 0.07391 
+DEAL::0.002794 -0.001328 -0.001397 0.0006638 -0.001397 0.0006638 0.0006984 -0.0003319 0.002794 -0.001328 0.002932 -0.001466 -0.001397 0.0006638 -0.001466 0.0007331 0.002794 -0.001328 -0.001397 0.0006638 0.002794 -0.001328 0.002932 -0.001466 0.002932 -0.001466 0.002932 
+DEAL::-0.001328 0.002516 0.0006638 -0.001258 0.0006638 -0.001258 -0.0003319 0.0006291 -0.001328 0.002516 0.002378 -0.001189 0.0006638 -0.001258 -0.001189 0.0005945 -0.001328 0.002516 0.0006638 -0.001258 -0.001328 0.002516 0.002378 -0.001189 0.002378 -0.001189 0.002378 
+DEAL::-0.001397 0.0006638 0.002794 -0.001328 0.0006984 -0.0003319 -0.001397 0.0006638 0.002794 -0.001328 -0.001466 0.002932 -0.001397 0.0006638 0.0007331 -0.001466 -0.001397 0.0006638 0.002794 -0.001328 0.002794 -0.001328 -0.001466 0.002932 0.002932 -0.001466 0.002932 
+DEAL::0.0006638 -0.001258 -0.001328 0.002516 -0.0003319 0.0006291 0.0006638 -0.001258 -0.001328 0.002516 -0.001189 0.002378 0.0006638 -0.001258 0.0005945 -0.001189 0.0006638 -0.001258 -0.001328 0.002516 -0.001328 0.002516 -0.001189 0.002378 0.002378 -0.001189 0.002378 
+DEAL::-0.001397 0.0006638 0.0006984 -0.0003319 0.002794 -0.001328 -0.001397 0.0006638 -0.001397 0.0006638 -0.001466 0.0007331 0.002794 -0.001328 0.002932 -0.001466 0.002794 -0.001328 -0.001397 0.0006638 0.002794 -0.001328 0.002932 -0.001466 -0.001466 0.002932 0.002932 
+DEAL::0.0006638 -0.001258 -0.0003319 0.0006291 -0.001328 0.002516 0.0006638 -0.001258 0.0006638 -0.001258 -0.001189 0.0005945 -0.001328 0.002516 0.002378 -0.001189 -0.001328 0.002516 0.0006638 -0.001258 -0.001328 0.002516 0.002378 -0.001189 -0.001189 0.002378 0.002378 
+DEAL::0.0006984 -0.0003319 -0.001397 0.0006638 -0.001397 0.0006638 0.002794 -0.001328 -0.001397 0.0006638 0.0007331 -0.001466 0.002794 -0.001328 -0.001466 0.002932 -0.001397 0.0006638 0.002794 -0.001328 0.002794 -0.001328 -0.001466 0.002932 -0.001466 0.002932 0.002932 
+DEAL::-0.0003319 0.0006291 0.0006638 -0.001258 0.0006638 -0.001258 -0.001328 0.002516 0.0006638 -0.001258 0.0005945 -0.001189 -0.001328 0.002516 -0.001189 0.002378 0.0006638 -0.001258 -0.001328 0.002516 -0.001328 0.002516 -0.001189 0.002378 -0.001189 0.002378 0.002378 
+DEAL::0.002794 -0.001328 0.002794 -0.001328 -0.001397 0.0006638 -0.001397 0.0006638 0.01117 -0.005310 0.002932 0.002932 -0.005587 0.002655 -0.001466 -0.001466 0.002794 -0.001328 0.002794 -0.001328 0.01117 -0.005310 0.002932 0.002932 0.01173 -0.005865 0.01173 
+DEAL::-0.001328 0.002516 -0.001328 0.002516 0.0006638 -0.001258 0.0006638 -0.001258 -0.005310 0.01007 0.002378 0.002378 0.002655 -0.005033 -0.001189 -0.001189 -0.001328 0.002516 -0.001328 0.002516 -0.005310 0.01007 0.002378 0.002378 0.009511 -0.004756 0.009511 
+DEAL::0.002932 0.002378 -0.001466 -0.001189 -0.001466 -0.001189 0.0007331 0.0005945 0.002932 0.002378 0.01062 -0.005310 -0.001466 -0.001189 -0.005310 0.002655 0.002932 0.002378 -0.001466 -0.001189 0.002932 0.002378 0.01062 -0.005310 0.01062 -0.005310 0.01062 
+DEAL::-0.001466 -0.001189 0.002932 0.002378 0.0007331 0.0005945 -0.001466 -0.001189 0.002932 0.002378 -0.005310 0.01062 -0.001466 -0.001189 0.002655 -0.005310 -0.001466 -0.001189 0.002932 0.002378 0.002932 0.002378 -0.005310 0.01062 0.01062 -0.005310 0.01062 
+DEAL::-0.001397 0.0006638 -0.001397 0.0006638 0.002794 -0.001328 0.002794 -0.001328 -0.005587 0.002655 -0.001466 -0.001466 0.01117 -0.005310 0.002932 0.002932 0.002794 -0.001328 0.002794 -0.001328 0.01117 -0.005310 0.002932 0.002932 -0.005865 0.01173 0.01173 
+DEAL::0.0006638 -0.001258 0.0006638 -0.001258 -0.001328 0.002516 -0.001328 0.002516 0.002655 -0.005033 -0.001189 -0.001189 -0.005310 0.01007 0.002378 0.002378 -0.001328 0.002516 -0.001328 0.002516 -0.005310 0.01007 0.002378 0.002378 -0.004756 0.009511 0.009511 
+DEAL::-0.001466 -0.001189 0.0007331 0.0005945 0.002932 0.002378 -0.001466 -0.001189 -0.001466 -0.001189 -0.005310 0.002655 0.002932 0.002378 0.01062 -0.005310 0.002932 0.002378 -0.001466 -0.001189 0.002932 0.002378 0.01062 -0.005310 -0.005310 0.01062 0.01062 
+DEAL::0.0007331 0.0005945 -0.001466 -0.001189 -0.001466 -0.001189 0.002932 0.002378 -0.001466 -0.001189 0.002655 -0.005310 0.002932 0.002378 -0.005310 0.01062 -0.001466 -0.001189 0.002932 0.002378 0.002932 0.002378 -0.005310 0.01062 -0.005310 0.01062 0.0106
+DEAL::0.002794 -0.001328 -0.001397 0.0006638 0.002794 -0.001328 -0.001397 0.0006638 0.002794 -0.001328 0.002932 -0.001466 0.002794 -0.001328 0.002932 -0.001466 0.01117 -0.005310 -0.005587 0.002655 0.01117 -0.005310 0.01173 -0.005865 0.002932 0.002932 0.01173 
+DEAL::-0.001328 0.002516 0.0006638 -0.001258 -0.001328 0.002516 0.0006638 -0.001258 -0.001328 0.002516 0.002378 -0.001189 -0.001328 0.002516 0.002378 -0.001189 -0.005310 0.01007 0.002655 -0.005033 -0.005310 0.01007 0.009511 -0.004756 0.002378 0.002378 0.009511 
+DEAL::-0.001397 0.0006638 0.002794 -0.001328 -0.001397 0.0006638 0.002794 -0.001328 0.002794 -0.001328 -0.001466 0.002932 0.002794 -0.001328 -0.001466 0.002932 -0.005587 0.002655 0.01117 -0.005310 0.01117 -0.005310 -0.005865 0.01173 0.002932 0.002932 0.01173 
+DEAL::0.0006638 -0.001258 -0.001328 0.002516 0.0006638 -0.001258 -0.001328 0.002516 -0.001328 0.002516 -0.001189 0.002378 -0.001328 0.002516 -0.001189 0.002378 0.002655 -0.005033 -0.005310 0.01007 -0.005310 0.01007 -0.004756 0.009511 0.002378 0.002378 0.009511 
+DEAL::0.002794 -0.001328 0.002794 -0.001328 0.002794 -0.001328 0.002794 -0.001328 0.01117 -0.005310 0.002932 0.002932 0.01117 -0.005310 0.002932 0.002932 0.01117 -0.005310 0.01117 -0.005310 0.04470 -0.02124 0.01173 0.01173 0.01173 0.01173 0.04692 
+DEAL::-0.001328 0.002516 -0.001328 0.002516 -0.001328 0.002516 -0.001328 0.002516 -0.005310 0.01007 0.002378 0.002378 -0.005310 0.01007 0.002378 0.002378 -0.005310 0.01007 -0.005310 0.01007 -0.02124 0.04026 0.009511 0.009511 0.009511 0.009511 0.03804 
+DEAL::0.002932 0.002378 -0.001466 -0.001189 0.002932 0.002378 -0.001466 -0.001189 0.002932 0.002378 0.01062 -0.005310 0.002932 0.002378 0.01062 -0.005310 0.01173 0.009511 -0.005865 -0.004756 0.01173 0.009511 0.04248 -0.02124 0.01062 0.01062 0.04248 
+DEAL::-0.001466 -0.001189 0.002932 0.002378 -0.001466 -0.001189 0.002932 0.002378 0.002932 0.002378 -0.005310 0.01062 0.002932 0.002378 -0.005310 0.01062 -0.005865 -0.004756 0.01173 0.009511 0.01173 0.009511 -0.02124 0.04248 0.01062 0.01062 0.04248 
+DEAL::0.002932 0.002378 0.002932 0.002378 -0.001466 -0.001189 -0.001466 -0.001189 0.01173 0.009511 0.01062 0.01062 -0.005865 -0.004756 -0.005310 -0.005310 0.002932 0.002378 0.002932 0.002378 0.01173 0.009511 0.01062 0.01062 0.04248 -0.02124 0.04248 
+DEAL::-0.001466 -0.001189 -0.001466 -0.001189 0.002932 0.002378 0.002932 0.002378 -0.005865 -0.004756 -0.005310 -0.005310 0.01173 0.009511 0.01062 0.01062 0.002932 0.002378 0.002932 0.002378 0.01173 0.009511 0.01062 0.01062 -0.02124 0.04248 0.04248 
+DEAL::0.002932 0.002378 0.002932 0.002378 0.002932 0.002378 0.002932 0.002378 0.01173 0.009511 0.01062 0.01062 0.01173 0.009511 0.01062 0.01062 0.01173 0.009511 0.01173 0.009511 0.04692 0.03804 0.04248 0.04248 0.04248 0.04248 0.1699 
 DEAL::laplace_matrix:
-DEAL::0.1223 -0.02510 -0.02189 -0.004440 -0.02189 -0.004440 -0.007037 0.008576 0.003618 0.02415 0.03501 0.03319 0.03974 -0.03097 0.03319 -0.03628 0.003618 0.02415 0.03974 -0.03097 -0.1124 0.07432 -0.09148 0.1105 -0.09148 0.1105 -0.184
-DEAL::-0.02510 0.07092 0.001985 -0.01598 0.001985 -0.01598 0.003376 0.0004581 0.01110 0.03551 -0.01101 0.005215 -0.01450 0.004161 0.005215 -0.006985 0.01110 0.03551 -0.01450 0.004161 0.05741 0.001033 -0.04566 0.01644 -0.04566 0.01644 -0.08662 
-DEAL::-0.02189 0.001985 0.04248 -0.01479 0.001524 0.002918 -0.01102 0.002675 -0.01389 0.007057 0.004571 0.01056 0.004751 -0.006585 -0.01162 0.006568 0.009976 -0.01115 0.01120 0.0002068 -0.02030 0.007905 0.04543 -0.02033 -0.01455 0.001842 -0.01554 
-DEAL::-0.004440 -0.01598 -0.01479 0.07667 0.002918 0.004126 0.006650 -0.02826 0.04541 -0.02951 0.02800 0.002700 -0.01581 0.01779 -0.01291 0.005728 -0.009973 -0.01478 -0.02203 0.08609 0.03381 -0.02704 0.03418 -0.01614 -0.07265 0.04785 -0.1076 
-DEAL::-0.02189 0.001985 0.001524 0.002918 0.04248 -0.01479 -0.01102 0.002675 0.009976 -0.01115 0.004571 -0.01162 0.01120 0.0002068 0.01056 0.006568 -0.01389 0.007057 0.004751 -0.006585 -0.02030 0.007905 -0.01455 0.001842 0.04543 -0.02033 -0.01554 
-DEAL::-0.004440 -0.01598 0.002918 0.004126 -0.01479 0.07667 0.006650 -0.02826 -0.009973 -0.01478 0.02800 -0.01291 -0.02203 0.08609 0.002700 0.005728 0.04541 -0.02951 -0.01581 0.01779 0.03381 -0.02704 -0.07265 0.04785 0.03418 -0.01614 -0.1076 
-DEAL::-0.007037 0.003376 -0.01102 0.006650 -0.01102 0.006650 0.06586 -0.03138 0.04348 -0.01612 -0.008973 -0.01943 -0.02107 0.01333 -0.01943 0.07482 0.04348 -0.01612 -0.02107 0.01333 -0.05751 0.04770 0.03004 -0.02514 0.03004 -0.02514 -0.0883
-DEAL::0.008576 0.0004581 0.002675 -0.02826 0.002675 -0.02826 -0.03138 0.07699 -0.03084 0.02819 -0.02128 0.01153 0.02290 -0.02252 0.01153 0.01321 -0.03084 0.02819 0.02290 -0.02252 0.03804 -0.06107 0.04139 -0.02988 0.04139 -0.02988 -0.01392 
-DEAL::0.003618 0.01110 -0.01389 0.04541 0.009976 -0.009973 0.04348 -0.03084 0.6490 -0.1225 -0.06429 -0.2163 -0.1058 -0.001276 0.03880 0.1274 -0.07986 0.04251 -0.2259 0.1347 0.3024 0.04412 -0.1845 -0.5367 0.3069 0.06781 -0.2355 
-DEAL::0.02415 0.03551 0.007057 -0.02951 -0.01115 -0.01478 -0.01612 0.02819 -0.1225 0.2567 -0.1229 0.03239 0.01952 -0.09283 0.03694 0.02095 0.04251 0.06188 0.05513 -0.1177 -0.02646 0.1383 -0.09797 0.02299 -0.06148 0.02457 -0.09341 
-DEAL::0.03501 -0.01101 0.004571 0.02800 0.004571 0.02800 -0.008973 -0.02128 -0.06429 -0.1229 0.4409 -0.06362 0.04406 0.08174 -0.06362 -0.01426 -0.06429 -0.1229 0.04406 0.08174 -0.2217 -0.2593 0.1539 0.08656 0.1539 0.08656 -0.2355 
-DEAL::0.03319 0.005215 0.01056 0.002700 -0.01162 -0.01291 -0.01943 0.01153 -0.2163 0.03239 -0.06362 0.2454 0.04984 0.02203 0.01020 -0.1216 0.03880 0.03694 0.09868 -0.07506 -0.1088 0.006450 -0.01085 0.2734 -0.2110 0.05071 -0.07694 
-DEAL::0.03974 -0.01450 0.004751 -0.01581 0.01120 -0.02203 -0.02107 0.02290 -0.1058 0.01952 0.04406 0.04984 0.2627 -0.1341 0.09868 -0.1141 -0.2259 0.05513 0.04674 0.01868 -0.1829 0.03635 -0.1136 0.03918 -0.04142 0.2624 -0.02070 
-DEAL::-0.03097 0.004161 -0.006585 0.01779 0.0002068 0.08609 0.01333 -0.02252 -0.001276 -0.09283 0.08174 0.02203 -0.1341 0.4592 -0.07506 -0.04238 0.1347 -0.1177 0.01868 -0.06856 0.1322 0.01819 -0.1906 -0.1486 0.08174 0.1070 -0.2460 
-DEAL::0.03319 0.005215 -0.01162 -0.01291 0.01056 0.002700 -0.01943 0.01153 0.03880 0.03694 -0.06362 0.01020 0.09868 -0.07506 0.2454 -0.1216 -0.2163 0.03239 0.04984 0.02203 -0.1088 0.006450 -0.2110 0.05071 -0.01085 0.2734 -0.07694 
-DEAL::-0.03628 -0.006985 0.006568 0.005728 0.006568 0.005728 0.07482 0.01321 0.1274 0.02095 -0.01426 -0.1216 -0.1141 -0.04238 -0.1216 0.4818 0.1274 0.02095 -0.1141 -0.04238 -0.1536 -0.1565 0.1490 0.006050 0.1490 0.006050 -0.2774 
-DEAL::0.003618 0.01110 0.009976 -0.009973 -0.01389 0.04541 0.04348 -0.03084 -0.07986 0.04251 -0.06429 0.03880 -0.2259 0.1347 -0.2163 0.1274 0.6490 -0.1225 -0.1058 -0.001276 0.3024 0.04412 0.3069 0.06781 -0.1845 -0.5367 -0.2355 
-DEAL::0.02415 0.03551 -0.01115 -0.01478 0.007057 -0.02951 -0.01612 0.02819 0.04251 0.06188 -0.1229 0.03694 0.05513 -0.1177 0.03239 0.02095 -0.1225 0.2567 0.01952 -0.09283 -0.02646 0.1383 -0.06148 0.02457 -0.09797 0.02299 -0.09341 
-DEAL::0.03974 -0.01450 0.01120 -0.02203 0.004751 -0.01581 -0.02107 0.02290 -0.2259 0.05513 0.04406 0.09868 0.04674 0.01868 0.04984 -0.1141 -0.1058 0.01952 0.2627 -0.1341 -0.1829 0.03635 -0.04142 0.2624 -0.1136 0.03918 -0.02070 
-DEAL::-0.03097 0.004161 0.0002068 0.08609 -0.006585 0.01779 0.01333 -0.02252 0.1347 -0.1177 0.08174 -0.07506 0.01868 -0.06856 0.02203 -0.04238 -0.001276 -0.09283 -0.1341 0.4592 0.1322 0.01819 0.08174 0.1070 -0.1906 -0.1486 -0.2460 
-DEAL::-0.1124 0.05741 -0.02030 0.03381 -0.02030 0.03381 -0.05751 0.03804 0.3024 -0.02646 -0.2217 -0.1088 -0.1829 0.1322 -0.1088 -0.1536 0.3024 -0.02646 -0.1829 0.1322 1.320 -0.2294 0.05689 -0.6022 0.05689 -0.6022 0.1898 
-DEAL::0.07432 0.001033 0.007905 -0.02704 0.007905 -0.02704 0.04770 -0.06107 0.04412 0.1383 -0.2593 0.006450 0.03635 0.01819 0.006450 -0.1565 0.04412 0.1383 0.03635 0.01819 -0.2294 1.083 -0.3656 -0.06094 -0.3656 -0.06094 -0.09484 
-DEAL::-0.09148 -0.04566 0.04543 0.03418 -0.01455 -0.07265 0.03004 0.04139 -0.1845 -0.09797 0.1539 -0.01085 -0.1136 -0.1906 -0.2110 0.1490 0.3069 -0.06148 -0.04142 0.08174 0.05689 -0.3656 1.274 -0.2473 -0.04871 -0.6694 0.2932 
-DEAL::0.1105 0.01644 -0.02033 -0.01614 0.001842 0.04785 -0.02514 -0.02988 -0.5367 0.02299 0.08656 0.2734 0.03918 -0.1486 0.05071 0.006050 0.06781 0.02457 0.2624 0.1070 -0.6022 -0.06094 -0.2473 1.664 -0.6694 -0.08390 -0.3407 
-DEAL::-0.09148 -0.04566 -0.01455 -0.07265 0.04543 0.03418 0.03004 0.04139 0.3069 -0.06148 0.1539 -0.2110 -0.04142 0.08174 -0.01085 0.1490 -0.1845 -0.09797 -0.1136 -0.1906 0.05689 -0.3656 -0.04871 -0.6694 1.274 -0.2473 0.2932 
-DEAL::0.1105 0.01644 0.001842 0.04785 -0.02033 -0.01614 -0.02514 -0.02988 0.06781 0.02457 0.08656 0.05071 0.2624 0.1070 0.2734 0.006050 -0.5367 0.02299 0.03918 -0.1486 -0.6022 -0.06094 -0.6694 -0.08390 -0.2473 1.664 -0.3407 
-DEAL::-0.1842 -0.08662 -0.01554 -0.1076 -0.01554 -0.1076 -0.08832 -0.01392 -0.2355 -0.09341 -0.2355 -0.07694 -0.02070 -0.2460 -0.07694 -0.2774 -0.2355 -0.09341 -0.02070 -0.2460 0.1898 -0.09484 0.2932 -0.3407 0.2932 -0.3407 2.477 
+DEAL::0.4768 -0.1162 -0.1394 0.02356 -0.1394 0.02356 0.02854 0.002027 0.2203 -0.01369 0.2466 -0.03212 -0.02780 -0.02077 -0.03212 -0.01977 0.2203 -0.01369 -0.02780 -0.02077 -0.1363 0.1301 -0.1371 0.1402 -0.1371 0.1402 -0.638
+DEAL::-0.1162 0.3135 0.02356 -0.09244 0.02356 -0.09244 0.002027 0.01953 -0.03040 0.1431 0.1131 -0.009685 -0.01242 -0.01816 -0.009685 -0.01457 -0.03040 0.1431 -0.01242 -0.01816 0.09672 -0.09273 -0.08163 0.07964 -0.08163 0.07964 -0.3245 
+DEAL::-0.1394 0.02356 0.4768 -0.1162 0.02854 0.002027 -0.1394 0.02356 0.2203 -0.01369 -0.03212 0.2466 -0.02780 -0.02077 -0.01977 -0.03212 -0.02780 -0.02077 0.2203 -0.01369 -0.1363 0.1301 0.1402 -0.1371 -0.1371 0.1402 -0.6382 
+DEAL::0.02356 -0.09244 -0.1162 0.3135 0.002027 0.01953 0.02356 -0.09244 -0.03040 0.1431 -0.009685 0.1131 -0.01242 -0.01816 -0.01457 -0.009685 -0.01242 -0.01816 -0.03040 0.1431 0.09672 -0.09273 0.07964 -0.08163 -0.08163 0.07964 -0.3245 
+DEAL::-0.1394 0.02356 0.02854 0.002027 0.4768 -0.1162 -0.1394 0.02356 -0.02780 -0.02077 -0.03212 -0.01977 0.2203 -0.01369 0.2466 -0.03212 0.2203 -0.01369 -0.02780 -0.02077 -0.1363 0.1301 -0.1371 0.1402 0.1402 -0.1371 -0.6382 
+DEAL::0.02356 -0.09244 0.002027 0.01953 -0.1162 0.3135 0.02356 -0.09244 -0.01242 -0.01816 -0.009685 -0.01457 -0.03040 0.1431 0.1131 -0.009685 -0.03040 0.1431 -0.01242 -0.01816 0.09672 -0.09273 -0.08163 0.07964 0.07964 -0.08163 -0.3245 
+DEAL::0.02854 0.002027 -0.1394 0.02356 -0.1394 0.02356 0.4768 -0.1162 -0.02780 -0.02077 -0.01977 -0.03212 0.2203 -0.01369 -0.03212 0.2466 -0.02780 -0.02077 0.2203 -0.01369 -0.1363 0.1301 0.1402 -0.1371 0.1402 -0.1371 -0.638
+DEAL::0.002027 0.01953 0.02356 -0.09244 0.02356 -0.09244 -0.1162 0.3135 -0.01242 -0.01816 -0.01457 -0.009685 -0.03040 0.1431 -0.009685 0.1131 -0.01242 -0.01816 -0.03040 0.1431 0.09672 -0.09273 0.07964 -0.08163 0.07964 -0.08163 -0.3245 
+DEAL::0.2203 -0.03040 0.2203 -0.03040 -0.02780 -0.01242 -0.02780 -0.01242 0.4679 -0.03763 0.1318 0.1318 -0.03798 -0.03650 0.005772 0.005772 0.1304 0.003119 0.1304 0.003119 -0.1582 0.2068 -0.05651 -0.05651 -0.1314 0.1959 -1.197 
+DEAL::-0.01369 0.1431 -0.01369 0.1431 -0.02077 -0.01816 -0.02077 -0.01816 -0.03763 0.3793 0.05357 0.05357 -0.03650 -0.06336 0.01204 0.01204 0.003119 0.08180 0.003119 0.08180 0.1399 -0.04042 -0.06086 -0.06086 -0.1377 0.1599 -0.7238 
+DEAL::0.2466 0.1131 -0.03212 -0.009685 -0.03212 -0.009685 -0.01977 -0.01457 0.1318 0.05357 0.4655 -0.05768 0.005772 0.01204 -0.05768 -0.03112 0.1318 0.05357 0.005772 0.01204 -0.1004 -0.05410 -0.08531 0.1626 -0.08531 0.1626 -0.9671 
+DEAL::-0.03212 -0.009685 0.2466 0.1131 -0.01977 -0.01457 -0.03212 -0.009685 0.1318 0.05357 -0.05768 0.4655 0.005772 0.01204 -0.03112 -0.05768 0.005772 0.01204 0.1318 0.05357 -0.1004 -0.05410 0.1626 -0.08531 -0.08531 0.1626 -0.9671 
+DEAL::-0.02780 -0.01242 -0.02780 -0.01242 0.2203 -0.03040 0.2203 -0.03040 -0.03798 -0.03650 0.005772 0.005772 0.4679 -0.03763 0.1318 0.1318 0.1304 0.003119 0.1304 0.003119 -0.1582 0.2068 -0.05651 -0.05651 0.1959 -0.1314 -1.197 
+DEAL::-0.02077 -0.01816 -0.02077 -0.01816 -0.01369 0.1431 -0.01369 0.1431 -0.03650 -0.06336 0.01204 0.01204 -0.03763 0.3793 0.05357 0.05357 0.003119 0.08180 0.003119 0.08180 0.1399 -0.04042 -0.06086 -0.06086 0.1599 -0.1377 -0.7238 
+DEAL::-0.03212 -0.009685 -0.01977 -0.01457 0.2466 0.1131 -0.03212 -0.009685 0.005772 0.01204 -0.05768 -0.03112 0.1318 0.05357 0.4655 -0.05768 0.1318 0.05357 0.005772 0.01204 -0.1004 -0.05410 -0.08531 0.1626 0.1626 -0.08531 -0.9671 
+DEAL::-0.01977 -0.01457 -0.03212 -0.009685 -0.03212 -0.009685 0.2466 0.1131 0.005772 0.01204 -0.03112 -0.05768 0.1318 0.05357 -0.05768 0.4655 0.005772 0.01204 0.1318 0.05357 -0.1004 -0.05410 0.1626 -0.08531 0.1626 -0.08531 -0.9671 
+DEAL::0.2203 -0.03040 -0.02780 -0.01242 0.2203 -0.03040 -0.02780 -0.01242 0.1304 0.003119 0.1318 0.005772 0.1304 0.003119 0.1318 0.005772 0.4679 -0.03763 -0.03798 -0.03650 -0.1582 0.2068 -0.1314 0.1959 -0.05651 -0.05651 -1.197 
+DEAL::-0.01369 0.1431 -0.02077 -0.01816 -0.01369 0.1431 -0.02077 -0.01816 0.003119 0.08180 0.05357 0.01204 0.003119 0.08180 0.05357 0.01204 -0.03763 0.3793 -0.03650 -0.06336 0.1399 -0.04042 -0.1377 0.1599 -0.06086 -0.06086 -0.7238 
+DEAL::-0.02780 -0.01242 0.2203 -0.03040 -0.02780 -0.01242 0.2203 -0.03040 0.1304 0.003119 0.005772 0.1318 0.1304 0.003119 0.005772 0.1318 -0.03798 -0.03650 0.4679 -0.03763 -0.1582 0.2068 0.1959 -0.1314 -0.05651 -0.05651 -1.197 
+DEAL::-0.02077 -0.01816 -0.01369 0.1431 -0.02077 -0.01816 -0.01369 0.1431 0.003119 0.08180 0.01204 0.05357 0.003119 0.08180 0.01204 0.05357 -0.03650 -0.06336 -0.03763 0.3793 0.1399 -0.04042 0.1599 -0.1377 -0.06086 -0.06086 -0.7238 
+DEAL::-0.1363 0.09672 -0.1363 0.09672 -0.1363 0.09672 -0.1363 0.09672 -0.1582 0.1399 -0.1004 -0.1004 -0.1582 0.1399 -0.1004 -0.1004 -0.1582 0.1399 -0.1582 0.1399 0.9151 -0.4283 -0.1217 -0.1217 -0.1217 -0.1217 0.6333 
+DEAL::0.1301 -0.09273 0.1301 -0.09273 0.1301 -0.09273 0.1301 -0.09273 0.2068 -0.04042 -0.05410 -0.05410 0.2068 -0.04042 -0.05410 -0.05410 0.2068 -0.04042 0.2068 -0.04042 -0.4283 1.160 -0.1830 -0.1830 -0.1830 -0.1830 -0.5986 
+DEAL::-0.1371 -0.08163 0.1402 0.07964 -0.1371 -0.08163 0.1402 0.07964 -0.05651 -0.06086 -0.08531 0.1626 -0.05651 -0.06086 -0.08531 0.1626 -0.1314 -0.1377 0.1959 0.1599 -0.1217 -0.1830 0.8708 -0.3970 -0.08456 -0.08456 -0.008681 
+DEAL::0.1402 0.07964 -0.1371 -0.08163 0.1402 0.07964 -0.1371 -0.08163 -0.05651 -0.06086 0.1626 -0.08531 -0.05651 -0.06086 0.1626 -0.08531 0.1959 0.1599 -0.1314 -0.1377 -0.1217 -0.1830 -0.3970 0.8708 -0.08456 -0.08456 -0.008681 
+DEAL::-0.1371 -0.08163 -0.1371 -0.08163 0.1402 0.07964 0.1402 0.07964 -0.1314 -0.1377 -0.08531 -0.08531 0.1959 0.1599 0.1626 0.1626 -0.05651 -0.06086 -0.05651 -0.06086 -0.1217 -0.1830 -0.08456 -0.08456 0.8708 -0.3970 -0.008681 
+DEAL::0.1402 0.07964 0.1402 0.07964 -0.1371 -0.08163 -0.1371 -0.08163 0.1959 0.1599 0.1626 0.1626 -0.1314 -0.1377 -0.08531 -0.08531 -0.05651 -0.06086 -0.05651 -0.06086 -0.1217 -0.1830 -0.08456 -0.08456 -0.3970 0.8708 -0.008681 
+DEAL::-0.6382 -0.3245 -0.6382 -0.3245 -0.6382 -0.3245 -0.6382 -0.3245 -1.197 -0.7238 -0.9671 -0.9671 -1.197 -0.7238 -0.9671 -0.9671 -1.197 -0.7238 -1.197 -0.7238 0.6333 -0.5986 -0.008681 -0.008681 -0.008681 -0.008681 15.40 
 DEAL::cell=1
 DEAL::mass_matrix:
 DEAL::0.001372 -0.0006859 -0.0006859 0.0003429 -0.0006859 0.0003429 0.0003429 -0.0001715 0.001372 -0.0006859 0.001372 -0.0006859 -0.0006859 0.0003429 -0.0006859 0.0003429 0.001372 -0.0006859 -0.0006859 0.0003429 0.001372 -0.0006859 0.001372 -0.0006859 0.001372 -0.0006859 0.001372 
index 414cb88be41a099a80ea4d29a7768fae0ad3f305..48be21c02c5523ebde727cde91ab1c6e72cce64d 100644 (file)
@@ -148,7 +148,7 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary;
+  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
   tr.set_boundary (1, boundary);
 
                                   // set boundary id on cell 1
index 5031deacca9ecb2be091e74771b4399f3adf7e19..d75bd98b1d3dd40f4df57f485e17ea4c60d89723 100644 (file)
@@ -5,15 +5,15 @@ DEAL::Mass matrices are not equal.
 DEAL::Laplace matrices are not equal.
 DEAL::cell=0
 DEAL::mass_matrix:
-DEAL::nan nan nan nan 
-DEAL::nan nan nan nan 
-DEAL::nan nan nan nan 
-DEAL::nan nan nan nan 
+DEAL::0.1804 0.08456 0.09018 0.04228 
+DEAL::0.08456 0.1579 0.04228 0.07894 
+DEAL::0.09018 0.04228 0.1804 0.08456 
+DEAL::0.04228 0.07894 0.08456 0.1579 
 DEAL::laplace_matrix:
-DEAL::nan nan nan nan 
-DEAL::nan nan nan nan 
-DEAL::nan nan nan nan 
-DEAL::nan nan nan nan 
+DEAL::0.8342 -0.2993 -0.2329 -0.3019 
+DEAL::-0.2993 0.7035 -0.3019 -0.1023 
+DEAL::-0.2329 -0.3019 0.8342 -0.2993 
+DEAL::-0.3019 -0.1023 -0.2993 0.7035 
 DEAL::cell=1
 DEAL::mass_matrix:
 DEAL::0.1111 0.05556 0.05556 0.02778 
index 3ced2f5828e456ee8ea5efd803bd6f34c1212651..19e140d417cae21aca0dbc9a5f28be2da0598870 100644 (file)
@@ -148,7 +148,7 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary;
+  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
   tr.set_boundary (1, boundary);
 
                                   // set boundary id on cell 1
index c9f1f3a18e3e58ed6ab40b2b785fbd9aef493435..55d00e98d0510b58a8197254dc970b9c43963614 100644 (file)
@@ -5,25 +5,25 @@ DEAL::Mass matrices are not equal.
 DEAL::Laplace matrices are not equal.
 DEAL::cell=0
 DEAL::mass_matrix:
-DEAL::nan nan nan nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan nan nan nan 
+DEAL::0.02004 -0.009396 -0.01002 0.004698 0.02004 -0.009396 0.02129 -0.01065 0.02129 
+DEAL::-0.009396 0.01754 0.004698 -0.008771 -0.009396 0.01754 0.01629 -0.008147 0.01629 
+DEAL::-0.01002 0.004698 0.02004 -0.009396 0.02004 -0.009396 -0.01065 0.02129 0.02129 
+DEAL::0.004698 -0.008771 -0.009396 0.01754 -0.009396 0.01754 -0.008147 0.01629 0.01629 
+DEAL::0.02004 -0.009396 0.02004 -0.009396 0.08016 -0.03758 0.02129 0.02129 0.08516 
+DEAL::-0.009396 0.01754 -0.009396 0.01754 -0.03758 0.07017 0.01629 0.01629 0.06517 
+DEAL::0.02129 0.01629 -0.01065 -0.008147 0.02129 0.01629 0.07517 -0.03758 0.07517 
+DEAL::-0.01065 -0.008147 0.02129 0.01629 0.02129 0.01629 -0.03758 0.07517 0.07517 
+DEAL::0.02129 0.01629 0.02129 0.01629 0.08516 0.06517 0.07517 0.07517 0.3007 
 DEAL::laplace_matrix:
-DEAL::nan nan nan nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan nan nan nan 
+DEAL::0.7717 -0.1078 -0.1509 -0.03354 0.07142 0.1934 0.01508 0.1997 -0.9591 
+DEAL::-0.1078 0.6443 -0.03354 -0.1525 0.1495 0.1179 -0.1185 0.2017 -0.7011 
+DEAL::-0.1509 -0.03354 0.7717 -0.1078 0.07142 0.1934 0.1997 0.01508 -0.9591 
+DEAL::-0.03354 -0.1525 -0.1078 0.6443 0.1495 0.1179 0.2017 -0.1185 -0.7011 
+DEAL::0.07142 0.1495 0.07142 0.1495 1.131 -0.002272 -0.3435 -0.3435 -0.8837 
+DEAL::0.1934 0.1179 0.1934 0.1179 -0.002272 1.498 -0.4364 -0.4364 -1.245 
+DEAL::0.01508 -0.1185 0.1997 0.2017 -0.3435 -0.4364 1.197 -0.1194 -0.5960 
+DEAL::0.1997 0.2017 0.01508 -0.1185 -0.3435 -0.4364 -0.1194 1.197 -0.5960 
+DEAL::-0.9591 -0.7011 -0.9591 -0.7011 -0.8837 -1.245 -0.5960 -0.5960 6.641 
 DEAL::cell=1
 DEAL::mass_matrix:
 DEAL::0.01235 -0.006173 -0.006173 0.003086 0.01235 -0.006173 0.01235 -0.006173 0.01235 
index 721f8107c1ff2ac1aba84e52f04b4b87968e51ab..1e07b7d18e2af64e21bc990307b8ed475c4a23f2 100644 (file)
@@ -148,7 +148,7 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary;
+  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
   tr.set_boundary (1, boundary);
 
                                   // set boundary id on cell 1
index 5a5cfcef464d6b5531286396f635fd3ed484098a..a9aa94651d791679b38cff2f22405a811dfa013c 100644 (file)
@@ -148,7 +148,7 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary;
+  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
   tr.set_boundary (1, boundary);
 
                                   // set boundary id on cell 1
index d6b9163249964faed0836923588562b2e1979dda..b091c11059793763dfa74193c524b82b4be8bf99 100644 (file)
@@ -148,7 +148,7 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary;
+  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
   tr.set_boundary (1, boundary);
 
                                   // set boundary id on cell 1
index 589b06efe1d96dd1be9d289ee8b962c4b42fd0e8..bad1cf83b0fe3da9e497547c2dba665833ff31ea 100644 (file)
@@ -148,7 +148,7 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary;
+  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
   tr.set_boundary (1, boundary);
 
                                   // set boundary id on cell 1
index 8e00250c5d5b402cb0ab9ae54245dcc4506d948d..752d4d9664285371e4cda311fe0fb4525f3dc11d 100644 (file)
@@ -148,7 +148,7 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary;
+  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
   tr.set_boundary (1, boundary);
 
                                   // set boundary id on cell 1
index 394b756ad1ab8de3fe40131f1ecef60976b97111..7ba94b1530068da588986e62ef59f4290e072579 100644 (file)
@@ -148,7 +148,7 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary;
+  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
   tr.set_boundary (1, boundary);
 
                                   // set boundary id on cell 1
index 5ea27062781476d47bbc83f0d49aa1f4f61d901f..0bef24ce1c1eee43ecbc4238b11f96f6434ea4e7 100644 (file)
@@ -148,7 +148,7 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary;
+  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
   tr.set_boundary (1, boundary);
 
                                   // set boundary id on cell 1
index 6ed6e3be786faeed17e9072a99964b04a3e5f0ff..0a53a5c6da84719dfc2e144bf04d4f2b93f546f8 100644 (file)
@@ -5,13 +5,13 @@ DEAL::Mass matrices are not equal.
 DEAL::Laplace matrices are not equal.
 DEAL::cell=0
 DEAL::mass_matrix:
-DEAL::0.9296 0.3715 0.3446 
-DEAL::0.3715 0.2166 0.1883 
-DEAL::0.3446 0.1883 0.189
+DEAL::2.000 0.9491 1.000 
+DEAL::0.9491 0.6157 0.4745 
+DEAL::1.000 0.4745 0.666
 DEAL::laplace_matrix:
 DEAL::0 0 0 
-DEAL::0 1.930 1.352 
-DEAL::0 1.352 4.80
+DEAL::0 1.503 0 
+DEAL::0 0 1.12
 DEAL::cell=1
 DEAL::mass_matrix:
 DEAL::1.000 0.5000 0.5000 
@@ -27,15 +27,15 @@ DEAL::Mass matrices are not equal.
 DEAL::Laplace matrices are not equal.
 DEAL::cell=0
 DEAL::mass_matrix:
-DEAL::0.8419 0.3942 0.3606 0.3606 
-DEAL::0.3942 0.2539 0.1892 0.1892 
-DEAL::0.3606 0.1892 0.2203 0.1794 
-DEAL::0.3606 0.1892 0.1794 0.2203 
+DEAL::1.936 0.9341 0.9678 0.9678 
+DEAL::0.9341 0.6115 0.4670 0.4670 
+DEAL::0.9678 0.4670 0.6452 0.4839 
+DEAL::0.9678 0.4670 0.4839 0.6452 
 DEAL::laplace_matrix:
 DEAL::0 0 0 0 
-DEAL::0 1.041 0.4212 0.4212 
-DEAL::0 0.4212 1.312 0.4953 
-DEAL::0 0.4212 0.4953 1.31
+DEAL::0 2.124 0 0 
+DEAL::0 0 2.032 0 
+DEAL::0 0 0 2.03
 DEAL::cell=1
 DEAL::mass_matrix:
 DEAL::1.000 0.5000 0.5000 0.5000 
index 0cf5ef131840441a7e87a054713b0e82b9ca68ea..c2071d883238dd57197accd6bf8cd1339d1bc509 100644 (file)
@@ -148,7 +148,7 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary;
+  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
   tr.set_boundary (1, boundary);
 
                                   // set boundary id on cell 1
index 2cc704bc169ad93684bae531074c65a6f61d2ed1..f0609c6893b3ec8e9c41d2993868e3f23512678c 100644 (file)
@@ -5,19 +5,19 @@ DEAL::Mass matrices are not equal.
 DEAL::Laplace matrices are not equal.
 DEAL::cell=0
 DEAL::mass_matrix:
-DEAL::0.9296 0.3715 0.3446 0.1883 0.2166 0.189
-DEAL::0.3715 0.2166 0.1883 0.1309 0.1547 0.1264 
-DEAL::0.3446 0.1883 0.1897 0.1264 0.1309 0.1323 
-DEAL::0.1883 0.1309 0.1264 0.09477 0.09948 0.09500 
-DEAL::0.2166 0.1547 0.1309 0.09948 0.1186 0.09477 
-DEAL::0.1897 0.1264 0.1323 0.09500 0.09477 0.1006 
+DEAL::2.000 0.9491 1.000 0.4745 0.6157 0.666
+DEAL::0.9491 0.6157 0.4745 0.3079 0.4576 0.3164 
+DEAL::1.000 0.4745 0.6667 0.3164 0.3079 0.5000 
+DEAL::0.4745 0.3079 0.3164 0.2052 0.2288 0.2373 
+DEAL::0.6157 0.4576 0.3079 0.2288 0.3549 0.2052 
+DEAL::0.6667 0.3164 0.5000 0.2373 0.2052 0.3889 
 DEAL::laplace_matrix:
 DEAL::0 0 0 0 0 0 
-DEAL::0 1.930 1.352 1.098 2.451 1.602 
-DEAL::0 1.352 4.800 2.053 0.9102 7.074 
-DEAL::0 1.098 2.053 1.205 1.216 2.949 
-DEAL::0 2.451 0.9102 1.216 3.615 0.863
-DEAL::0 1.602 7.074 2.949 0.8631 10.95 
+DEAL::0 1.503 0 0.7515 1.515 -0.1120 
+DEAL::0 0 1.120 0.4402 0 1.120 
+DEAL::0 0.7515 0.4402 0.8880 0.7573 0.3842 
+DEAL::0 1.515 0 0.7573 2.027 0.155
+DEAL::0 -0.1120 1.120 0.3842 0.1551 1.494 
 DEAL::cell=1
 DEAL::mass_matrix:
 DEAL::1.000 0.5000 0.5000 0.2500 0.3333 0.3333 
@@ -39,27 +39,27 @@ DEAL::Mass matrices are not equal.
 DEAL::Laplace matrices are not equal.
 DEAL::cell=0
 DEAL::mass_matrix:
-DEAL::0.8419 0.3942 0.3606 0.3606 0.1892 0.1892 0.1794 0.2539 0.2203 0.2203 
-DEAL::0.3942 0.2539 0.1892 0.1892 0.1291 0.1291 0.1017 0.1882 0.1235 0.1235 
-DEAL::0.3606 0.1892 0.2203 0.1794 0.1235 0.1017 0.1193 0.1291 0.1602 0.1193 
-DEAL::0.3606 0.1892 0.1794 0.2203 0.1017 0.1235 0.1193 0.1291 0.1193 0.1602 
-DEAL::0.1892 0.1291 0.1235 0.1017 0.08679 0.07184 0.07021 0.09758 0.09198 0.07021 
-DEAL::0.1892 0.1291 0.1017 0.1235 0.07184 0.08679 0.07021 0.09758 0.07021 0.09198 
-DEAL::0.1794 0.1017 0.1193 0.1193 0.07021 0.07021 0.08261 0.07184 0.08942 0.08942 
-DEAL::0.2539 0.1882 0.1291 0.1291 0.09758 0.09758 0.07184 0.1459 0.08679 0.08679 
-DEAL::0.2203 0.1235 0.1602 0.1193 0.09198 0.07021 0.08942 0.08679 0.1235 0.0826
-DEAL::0.2203 0.1235 0.1193 0.1602 0.07021 0.09198 0.08942 0.08679 0.08261 0.1235 
+DEAL::1.936 0.9341 0.9678 0.9678 0.4670 0.4670 0.4839 0.6115 0.6452 0.6452 
+DEAL::0.9341 0.6115 0.4670 0.4670 0.3057 0.3057 0.2335 0.4558 0.3114 0.3114 
+DEAL::0.9678 0.4670 0.6452 0.4839 0.3114 0.2335 0.3226 0.3057 0.4839 0.3226 
+DEAL::0.9678 0.4670 0.4839 0.6452 0.2335 0.3114 0.3226 0.3057 0.3226 0.4839 
+DEAL::0.4670 0.3057 0.3114 0.2335 0.2038 0.1529 0.1557 0.2279 0.2335 0.1557 
+DEAL::0.4670 0.3057 0.2335 0.3114 0.1529 0.2038 0.1557 0.2279 0.1557 0.2335 
+DEAL::0.4839 0.2335 0.3226 0.3226 0.1557 0.1557 0.2151 0.1529 0.2419 0.2419 
+DEAL::0.6115 0.4558 0.3057 0.3057 0.2279 0.2279 0.1529 0.3539 0.2038 0.2038 
+DEAL::0.6452 0.3114 0.4839 0.3226 0.2335 0.1557 0.2419 0.2038 0.3764 0.215
+DEAL::0.6452 0.3114 0.3226 0.4839 0.1557 0.2335 0.2419 0.2038 0.2151 0.3764 
 DEAL::laplace_matrix:
 DEAL::0 0 0 0 0 0 0 0 0 0 
-DEAL::0 1.041 0.4212 0.4212 0.6216 0.6216 0.3551 1.086 0.4638 0.4638 
-DEAL::0 0.4212 1.312 0.4953 0.8130 0.3591 0.8161 0.3333 1.524 0.4967 
-DEAL::0 0.4212 0.4953 1.312 0.3591 0.8130 0.8161 0.3333 0.4967 1.524 
-DEAL::0 0.6216 0.8130 0.3591 0.8054 0.4236 0.5330 0.6581 0.9996 0.3187 
-DEAL::0 0.6216 0.3591 0.8130 0.4236 0.8054 0.5330 0.6581 0.3187 0.9996 
-DEAL::0 0.3551 0.8161 0.8161 0.5330 0.5330 0.8905 0.2349 0.9729 0.9729 
-DEAL::0 1.086 0.3333 0.3333 0.6581 0.6581 0.2349 1.478 0.3230 0.323
-DEAL::0 0.4638 1.524 0.4967 0.9996 0.3187 0.9729 0.3230 2.173 0.3850 
-DEAL::0 0.4638 0.4967 1.524 0.3187 0.9996 0.9729 0.3230 0.3850 2.173 
+DEAL::0 2.124 0 0 1.062 1.062 0 2.036 -0.2256 -0.2256 
+DEAL::0 0 2.032 0 0.7899 0 1.016 0 2.032 0 
+DEAL::0 0 0 2.032 0 0.7899 1.016 0 0 2.032 
+DEAL::0 1.062 0.7899 0 1.387 0.5310 0.3949 1.018 0.6771 -0.1128 
+DEAL::0 1.062 0 0.7899 0.5310 1.387 0.3949 1.018 -0.1128 0.6771 
+DEAL::0 0 1.016 1.016 0.3949 0.3949 1.603 0 1.016 1.016 
+DEAL::0 2.036 0 0 1.018 1.018 0 2.656 0.2290 0.229
+DEAL::0 -0.2256 2.032 0 0.6771 -0.1128 1.016 0.2290 2.709 0.4964 
+DEAL::0 -0.2256 0 2.032 -0.1128 0.6771 1.016 0.2290 0.4964 2.709 
 DEAL::cell=1
 DEAL::mass_matrix:
 DEAL::1.000 0.5000 0.5000 0.5000 0.2500 0.2500 0.2500 0.3333 0.3333 0.3333 
index 6b3dc96c4e4ad6b859142585298adcc7af096fdf..0afe1e8da51d68b1baebc142909ce3fcba8e5341 100644 (file)
@@ -148,7 +148,7 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary;
+  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
   tr.set_boundary (1, boundary);
 
                                   // set boundary id on cell 1
index 74460e5960b20156e73419d46a8c8611449389fe..ecfb2838833f72d72cf29c7b6a2cc686978c8d15 100644 (file)
@@ -5,13 +5,13 @@ DEAL::Mass matrices are not equal.
 DEAL::Laplace matrices are not equal.
 DEAL::cell=0
 DEAL::mass_matrix:
-DEAL::nan nan nan 
-DEAL::nan nan nan 
-DEAL::nan nan nan 
+DEAL::1.522 0.7273 0.7611 
+DEAL::0.7273 0.4736 0.3637 
+DEAL::0.7611 0.3637 0.5074 
 DEAL::laplace_matrix:
-DEAL::nan nan nan 
-DEAL::nan nan nan 
-DEAL::nan nan nan 
+DEAL::0 0 0 
+DEAL::0 1.202 0 
+DEAL::0 0 0.9390 
 DEAL::cell=1
 DEAL::mass_matrix:
 DEAL::1.000 0.5000 0.5000 
index c9da1735576920b0f640b88664742e18a6d8b442..0ed3ed52866192a713cc11fa09f2d5e2bb2e7abe 100644 (file)
@@ -148,7 +148,7 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary;
+  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
   tr.set_boundary (1, boundary);
 
                                   // set boundary id on cell 1
index ef93432ba51012841a8d74ca274808e55b0da3fb..cfc564db3c3a6574f4a3a443c34dfe6e71206984 100644 (file)
@@ -5,19 +5,19 @@ DEAL::Mass matrices are not equal.
 DEAL::Laplace matrices are not equal.
 DEAL::cell=0
 DEAL::mass_matrix:
-DEAL::nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan 
+DEAL::1.522 0.7273 0.7611 0.3637 0.4736 0.5074 
+DEAL::0.7273 0.4736 0.3637 0.2368 0.3524 0.2424 
+DEAL::0.7611 0.3637 0.5074 0.2424 0.2368 0.3805 
+DEAL::0.3637 0.2368 0.2424 0.1579 0.1762 0.1818 
+DEAL::0.4736 0.3524 0.2368 0.1762 0.2735 0.1579 
+DEAL::0.5074 0.2424 0.3805 0.1818 0.1579 0.2960 
 DEAL::laplace_matrix:
-DEAL::nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan 
+DEAL::0 0 0 0 0 0 
+DEAL::0 1.202 0 0.6012 1.239 -0.06579 
+DEAL::0 0 0.9390 0.4042 0 0.9390 
+DEAL::0 0.6012 0.4042 0.7035 0.6196 0.3713 
+DEAL::0 1.239 0 0.6196 1.677 0.04426 
+DEAL::0 -0.06579 0.9390 0.3713 0.04426 1.252 
 DEAL::cell=1
 DEAL::mass_matrix:
 DEAL::1.000 0.5000 0.5000 0.2500 0.3333 0.3333 
index 6e35fc175bd9ddb24383305c4feafe30e14d9bff..fc16a101bcb7d682e58154d0fccc16716c1bd46a 100644 (file)
@@ -148,7 +148,7 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary;
+  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
   tr.set_boundary (1, boundary);
 
                                   // set boundary id on cell 1
index d6adcfa9e9bf699cad40a34d55345fde61339366..55366580229341f21eae777ebe0977dfc0f0654c 100644 (file)
@@ -148,7 +148,7 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary;
+  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
   tr.set_boundary (1, boundary);
 
                                   // set boundary id on cell 1
index cc08baaaa47e11611ca16c29203b8d9ee5bcfb5f..971e8091d00c28b889cbacb70c1d84d5fddd33dc 100644 (file)
@@ -148,7 +148,7 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary;
+  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
   tr.set_boundary (1, boundary);
 
                                   // set boundary id on cell 1
index cb08263629b67544ffd4456f376eb38600112ee4..2d529983b9bd4339d00373e02d894ca3a654a0be 100644 (file)
@@ -148,7 +148,7 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary;
+  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
   tr.set_boundary (1, boundary);
 
                                   // set boundary id on cell 1
index 0f0a01e193f2a1dfff8ea29e233d2b1d6c1ae9ae..78b8e8c232ae2e932a63ff42427d54d84a789e82 100644 (file)
@@ -148,7 +148,7 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary;
+  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
   tr.set_boundary (1, boundary);
 
                                   // set boundary id on cell 1
index 3c7adeab3e61ffdd0a301d50de1a013af14dda98..f63bb107e2a90c24dae13f69c5a70641ccf490f7 100644 (file)
@@ -148,7 +148,7 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary;
+  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
   tr.set_boundary (1, boundary);
 
                                   // set boundary id on cell 1
index ad71c2b4bdaad6dca39b24e03ef55a94eb498e5b..f02787f13444ae81db59e26ca330d130858f3f94 100644 (file)
@@ -148,7 +148,7 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary;
+  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
   tr.set_boundary (1, boundary);
 
                                   // set boundary id on cell 1
index 2930838742c0ea167e2f566303083087f6ad6b09..1f829bc776558e98b892bd45ffde26216226b4ce 100644 (file)
@@ -5,13 +5,13 @@ DEAL::Mass matrices are not equal.
 DEAL::Laplace matrices are not equal.
 DEAL::cell=0
 DEAL::mass_matrix:
-DEAL::0.9296 0.1103 -0.03228 
-DEAL::0.1103 0.4840 0.1078 
-DEAL::-0.03228 0.1078 0.3581 
+DEAL::2.000 -0.9088 0 
+DEAL::-0.9088 3.941 0 
+DEAL::0 0 5.556 
 DEAL::laplace_matrix:
 DEAL::0 0 0 
-DEAL::0 11.16 0 
-DEAL::0 0 11.16 
+DEAL::0 24.00 0 
+DEAL::0 0 24.00 
 DEAL::cell=1
 DEAL::mass_matrix:
 DEAL::1.000 3.464 0 
@@ -27,15 +27,15 @@ DEAL::Mass matrices are not equal.
 DEAL::Laplace matrices are not equal.
 DEAL::cell=0
 DEAL::mass_matrix:
-DEAL::0.8419 0.02942 -0.1157 -0.1157 
-DEAL::0.02942 0.7414 0.003550 0.00355
-DEAL::-0.1157 0.003550 0.5949 0.05282 
-DEAL::-0.1157 0.003550 0.05282 0.5949 
+DEAL::1.936 -0.4737 0 0 
+DEAL::-0.4737 2.997 0 
+DEAL::0 0 4.038 0 
+DEAL::0 0 0 4.038 
 DEAL::laplace_matrix:
 DEAL::0 0 0 0 
-DEAL::0 10.10 0 0 
-DEAL::0 0 10.10 0 
-DEAL::0 0 0 10.10 
+DEAL::0 23.23 0 0 
+DEAL::0 0 23.23 0 
+DEAL::0 0 0 23.23 
 DEAL::cell=1
 DEAL::mass_matrix:
 DEAL::1.000 3.464 0 0 
index ed86fe8590657256ab59364bb9f7369419d28d8a..ea28809c5a67a01d593cb968bf8f4ec7951094e4 100644 (file)
@@ -148,7 +148,7 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary;
+  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
   tr.set_boundary (1, boundary);
 
                                   // set boundary id on cell 1
index ab3cd27ebab493e8aadbcac4d5d3c46ef9fa5951..240689c9b4d461b5304cbcf331d752ff2e047ebe 100644 (file)
@@ -5,19 +5,19 @@ DEAL::Mass matrices are not equal.
 DEAL::Laplace matrices are not equal.
 DEAL::cell=0
 DEAL::mass_matrix:
-DEAL::0.9296 0.1103 -0.4982 -0.03228 0.1078 -0.6390 
-DEAL::0.1103 0.4840 0.3899 0.1078 0.003193 0.1013 
-DEAL::-0.4982 0.3899 0.8769 0.03966 -0.3244 0.6142 
-DEAL::-0.03228 0.1078 0.03966 0.3581 0.2009 0.1898 
-DEAL::0.1078 0.003193 -0.3244 0.2009 0.4039 -0.1713 
-DEAL::-0.6390 0.1013 0.6142 0.1898 -0.1713 0.6670 
+DEAL::2.000 -0.9088 2.170 0 0 3.975 
+DEAL::-0.9088 3.941 -3.643 0 0 -1.806 
+DEAL::2.170 -3.643 4.355 0 0 4.313 
+DEAL::0 0 0 5.556 -2.524 0 
+DEAL::0 0 0 -2.524 10.95 0 
+DEAL::3.975 -1.806 4.313 0 0 7.901 
 DEAL::laplace_matrix:
 DEAL::0 0 0 0 0 0 
-DEAL::0 11.16 2.959 0 -0.3873 0 
-DEAL::0 2.959 29.04 0 2.893 0 
-DEAL::0 0 0 11.16 1.323 -0.8661 
-DEAL::0 -0.3873 2.893 1.323 10.11 2.893 
-DEAL::0 0 0 -0.8661 2.893 21.49 
+DEAL::0 24.00 -24.39 0 0 0 
+DEAL::0 -24.39 236.5 0 0 0 
+DEAL::0 0 0 24.00 -10.91 0 
+DEAL::0 0 0 -10.91 114.0 0 
+DEAL::0 0 0 0 0 333.3 
 DEAL::cell=1
 DEAL::mass_matrix:
 DEAL::1.000 3.464 13.42 0 0 0 
@@ -39,27 +39,27 @@ DEAL::Mass matrices are not equal.
 DEAL::Laplace matrices are not equal.
 DEAL::cell=0
 DEAL::mass_matrix:
-DEAL::0.8419 0.02942 -0.1123 -0.1157 0.003550 -0.2762 -0.1157 0.003550 0.05282 -0.2762 
-DEAL::0.02942 0.7414 0.2219 0.003550 -0.08994 0.03889 0.003550 -0.08994 -0.03985 0.03889 
-DEAL::-0.1123 0.2219 0.2460 0.02874 -0.2454 0.1126 0.02874 -0.2454 -0.02620 0.1126 
-DEAL::-0.1157 0.003550 0.02874 0.5949 0.06421 0.1207 0.05282 -0.03985 -0.04795 0.0757
-DEAL::0.003550 -0.08994 -0.2454 0.06421 0.5845 -0.1252 -0.03985 0.02939 -0.01761 -0.02365 
-DEAL::-0.2762 0.03889 0.1126 0.1207 -0.1252 0.1920 0.07570 -0.02365 -0.1859 0.1656 
-DEAL::-0.1157 0.003550 0.02874 0.05282 -0.03985 0.07570 0.5949 0.06421 -0.04795 0.1207 
-DEAL::0.003550 -0.08994 -0.2454 -0.03985 0.02939 -0.02365 0.06421 0.5845 -0.01761 -0.1252 
-DEAL::0.05282 -0.03985 -0.02620 -0.04795 -0.01761 -0.1859 -0.04795 -0.01761 0.4803 -0.1859 
-DEAL::-0.2762 0.03889 0.1126 0.07570 -0.02365 0.1656 0.1207 -0.1252 -0.1859 0.1920 
+DEAL::1.936 -0.4737 1.187 0 0 2.351 0 0 0 2.351 
+DEAL::-0.4737 2.997 -1.392 0 0 -0.5754 0 0 0 -0.5754 
+DEAL::1.187 -1.392 1.149 0 0 1.441 0 0 0 1.441 
+DEAL::0 0 0 4.038 -0.9884 0 0 0 0 
+DEAL::0 0 0 -0.9884 6.253 0 0 0 0 0 
+DEAL::2.351 -0.5754 1.441 0 0 2.856 0 0 0 2.856 
+DEAL::0 0 0 0 0 0 4.038 -0.9884 0 0 
+DEAL::0 0 0 0 0 0 -0.9884 6.253 0 0 
+DEAL::0 0 0 0 0 0 0 0 8.426 0 
+DEAL::2.351 -0.5754 1.441 0 0 2.856 0 0 0 2.856 
 DEAL::laplace_matrix:
 DEAL::0 0 0 0 0 0 0 0 0 0 
-DEAL::0 10.10 0.7895 0 -1.388 0 0 -1.388 0 0 
-DEAL::0 0.7895 44.49 0 0.09526 0 0 0.09526 0 0 
-DEAL::0 0 0 10.10 0.3531 -3.103 0 0 -1.388 0 
-DEAL::0 -1.388 0.09526 0.3531 16.04 0.09526 0 0.6338 0.04260 0 
-DEAL::0 0 0 -3.103 0.09526 35.69 0 0 1.417 0 
-DEAL::0 0 0 0 0 0 10.10 0.3531 -1.388 -3.103 
-DEAL::0 -1.388 0.09526 0 0.6338 0 0.3531 16.04 0.04260 0.09526 
-DEAL::0 0 0 -1.388 0.04260 1.417 -1.388 0.04260 14.28 1.417 
-DEAL::0 0 0 0 0 0 -3.103 0.09526 1.417 35.69 
+DEAL::0 23.23 -12.71 0 0 0 0 0 0 0 
+DEAL::0 -12.71 179.8 0 0 0 0 0 0 0 
+DEAL::0 0 0 23.23 -5.685 0 0 0 0 0 
+DEAL::0 0 0 -5.685 84.42 0 0 0 0 0 
+DEAL::0 0 0 0 0 242.3 0 0 0 0 
+DEAL::0 0 0 0 0 0 23.23 -5.685 0 0 
+DEAL::0 0 0 0 0 0 -5.685 84.42 0 0 
+DEAL::0 0 0 0 0 0 0 0 96.92 0 
+DEAL::0 0 0 0 0 0 0 0 0 242.3 
 DEAL::cell=1
 DEAL::mass_matrix:
 DEAL::1.000 3.464 13.42 0 0 0 0 0 0 0 
index 14154bbb5a9fb5a4656de7771aa39e912915cbd7..7b3d7112b01d2913ea59402d8b8fde5f34521068 100644 (file)
@@ -148,7 +148,7 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary;
+  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
   tr.set_boundary (1, boundary);
 
                                   // set boundary id on cell 1
index f6a9046ff08ea30574d218940001da082de3217d..eb5e1fad2db8b67c62e56ad39883465289a2de1b 100644 (file)
@@ -5,13 +5,13 @@ DEAL::Mass matrices are not equal.
 DEAL::Laplace matrices are not equal.
 DEAL::cell=0
 DEAL::mass_matrix:
-DEAL::nan nan nan 
-DEAL::nan nan nan 
-DEAL::nan nan nan 
+DEAL::1.522 -0.3883 0 
+DEAL::-0.3883 2.029 0 
+DEAL::0 0 2.411 
 DEAL::laplace_matrix:
-DEAL::nan nan nan 
-DEAL::nan nan nan 
-DEAL::nan nan nan 
+DEAL::0 0 0 
+DEAL::0 18.27 0 
+DEAL::0 0 18.27 
 DEAL::cell=1
 DEAL::mass_matrix:
 DEAL::1.000 3.464 0 
index c5294b5ca20bf4f6a7a9a4418b693808d8c7545c..b1391248eafd7091482676ce1ac783f03653b7a9 100644 (file)
@@ -148,7 +148,7 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary;
+  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
   tr.set_boundary (1, boundary);
 
                                   // set boundary id on cell 1
index 7365b3baedfdb8afcd620d3e4813815ce6d39511..275d2528d3b6946caca73d6b697c2184b11be99e 100644 (file)
@@ -5,19 +5,19 @@ DEAL::Mass matrices are not equal.
 DEAL::Laplace matrices are not equal.
 DEAL::cell=0
 DEAL::mass_matrix:
-DEAL::nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan 
+DEAL::1.522 -0.3883 0.5664 0 0 0.9937 
+DEAL::-0.3883 2.029 -0.8711 0 0 -0.2535 
+DEAL::0.5664 -0.8711 0.4843 0 0 0.3698 
+DEAL::0 0 0 2.411 -0.6150 0 
+DEAL::0 0 0 -0.6150 3.213 0 
+DEAL::0.9937 -0.2535 0.3698 0 0 0.6488 
 DEAL::laplace_matrix:
-DEAL::nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan 
-DEAL::nan nan nan nan nan nan 
+DEAL::0 0 0 0 0 0 
+DEAL::0 18.27 -10.42 0 0 0 
+DEAL::0 -10.42 121.7 0 0 0 
+DEAL::0 0 0 18.27 -4.659 0 
+DEAL::0 0 0 -4.659 53.28 0 
+DEAL::0 0 0 0 0 144.7 
 DEAL::cell=1
 DEAL::mass_matrix:
 DEAL::1.000 3.464 13.42 0 0 0 
index 3ce4f84f568b86d71f2b82c3d95f96a741cea7bf..40cb33556d558492f9a1e844ce32ce37e7a32a50 100644 (file)
@@ -148,7 +148,7 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary;
+  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
   tr.set_boundary (1, boundary);
 
                                   // set boundary id on cell 1
index f7d0266facd587e0248064f90fe9827019f9f135..d9acc58a0c05d02a78e7ee36a50440ff952720f7 100644 (file)
@@ -148,7 +148,7 @@ void test()
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
 
-  static const HyperBallBoundary<dim> boundary;
+  static const HyperBallBoundary<dim> boundary(tr.begin_active()->center());
   tr.set_boundary (1, boundary);
 
                                   // set boundary id on cell 1

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.