]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
new function tensor_degree in FiniteElementData
authorguido <guido@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 28 Nov 2003 11:52:35 +0000 (11:52 +0000)
committerguido <guido@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 28 Nov 2003 11:52:35 +0000 (11:52 +0000)
git-svn-id: https://svn.dealii.org/trunk@8207 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/include/fe/fe_base.h
deal.II/deal.II/source/fe/fe_data.cc
deal.II/deal.II/source/fe/fe_dgp.cc
deal.II/deal.II/source/fe/fe_dgq.cc
deal.II/deal.II/source/fe/fe_q.cc
deal.II/deal.II/source/fe/fe_q_hierarchical.cc
deal.II/deal.II/source/fe/fe_raviart_thomas.cc
deal.II/deal.II/source/fe/fe_system.cc
deal.II/doc/news/c-4-0.html
tests/fe/fe_data_test.cc
tests/results/i686-pc-linux-gnu+gcc3.2/fe/fe_data_test.output

index b3561e90b72c18771057f50a7e773263b124d608..e3e6c43a76028daba7d5ffce59ffebeff47c9662 100644 (file)
@@ -45,7 +45,7 @@ template<int dim> class FESystem;
  * @ref{FiniteElementBase} should be private base classes of
  * @ref{FiniteElement}.
  *
- * @author Wolfgang Bangerth, Guido Kanschat, 1998, 1999, 2000, 2001
+ * @author Wolfgang Bangerth, Guido Kanschat, 1998, 1999, 2000, 2001, 2003
  */
 template<int dim>
 class FiniteElementData
@@ -137,6 +137,13 @@ class FiniteElementData
                                      */
     const unsigned int components;
 
+                                    /**
+                                     * Maximal polynomial degree of a
+                                     * shape function in a single
+                                     * coordinate direction.
+                                     */
+    const unsigned int degree;
+    
                                     /**
                                      * Default
                                      * constructor. Constructs an
@@ -162,9 +169,17 @@ class FiniteElementData
                                      *
                                      * Hence this constructor requires
                                      * @p{dofs_per_object.size()==dim+1}.
+                                     *
+                                     * @param n_components Number of
+                                     * vector components of the
+                                     * element.
+                                     * @param degree
+                                     * Maximal polynomial degree in a
+                                     * single direction.
                                      */
     FiniteElementData (const std::vector<unsigned int> &dofs_per_object,
-                      const unsigned int n_components);
+                      const unsigned int n_components,
+                      const unsigned int degree = deal_II_numbers::invalid_unsigned_int);
 
                                     /**
                                      * Return the @p{dofs_per_vertex}.
@@ -201,6 +216,17 @@ class FiniteElementData
                                      */
     unsigned int n_components () const;
 
+                                    /**
+                                     * Maximal polynomial degree of a
+                                     * shape function in a single
+                                     * coordinate direction.
+                                     *
+                                     * This function can be used to
+                                     * determine the optimal
+                                     * quadrature rule.
+                                     */
+    unsigned int tensor_degree () const;
+    
                                     /**
                                      * Comparison operator.
                                      */
@@ -1778,6 +1804,16 @@ FiniteElementData<dim>::n_components () const
 
 
 
+template <int dim>
+inline
+unsigned int 
+FiniteElementData<dim>::tensor_degree () const
+{
+  return degree;
+}
+
+
+
 template <int dim>  
 inline
 std::pair<unsigned int,unsigned int>
index bd13f0ee2585cc60783ccd836fcf80703eea00ac..6ecadf0b5aa8543c39e3902f47eb61c3490eadfa 100644 (file)
@@ -18,7 +18,8 @@
 
 template <int dim>
 FiniteElementData<dim>::FiniteElementData (const std::vector<unsigned int> &dofs_per_object,
-                                          const unsigned int n_components) :
+                                          const unsigned int n_components,
+                                          const unsigned int degree) :
                dofs_per_vertex(dofs_per_object[0]),
                dofs_per_line(dofs_per_object[1]),
                dofs_per_quad(dim>1? dofs_per_object[2]:0),
@@ -45,7 +46,8 @@ FiniteElementData<dim>::FiniteElementData (const std::vector<unsigned int> &dofs
                               GeometryInfo<dim>::lines_per_cell * dofs_per_line +
                               GeometryInfo<dim>::quads_per_cell * dofs_per_quad +
                               GeometryInfo<dim>::hexes_per_cell * dofs_per_hex),
-               components(n_components)
+               components(n_components),
+               degree(degree)
 {
   Assert(dofs_per_object.size()==dim+1, ExcDimensionMismatch(dofs_per_object.size()-1,dim));
 }
index 90d06302b15b1809e69359cf6c931365d5ff3290..9230f796cf5dbe0dfbafc5d0d1433d5a93d0f157 100644 (file)
@@ -30,9 +30,9 @@
 template <int dim>
 FE_DGP<dim>::FE_DGP (const unsigned int degree)
                :
