More comments. One more test

git-svn-id: https://svn.dealii.org/trunk@30766 0785d39b-7218-0410-832d-ea1e28bc413d

diff --git a/deal.II/include/deal.II/base/polynomial_space.h b/deal.II/include/deal.II/base/polynomial_space.h
index 7865f3f5818e09e5f05b740fd90fbf47623e4e76..56a9f26e9e0d6ad7045273f4a3a1a6220e798741 100644 (file)
--- 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:
    * p<sub>i</sub>(x)p<sub>j</sub>(y)p<sub>k</sub>(z)</i>.
    */
   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 e6464118ba83f0fdf89d31dba10e84c855fa23be..64bb4e4f2b085467a5486154d9efa97ecd52da4e 100644 (file)
--- 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 a6d48f1b3fa76f391f636d3c66aafdba569aadc6..cecc03e951888d121c0769fecf32ee6545348009 100644 (file)
--- 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<dim> (0));
       deallog << (*fe_datas.rbegin())->get_name() << std::endl;
       fe_datas.push_back(new FE_FaceQ<dim> (1));
       deallog << (*fe_datas.rbegin())->get_name() << std::endl;
       fe_datas.push_back(new FE_FaceQ<dim> (3));
       deallog << (*fe_datas.rbegin())->get_name() << std::endl;
+      // Face P elements
+      fe_datas.push_back(new FE_FaceP<dim> (0));
+      deallog << (*fe_datas.rbegin())->get_name() << std::endl;
+      fe_datas.push_back(new FE_FaceP<dim> (1));
+      deallog << (*fe_datas.rbegin())->get_name() << std::endl;
+      fe_datas.push_back(new FE_FaceP<dim> (3));
+      deallog << (*fe_datas.rbegin())->get_name() << std::endl;
       // Vector DG elements
       fe_datas.push_back(
         new FE_DGRaviartThomas<dim>(0));

diff --git a/tests/fe/fe_data_test/cmp/generic b/tests/fe/fe_data_test/cmp/generic
index 5387c99b32951b955968bf59b1726338b80057e4..11757b20eaab56162107421907385560b4406b68 100644 (file)
--- 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