From 3ea00be96dc39277605edbd802369a900da49f47 Mon Sep 17 00:00:00 2001 From: kronbichler Date: Tue, 17 Sep 2013 14:38:55 +0000 Subject: [PATCH] More comments. One more test git-svn-id: https://svn.dealii.org/trunk@30766 0785d39b-7218-0410-832d-ea1e28bc413d --- .../include/deal.II/base/polynomial_space.h | 2 +- deal.II/include/deal.II/fe/fe_face.h | 8 +- tests/fe/fe_data_test.cc | 9 +- tests/fe/fe_data_test/cmp/generic | 286 +++++++++++++++--- 4 files changed, 260 insertions(+), 45 deletions(-) diff --git a/deal.II/include/deal.II/base/polynomial_space.h b/deal.II/include/deal.II/base/polynomial_space.h index 7865f3f581..56a9f26e9e 100644 --- a/deal.II/include/deal.II/base/polynomial_space.h +++ b/deal.II/include/deal.II/base/polynomial_space.h @@ -236,7 +236,7 @@ protected: * pi(x)pj(y)pk(z). */ void compute_index (const unsigned int n, - unsigned int (&index)[dim]) const; + unsigned int (&index)[dim>0?dim:1]) const; private: /** diff --git a/deal.II/include/deal.II/fe/fe_face.h b/deal.II/include/deal.II/fe/fe_face.h index e6464118ba..64bb4e4f2b 100644 --- a/deal.II/include/deal.II/fe/fe_face.h +++ b/deal.II/include/deal.II/fe/fe_face.h @@ -33,8 +33,9 @@ DEAL_II_NAMESPACE_OPEN * and two, the polynomials hence correspond to the usual Lagrange polynomials * on equidistant points. * - * This finite element is the trace space of FE_RaviartThomas on the - * faces and serves in hybridized methods. + * This finite element is the trace space of FE_RaviartThomas on the faces and + * serves in hybridized methods, e.g. in combination with the FE_DGQ + * element. Its use is demonstrated in the step-51 tutorial program. * * @note Since these are only finite elements on faces, only * FEFaceValues and FESubfaceValues will be able to extract reasonable @@ -136,6 +137,9 @@ private: * and undefined in the interior of the cells. The basis functions on * the faces are from Polynomials::Legendre. * + * This element is used in a hybridized method together with the FE_DGP + * element for the interior degrees of freedom. + * * @note Since these are only finite elements on faces, only * FEFaceValues and FESubfaceValues will be able to extract reasonable * values from any face polynomial. In order to make the use of diff --git a/tests/fe/fe_data_test.cc b/tests/fe/fe_data_test.cc index a6d48f1b3f..cecc03e951 100644 --- a/tests/fe/fe_data_test.cc +++ b/tests/fe/fe_data_test.cc @@ -76,13 +76,20 @@ void test_fe_datas() // Check vector elements in 2d and higher only if (dim>1) { - // Face elements + // Face Q elements fe_datas.push_back(new FE_FaceQ (0)); deallog << (*fe_datas.rbegin())->get_name() << std::endl; fe_datas.push_back(new FE_FaceQ (1)); deallog << (*fe_datas.rbegin())->get_name() << std::endl; fe_datas.push_back(new FE_FaceQ (3)); deallog << (*fe_datas.rbegin())->get_name() << std::endl; + // Face P elements + fe_datas.push_back(new FE_FaceP (0)); + deallog << (*fe_datas.rbegin())->get_name() << std::endl; + fe_datas.push_back(new FE_FaceP (1)); + deallog << (*fe_datas.rbegin())->get_name() << std::endl; + fe_datas.push_back(new FE_FaceP (3)); + deallog << (*fe_datas.rbegin())->get_name() << std::endl; // Vector DG elements fe_datas.push_back( new FE_DGRaviartThomas(0)); diff --git a/tests/fe/fe_data_test/cmp/generic b/tests/fe/fe_data_test/cmp/generic index 5387c99b32..11757b20ea 100644 --- a/tests/fe/fe_data_test/cmp/generic +++ b/tests/fe/fe_data_test/cmp/generic @@ -577,6 +577,9 @@ DEAL::FESystem<2>[FE_Q<2>(1)^2-FE_Q<2>(2)] DEAL::FE_FaceQ<2>(0) DEAL::FE_FaceQ<2>(1) DEAL::FE_FaceQ<2>(3) +DEAL::FE_FaceP<2>(0) +DEAL::FE_FaceP<2>(1) +DEAL::FE_FaceP<2>(3) DEAL::FE_DGRaviartThomas<2,2>(0) DEAL::FE_DGRaviartThomas<2,2>(1) DEAL::FE_DGNedelec<2,2>(0) @@ -1099,7 +1102,100 @@ DEAL::support on face 1: 4 5 6 7 DEAL::support on face 2: 8 9 10 11 DEAL::support on face 3: 12 13 14 15 DEAL:: -DEAL::fe_data[16]:FE_DGRaviartThomas<2,2>(0) +DEAL::fe_data[16]:FE_FaceP<2>(0) +DEAL::dofs_per_vertex=0 +DEAL::dofs_per_line=1 +DEAL::dofs_per_quad=0 +DEAL::dofs_per_hex=0 +DEAL::first_line_index=0 +DEAL::first_quad_index=4 +DEAL::first_hex_index=4 +DEAL::first_face_line_index=0 +DEAL::first_face_quad_index=1 +DEAL::dofs_per_face=1 +DEAL::dofs_per_cell=4 +DEAL::primitive=yes +DEAL::components=1 +DEAL::blocks=1:[4]->4 +DEAL::degree=0 +DEAL::conformity= L2 +DEAL::unit_support_points=0 +DEAL::unit_face_support_points=0 +DEAL::generalized_support_points=0 +DEAL::generalized_face_support_points=0 +DEAL::face_to_equivalent_cell_index: 0 +DEAL::face_to_cell_index: 0 +DEAL::face_to_cell_index: 1 +DEAL::face_to_cell_index: 2 +DEAL::face_to_cell_index: 3 +DEAL::support on face 0: 0 +DEAL::support on face 1: 1 +DEAL::support on face 2: 2 +DEAL::support on face 3: 3 +DEAL:: +DEAL::fe_data[17]:FE_FaceP<2>(1) +DEAL::dofs_per_vertex=0 +DEAL::dofs_per_line=2 +DEAL::dofs_per_quad=0 +DEAL::dofs_per_hex=0 +DEAL::first_line_index=0 +DEAL::first_quad_index=8 +DEAL::first_hex_index=8 +DEAL::first_face_line_index=0 +DEAL::first_face_quad_index=2 +DEAL::dofs_per_face=2 +DEAL::dofs_per_cell=8 +DEAL::primitive=yes +DEAL::components=1 +DEAL::blocks=1:[8]->8 +DEAL::degree=1 +DEAL::conformity= L2 +DEAL::unit_support_points=0 +DEAL::unit_face_support_points=0 +DEAL::generalized_support_points=0 +DEAL::generalized_face_support_points=0 +DEAL::face_to_equivalent_cell_index: 0 1 +DEAL::face_to_cell_index: 0 1 +DEAL::face_to_cell_index: 2 3 +DEAL::face_to_cell_index: 4 5 +DEAL::face_to_cell_index: 6 7 +DEAL::support on face 0: 0 1 +DEAL::support on face 1: 2 3 +DEAL::support on face 2: 4 5 +DEAL::support on face 3: 6 7 +DEAL:: +DEAL::fe_data[18]:FE_FaceP<2>(3) +DEAL::dofs_per_vertex=0 +DEAL::dofs_per_line=4 +DEAL::dofs_per_quad=0 +DEAL::dofs_per_hex=0 +DEAL::first_line_index=0 +DEAL::first_quad_index=16 +DEAL::first_hex_index=16 +DEAL::first_face_line_index=0 +DEAL::first_face_quad_index=4 +DEAL::dofs_per_face=4 +DEAL::dofs_per_cell=16 +DEAL::primitive=yes +DEAL::components=1 +DEAL::blocks=1:[16]->16 +DEAL::degree=3 +DEAL::conformity= L2 +DEAL::unit_support_points=0 +DEAL::unit_face_support_points=0 +DEAL::generalized_support_points=0 +DEAL::generalized_face_support_points=0 +DEAL::face_to_equivalent_cell_index: 0 1 2 3 +DEAL::face_to_cell_index: 0 1 2 3 +DEAL::face_to_cell_index: 4 5 6 7 +DEAL::face_to_cell_index: 8 9 10 11 +DEAL::face_to_cell_index: 12 13 14 15 +DEAL::support on face 0: 0 1 2 3 +DEAL::support on face 1: 4 5 6 7 +DEAL::support on face 2: 8 9 10 11 +DEAL::support on face 3: 12 13 14 15 +DEAL:: +DEAL::fe_data[19]:FE_DGRaviartThomas<2,2>(0) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=4 @@ -1130,7 +1226,7 @@ DEAL::support on face 1: 0 1 2 3 DEAL::support on face 2: 0 1 2 3 DEAL::support on face 3: 0 1 2 3 DEAL:: -DEAL::fe_data[17]:FE_DGRaviartThomas<2,2>(1) +DEAL::fe_data[20]:FE_DGRaviartThomas<2,2>(1) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=12 @@ -1161,7 +1257,7 @@ DEAL::support on face 1: 0 1 2 3 4 5 6 7 8 9 10 11 DEAL::support on face 2: 0 1 2 3 4 5 6 7 8 9 10 11 DEAL::support on face 3: 0 1 2 3 4 5 6 7 8 9 10 11 DEAL:: -DEAL::fe_data[18]:FE_DGNedelec<2,2>(0) +DEAL::fe_data[21]:FE_DGNedelec<2,2>(0) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=4 @@ -1192,7 +1288,7 @@ DEAL::support on face 1: 0 1 2 3 DEAL::support on face 2: 0 1 2 3 DEAL::support on face 3: 0 1 2 3 DEAL:: -DEAL::fe_data[19]:FE_DGNedelec<2,2>(1) +DEAL::fe_data[22]:FE_DGNedelec<2,2>(1) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=12 @@ -1223,7 +1319,7 @@ DEAL::support on face 1: 0 1 2 3 4 5 6 7 8 9 10 11 DEAL::support on face 2: 0 1 2 3 4 5 6 7 8 9 10 11 DEAL::support on face 3: 0 1 2 3 4 5 6 7 8 9 10 11 DEAL:: -DEAL::fe_data[20]:FE_RaviartThomas<2>(0) +DEAL::fe_data[23]:FE_RaviartThomas<2>(0) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=1 DEAL::dofs_per_quad=0 @@ -1254,7 +1350,7 @@ DEAL::support on face 1: 1 2 3 DEAL::support on face 2: 0 1 2 DEAL::support on face 3: 0 1 3 DEAL:: -DEAL::fe_data[21]:FE_RaviartThomas<2>(1) +DEAL::fe_data[24]:FE_RaviartThomas<2>(1) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=2 DEAL::dofs_per_quad=4 @@ -1285,7 +1381,7 @@ DEAL::support on face 1: 0 1 2 3 4 5 6 7 8 9 10 11 DEAL::support on face 2: 0 1 2 3 4 5 6 7 8 9 10 11 DEAL::support on face 3: 0 1 2 3 4 5 6 7 8 9 10 11 DEAL:: -DEAL::fe_data[22]:FE_RaviartThomas<2>(2) +DEAL::fe_data[25]:FE_RaviartThomas<2>(2) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=3 DEAL::dofs_per_quad=12 @@ -1316,7 +1412,7 @@ DEAL::support on face 1: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 DEAL::support on face 2: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 DEAL::support on face 3: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 DEAL:: -DEAL::fe_data[23]:FESystem<2>[FE_RaviartThomas<2>(1)-FE_DGQ<2>(1)] +DEAL::fe_data[26]:FESystem<2>[FE_RaviartThomas<2>(1)-FE_DGQ<2>(1)] DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=2 DEAL::dofs_per_quad=8 @@ -1347,7 +1443,7 @@ DEAL::support on face 1: 0 1 2 3 4 5 6 7 8 9 10 11 13 15 DEAL::support on face 2: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 DEAL::support on face 3: 0 1 2 3 4 5 6 7 8 9 10 11 14 15 DEAL:: -DEAL::fe_data[24]:FE_Nedelec<2>(0) +DEAL::fe_data[27]:FE_Nedelec<2>(0) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=1 DEAL::dofs_per_quad=0 @@ -1378,7 +1474,7 @@ DEAL::support on face 1: 1 2 3 DEAL::support on face 2: 0 1 2 DEAL::support on face 3: 0 1 3 DEAL:: -DEAL::fe_data[25]:FE_Nedelec<2>(1) +DEAL::fe_data[28]:FE_Nedelec<2>(1) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=2 DEAL::dofs_per_quad=4 @@ -1409,7 +1505,7 @@ DEAL::support on face 1: 2 3 4 5 6 7 DEAL::support on face 2: 0 1 2 3 4 5 DEAL::support on face 3: 0 1 2 3 6 7 8 9 10 11 DEAL:: -DEAL::fe_data[26]:FE_DGBDM<2,2>(1) +DEAL::fe_data[29]:FE_DGBDM<2,2>(1) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=8 @@ -1440,7 +1536,7 @@ DEAL::support on face 1: 0 1 2 3 4 5 6 7 DEAL::support on face 2: 0 1 2 3 4 5 6 7 DEAL::support on face 3: 0 1 2 3 4 5 6 7 DEAL:: -DEAL::fe_data[27]:FE_DGBDM<2,2>(2) +DEAL::fe_data[30]:FE_DGBDM<2,2>(2) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=14 @@ -1471,7 +1567,7 @@ DEAL::support on face 1: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 DEAL::support on face 2: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 DEAL::support on face 3: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 DEAL:: -DEAL::fe_data[28]:FE_BDM<2>(1) +DEAL::fe_data[31]:FE_BDM<2>(1) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=2 DEAL::dofs_per_quad=0 @@ -1502,7 +1598,7 @@ DEAL::support on face 1: 0 1 2 3 4 5 6 7 DEAL::support on face 2: 0 1 2 3 4 5 6 7 DEAL::support on face 3: 0 1 2 3 4 5 6 7 DEAL:: -DEAL::fe_data[29]:FE_BDM<2>(2) +DEAL::fe_data[32]:FE_BDM<2>(2) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=3 DEAL::dofs_per_quad=2 @@ -1533,7 +1629,7 @@ DEAL::support on face 1: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 DEAL::support on face 2: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 DEAL::support on face 3: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 DEAL:: -DEAL::fe_data[30]:FE_RaviartThomasNodal<2>(0) +DEAL::fe_data[33]:FE_RaviartThomasNodal<2>(0) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=1 DEAL::dofs_per_quad=0 @@ -1564,7 +1660,7 @@ DEAL::support on face 1: 1 2 3 DEAL::support on face 2: 0 1 2 DEAL::support on face 3: 0 1 3 DEAL:: -DEAL::fe_data[31]:FE_RaviartThomasNodal<2>(1) +DEAL::fe_data[34]:FE_RaviartThomasNodal<2>(1) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=2 DEAL::dofs_per_quad=4 @@ -1595,7 +1691,7 @@ DEAL::support on face 1: 2 3 4 5 6 7 8 9 10 11 DEAL::support on face 2: 0 1 2 3 4 5 8 9 10 11 DEAL::support on face 3: 0 1 2 3 6 7 8 9 10 11 DEAL:: -DEAL::fe_data[32]:FESystem<2>[FE_RaviartThomasNodal<2>(1)-FE_DGQ<2>(1)] +DEAL::fe_data[35]:FESystem<2>[FE_RaviartThomasNodal<2>(1)-FE_DGQ<2>(1)] DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=2 DEAL::dofs_per_quad=8 @@ -1626,7 +1722,7 @@ DEAL::support on face 1: 2 3 4 5 6 7 8 9 10 11 13 15 DEAL::support on face 2: 0 1 2 3 4 5 8 9 10 11 12 13 DEAL::support on face 3: 0 1 2 3 6 7 8 9 10 11 14 15 DEAL:: -DEAL::fe_data[33]:FESystem<2>[FE_Q<2>(3)^2] +DEAL::fe_data[36]:FESystem<2>[FE_Q<2>(3)^2] DEAL::dofs_per_vertex=2 DEAL::dofs_per_line=4 DEAL::dofs_per_quad=8 @@ -1657,7 +1753,7 @@ DEAL::support on face 1: 2 3 6 7 12 13 14 15 DEAL::support on face 2: 0 1 2 3 16 17 18 19 DEAL::support on face 3: 4 5 6 7 20 21 22 23 DEAL:: -DEAL::fe_data[34]:FESystem<2>[FE_Q<2>(1)^2-FE_Q<2>(3)] +DEAL::fe_data[37]:FESystem<2>[FE_Q<2>(1)^2-FE_Q<2>(3)] DEAL::dofs_per_vertex=3 DEAL::dofs_per_line=2 DEAL::dofs_per_quad=4 @@ -1688,7 +1784,7 @@ DEAL::support on face 1: 3 4 5 9 10 11 14 15 DEAL::support on face 2: 0 1 2 3 4 5 16 17 DEAL::support on face 3: 6 7 8 9 10 11 18 19 DEAL:: -DEAL::fe_data[35]:FESystem<2>[FE_Q<2>(4)^2] +DEAL::fe_data[38]:FESystem<2>[FE_Q<2>(4)^2] DEAL::dofs_per_vertex=2 DEAL::dofs_per_line=6 DEAL::dofs_per_quad=18 @@ -1719,7 +1815,7 @@ DEAL::support on face 1: 2 3 6 7 14 15 16 17 18 19 DEAL::support on face 2: 0 1 2 3 20 21 22 23 24 25 DEAL::support on face 3: 4 5 6 7 26 27 28 29 30 31 DEAL:: -DEAL::fe_data[36]:FESystem<2>[FESystem<2>[FE_Q<2>(1)^2]^2] +DEAL::fe_data[39]:FESystem<2>[FESystem<2>[FE_Q<2>(1)^2]^2] DEAL::dofs_per_vertex=4 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=0 @@ -1750,7 +1846,7 @@ DEAL::support on face 1: 4 5 6 7 12 13 14 15 DEAL::support on face 2: 0 1 2 3 4 5 6 7 DEAL::support on face 3: 8 9 10 11 12 13 14 15 DEAL:: -DEAL::fe_data[37]:FESystem<2>[FESystem<2>[FE_Q<2>(1)^2]-FESystem<2>[FE_DGQ<2>(1)^2]] +DEAL::fe_data[40]:FESystem<2>[FESystem<2>[FE_Q<2>(1)^2]-FESystem<2>[FE_DGQ<2>(1)^2]] DEAL::dofs_per_vertex=2 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=8 @@ -1781,7 +1877,7 @@ DEAL::support on face 1: 2 3 6 7 9 11 13 15 DEAL::support on face 2: 0 1 2 3 8 9 12 13 DEAL::support on face 3: 4 5 6 7 10 11 14 15 DEAL:: -DEAL::fe_data[38]:FESystem<2>[FESystem<2>[FE_Q<2>(1)-FE_Q<2>(2)]-FESystem<2>[FE_Q<2>(2)^2]-FESystem<2>[FE_DGQ<2>(2)^2]] +DEAL::fe_data[41]:FESystem<2>[FESystem<2>[FE_Q<2>(1)-FE_Q<2>(2)]-FESystem<2>[FE_Q<2>(2)^2]-FESystem<2>[FE_DGQ<2>(2)^2]] DEAL::dofs_per_vertex=4 DEAL::dofs_per_line=3 DEAL::dofs_per_quad=21 @@ -1812,7 +1908,7 @@ DEAL::support on face 1: 4 5 6 7 12 13 14 15 19 20 21 33 36 39 42 45 48 DEAL::support on face 2: 0 1 2 3 4 5 6 7 22 23 24 31 32 33 40 41 42 DEAL::support on face 3: 8 9 10 11 12 13 14 15 25 26 27 37 38 39 46 47 48 DEAL:: -DEAL::fe_data[39]:FESystem<2>[FESystem<2>[FESystem<2>[FE_Q<2>(1)^2]^2]^2-FESystem<2>[FESystem<2>[FE_Q<2>(1)^2]-FESystem<2>[FE_DGQ<2>(1)^2]]-FESystem<2>[FESystem<2>[FE_Q<2>(1)-FE_Q<2>(2)]-FESystem<2>[FE_Q<2>(2)^2]-FESystem<2>[FE_DGQ<2>(2)^2]]^2] +DEAL::fe_data[42]:FESystem<2>[FESystem<2>[FESystem<2>[FE_Q<2>(1)^2]^2]^2-FESystem<2>[FESystem<2>[FE_Q<2>(1)^2]-FESystem<2>[FE_DGQ<2>(1)^2]]-FESystem<2>[FESystem<2>[FE_Q<2>(1)-FE_Q<2>(2)]-FESystem<2>[FE_Q<2>(2)^2]-FESystem<2>[FE_DGQ<2>(2)^2]]^2] DEAL::dofs_per_vertex=18 DEAL::dofs_per_line=6 DEAL::dofs_per_quad=50 @@ -1859,6 +1955,9 @@ DEAL::FESystem<3>[FE_Q<3>(1)^2-FE_Q<3>(2)] DEAL::FE_FaceQ<3>(0) DEAL::FE_FaceQ<3>(1) DEAL::FE_FaceQ<3>(3) +DEAL::FE_FaceP<3>(0) +DEAL::FE_FaceP<3>(1) +DEAL::FE_FaceP<3>(3) DEAL::FE_DGRaviartThomas<3,3>(0) DEAL::FE_DGRaviartThomas<3,3>(1) DEAL::FE_DGNedelec<3,3>(0) @@ -2438,7 +2537,112 @@ DEAL::support on face 3: 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 DEAL::support on face 4: 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 DEAL::support on face 5: 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 DEAL:: -DEAL::fe_data[16]:FE_DGRaviartThomas<3,3>(0) +DEAL::fe_data[16]:FE_FaceP<3>(0) +DEAL::dofs_per_vertex=0 +DEAL::dofs_per_line=0 +DEAL::dofs_per_quad=1 +DEAL::dofs_per_hex=0 +DEAL::first_line_index=0 +DEAL::first_quad_index=0 +DEAL::first_hex_index=6 +DEAL::first_face_line_index=0 +DEAL::first_face_quad_index=0 +DEAL::dofs_per_face=1 +DEAL::dofs_per_cell=6 +DEAL::primitive=yes +DEAL::components=1 +DEAL::blocks=1:[6]->6 +DEAL::degree=0 +DEAL::conformity= L2 +DEAL::unit_support_points=0 +DEAL::unit_face_support_points=0 +DEAL::generalized_support_points=0 +DEAL::generalized_face_support_points=0 +DEAL::face_to_equivalent_cell_index: 0 +DEAL::face_to_cell_index: 0 +DEAL::face_to_cell_index: 1 +DEAL::face_to_cell_index: 2 +DEAL::face_to_cell_index: 3 +DEAL::face_to_cell_index: 4 +DEAL::face_to_cell_index: 5 +DEAL::support on face 0: 0 +DEAL::support on face 1: 1 +DEAL::support on face 2: 2 +DEAL::support on face 3: 3 +DEAL::support on face 4: 4 +DEAL::support on face 5: 5 +DEAL:: +DEAL::fe_data[17]:FE_FaceP<3>(1) +DEAL::dofs_per_vertex=0 +DEAL::dofs_per_line=0 +DEAL::dofs_per_quad=3 +DEAL::dofs_per_hex=0 +DEAL::first_line_index=0 +DEAL::first_quad_index=0 +DEAL::first_hex_index=18 +DEAL::first_face_line_index=0 +DEAL::first_face_quad_index=0 +DEAL::dofs_per_face=3 +DEAL::dofs_per_cell=18 +DEAL::primitive=yes +DEAL::components=1 +DEAL::blocks=1:[18]->18 +DEAL::degree=1 +DEAL::conformity= L2 +DEAL::unit_support_points=0 +DEAL::unit_face_support_points=0 +DEAL::generalized_support_points=0 +DEAL::generalized_face_support_points=0 +DEAL::face_to_equivalent_cell_index: 0 1 2 +DEAL::face_to_cell_index: 0 1 2 +DEAL::face_to_cell_index: 3 4 5 +DEAL::face_to_cell_index: 6 7 8 +DEAL::face_to_cell_index: 9 10 11 +DEAL::face_to_cell_index: 12 13 14 +DEAL::face_to_cell_index: 15 16 17 +DEAL::support on face 0: 0 1 2 +DEAL::support on face 1: 3 4 5 +DEAL::support on face 2: 6 7 8 +DEAL::support on face 3: 9 10 11 +DEAL::support on face 4: 12 13 14 +DEAL::support on face 5: 15 16 17 +DEAL:: +DEAL::fe_data[18]:FE_FaceP<3>(3) +DEAL::dofs_per_vertex=0 +DEAL::dofs_per_line=0 +DEAL::dofs_per_quad=10 +DEAL::dofs_per_hex=0 +DEAL::first_line_index=0 +DEAL::first_quad_index=0 +DEAL::first_hex_index=60 +DEAL::first_face_line_index=0 +DEAL::first_face_quad_index=0 +DEAL::dofs_per_face=10 +DEAL::dofs_per_cell=60 +DEAL::primitive=yes +DEAL::components=1 +DEAL::blocks=1:[60]->60 +DEAL::degree=3 +DEAL::conformity= L2 +DEAL::unit_support_points=0 +DEAL::unit_face_support_points=0 +DEAL::generalized_support_points=0 +DEAL::generalized_face_support_points=0 +DEAL::face_to_equivalent_cell_index: 0 1 2 3 4 5 6 7 8 9 +DEAL::face_to_cell_index: 0 1 2 3 4 5 6 7 8 9 +DEAL::face_to_cell_index: 10 11 12 13 14 15 16 17 18 19 +DEAL::face_to_cell_index: 20 21 22 23 24 25 26 27 28 29 +DEAL::face_to_cell_index: 30 31 32 33 34 35 36 37 38 39 +DEAL::face_to_cell_index: 40 41 42 43 44 45 46 47 48 49 +DEAL::face_to_cell_index: 50 51 52 53 54 55 56 57 58 59 +DEAL::support on face 0: 0 1 2 3 4 5 6 7 8 9 +DEAL::support on face 1: 10 11 12 13 14 15 16 17 18 19 +DEAL::support on face 2: 20 21 22 23 24 25 26 27 28 29 +DEAL::support on face 3: 30 31 32 33 34 35 36 37 38 39 +DEAL::support on face 4: 40 41 42 43 44 45 46 47 48 49 +DEAL::support on face 5: 50 51 52 53 54 55 56 57 58 59 +DEAL:: +DEAL::fe_data[19]:FE_DGRaviartThomas<3,3>(0) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=0 @@ -2473,7 +2677,7 @@ DEAL::support on face 3: 0 1 2 3 4 5 DEAL::support on face 4: 0 1 2 3 4 5 DEAL::support on face 5: 0 1 2 3 4 5 DEAL:: -DEAL::fe_data[17]:FE_DGRaviartThomas<3,3>(1) +DEAL::fe_data[20]:FE_DGRaviartThomas<3,3>(1) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=0 @@ -2508,7 +2712,7 @@ DEAL::support on face 3: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 DEAL::support on face 4: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 DEAL::support on face 5: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 DEAL:: -DEAL::fe_data[18]:FE_DGNedelec<3,3>(0) +DEAL::fe_data[21]:FE_DGNedelec<3,3>(0) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=0 @@ -2543,7 +2747,7 @@ DEAL::support on face 3: 0 1 2 3 4 5 6 7 8 9 10 11 DEAL::support on face 4: 0 1 2 3 4 5 6 7 8 9 10 11 DEAL::support on face 5: 0 1 2 3 4 5 6 7 8 9 10 11 DEAL:: -DEAL::fe_data[19]:FE_DGNedelec<3,3>(1) +DEAL::fe_data[22]:FE_DGNedelec<3,3>(1) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=0 @@ -2578,7 +2782,7 @@ DEAL::support on face 3: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 DEAL::support on face 4: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 DEAL::support on face 5: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 DEAL:: -DEAL::fe_data[20]:FE_RaviartThomas<3>(0) +DEAL::fe_data[23]:FE_RaviartThomas<3>(0) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=1 @@ -2613,7 +2817,7 @@ DEAL::support on face 3: 0 1 2 3 4 5 DEAL::support on face 4: 0 1 2 3 4 5 DEAL::support on face 5: 0 1 2 3 4 5 DEAL:: -DEAL::fe_data[21]:FE_RaviartThomas<3>(1) +DEAL::fe_data[24]:FE_RaviartThomas<3>(1) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=4 @@ -2648,7 +2852,7 @@ DEAL::support on face 3: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 DEAL::support on face 4: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 DEAL::support on face 5: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 DEAL:: -DEAL::fe_data[22]:FE_RaviartThomas<3>(2) +DEAL::fe_data[25]:FE_RaviartThomas<3>(2) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=9 @@ -2683,7 +2887,7 @@ DEAL::support on face 3: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 DEAL::support on face 4: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 DEAL::support on face 5: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 DEAL:: -DEAL::fe_data[23]:FESystem<3>[FE_RaviartThomas<3>(1)-FE_DGQ<3>(1)] +DEAL::fe_data[26]:FESystem<3>[FE_RaviartThomas<3>(1)-FE_DGQ<3>(1)] DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=4 @@ -2718,7 +2922,7 @@ DEAL::support on face 3: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 DEAL::support on face 4: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 DEAL::support on face 5: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 40 41 42 43 DEAL:: -DEAL::fe_data[24]:FE_Nedelec<3>(0) +DEAL::fe_data[27]:FE_Nedelec<3>(0) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=1 DEAL::dofs_per_quad=0 @@ -2753,7 +2957,7 @@ DEAL::support on face 3: 0 1 3 4 5 7 10 11 DEAL::support on face 4: 0 1 2 3 8 9 10 11 DEAL::support on face 5: 4 5 6 7 8 9 10 11 DEAL:: -DEAL::fe_data[25]:FE_Nedelec<3>(1) +DEAL::fe_data[28]:FE_Nedelec<3>(1) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=2 DEAL::dofs_per_quad=4 @@ -2788,7 +2992,7 @@ DEAL::support on face 3: 0 1 2 3 6 7 8 9 10 11 14 15 20 21 22 23 26 27 30 31 36 DEAL::support on face 4: 0 1 2 3 4 5 6 7 16 17 18 19 20 21 22 23 24 25 28 29 34 35 38 39 40 41 42 43 DEAL::support on face 5: 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 28 29 34 35 38 39 44 45 46 47 DEAL:: -DEAL::fe_data[26]:FE_RaviartThomasNodal<3>(0) +DEAL::fe_data[29]:FE_RaviartThomasNodal<3>(0) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=1 @@ -2823,7 +3027,7 @@ DEAL::support on face 3: 0 1 3 4 5 DEAL::support on face 4: 0 1 2 3 4 DEAL::support on face 5: 0 1 2 3 5 DEAL:: -DEAL::fe_data[27]:FE_RaviartThomasNodal<3>(1) +DEAL::fe_data[30]:FE_RaviartThomasNodal<3>(1) DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=4 @@ -2858,7 +3062,7 @@ DEAL::support on face 3: 0 1 2 3 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 DEAL::support on face 4: 0 1 2 3 4 5 6 7 8 9 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 DEAL::support on face 5: 0 1 2 3 4 5 6 7 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 DEAL:: -DEAL::fe_data[28]:FESystem<3>[FE_RaviartThomasNodal<3>(1)-FE_DGQ<3>(1)] +DEAL::fe_data[31]:FESystem<3>[FE_RaviartThomasNodal<3>(1)-FE_DGQ<3>(1)] DEAL::dofs_per_vertex=0 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=4 @@ -2893,7 +3097,7 @@ DEAL::support on face 3: 0 1 2 3 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 DEAL::support on face 4: 0 1 2 3 4 5 6 7 8 9 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 DEAL::support on face 5: 0 1 2 3 4 5 6 7 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 40 41 42 43 DEAL:: -DEAL::fe_data[29]:FESystem<3>[FESystem<3>[FE_Q<3>(1)^2]^2] +DEAL::fe_data[32]:FESystem<3>[FESystem<3>[FE_Q<3>(1)^2]^2] DEAL::dofs_per_vertex=4 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=0 @@ -2928,7 +3132,7 @@ DEAL::support on face 3: 8 9 10 11 12 13 14 15 24 25 26 27 28 29 30 31 DEAL::support on face 4: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 DEAL::support on face 5: 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 DEAL:: -DEAL::fe_data[30]:FESystem<3>[FESystem<3>[FE_Q<3>(1)^2]-FESystem<3>[FE_DGQ<3>(1)^2]] +DEAL::fe_data[33]:FESystem<3>[FESystem<3>[FE_Q<3>(1)^2]-FESystem<3>[FE_DGQ<3>(1)^2]] DEAL::dofs_per_vertex=2 DEAL::dofs_per_line=0 DEAL::dofs_per_quad=0 @@ -2963,7 +3167,7 @@ DEAL::support on face 3: 4 5 6 7 12 13 14 15 18 19 22 23 26 27 30 31 DEAL::support on face 4: 0 1 2 3 4 5 6 7 16 17 18 19 24 25 26 27 DEAL::support on face 5: 8 9 10 11 12 13 14 15 20 21 22 23 28 29 30 31 DEAL:: -DEAL::fe_data[31]:FESystem<3>[FESystem<3>[FE_Q<3>(1)-FE_Q<3>(2)]-FESystem<3>[FE_Q<3>(2)^2]-FESystem<3>[FE_DGQ<3>(2)^2]] +DEAL::fe_data[34]:FESystem<3>[FESystem<3>[FE_Q<3>(1)-FE_Q<3>(2)]-FESystem<3>[FE_Q<3>(2)^2]-FESystem<3>[FE_DGQ<3>(2)^2]] DEAL::dofs_per_vertex=4 DEAL::dofs_per_line=3 DEAL::dofs_per_quad=3 @@ -2998,7 +3202,7 @@ DEAL::support on face 3: 8 9 10 11 12 13 14 15 24 25 26 27 28 29 30 31 41 42 43 DEAL::support on face 4: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 32 33 34 35 36 37 38 39 40 41 42 43 80 81 82 89 90 91 92 93 94 95 96 97 116 117 118 119 120 121 122 123 124 DEAL::support on face 5: 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 44 45 46 47 48 49 50 51 52 53 54 55 83 84 85 107 108 109 110 111 112 113 114 115 134 135 136 137 138 139 140 141 142 DEAL:: -DEAL::fe_data[32]:FESystem<3>[FESystem<3>[FESystem<3>[FE_Q<3>(1)^2]^2]^2-FESystem<3>[FESystem<3>[FE_Q<3>(1)^2]-FESystem<3>[FE_DGQ<3>(1)^2]]-FESystem<3>[FESystem<3>[FE_Q<3>(1)-FE_Q<3>(2)]-FESystem<3>[FE_Q<3>(2)^2]-FESystem<3>[FE_DGQ<3>(2)^2]]^2] +DEAL::fe_data[35]:FESystem<3>[FESystem<3>[FESystem<3>[FE_Q<3>(1)^2]^2]^2-FESystem<3>[FESystem<3>[FE_Q<3>(1)^2]-FESystem<3>[FE_DGQ<3>(1)^2]]-FESystem<3>[FESystem<3>[FE_Q<3>(1)-FE_Q<3>(2)]-FESystem<3>[FE_Q<3>(2)^2]-FESystem<3>[FE_DGQ<3>(2)^2]]^2] DEAL::dofs_per_vertex=18 DEAL::dofs_per_line=6 DEAL::dofs_per_quad=6 -- 2.39.5