-               FiniteElement<dim> (FiniteElementData<dim>(get_dpo_vector(degree),1),
-                                   std::vector<bool>(FiniteElementData<dim>(get_dpo_vector(degree),1).dofs_per_cell,true),
-                                   std::vector<std::vector<bool> >(FiniteElementData<dim>(get_dpo_vector(degree),1).dofs_per_cell,
+               FiniteElement<dim> (FiniteElementData<dim>(get_dpo_vector(degree), 1, degree),
+                                   std::vector<bool>(FiniteElementData<dim>(get_dpo_vector(degree), 1, degree).dofs_per_cell,true),
+                                   std::vector<std::vector<bool> >(FiniteElementData<dim>(get_dpo_vector(degree), 1, degree).dofs_per_cell,
                                                                    std::vector<bool>(1,true))),
                 degree(degree),
                 polynomial_space (Polynomials::Legendre::generate_complete_basis(degree))
index 7e1d50b4c4717ff66e70d7f47aea531131ae205f..f1286405567010692ef6fb28866f216a1b0e08f0 100644 (file)
@@ -45,7 +45,7 @@ namespace
                                   // given an integer N, compute its
                                   // integer square root (if it
                                   // exists, otherwise give up)
-  unsigned int int_sqrt (const unsigned int N)
+  inline unsigned int int_sqrt (const unsigned int N)
   {
     for (unsigned int i=0; i<=N; ++i)
       if (i*i == N)
@@ -58,7 +58,7 @@ namespace
                                   // given an integer N, compute its
                                   // integer cube root (if it
                                   // exists, otherwise give up)
-  unsigned int int_cuberoot (const unsigned int N)
+  inline unsigned int int_cuberoot (const unsigned int N)
   {
     for (unsigned int i=0; i<=N; ++i)
       if (i*i*i == N)
@@ -71,7 +71,7 @@ namespace
                                   // given N, generate i=0...N-1
                                   // equidistant points in the
                                   // interior of the interval [0,1]
-  Point<1>
+  inline Point<1>
   generate_unit_point (const unsigned int i,
                       const unsigned int N,
                       const int2type<1>  )
@@ -90,7 +90,7 @@ namespace
                                   // given N, generate i=0...N-1
                                   // equidistant points in the domain
                                   // [0,1]^2
-  Point<2>
+  inline Point<2>
   generate_unit_point (const unsigned int i,
                       const unsigned int N,
                       const int2type<2>  )
@@ -116,7 +116,7 @@ namespace
                                   // given N, generate i=0...N-1
                                   // equidistant points in the domain
                                   // [0,1]^3
-  Point<3>
+  inline Point<3>
   generate_unit_point (const unsigned int i,
                       const unsigned int N,
                       const int2type<3>  )
@@ -144,10 +144,10 @@ namespace
 template <int dim>
 FE_DGQ<dim>::FE_DGQ (const unsigned int degree)
                :
-               FiniteElement<dim> (FiniteElementData<dim>(get_dpo_vector(degree),1),
-                                   std::vector<bool>(FiniteElementData<dim>(get_dpo_vector(degree),1).dofs_per_cell,
+               FiniteElement<dim> (FiniteElementData<dim>(get_dpo_vector(degree), 1, degree),
+                                   std::vector<bool>(FiniteElementData<dim>(get_dpo_vector(degree),1, degree).dofs_per_cell,
                                                       true),
-                                   std::vector<std::vector<bool> >(FiniteElementData<dim>(get_dpo_vector(degree),1).dofs_per_cell,
+                                   std::vector<std::vector<bool> >(FiniteElementData<dim>(get_dpo_vector(degree),1, degree).dofs_per_cell,
                                                                    std::vector<bool>(1,true))),
                 degree(degree),
                 polynomial_space (Polynomials::LagrangeEquidistant::generate_complete_basis(degree))
index db4bfa1b1c58015344ba34db4586bbdd592564d1..53e7be3370854b3bc20f4ce0e98a8369a95a35b6 100644 (file)
@@ -431,10 +431,10 @@ namespace FE_Q_Helper
 template <int dim>
 FE_Q<dim>::FE_Q (const unsigned int degree)
                :
-               FiniteElement<dim> (FiniteElementData<dim>(get_dpo_vector(degree),1),
-                                   std::vector<bool> (FiniteElementData<dim>(get_dpo_vector(degree),1).dofs_per_cell,
+               FiniteElement<dim> (FiniteElementData<dim>(get_dpo_vector(degree),1, degree),
+                                   std::vector<bool> (FiniteElementData<dim>(get_dpo_vector(degree),1, degree).dofs_per_cell,
                                                        false),
-                                   std::vector<std::vector<bool> >(FiniteElementData<dim>(get_dpo_vector(degree),1).dofs_per_cell,
+                                   std::vector<std::vector<bool> >(FiniteElementData<dim>(get_dpo_vector(degree),1, degree).dofs_per_cell,
                                                                    std::vector<bool>(1,true))),
                 degree(degree),
                 renumber(lexicographic_to_hierarchic_numbering (*this, degree)),
index 7a3413c53ac75aff4a6ae8ef83ecd591bfa4ad3c..edf384f4091e782b0fe5e31c2637165d7f0c2c8f 100644 (file)
@@ -47,10 +47,10 @@ namespace
 template <int dim>
 FE_Q_Hierarchical<dim>::FE_Q_Hierarchical (const unsigned int degree)
                :
-               FiniteElement<dim> (FiniteElementData<dim>(get_dpo_vector(degree),1),
-                                   std::vector<bool> (FiniteElementData<dim>(get_dpo_vector(degree),1).dofs_per_cell,
+               FiniteElement<dim> (FiniteElementData<dim>(get_dpo_vector(degree),1, degree),
+                                   std::vector<bool> (FiniteElementData<dim>(get_dpo_vector(degree),1, degree).dofs_per_cell,
                                                        false),
-                                   std::vector<std::vector<bool> >(FiniteElementData<dim>(get_dpo_vector(degree),1).dofs_per_cell,
+                                   std::vector<std::vector<bool> >(FiniteElementData<dim>(get_dpo_vector(degree),1, degree).dofs_per_cell,
                                                                    std::vector<bool>(1,true))),
                degree(degree),
                renumber(lexicographic_to_hierarchic_numbering (*this, degree)),
index b63f5b854b2d40922a2097157035fdae2a115dd7..ddcc7c65d93de1484a9be4d0cd3e4ed6c5527da2 100644 (file)
@@ -128,9 +128,9 @@ template <int dim>
 FE_RaviartThomas<dim>::FE_RaviartThomas (const unsigned int degree)
                :
                FiniteElement<dim> (FiniteElementData<dim>(get_dpo_vector(degree),
-                                                          dim),
+                                                          dim, degree+1),
                                    get_ria_vector (degree),
-                                   std::vector<std::vector<bool> >(FiniteElementData<dim>(get_dpo_vector(degree),dim).dofs_per_cell,
+                                   std::vector<std::vector<bool> >(FiniteElementData<dim>(get_dpo_vector(degree),dim,degree+1).dofs_per_cell,
                                                                    std::vector<bool>(dim,true))),
                degree(degree),
                 polynomials (create_polynomials(degree)),
index 6aac1b39ff3c5b83c871b71580b64385d0ea2cac..9ae5a5cbb06d494cc89b58c9487e47b1a7baf4c0 100644 (file)
@@ -20,7 +20,7 @@
 #include <fe/mapping.h>
 #include <fe/fe_system.h>
 #include <fe/fe_values.h>
-
+#include <iostream>
 #ifdef HAVE_STD_STRINGSTREAM
 #  include <sstream>
 #else
@@ -1732,7 +1732,7 @@ FESystem<dim>::multiply_dof_numbers (const FiniteElementData<dim> &fe_data,
   if (dim>1) dpo.push_back(fe_data.dofs_per_quad * N);
   if (dim>2) dpo.push_back(fe_data.dofs_per_hex * N);
   
-  return FiniteElementData<dim> (dpo, fe_data.n_components() * N);
+  return FiniteElementData<dim> (dpo, fe_data.n_components() * N, fe_data.tensor_degree());
 }
 
 
@@ -1843,10 +1843,17 @@ FESystem<dim>::multiply_dof_numbers (const FiniteElementData<dim> &fe1,
   dpo.push_back(fe1.dofs_per_line * N1 + fe2.dofs_per_line * N2);
   if (dim>1) dpo.push_back(fe1.dofs_per_quad * N1 + fe2.dofs_per_quad * N2);
   if (dim>2) dpo.push_back(fe1.dofs_per_hex * N1 + fe2.dofs_per_hex * N2);
-  
+
+                                  // degree is the maximal degree of
+                                  // the components.  max also makes
+                                  // sure that one unknown degree
+                                  // makes the degree of the system
+                                  // unknown.
+  unsigned int degree = std::max(fe1.tensor_degree(), fe2.tensor_degree());
   return FiniteElementData<dim> (dpo,
                                 fe1.n_components() * N1 +
-                                fe2.n_components() * N2);
+                                fe2.n_components() * N2,
+                                degree);
 }
 
 
@@ -1873,11 +1880,13 @@ FESystem<dim>::multiply_dof_numbers (const FiniteElementData<dim> &fe1,
   if (dim>2) dpo.push_back(fe1.dofs_per_hex * N1 +
                           fe2.dofs_per_hex * N2 +
                           fe3.dofs_per_hex * N3);
-  
+                                  // degree is the maximal degree of the components
+  unsigned int degree = std::max(fe1.tensor_degree(), fe2.tensor_degree());
+  degree = std::max(degree, fe3.tensor_degree());
   return FiniteElementData<dim> (dpo,
                                 fe1.n_components() * N1 +
                                 fe2.n_components() * N2 +
-                                fe3.n_components() * N3);
+                                fe3.n_components() * N3, degree);
 }
 
 
index 7de7e4217578385caaf5af20c821a67ba66ef525..22a5c5f8b9a5baee1181e43a9fb29034e371cf6b 100644 (file)
@@ -303,10 +303,19 @@ contributor's names are abbreviated by WB (Wolfgang Bangerth), GK
 <h3>deal.II</h3>
 
 <ol>
+
+  <li> <p>
+  Improved: <code class="class">FiniteElementData</code> has a function
+  <class="member">tensor_degree()</code>, returning the degree of the
+  polynomial space suitable for choosing a tensor product quadrature formula.
+  <br>
+  (GK 2003/11/28)
+  </p>
+
   <li> <p>
        New: Long requested but never implemented before in the
        library: there is now a function <code
-       <class="member">GridTool::find_active_cell_around_point</code>
+       <class="class">GridTool</code>::<class="member">find_active_cell_around_point</code>
        that, given a point, finds the active cell in which this point
        lies.
        <br>
index 04240d56eaa21e9ae0052dfe2c82e6f937e87cf3..e542e044a0ebeb8bfa2d97a2c1dcabe6b3b32de2 100644 (file)
@@ -95,7 +95,8 @@ void test_fe_datas()
       deallog << "first_face_quad_index=" << fe_data->first_face_quad_index << std::endl;
       deallog << "dofs_per_face=" << fe_data->dofs_per_face << std::endl;
       deallog << "dofs_per_cell=" << fe_data->dofs_per_cell << std::endl;
-      deallog << "components=" << fe_data->components << std::endl;
+      deallog << "components=" << fe_data->components << std::endl
+             << "degree=" << fe_data->tensor_degree() << std::endl;
 
       for (unsigned int f=0;f<GeometryInfo<dim>::faces_per_cell;++f)
        {
@@ -105,6 +106,7 @@ void test_fe_datas()
              deallog << '\t' << s;
          deallog << std::endl;
        }
+      deallog << std::endl;
     }
 
                                   // delete all FiniteElementDatas
index f2ab5f76969d9a71ba059753bf0d511b7be6594d..c481055b3781653ec8d48f09566c93e17673a26a 100644 (file)
@@ -14,8 +14,10 @@ DEAL::first_face_quad_index=2
 DEAL::dofs_per_face=1
 DEAL::dofs_per_cell=2
 DEAL::components=1
+DEAL::degree=1
 DEAL::support on face 0:       0
 DEAL::support on face 1:       1
+DEAL::
 DEAL::fe_data[1]:FE_Q<1>(2)
 DEAL::dofs_per_vertex=1
 DEAL::dofs_per_line=1
@@ -29,8 +31,10 @@ DEAL::first_face_quad_index=2
 DEAL::dofs_per_face=1
 DEAL::dofs_per_cell=3
 DEAL::components=1
+DEAL::degree=2
 DEAL::support on face 0:       0
 DEAL::support on face 1:       1
+DEAL::
 DEAL::fe_data[2]:FE_Q<1>(4)
 DEAL::dofs_per_vertex=1
 DEAL::dofs_per_line=3
@@ -44,8 +48,10 @@ DEAL::first_face_quad_index=2
 DEAL::dofs_per_face=1
 DEAL::dofs_per_cell=5
 DEAL::components=1
+DEAL::degree=4
 DEAL::support on face 0:       0
 DEAL::support on face 1:       1
+DEAL::
 DEAL::fe_data[3]:FE_DGQ<1>(1)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=2
@@ -59,8 +65,10 @@ DEAL::first_face_quad_index=0
 DEAL::dofs_per_face=0
 DEAL::dofs_per_cell=2
 DEAL::components=1
+DEAL::degree=1
 DEAL::support on face 0:       0
 DEAL::support on face 1:       1
+DEAL::
 DEAL::fe_data[4]:FE_DGQ<1>(2)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=3
@@ -74,8 +82,10 @@ DEAL::first_face_quad_index=0
 DEAL::dofs_per_face=0
 DEAL::dofs_per_cell=3
 DEAL::components=1
+DEAL::degree=2
 DEAL::support on face 0:       0
 DEAL::support on face 1:       1
+DEAL::
 DEAL::fe_data[5]:FE_DGQ<1>(4)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=5
@@ -89,8 +99,10 @@ DEAL::first_face_quad_index=0
 DEAL::dofs_per_face=0
 DEAL::dofs_per_cell=5
 DEAL::components=1
+DEAL::degree=4
 DEAL::support on face 0:       0
 DEAL::support on face 1:       1
+DEAL::
 DEAL::fe_data[6]:FE_DGP<1>(1)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=2
@@ -104,8 +116,10 @@ DEAL::first_face_quad_index=0
 DEAL::dofs_per_face=0
 DEAL::dofs_per_cell=2
 DEAL::components=1
+DEAL::degree=1
 DEAL::support on face 0:       0       1
 DEAL::support on face 1:       0       1
+DEAL::
 DEAL::fe_data[7]:FE_DGP<1>(2)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=3
@@ -119,8 +133,10 @@ DEAL::first_face_quad_index=0
 DEAL::dofs_per_face=0
 DEAL::dofs_per_cell=3
 DEAL::components=1
+DEAL::degree=2
 DEAL::support on face 0:       0       1       2
 DEAL::support on face 1:       0       1       2
+DEAL::
 DEAL::fe_data[8]:FE_DGP<1>(4)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=5
@@ -134,8 +150,10 @@ DEAL::first_face_quad_index=0
 DEAL::dofs_per_face=0
 DEAL::dofs_per_cell=5
 DEAL::components=1
+DEAL::degree=4
 DEAL::support on face 0:       0       1       2       3       4
 DEAL::support on face 1:       0       1       2       3       4
+DEAL::
 DEAL::fe_data[9]:FESystem<1>[FE_Q<1>(2)^2]
 DEAL::dofs_per_vertex=2
 DEAL::dofs_per_line=2
@@ -149,8 +167,10 @@ DEAL::first_face_quad_index=4
 DEAL::dofs_per_face=2
 DEAL::dofs_per_cell=6
 DEAL::components=2
+DEAL::degree=2
 DEAL::support on face 0:       0       1
 DEAL::support on face 1:       2       3
+DEAL::
 DEAL::fe_data[10]:FESystem<1>[FE_Q<1>(1)^2-FE_Q<1>(2)]
 DEAL::dofs_per_vertex=3
 DEAL::dofs_per_line=1
@@ -164,8 +184,10 @@ DEAL::first_face_quad_index=6
 DEAL::dofs_per_face=3
 DEAL::dofs_per_cell=7
 DEAL::components=3
+DEAL::degree=2
 DEAL::support on face 0:       0       1       2
 DEAL::support on face 1:       3       4       5
+DEAL::
 DEAL::fe_data[11]:FESystem<1>[FE_Q<1>(3)^2]
 DEAL::dofs_per_vertex=2
 DEAL::dofs_per_line=4
@@ -179,8 +201,10 @@ DEAL::first_face_quad_index=4
 DEAL::dofs_per_face=2
 DEAL::dofs_per_cell=8
 DEAL::components=2
+DEAL::degree=3
 DEAL::support on face 0:       0       1
 DEAL::support on face 1:       2       3
+DEAL::
 DEAL::fe_data[12]:FESystem<1>[FE_Q<1>(1)^2-FE_Q<1>(3)]
 DEAL::dofs_per_vertex=3
 DEAL::dofs_per_line=2
@@ -194,8 +218,10 @@ DEAL::first_face_quad_index=6
 DEAL::dofs_per_face=3
 DEAL::dofs_per_cell=8
 DEAL::components=3
+DEAL::degree=3
 DEAL::support on face 0:       0       1       2
 DEAL::support on face 1:       3       4       5
+DEAL::
 DEAL::fe_data[13]:FESystem<1>[FE_Q<1>(4)^2]
 DEAL::dofs_per_vertex=2
 DEAL::dofs_per_line=6
@@ -209,8 +235,10 @@ DEAL::first_face_quad_index=4
 DEAL::dofs_per_face=2
 DEAL::dofs_per_cell=10
 DEAL::components=2
+DEAL::degree=4
 DEAL::support on face 0:       0       1
 DEAL::support on face 1:       2       3
+DEAL::
 DEAL::fe_data[14]:FESystem<1>[FESystem<1>[FE_Q<1>(1)^2]^2]
 DEAL::dofs_per_vertex=4
 DEAL::dofs_per_line=0
@@ -224,8 +252,10 @@ DEAL::first_face_quad_index=8
 DEAL::dofs_per_face=4
 DEAL::dofs_per_cell=8
 DEAL::components=4
+DEAL::degree=1
 DEAL::support on face 0:       0       1       2       3
 DEAL::support on face 1:       4       5       6       7
+DEAL::
 DEAL::fe_data[15]:FESystem<1>[FESystem<1>[FE_Q<1>(1)^2]-FESystem<1>[FE_DGQ<1>(1)^2]]
 DEAL::dofs_per_vertex=2
 DEAL::dofs_per_line=4
@@ -239,8 +269,10 @@ DEAL::first_face_quad_index=4
 DEAL::dofs_per_face=2
 DEAL::dofs_per_cell=8
 DEAL::components=4
+DEAL::degree=1
 DEAL::support on face 0:       0       1       4       6
 DEAL::support on face 1:       2       3       5       7
+DEAL::
 DEAL::fe_data[16]:FESystem<1>[FESystem<1>[FE_Q<1>(1)-FE_Q<1>(2)]-FESystem<1>[FE_Q<1>(2)^2]-FESystem<1>[FE_DGQ<1>(2)^2]]
 DEAL::dofs_per_vertex=4
 DEAL::dofs_per_line=9
@@ -254,8 +286,10 @@ DEAL::first_face_quad_index=8
 DEAL::dofs_per_face=4
 DEAL::dofs_per_cell=17
 DEAL::components=6
+DEAL::degree=2
 DEAL::support on face 0:       0       1       2       3       11      14
 DEAL::support on face 1:       4       5       6       7       12      15
+DEAL::
 DEAL::fe_data[17]:FESystem<1>[FESystem<1>[FESystem<1>[FE_Q<1>(1)^2]^2]^2-FESystem<1>[FESystem<1>[FE_Q<1>(1)^2]-FESystem<1>[FE_DGQ<1>(1)^2]]-FESystem<1>[FESystem<1>[FE_Q<1>(1)-FE_Q<1>(2)]-FESystem<1>[FE_Q<1>(2)^2]-FESystem<1>[FE_DGQ<1>(2)^2]]^2]
 DEAL::dofs_per_vertex=18
 DEAL::dofs_per_line=22
@@ -269,9 +303,11 @@ DEAL::first_face_quad_index=36
 DEAL::dofs_per_face=18
 DEAL::dofs_per_cell=58
 DEAL::components=24
+DEAL::degree=2
 DEAL::support on face 0:       0       1       2       3       4       5       6       7       8       9       10      11      12      13      14      15      16      17      36      38      43      46      52      55
 DEAL::support on face 1:       18      19      20      21      22      23      24      25      26      27      28      29      30      31      32      33      34      35      37      39      44      47      53      56
 DEAL::
+DEAL::
 DEAL::dim=2
 DEAL::fe_data[0]:FE_Q<2>(1)
 DEAL::dofs_per_vertex=1
@@ -286,10 +322,12 @@ DEAL::first_face_quad_index=4
 DEAL::dofs_per_face=2
 DEAL::dofs_per_cell=4
 DEAL::components=1
+DEAL::degree=1
 DEAL::support on face 0:       0       1
 DEAL::support on face 1:       1       2
 DEAL::support on face 2:       2       3
 DEAL::support on face 3:       0       3
+DEAL::
 DEAL::fe_data[1]:FE_Q<2>(2)
 DEAL::dofs_per_vertex=1
 DEAL::dofs_per_line=1
@@ -303,10 +341,12 @@ DEAL::first_face_quad_index=5
 DEAL::dofs_per_face=3
 DEAL::dofs_per_cell=9
 DEAL::components=1
+DEAL::degree=2
 DEAL::support on face 0:       0       1       4
 DEAL::support on face 1:       1       2       5
 DEAL::support on face 2:       2       3       6
 DEAL::support on face 3:       0       3       7
+DEAL::
 DEAL::fe_data[2]:FE_Q<2>(4)
 DEAL::dofs_per_vertex=1
 DEAL::dofs_per_line=3
@@ -320,10 +360,12 @@ DEAL::first_face_quad_index=7
 DEAL::dofs_per_face=5
 DEAL::dofs_per_cell=25
 DEAL::components=1
+DEAL::degree=4
 DEAL::support on face 0:       0       1       4       5       6
 DEAL::support on face 1:       1       2       7       8       9
 DEAL::support on face 2:       2       3       10      11      12
 DEAL::support on face 3:       0       3       13      14      15
+DEAL::
 DEAL::fe_data[3]:FE_DGQ<2>(1)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
@@ -337,10 +379,12 @@ DEAL::first_face_quad_index=0
 DEAL::dofs_per_face=0
 DEAL::dofs_per_cell=4
 DEAL::components=1
+DEAL::degree=1
 DEAL::support on face 0:       0       1
 DEAL::support on face 1:       1       3
 DEAL::support on face 2:       2       3
 DEAL::support on face 3:       0       2
+DEAL::
 DEAL::fe_data[4]:FE_DGQ<2>(2)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
@@ -354,10 +398,12 @@ DEAL::first_face_quad_index=0
 DEAL::dofs_per_face=0
 DEAL::dofs_per_cell=9
 DEAL::components=1
+DEAL::degree=2
 DEAL::support on face 0:       0       1       2
 DEAL::support on face 1:       2       5       8
 DEAL::support on face 2:       6       7       8
 DEAL::support on face 3:       0       3       6
+DEAL::
 DEAL::fe_data[5]:FE_DGQ<2>(4)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
@@ -371,10 +417,12 @@ DEAL::first_face_quad_index=0
 DEAL::dofs_per_face=0
 DEAL::dofs_per_cell=25
 DEAL::components=1
+DEAL::degree=4
 DEAL::support on face 0:       0       1       2       3       4
 DEAL::support on face 1:       4       9       14      19      24
 DEAL::support on face 2:       20      21      22      23      24
 DEAL::support on face 3:       0       5       10      15      20
+DEAL::
 DEAL::fe_data[6]:FE_DGP<2>(1)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
@@ -388,10 +436,12 @@ DEAL::first_face_quad_index=0
 DEAL::dofs_per_face=0
 DEAL::dofs_per_cell=3
 DEAL::components=1
+DEAL::degree=1
 DEAL::support on face 0:       0       1       2
 DEAL::support on face 1:       0       1       2
 DEAL::support on face 2:       0       1       2
 DEAL::support on face 3:       0       1       2
+DEAL::
 DEAL::fe_data[7]:FE_DGP<2>(2)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
@@ -405,10 +455,12 @@ DEAL::first_face_quad_index=0
 DEAL::dofs_per_face=0
 DEAL::dofs_per_cell=6
 DEAL::components=1
+DEAL::degree=2
 DEAL::support on face 0:       0       1       2       3       4       5
 DEAL::support on face 1:       0       1       2       3       4       5
 DEAL::support on face 2:       0       1       2       3       4       5
 DEAL::support on face 3:       0       1       2       3       4       5
+DEAL::
 DEAL::fe_data[8]:FE_DGP<2>(4)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
@@ -422,10 +474,12 @@ DEAL::first_face_quad_index=0
 DEAL::dofs_per_face=0
 DEAL::dofs_per_cell=15
 DEAL::components=1
+DEAL::degree=4
 DEAL::support on face 0:       0       1       2       3       4       5       6       7       8       9       10      11      12      13      14
 DEAL::support on face 1:       0       1       2       3       4       5       6       7       8       9       10      11      12      13      14
 DEAL::support on face 2:       0       1       2       3       4       5       6       7       8       9       10      11      12      13      14
 DEAL::support on face 3:       0       1       2       3       4       5       6       7       8       9       10      11      12      13      14
+DEAL::
 DEAL::fe_data[9]:FESystem<2>[FE_Q<2>(2)^2]
 DEAL::dofs_per_vertex=2
 DEAL::dofs_per_line=2
@@ -439,10 +493,12 @@ DEAL::first_face_quad_index=10
 DEAL::dofs_per_face=6
 DEAL::dofs_per_cell=18
 DEAL::components=2
+DEAL::degree=2
 DEAL::support on face 0:       0       1       2       3       8       9
 DEAL::support on face 1:       2       3       4       5       10      11
 DEAL::support on face 2:       4       5       6       7       12      13
 DEAL::support on face 3:       0       1       6       7       14      15
+DEAL::
 DEAL::fe_data[10]:FESystem<2>[FE_Q<2>(1)^2-FE_Q<2>(2)]
 DEAL::dofs_per_vertex=3
 DEAL::dofs_per_line=1
@@ -456,10 +512,12 @@ DEAL::first_face_quad_index=13
 DEAL::dofs_per_face=7
 DEAL::dofs_per_cell=17
 DEAL::components=3
+DEAL::degree=2
 DEAL::support on face 0:       0       1       2       3       4       5       12
 DEAL::support on face 1:       3       4       5       6       7       8       13
 DEAL::support on face 2:       6       7       8       9       10      11      14
 DEAL::support on face 3:       0       1       2       9       10      11      15
+DEAL::
 DEAL::fe_data[11]:FE_RaviartThomas<2>(0)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=1
@@ -473,10 +531,12 @@ DEAL::first_face_quad_index=1
 DEAL::dofs_per_face=1
 DEAL::dofs_per_cell=4
 DEAL::components=2
+DEAL::degree=1
 DEAL::support on face 0:       0       1       3
 DEAL::support on face 1:       0       1       2
 DEAL::support on face 2:       1       2       3
 DEAL::support on face 3:       0       2       3
+DEAL::
 DEAL::fe_data[12]:FE_RaviartThomas<2>(1)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=2
@@ -490,10 +550,12 @@ DEAL::first_face_quad_index=2
 DEAL::dofs_per_face=2
 DEAL::dofs_per_cell=12
 DEAL::components=2
+DEAL::degree=2
 DEAL::support on face 0:       0       1       2       3       4       5       6       7       8       9       10      11
 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[13]:FE_RaviartThomas<2>(2)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=3
@@ -507,10 +569,12 @@ DEAL::first_face_quad_index=3
 DEAL::dofs_per_face=3
 DEAL::dofs_per_cell=24
 DEAL::components=2
+DEAL::degree=3
 DEAL::support on face 0:       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 1:       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 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[14]:FESystem<2>[FE_Q<2>(3)^2]
 DEAL::dofs_per_vertex=2
 DEAL::dofs_per_line=4
@@ -524,10 +588,12 @@ DEAL::first_face_quad_index=12
 DEAL::dofs_per_face=8
 DEAL::dofs_per_cell=32
 DEAL::components=2
+DEAL::degree=3
 DEAL::support on face 0:       0       1       2       3       8       9       10      11
 DEAL::support on face 1:       2       3       4       5       12      13      14      15
 DEAL::support on face 2:       4       5       6       7       16      17      18      19
 DEAL::support on face 3:       0       1       6       7       20      21      22      23
+DEAL::
 DEAL::fe_data[15]:FESystem<2>[FE_Q<2>(1)^2-FE_Q<2>(3)]
 DEAL::dofs_per_vertex=3
 DEAL::dofs_per_line=2
@@ -541,10 +607,12 @@ DEAL::first_face_quad_index=14
 DEAL::dofs_per_face=8
 DEAL::dofs_per_cell=24
 DEAL::components=3
+DEAL::degree=3
 DEAL::support on face 0:       0       1       2       3       4       5       12      13
 DEAL::support on face 1:       3       4       5       6       7       8       14      15
 DEAL::support on face 2:       6       7       8       9       10      11      16      17
 DEAL::support on face 3:       0       1       2       9       10      11      18      19
+DEAL::
 DEAL::fe_data[16]:FESystem<2>[FE_Q<2>(4)^2]
 DEAL::dofs_per_vertex=2
 DEAL::dofs_per_line=6
@@ -558,10 +626,12 @@ DEAL::first_face_quad_index=14
 DEAL::dofs_per_face=10
 DEAL::dofs_per_cell=50
 DEAL::components=2
+DEAL::degree=4
 DEAL::support on face 0:       0       1       2       3       8       9       10      11      12      13
 DEAL::support on face 1:       2       3       4       5       14      15      16      17      18      19
 DEAL::support on face 2:       4       5       6       7       20      21      22      23      24      25
 DEAL::support on face 3:       0       1       6       7       26      27      28      29      30      31
+DEAL::
 DEAL::fe_data[17]:FESystem<2>[FESystem<2>[FE_Q<2>(1)^2]^2]
 DEAL::dofs_per_vertex=4
 DEAL::dofs_per_line=0
@@ -575,10 +645,12 @@ DEAL::first_face_quad_index=16
 DEAL::dofs_per_face=8
 DEAL::dofs_per_cell=16
 DEAL::components=4
+DEAL::degree=1
 DEAL::support on face 0:       0       1       2       3       4       5       6       7
 DEAL::support on face 1:       4       5       6       7       8       9       10      11
 DEAL::support on face 2:       8       9       10      11      12      13      14      15
 DEAL::support on face 3:       0       1       2       3       12      13      14      15
+DEAL::
 DEAL::fe_data[18]: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
@@ -592,10 +664,12 @@ DEAL::first_face_quad_index=8
 DEAL::dofs_per_face=4
 DEAL::dofs_per_cell=16
 DEAL::components=4
+DEAL::degree=1
 DEAL::support on face 0:       0       1       2       3       8       9       12      13
 DEAL::support on face 1:       2       3       4       5       9       11      13      15
 DEAL::support on face 2:       4       5       6       7       10      11      14      15
 DEAL::support on face 3:       0       1       6       7       8       10      12      14
+DEAL::
 DEAL::fe_data[19]: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
@@ -609,10 +683,12 @@ DEAL::first_face_quad_index=19
 DEAL::dofs_per_face=11
 DEAL::dofs_per_cell=49
 DEAL::components=6
+DEAL::degree=2
 DEAL::support on face 0:       0       1       2       3       4       5       6       7       16      17      18      31      32      33      40      41      42
 DEAL::support on face 1:       4       5       6       7       8       9       10      11      19      20      21      33      36      39      42      45      48
 DEAL::support on face 2:       8       9       10      11      12      13      14      15      22      23      24      37      38      39      46      47      48
 DEAL::support on face 3:       0       1       2       3       12      13      14      15      25      26      27      31      34      37      40      43      46
+DEAL::
 DEAL::fe_data[20]: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
@@ -626,11 +702,13 @@ DEAL::first_face_quad_index=78
 DEAL::dofs_per_face=42
 DEAL::dofs_per_cell=146
 DEAL::components=24
+DEAL::degree=2
 DEAL::support on face 0:       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      72      73      74      75      76      77      96      97      100     101     107     108     109     116     117     118     128     129     130     137     138     139
 DEAL::support on face 1:       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      78      79      80      81      82      83      97      99      101     103     109     112     115     118     121     124     130     133     136     139     142     145
 DEAL::support on face 2:       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      84      85      86      87      88      89      98      99      102     103     113     114     115     122     123     124     134     135     136     143     144     145
 DEAL::support on face 3:       0       1       2       3       4       5       6       7       8       9       10      11      12      13      14      15      16      17      54      55      56      57      58      59      60      61      62      63      64      65      66      67      68      69      70      71      90      91      92      93      94      95      96      98      100     102     107     110     113     116     119     122     128     131     134     137     140     143
 DEAL::
+DEAL::
 DEAL::dim=3
 DEAL::fe_data[0]:FE_Q<3>(1)
 DEAL::dofs_per_vertex=1
@@ -645,12 +723,14 @@ DEAL::first_face_quad_index=4
 DEAL::dofs_per_face=4
 DEAL::dofs_per_cell=8
 DEAL::components=1
+DEAL::degree=1
 DEAL::support on face 0:       0       1       2       3
 DEAL::support on face 1:       4       5       6       7
 DEAL::support on face 2:       0       1       4       5
 DEAL::support on face 3:       1       2       5       6
 DEAL::support on face 4:       2       3       6       7
 DEAL::support on face 5:       0       3       4       7
+DEAL::
 DEAL::fe_data[1]:FE_Q<3>(2)
 DEAL::dofs_per_vertex=1
 DEAL::dofs_per_line=1
@@ -664,12 +744,14 @@ DEAL::first_face_quad_index=8
 DEAL::dofs_per_face=9
 DEAL::dofs_per_cell=27
 DEAL::components=1
+DEAL::degree=2
 DEAL::support on face 0:       0       1       2       3       8       9       10      11      20
 DEAL::support on face 1:       4       5       6       7       12      13      14      15      21
 DEAL::support on face 2:       0       1       4       5       8       12      16      17      22
 DEAL::support on face 3:       1       2       5       6       9       13      17      18      23
 DEAL::support on face 4:       2       3       6       7       10      14      18      19      24
 DEAL::support on face 5:       0       3       4       7       11      15      16      19      25
+DEAL::
 DEAL::fe_data[2]:FE_Q<3>(4)
 DEAL::dofs_per_vertex=1
 DEAL::dofs_per_line=3
@@ -683,12 +765,14 @@ DEAL::first_face_quad_index=16
 DEAL::dofs_per_face=25
 DEAL::dofs_per_cell=125
 DEAL::components=1
+DEAL::degree=4
 DEAL::support on face 0:       0       1       2       3       8       9       10      11      12      13      14      15      16      17      18      19      44      45      46      47      48      49      50      51      52
 DEAL::support on face 1:       4       5       6       7       20      21      22      23      24      25      26      27      28      29      30      31      53      54      55      56      57      58      59      60      61
 DEAL::support on face 2:       0       1       4       5       8       9       10      20      21      22      32      33      34      35      36      37      62      63      64      65      66      67      68      69      70
 DEAL::support on face 3:       1       2       5       6       11      12      13      23      24      25      35      36      37      38      39      40      71      72      73      74      75      76      77      78      79
 DEAL::support on face 4:       2       3       6       7       14      15      16      26      27      28      38      39      40      41      42      43      80      81      82      83      84      85      86      87      88
 DEAL::support on face 5:       0       3       4       7       17      18      19      29      30      31      32      33      34      41      42      43      89      90      91      92      93      94      95      96      97
+DEAL::
 DEAL::fe_data[3]:FE_DGQ<3>(1)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
@@ -702,12 +786,14 @@ DEAL::first_face_quad_index=0
 DEAL::dofs_per_face=0
 DEAL::dofs_per_cell=8
 DEAL::components=1
+DEAL::degree=1
 DEAL::support on face 0:       0       1       4       5
 DEAL::support on face 1:       2       3       6       7
 DEAL::support on face 2:       0       1       2       3
 DEAL::support on face 3:       1       3       5       7
 DEAL::support on face 4:       4       5       6       7
 DEAL::support on face 5:       0       2       4       6
+DEAL::
 DEAL::fe_data[4]:FE_DGQ<3>(2)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
@@ -721,12 +807,14 @@ DEAL::first_face_quad_index=0
 DEAL::dofs_per_face=0
 DEAL::dofs_per_cell=27
 DEAL::components=1
+DEAL::degree=2
 DEAL::support on face 0:       0       1       2       9       10      11      18      19      20
 DEAL::support on face 1:       6       7       8       15      16      17      24      25      26
 DEAL::support on face 2:       0       1       2       3       4       5       6       7       8
 DEAL::support on face 3:       2       5       8       11      14      17      20      23      26
 DEAL::support on face 4:       18      19      20      21      22      23      24      25      26
 DEAL::support on face 5:       0       3       6       9       12      15      18      21      24
+DEAL::
 DEAL::fe_data[5]:FE_DGQ<3>(4)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
@@ -740,12 +828,14 @@ DEAL::first_face_quad_index=0
 DEAL::dofs_per_face=0
 DEAL::dofs_per_cell=125
 DEAL::components=1
+DEAL::degree=4
 DEAL::support on face 0:       0       1       2       3       4       25      26      27      28      29      50      51      52      53      54      75      76      77      78      79      100     101     102     103     104
 DEAL::support on face 1:       20      21      22      23      24      45      46      47      48      49      70      71      72      73      74      95      96      97      98      99      120     121     122     123     124
 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      24
 DEAL::support on face 3:       4       9       14      19      24      29      34      39      44      49      54      59      64      69      74      79      84      89      94      99      104     109     114     119     124
 DEAL::support on face 4:       100     101     102     103     104     105     106     107     108     109     110     111     112     113     114     115     116     117     118     119     120     121     122     123     124
 DEAL::support on face 5:       0       5       10      15      20      25      30      35      40      45      50      55      60      65      70      75      80      85      90      95      100     105     110     115     120
+DEAL::
 DEAL::fe_data[6]:FE_DGP<3>(1)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
@@ -759,12 +849,14 @@ DEAL::first_face_quad_index=0
 DEAL::dofs_per_face=0
 DEAL::dofs_per_cell=4
 DEAL::components=1
+DEAL::degree=1
 DEAL::support on face 0:       0       1       2       3
 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::support on face 4:       0       1       2       3
 DEAL::support on face 5:       0       1       2       3
+DEAL::
 DEAL::fe_data[7]:FE_DGP<3>(2)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
@@ -778,12 +870,14 @@ DEAL::first_face_quad_index=0
 DEAL::dofs_per_face=0
 DEAL::dofs_per_cell=10
 DEAL::components=1
+DEAL::degree=2
 DEAL::support on face 0:       0       1       2       3       4       5       6       7       8       9
 DEAL::support on face 1:       0       1       2       3       4       5       6       7       8       9
 DEAL::support on face 2:       0       1       2       3       4       5       6       7       8       9
 DEAL::support on face 3:       0       1       2       3       4       5       6       7       8       9
 DEAL::support on face 4:       0       1       2       3       4       5       6       7       8       9
 DEAL::support on face 5:       0       1       2       3       4       5       6       7       8       9
+DEAL::
 DEAL::fe_data[8]:FE_DGP<3>(4)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
@@ -797,12 +891,14 @@ DEAL::first_face_quad_index=0
 DEAL::dofs_per_face=0
 DEAL::dofs_per_cell=35
 DEAL::components=1
+DEAL::degree=4
 DEAL::support on face 0:       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
 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      22      23      24      25      26      27      28      29      30      31      32      33      34
 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      24      25      26      27      28      29      30      31      32      33      34
 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      24      25      26      27      28      29      30      31      32      33      34
 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
 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
+DEAL::
 DEAL::fe_data[9]:FESystem<3>[FE_Q<3>(2)^2]
 DEAL::dofs_per_vertex=2
 DEAL::dofs_per_line=2
@@ -816,12 +912,14 @@ DEAL::first_face_quad_index=16
 DEAL::dofs_per_face=18
 DEAL::dofs_per_cell=54
 DEAL::components=2
+DEAL::degree=2
 DEAL::support on face 0:       0       1       2       3       4       5       6       7       16      17      18      19      20      21      22      23      40      41
 DEAL::support on face 1:       8       9       10      11      12      13      14      15      24      25      26      27      28      29      30      31      42      43
 DEAL::support on face 2:       0       1       2       3       8       9       10      11      16      17      24      25      32      33      34      35      44      45
 DEAL::support on face 3:       2       3       4       5       10      11      12      13      18      19      26      27      34      35      36      37      46      47
 DEAL::support on face 4:       4       5       6       7       12      13      14      15      20      21      28      29      36      37      38      39      48      49
 DEAL::support on face 5:       0       1       6       7       8       9       14      15      22      23      30      31      32      33      38      39      50      51
+DEAL::
 DEAL::fe_data[10]:FESystem<3>[FE_Q<3>(1)^2-FE_Q<3>(2)]
 DEAL::dofs_per_vertex=3
 DEAL::dofs_per_line=1
@@ -835,12 +933,14 @@ DEAL::first_face_quad_index=16
 DEAL::dofs_per_face=17
 DEAL::dofs_per_cell=43
 DEAL::components=3
+DEAL::degree=2
 DEAL::support on face 0:       0       1       2       3       4       5       6       7       8       9       10      11      24      25      26      27      36
 DEAL::support on face 1:       12      13      14      15      16      17      18      19      20      21      22      23      28      29      30      31      37
 DEAL::support on face 2:       0       1       2       3       4       5       12      13      14      15      16      17      24      28      32      33      38
 DEAL::support on face 3:       3       4       5       6       7       8       15      16      17      18      19      20      25      29      33      34      39
 DEAL::support on face 4:       6       7       8       9       10      11      18      19      20      21      22      23      26      30      34      35      40
 DEAL::support on face 5:       0       1       2       9       10      11      12      13      14      21      22      23      27      31      32      35      41
+DEAL::
 DEAL::fe_data[11]:FESystem<3>[FESystem<3>[FE_Q<3>(1)^2]^2]
 DEAL::dofs_per_vertex=4
 DEAL::dofs_per_line=0
@@ -854,12 +954,14 @@ DEAL::first_face_quad_index=16
 DEAL::dofs_per_face=16
 DEAL::dofs_per_cell=32
 DEAL::components=4
+DEAL::degree=1
 DEAL::support on face 0:       0       1       2       3       4       5       6       7       8       9       10      11      12      13      14      15
 DEAL::support on face 1:       16      17      18      19      20      21      22      23      24      25      26      27      28      29      30      31
 DEAL::support on face 2:       0       1       2       3       4       5       6       7       16      17      18      19      20      21      22      23
 DEAL::support on face 3:       4       5       6       7       8       9       10      11      20      21      22      23      24      25      26      27
 DEAL::support on face 4:       8       9       10      11      12      13      14      15      24      25      26      27      28      29      30      31
 DEAL::support on face 5:       0       1       2       3       12      13      14      15      16      17      18      19      28      29      30      31
+DEAL::
 DEAL::fe_data[12]: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
@@ -873,12 +975,14 @@ DEAL::first_face_quad_index=8
 DEAL::dofs_per_face=8
 DEAL::dofs_per_cell=32
 DEAL::components=4
+DEAL::degree=1
 DEAL::support on face 0:       0       1       2       3       4       5       6       7       16      17      20      21      24      25      28      29
 DEAL::support on face 1:       8       9       10      11      12      13      14      15      18      19      22      23      26      27      30      31
 DEAL::support on face 2:       0       1       2       3       8       9       10      11      16      17      18      19      24      25      26      27
 DEAL::support on face 3:       2       3       4       5       10      11      12      13      17      19      21      23      25      27      29      31
 DEAL::support on face 4:       4       5       6       7       12      13      14      15      20      21      22      23      28      29      30      31
 DEAL::support on face 5:       0       1       6       7       8       9       14      15      16      18      20      22      24      26      28      30
+DEAL::
 DEAL::fe_data[13]: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
@@ -892,12 +996,14 @@ DEAL::first_face_quad_index=28
 DEAL::dofs_per_face=31
 DEAL::dofs_per_cell=143
 DEAL::components=6
+DEAL::degree=2
 DEAL::support on face 0:       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      68      69      70      89      90      91      98      99      100     107     108     109     116     117     118     125     126     127     134     135     136
 DEAL::support on face 1:       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      71      72      73      95      96      97      104     105     106     113     114     115     122     123     124     131     132     133     140     141     142
 DEAL::support on face 2:       0       1       2       3       4       5       6       7       16      17      18      19      20      21      22      23      32      33      34      44      45      46      56      57      58      59      60      61      74      75      76      89      90      91      92      93      94      95      96      97      116     117     118     119     120     121     122     123     124
 DEAL::support on face 3:       4       5       6       7       8       9       10      11      20      21      22      23      24      25      26      27      35      36      37      47      48      49      59      60      61      62      63      64      77      78      79      91      94      97      100     103     106     109     112     115     118     121     124     127     130     133     136     139     142
 DEAL::support on face 4:       8       9       10      11      12      13      14      15      24      25      26      27      28      29      30      31      38      39      40      50      51      52      62      63      64      65      66      67      80      81      82      107     108     109     110     111     112     113     114     115     134     135     136     137     138     139     140     141     142
 DEAL::support on face 5:       0       1       2       3       12      13      14      15      16      17      18      19      28      29      30      31      41      42      43      53      54      55      56      57      58      65      66      67      83      84      85      89      92      95      98      101     104     107     110     113     116     119     122     125     128     131     134     137     140
+DEAL::
 DEAL::fe_data[14]: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
@@ -911,9 +1017,11 @@ DEAL::first_face_quad_index=96
 DEAL::dofs_per_face=102
 DEAL::dofs_per_cell=382
 DEAL::components=24
+DEAL::degree=2
 DEAL::support on face 0:       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      144     145     146     147     148     149     150     151     152     153     154     155     156     157     158     159     160     161     162     163     164     165     166     167     216     217     218     219     220     221     252     253     256     257     260     261     264     265     271     272     273     280     281     282     289     290     291     298     299     300     307     308     309     316     317     318     328     329     330     337     338     339     346     347     348     355     356     357     364     365     366     373     374     375
 DEAL::support on face 1:       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     108     109     110     111     112     113     114     115     116     117     118     119     120     121     122     123     124     125     126     127     128     129     130     131     132     133     134     135     136     137     138     139     140     141     142     143     168     169     170     171     172     173     174     175     176     177     178     179     180     181     182     183     184     185     186     187     188     189     190     191     222     223     224     225     226     227     254     255     258     259     262     263     266     267     277     278     279     286     287     288     295     296     297     304     305     306     313     314     315     322     323     324     334     335     336     343     344     345     352     353     354     361     362     363     370     371     372     379     380     381
 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      24      25      26      27      28      29      30      31      32      33      34      35      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     144     145     146     147     148     149     168     169     170     171     172     173     192     193     194     195     196     197     198     199     200     201     202     203     228     229     230     231     232     233     252     253     254     255     260     261     262     263     271     272     273     274     275     276     277     278     279     298     299     300     301     302     303     304     305     306     328     329     330     331     332     333     334     335     336     355     356     357     358     359     360     361     362     363
 DEAL::support on face 3:       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      90      91      92      93      94      95      96      97      98      99      100     101     102     103     104     105     106     107     108     109     110     111     112     113     114     115     116     117     118     119     120     121     122     123     124     125     150     151     152     153     154     155     174     175     176     177     178     179     198     199     200     201     202     203     204     205     206     207     208     209     234     235     236     237     238     239     253     255     257     259     261     263     265     267     273     276     279     282     285     288     291     294     297     300     303     306     309     312     315     318     321     324     330     333     336     339     342     345     348     351     354     357     360     363     366     369     372     375     378     381
 DEAL::support on face 4:       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      108     109     110     111     112     113     114     115     116     117     118     119     120     121     122     123     124     125     126     127     128     129     130     131     132     133     134     135     136     137     138     139     140     141     142     143     156     157     158     159     160     161     180     181     182     183     184     185     204     205     206     207     208     209     210     211     212     213     214     215     240     241     242     243     244     245     256     257     258     259     264     265     266     267     289     290     291     292     293     294     295     296     297     316     317     318     319     320     321     322     323     324     346     347     348     349     350     351     352     353     354     373     374     375     376     377     378     379     380     381
 DEAL::support on face 5:       0       1       2       3       4       5       6       7       8       9       10      11      12      13      14      15      16      17      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      126     127     128     129     130     131     132     133     134     135     136     137     138     139     140     141     142     143     162     163     164     165     166     167     186     187     188     189     190     191     192     193     194     195     196     197     210     211     212     213     214     215     246     247     248     249     250     251     252     254     256     258     260     262     264     266     271     274     277     280     283     286     289     292     295     298     301     304     307     310     313     316     319     322     328     331     334     337     340     343     346     349     352     355     358     361     364     367     370     373     376     379
+DEAL::

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.