]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Fix FE_DGQArbitraryNodes::has_support_on_face
authorMartin Kronbichler <kronbichler@lnm.mw.tum.de>
Wed, 30 Apr 2014 09:15:00 +0000 (09:15 +0000)
committerMartin Kronbichler <kronbichler@lnm.mw.tum.de>
Wed, 30 Apr 2014 09:15:00 +0000 (09:15 +0000)
git-svn-id: https://svn.dealii.org/trunk@32863 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/doc/news/changes.h
deal.II/source/fe/fe_dgq.cc
tests/fe/fe_data_test.cc
tests/fe/fe_data_test.output

index c442348d1bbe73d8e99fe497f33e687212a1d7b3..d4be862392491b7bae1adc2eea3547bf993f0f78 100644 (file)
@@ -148,6 +148,13 @@ inconvenience this causes.
 <h3>Specific improvements</h3>
 
 <ol>
+  <li> Fixed: FE_DGQArbitraryNodes::has_support_on_face was broken when
+  polynomials with support points not on the element boundary were used.
+  This is now fixed.
+  <br>
+  (Martin Kronbichler, 2014/04/30)
+  </li>
+
   <li> Fixed: parallel::distributed::Triangulation::load now has an
   additional parameter <code>autopartition</code> to control p4est's behavior of
   rebalancing triangulations between MPI nodes upon reading. It is
index 8fe34b6c6028adbc68802c063d8495bf3ee56dd6..0ca4d1c021475f1d38f533bc03261857888f411d 100644 (file)
@@ -696,10 +696,17 @@ FE_DGQ<dim, spacedim>::has_support_on_face (const unsigned int shape_index,
 
   unsigned int n = this->degree+1;
 
-  // for DGQ(0) elements, the single
-  // shape functions is constant and
+  // for DGQ(0) elements or arbitrary node DGQ with support points not located
+  // at the element boundary, the single shape functions is constant and
   // therefore lives on the boundary
-  if (this->degree == 0)
+  bool support_points_on_boundary = true;
+  for (unsigned int d=0; d<dim; ++d)
+    if (std::abs(this->unit_support_points[0][d]) > 1e-13)
+      support_points_on_boundary = false;
+  for (unsigned int d=0; d<dim; ++d)
+    if (std::abs(this->unit_support_points.back()[d]-1.) > 1e-13)
+      support_points_on_boundary = false;
+  if (support_points_on_boundary == false)
     return true;
 
   unsigned int n2 = n*n;
index a9e16506a10e6917395fad236f855145929069ba..d554c5cb35cbfb1ce917e2fd85de6230d31dcd26 100644 (file)
@@ -21,6 +21,7 @@
 #include <fstream>
 
 #include <deal.II/base/logstream.h>
+#include <deal.II/base/quadrature_lib.h>
 #include <deal.II/base/polynomials_raviart_thomas.h>
 
 #include <deal.II/grid/tria_iterator.h>
@@ -112,6 +113,10 @@ void test_fe_datas()
   deallog << (*fe_datas.rbegin())->get_name() << std::endl;
   fe_datas.push_back(new FE_DGQ<dim> (4));
   deallog << (*fe_datas.rbegin())->get_name() << std::endl;
+  fe_datas.push_back(new FE_DGQArbitraryNodes<dim> (QGaussLobatto<1>(5)));
+  deallog << (*fe_datas.rbegin())->get_name() << std::endl;
+  fe_datas.push_back(new FE_DGQArbitraryNodes<dim> (QGauss<1>(3)));
+  deallog << (*fe_datas.rbegin())->get_name() << std::endl;
   fe_datas.push_back(new FE_DGP<dim> (1));
   deallog << (*fe_datas.rbegin())->get_name() << std::endl;
   fe_datas.push_back(new FE_DGP<dim> (2));
index f99beab9f4b408bd69fa1112692df27e4033daf7..b3cdb29e7c1dd9a8cf795785d62bc185c1662c97 100644 (file)
@@ -8,6 +8,8 @@ DEAL::FE_Q_Hierarchical<1>(4)
 DEAL::FE_DGQ<1>(1)
 DEAL::FE_DGQ<1>(2)
 DEAL::FE_DGQ<1>(4)
+DEAL::FE_DGQArbitraryNodes<1>(QGaussLobatto(5))
+DEAL::FE_DGQArbitraryNodes<1>(QUnknownNodes(2))
 DEAL::FE_DGP<1>(1)
 DEAL::FE_DGP<1>(2)
 DEAL::FESystem<1>[FE_Q<1>(2)^2]
@@ -270,7 +272,61 @@ DEAL::face_to_cell_index:
 DEAL::support on face 0:       0
 DEAL::support on face 1:       4
 DEAL::
-DEAL::fe_data[9]:FE_DGP<1>(1)
+DEAL::fe_data[9]:FE_DGQArbitraryNodes<1>(QGaussLobatto(5))
+DEAL::dofs_per_vertex=0
+DEAL::dofs_per_line=5
+DEAL::dofs_per_quad=0
+DEAL::dofs_per_hex=0
+DEAL::first_line_index=0
+DEAL::first_quad_index=5
+DEAL::first_hex_index=5
+DEAL::first_face_line_index=0
+DEAL::first_face_quad_index=0
+DEAL::dofs_per_face=0
+DEAL::dofs_per_cell=5
+DEAL::primitive=yes
+DEAL::components=1
+DEAL::blocks=1:[5]->5
+DEAL::degree=4
+DEAL::conformity= L2
+DEAL::unit_support_points=5
+DEAL::unit_face_support_points=0
+DEAL::generalized_support_points=5
+DEAL::generalized_face_support_points=0
+DEAL::face_to_equivalent_cell_index:
+DEAL::face_to_cell_index:
+DEAL::face_to_cell_index:
+DEAL::support on face 0:       0
+DEAL::support on face 1:       4
+DEAL::
+DEAL::fe_data[10]:FE_DGQArbitraryNodes<1>(QUnknownNodes(2))
+DEAL::dofs_per_vertex=0
+DEAL::dofs_per_line=3
+DEAL::dofs_per_quad=0
+DEAL::dofs_per_hex=0
+DEAL::first_line_index=0
+DEAL::first_quad_index=3
+DEAL::first_hex_index=3
+DEAL::first_face_line_index=0
+DEAL::first_face_quad_index=0
+DEAL::dofs_per_face=0
+DEAL::dofs_per_cell=3
+DEAL::primitive=yes
+DEAL::components=1
+DEAL::blocks=1:[3]->3
+DEAL::degree=2
+DEAL::conformity= L2
+DEAL::unit_support_points=3
+DEAL::unit_face_support_points=0
+DEAL::generalized_support_points=3
+DEAL::generalized_face_support_points=0
+DEAL::face_to_equivalent_cell_index:
+DEAL::face_to_cell_index:
+DEAL::face_to_cell_index:
+DEAL::support on face 0:       0       1       2
+DEAL::support on face 1:       0       1       2
+DEAL::
+DEAL::fe_data[11]:FE_DGP<1>(1)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=2
 DEAL::dofs_per_quad=0
@@ -297,7 +353,7 @@ DEAL::face_to_cell_index:
 DEAL::support on face 0:       0       1
 DEAL::support on face 1:       0       1
 DEAL::
-DEAL::fe_data[10]:FE_DGP<1>(2)
+DEAL::fe_data[12]:FE_DGP<1>(2)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=3
 DEAL::dofs_per_quad=0
@@ -324,7 +380,7 @@ DEAL::face_to_cell_index:
 DEAL::support on face 0:       0       1       2
 DEAL::support on face 1:       0       1       2
 DEAL::
-DEAL::fe_data[11]:FESystem<1>[FE_Q<1>(2)^2]
+DEAL::fe_data[13]:FESystem<1>[FE_Q<1>(2)^2]
 DEAL::dofs_per_vertex=2
 DEAL::dofs_per_line=2
 DEAL::dofs_per_quad=0
@@ -351,7 +407,7 @@ DEAL::face_to_cell_index: 2 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>(2)]
+DEAL::fe_data[14]:FESystem<1>[FE_Q<1>(1)^2-FE_Q<1>(2)]
 DEAL::dofs_per_vertex=3
 DEAL::dofs_per_line=1
 DEAL::dofs_per_quad=0
@@ -378,7 +434,7 @@ DEAL::face_to_cell_index: 3 4 5
 DEAL::support on face 0:       0       1       2
 DEAL::support on face 1:       3       4       5
 DEAL::
-DEAL::fe_data[13]:FE_FaceQ<1>(0)
+DEAL::fe_data[15]:FE_FaceQ<1>(0)
 DEAL::dofs_per_vertex=1
 DEAL::dofs_per_line=0
 DEAL::dofs_per_quad=0
@@ -405,7 +461,7 @@ DEAL::face_to_cell_index: 1
 DEAL::support on face 0:       0
 DEAL::support on face 1:       1
 DEAL::
-DEAL::fe_data[14]:FE_FaceQ<1>(1)
+DEAL::fe_data[16]:FE_FaceQ<1>(1)
 DEAL::dofs_per_vertex=1
 DEAL::dofs_per_line=0
 DEAL::dofs_per_quad=0
@@ -432,7 +488,7 @@ DEAL::face_to_cell_index: 1
 DEAL::support on face 0:       0
 DEAL::support on face 1:       1
 DEAL::
-DEAL::fe_data[15]:FE_FaceQ<1>(3)
+DEAL::fe_data[17]:FE_FaceQ<1>(3)
 DEAL::dofs_per_vertex=1
 DEAL::dofs_per_line=0
 DEAL::dofs_per_quad=0
@@ -459,7 +515,7 @@ DEAL::face_to_cell_index: 1
 DEAL::support on face 0:       0
 DEAL::support on face 1:       1
 DEAL::
-DEAL::fe_data[16]:FE_FaceP<1>(0)
+DEAL::fe_data[18]:FE_FaceP<1>(0)
 DEAL::dofs_per_vertex=1
 DEAL::dofs_per_line=0
 DEAL::dofs_per_quad=0
@@ -486,7 +542,7 @@ DEAL::face_to_cell_index: 1
 DEAL::support on face 0:       0
 DEAL::support on face 1:       1
 DEAL::
-DEAL::fe_data[17]:FE_FaceP<1>(1)
+DEAL::fe_data[19]:FE_FaceP<1>(1)
 DEAL::dofs_per_vertex=1
 DEAL::dofs_per_line=0
 DEAL::dofs_per_quad=0
@@ -513,7 +569,7 @@ DEAL::face_to_cell_index: 1
 DEAL::support on face 0:       0
 DEAL::support on face 1:       1
 DEAL::
-DEAL::fe_data[18]:FE_FaceP<1>(3)
+DEAL::fe_data[20]:FE_FaceP<1>(3)
 DEAL::dofs_per_vertex=1
 DEAL::dofs_per_line=0
 DEAL::dofs_per_quad=0
@@ -540,7 +596,7 @@ DEAL::face_to_cell_index: 1
 DEAL::support on face 0:       0
 DEAL::support on face 1:       1
 DEAL::
-DEAL::fe_data[19]:FESystem<1>[FE_Q<1>(3)^2]
+DEAL::fe_data[21]:FESystem<1>[FE_Q<1>(3)^2]
 DEAL::dofs_per_vertex=2
 DEAL::dofs_per_line=4
 DEAL::dofs_per_quad=0
@@ -567,7 +623,7 @@ DEAL::face_to_cell_index: 2 3
 DEAL::support on face 0:       0       1
 DEAL::support on face 1:       2       3
 DEAL::
-DEAL::fe_data[20]:FESystem<1>[FE_Q<1>(1)^2-FE_Q<1>(3)]
+DEAL::fe_data[22]:FESystem<1>[FE_Q<1>(1)^2-FE_Q<1>(3)]
 DEAL::dofs_per_vertex=3
 DEAL::dofs_per_line=2
 DEAL::dofs_per_quad=0
@@ -594,7 +650,7 @@ DEAL::face_to_cell_index: 3 4 5
 DEAL::support on face 0:       0       1       2
 DEAL::support on face 1:       3       4       5
 DEAL::
-DEAL::fe_data[21]:FESystem<1>[FE_Q<1>(4)^2]
+DEAL::fe_data[23]:FESystem<1>[FE_Q<1>(4)^2]
 DEAL::dofs_per_vertex=2
 DEAL::dofs_per_line=6
 DEAL::dofs_per_quad=0
@@ -621,7 +677,7 @@ DEAL::face_to_cell_index: 2 3
 DEAL::support on face 0:       0       1
 DEAL::support on face 1:       2       3
 DEAL::
-DEAL::fe_data[22]:FESystem<1>[FESystem<1>[FE_Q<1>(1)^2]^2]
+DEAL::fe_data[24]:FESystem<1>[FESystem<1>[FE_Q<1>(1)^2]^2]
 DEAL::dofs_per_vertex=4
 DEAL::dofs_per_line=0
 DEAL::dofs_per_quad=0
@@ -648,7 +704,7 @@ DEAL::face_to_cell_index: 4 5 6 7
 DEAL::support on face 0:       0       1       2       3
 DEAL::support on face 1:       4       5       6       7
 DEAL::
-DEAL::fe_data[23]:FESystem<1>[FESystem<1>[FE_Q<1>(1)^2]-FESystem<1>[FE_DGQ<1>(1)^2]]
+DEAL::fe_data[25]: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
 DEAL::dofs_per_quad=0
@@ -675,7 +731,7 @@ DEAL::face_to_cell_index: 2 3
 DEAL::support on face 0:       0       1       4       6
 DEAL::support on face 1:       2       3       5       7
 DEAL::
-DEAL::fe_data[24]: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::fe_data[26]: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
 DEAL::dofs_per_quad=0
@@ -702,7 +758,7 @@ DEAL::face_to_cell_index: 4 5 6 7
 DEAL::support on face 0:       0       1       2       3       11      14
 DEAL::support on face 1:       4       5       6       7       13      16
 DEAL::
-DEAL::fe_data[25]: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::fe_data[27]: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
 DEAL::dofs_per_quad=0
@@ -738,6 +794,8 @@ DEAL::FE_Q_Hierarchical<2>(4)
 DEAL::FE_DGQ<2>(1)
 DEAL::FE_DGQ<2>(2)
 DEAL::FE_DGQ<2>(4)
+DEAL::FE_DGQArbitraryNodes<2>(QGaussLobatto(5))
+DEAL::FE_DGQArbitraryNodes<2>(QUnknownNodes(2))
 DEAL::FE_DGP<2>(1)
 DEAL::FE_DGP<2>(2)
 DEAL::FESystem<2>[FE_Q<2>(2)^2]
@@ -1053,7 +1111,69 @@ DEAL::support on face 1: 4       9       14      19      24
 DEAL::support on face 2:       0       1       2       3       4
 DEAL::support on face 3:       20      21      22      23      24
 DEAL::
-DEAL::fe_data[9]:FE_DGP<2>(1)
+DEAL::fe_data[9]:FE_DGQArbitraryNodes<2>(QGaussLobatto(5))
+DEAL::dofs_per_vertex=0
+DEAL::dofs_per_line=0
+DEAL::dofs_per_quad=25
+DEAL::dofs_per_hex=0
+DEAL::first_line_index=0
+DEAL::first_quad_index=0
+DEAL::first_hex_index=25
+DEAL::first_face_line_index=0
+DEAL::first_face_quad_index=0
+DEAL::dofs_per_face=0
+DEAL::dofs_per_cell=25
+DEAL::primitive=yes
+DEAL::components=1
+DEAL::blocks=1:[25]->25
+DEAL::degree=4
+DEAL::conformity= L2
+DEAL::unit_support_points=25
+DEAL::unit_face_support_points=0
+DEAL::generalized_support_points=25
+DEAL::generalized_face_support_points=0
+DEAL::face_to_equivalent_cell_index:
+DEAL::face_to_cell_index:
+DEAL::face_to_cell_index:
+DEAL::face_to_cell_index:
+DEAL::face_to_cell_index:
+DEAL::support on face 0:       0       5       10      15      20
+DEAL::support on face 1:       4       9       14      19      24
+DEAL::support on face 2:       0       1       2       3       4
+DEAL::support on face 3:       20      21      22      23      24
+DEAL::
+DEAL::fe_data[10]:FE_DGQArbitraryNodes<2>(QUnknownNodes(2))
+DEAL::dofs_per_vertex=0
+DEAL::dofs_per_line=0
+DEAL::dofs_per_quad=9
+DEAL::dofs_per_hex=0
+DEAL::first_line_index=0
+DEAL::first_quad_index=0
+DEAL::first_hex_index=9
+DEAL::first_face_line_index=0
+DEAL::first_face_quad_index=0
+DEAL::dofs_per_face=0
+DEAL::dofs_per_cell=9
+DEAL::primitive=yes
+DEAL::components=1
+DEAL::blocks=1:[9]->9
+DEAL::degree=2
+DEAL::conformity= L2
+DEAL::unit_support_points=9
+DEAL::unit_face_support_points=0
+DEAL::generalized_support_points=9
+DEAL::generalized_face_support_points=0
+DEAL::face_to_equivalent_cell_index:
+DEAL::face_to_cell_index:
+DEAL::face_to_cell_index:
+DEAL::face_to_cell_index:
+DEAL::face_to_cell_index:
+DEAL::support on face 0:       0       1       2       3       4       5       6       7       8
+DEAL::support on face 1:       0       1       2       3       4       5       6       7       8
+DEAL::support on face 2:       0       1       2       3       4       5       6       7       8
+DEAL::support on face 3:       0       1       2       3       4       5       6       7       8
+DEAL::
+DEAL::fe_data[11]:FE_DGP<2>(1)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
 DEAL::dofs_per_quad=3
@@ -1084,7 +1204,7 @@ 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[10]:FE_DGP<2>(2)
+DEAL::fe_data[12]:FE_DGP<2>(2)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
 DEAL::dofs_per_quad=6
@@ -1115,7 +1235,7 @@ 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[11]:FESystem<2>[FE_Q<2>(2)^2]
+DEAL::fe_data[13]:FESystem<2>[FE_Q<2>(2)^2]
 DEAL::dofs_per_vertex=2
 DEAL::dofs_per_line=2
 DEAL::dofs_per_quad=2
@@ -1146,7 +1266,7 @@ DEAL::support on face 1:  2       3       6       7       10      11
 DEAL::support on face 2:       0       1       2       3       12      13
 DEAL::support on face 3:       4       5       6       7       14      15
 DEAL::
-DEAL::fe_data[12]:FESystem<2>[FE_Q<2>(1)^2-FE_Q<2>(2)]
+DEAL::fe_data[14]:FESystem<2>[FE_Q<2>(1)^2-FE_Q<2>(2)]
 DEAL::dofs_per_vertex=3
 DEAL::dofs_per_line=1
 DEAL::dofs_per_quad=1
@@ -1177,7 +1297,7 @@ DEAL::support on face 1:  3       4       5       9       10      11      13
 DEAL::support on face 2:       0       1       2       3       4       5       14
 DEAL::support on face 3:       6       7       8       9       10      11      15
 DEAL::
-DEAL::fe_data[13]:FE_FaceQ<2>(0)
+DEAL::fe_data[15]:FE_FaceQ<2>(0)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=1
 DEAL::dofs_per_quad=0
@@ -1208,7 +1328,7 @@ DEAL::support on face 1:  1
 DEAL::support on face 2:       2
 DEAL::support on face 3:       3
 DEAL::
-DEAL::fe_data[14]:FE_FaceQ<2>(1)
+DEAL::fe_data[16]:FE_FaceQ<2>(1)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=2
 DEAL::dofs_per_quad=0
@@ -1239,7 +1359,7 @@ 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[15]:FE_FaceQ<2>(3)
+DEAL::fe_data[17]:FE_FaceQ<2>(3)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=4
 DEAL::dofs_per_quad=0
@@ -1270,7 +1390,7 @@ 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_FaceP<2>(0)
+DEAL::fe_data[18]:FE_FaceP<2>(0)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=1
 DEAL::dofs_per_quad=0
@@ -1301,7 +1421,7 @@ 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::fe_data[19]:FE_FaceP<2>(1)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=2
 DEAL::dofs_per_quad=0
@@ -1332,7 +1452,7 @@ 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::fe_data[20]:FE_FaceP<2>(3)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=4
 DEAL::dofs_per_quad=0
@@ -1363,7 +1483,7 @@ 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::fe_data[21]:FE_DGRaviartThomas<2,2>(0)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
 DEAL::dofs_per_quad=4
@@ -1394,7 +1514,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[20]:FE_DGRaviartThomas<2,2>(1)
+DEAL::fe_data[22]:FE_DGRaviartThomas<2,2>(1)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
 DEAL::dofs_per_quad=12
@@ -1425,7 +1545,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[21]:FE_DGNedelec<2,2>(0)
+DEAL::fe_data[23]:FE_DGNedelec<2,2>(0)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
 DEAL::dofs_per_quad=4
@@ -1456,7 +1576,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[22]:FE_DGNedelec<2,2>(1)
+DEAL::fe_data[24]:FE_DGNedelec<2,2>(1)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
 DEAL::dofs_per_quad=12
@@ -1487,7 +1607,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[23]:FE_RaviartThomas<2>(0)
+DEAL::fe_data[25]:FE_RaviartThomas<2>(0)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=1
 DEAL::dofs_per_quad=0
@@ -1518,7 +1638,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[24]:FE_RaviartThomas<2>(1)
+DEAL::fe_data[26]:FE_RaviartThomas<2>(1)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=2
 DEAL::dofs_per_quad=4
@@ -1549,7 +1669,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[25]:FE_RaviartThomas<2>(2)
+DEAL::fe_data[27]:FE_RaviartThomas<2>(2)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=3
 DEAL::dofs_per_quad=12
@@ -1580,7 +1700,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[26]:FESystem<2>[FE_RaviartThomas<2>(1)-FE_DGQ<2>(1)]
+DEAL::fe_data[28]: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
@@ -1611,7 +1731,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[27]:FE_Nedelec<2>(0)
+DEAL::fe_data[29]:FE_Nedelec<2>(0)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=1
 DEAL::dofs_per_quad=0
@@ -1642,7 +1762,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[28]:FE_Nedelec<2>(1)
+DEAL::fe_data[30]:FE_Nedelec<2>(1)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=2
 DEAL::dofs_per_quad=4
@@ -1673,7 +1793,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[29]:FE_DGBDM<2,2>(1)
+DEAL::fe_data[31]:FE_DGBDM<2,2>(1)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
 DEAL::dofs_per_quad=8
@@ -1704,7 +1824,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[30]:FE_DGBDM<2,2>(2)
+DEAL::fe_data[32]:FE_DGBDM<2,2>(2)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
 DEAL::dofs_per_quad=14
@@ -1735,7 +1855,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[31]:FE_BDM<2>(1)
+DEAL::fe_data[33]:FE_BDM<2>(1)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=2
 DEAL::dofs_per_quad=0
@@ -1766,7 +1886,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[32]:FE_BDM<2>(2)
+DEAL::fe_data[34]:FE_BDM<2>(2)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=3
 DEAL::dofs_per_quad=2
@@ -1797,7 +1917,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[33]:FE_RaviartThomasNodal<2>(0)
+DEAL::fe_data[35]:FE_RaviartThomasNodal<2>(0)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=1
 DEAL::dofs_per_quad=0
@@ -1828,7 +1948,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[34]:FE_RaviartThomasNodal<2>(1)
+DEAL::fe_data[36]:FE_RaviartThomasNodal<2>(1)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=2
 DEAL::dofs_per_quad=4
@@ -1859,7 +1979,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[35]:FESystem<2>[FE_RaviartThomasNodal<2>(1)-FE_DGQ<2>(1)]
+DEAL::fe_data[37]: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
@@ -1890,7 +2010,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[36]:FESystem<2>[FE_Q<2>(3)^2]
+DEAL::fe_data[38]:FESystem<2>[FE_Q<2>(3)^2]
 DEAL::dofs_per_vertex=2
 DEAL::dofs_per_line=4
 DEAL::dofs_per_quad=8
@@ -1921,7 +2041,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[37]:FESystem<2>[FE_Q<2>(1)^2-FE_Q<2>(3)]
+DEAL::fe_data[39]: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
@@ -1952,7 +2072,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[38]:FESystem<2>[FE_Q<2>(4)^2]
+DEAL::fe_data[40]:FESystem<2>[FE_Q<2>(4)^2]
 DEAL::dofs_per_vertex=2
 DEAL::dofs_per_line=6
 DEAL::dofs_per_quad=18
@@ -1983,7 +2103,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[39]:FESystem<2>[FESystem<2>[FE_Q<2>(1)^2]^2]
+DEAL::fe_data[41]: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
@@ -2014,7 +2134,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[40]:FESystem<2>[FESystem<2>[FE_Q<2>(1)^2]-FESystem<2>[FE_DGQ<2>(1)^2]]
+DEAL::fe_data[42]: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
@@ -2045,7 +2165,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[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::fe_data[43]: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
@@ -2076,7 +2196,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[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::fe_data[44]: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
@@ -2116,6 +2236,8 @@ DEAL::FE_Q_Hierarchical<3>(4)
 DEAL::FE_DGQ<3>(1)
 DEAL::FE_DGQ<3>(2)
 DEAL::FE_DGQ<3>(4)
+DEAL::FE_DGQArbitraryNodes<3>(QGaussLobatto(5))
+DEAL::FE_DGQArbitraryNodes<3>(QUnknownNodes(2))
 DEAL::FE_DGP<3>(1)
 DEAL::FE_DGP<3>(2)
 DEAL::FESystem<3>[FE_Q<3>(2)^2]
@@ -2460,7 +2582,77 @@ DEAL::support on face 3: 20      21      22      23      24      45      46      47      48      49      70      71      72      73      74      95      96      97      9
 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
 DEAL::support on face 5:       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::
-DEAL::fe_data[9]:FE_DGP<3>(1)
+DEAL::fe_data[9]:FE_DGQArbitraryNodes<3>(QGaussLobatto(5))
+DEAL::dofs_per_vertex=0
+DEAL::dofs_per_line=0
+DEAL::dofs_per_quad=0
+DEAL::dofs_per_hex=125
+DEAL::first_line_index=0
+DEAL::first_quad_index=0
+DEAL::first_hex_index=0
+DEAL::first_face_line_index=0
+DEAL::first_face_quad_index=0
+DEAL::dofs_per_face=0
+DEAL::dofs_per_cell=125
+DEAL::primitive=yes
+DEAL::components=1
+DEAL::blocks=1:[125]->125
+DEAL::degree=4
+DEAL::conformity= L2
+DEAL::unit_support_points=125
+DEAL::unit_face_support_points=0
+DEAL::generalized_support_points=125
+DEAL::generalized_face_support_points=0
+DEAL::face_to_equivalent_cell_index:
+DEAL::face_to_cell_index:
+DEAL::face_to_cell_index:
+DEAL::face_to_cell_index:
+DEAL::face_to_cell_index:
+DEAL::face_to_cell_index:
+DEAL::face_to_cell_index:
+DEAL::support on face 0:       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::support on face 1:       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 2:       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 3:       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 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
+DEAL::support on face 5:       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::
+DEAL::fe_data[10]:FE_DGQArbitraryNodes<3>(QUnknownNodes(2))
+DEAL::dofs_per_vertex=0
+DEAL::dofs_per_line=0
+DEAL::dofs_per_quad=0
+DEAL::dofs_per_hex=27
+DEAL::first_line_index=0
+DEAL::first_quad_index=0
+DEAL::first_hex_index=0
+DEAL::first_face_line_index=0
+DEAL::first_face_quad_index=0
+DEAL::dofs_per_face=0
+DEAL::dofs_per_cell=27
+DEAL::primitive=yes
+DEAL::components=1
+DEAL::blocks=1:[27]->27
+DEAL::degree=2
+DEAL::conformity= L2
+DEAL::unit_support_points=27
+DEAL::unit_face_support_points=0
+DEAL::generalized_support_points=27
+DEAL::generalized_face_support_points=0
+DEAL::face_to_equivalent_cell_index:
+DEAL::face_to_cell_index:
+DEAL::face_to_cell_index:
+DEAL::face_to_cell_index:
+DEAL::face_to_cell_index:
+DEAL::face_to_cell_index:
+DEAL::face_to_cell_index:
+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
+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
+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
+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
+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
+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
+DEAL::
+DEAL::fe_data[11]:FE_DGP<3>(1)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
 DEAL::dofs_per_quad=0
@@ -2495,7 +2687,7 @@ 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[10]:FE_DGP<3>(2)
+DEAL::fe_data[12]:FE_DGP<3>(2)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
 DEAL::dofs_per_quad=0
@@ -2530,7 +2722,7 @@ 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[11]:FESystem<3>[FE_Q<3>(2)^2]
+DEAL::fe_data[13]:FESystem<3>[FE_Q<3>(2)^2]
 DEAL::dofs_per_vertex=2
 DEAL::dofs_per_line=2
 DEAL::dofs_per_quad=2
@@ -2565,7 +2757,7 @@ DEAL::support on face 3:  4       5       6       7       12      13      14      15      22      23      30      31      36      37      38      39      46      47
 DEAL::support on face 4:       0       1       2       3       4       5       6       7       16      17      18      19      20      21      22      23      48      49
 DEAL::support on face 5:       8       9       10      11      12      13      14      15      24      25      26      27      28      29      30      31      50      51
 DEAL::
-DEAL::fe_data[12]:FESystem<3>[FE_Q<3>(1)^2-FE_Q<3>(2)]
+DEAL::fe_data[14]:FESystem<3>[FE_Q<3>(1)^2-FE_Q<3>(2)]
 DEAL::dofs_per_vertex=3
 DEAL::dofs_per_line=1
 DEAL::dofs_per_quad=1
@@ -2600,7 +2792,7 @@ DEAL::support on face 3:  6       7       8       9       10      11      18      19      20      21      22      23      27      31      34      35      39
 DEAL::support on face 4:       0       1       2       3       4       5       6       7       8       9       10      11      24      25      26      27      40
 DEAL::support on face 5:       12      13      14      15      16      17      18      19      20      21      22      23      28      29      30      31      41
 DEAL::
-DEAL::fe_data[13]:FE_FaceQ<3>(0)
+DEAL::fe_data[15]:FE_FaceQ<3>(0)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
 DEAL::dofs_per_quad=1
@@ -2635,7 +2827,7 @@ DEAL::support on face 3:  3
 DEAL::support on face 4:       4
 DEAL::support on face 5:       5
 DEAL::
-DEAL::fe_data[14]:FE_FaceQ<3>(1)
+DEAL::fe_data[16]:FE_FaceQ<3>(1)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
 DEAL::dofs_per_quad=4
@@ -2670,7 +2862,7 @@ DEAL::support on face 3:  12      13      14      15
 DEAL::support on face 4:       16      17      18      19
 DEAL::support on face 5:       20      21      22      23
 DEAL::
-DEAL::fe_data[15]:FE_FaceQ<3>(3)
+DEAL::fe_data[17]:FE_FaceQ<3>(3)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
 DEAL::dofs_per_quad=16
@@ -2705,7 +2897,7 @@ 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_FaceP<3>(0)
+DEAL::fe_data[18]:FE_FaceP<3>(0)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
 DEAL::dofs_per_quad=1
@@ -2740,7 +2932,7 @@ 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::fe_data[19]:FE_FaceP<3>(1)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
 DEAL::dofs_per_quad=3
@@ -2775,7 +2967,7 @@ 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::fe_data[20]:FE_FaceP<3>(3)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
 DEAL::dofs_per_quad=10
@@ -2810,7 +3002,7 @@ 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::fe_data[21]:FE_DGRaviartThomas<3,3>(0)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
 DEAL::dofs_per_quad=0
@@ -2845,7 +3037,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[20]:FE_DGRaviartThomas<3,3>(1)
+DEAL::fe_data[22]:FE_DGRaviartThomas<3,3>(1)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
 DEAL::dofs_per_quad=0
@@ -2880,7 +3072,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[21]:FE_DGNedelec<3,3>(0)
+DEAL::fe_data[23]:FE_DGNedelec<3,3>(0)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
 DEAL::dofs_per_quad=0
@@ -2915,7 +3107,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[22]:FE_DGNedelec<3,3>(1)
+DEAL::fe_data[24]:FE_DGNedelec<3,3>(1)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
 DEAL::dofs_per_quad=0
@@ -2950,7 +3142,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[23]:FE_RaviartThomas<3>(0)
+DEAL::fe_data[25]:FE_RaviartThomas<3>(0)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
 DEAL::dofs_per_quad=1
@@ -2985,7 +3177,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[24]:FE_RaviartThomas<3>(1)
+DEAL::fe_data[26]:FE_RaviartThomas<3>(1)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
 DEAL::dofs_per_quad=4
@@ -3020,7 +3212,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[25]:FE_RaviartThomas<3>(2)
+DEAL::fe_data[27]:FE_RaviartThomas<3>(2)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
 DEAL::dofs_per_quad=9
@@ -3055,7 +3247,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[26]:FESystem<3>[FE_RaviartThomas<3>(1)-FE_DGQ<3>(1)]
+DEAL::fe_data[28]: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
@@ -3090,7 +3282,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[27]:FE_Nedelec<3>(0)
+DEAL::fe_data[29]:FE_Nedelec<3>(0)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=1
 DEAL::dofs_per_quad=0
@@ -3125,7 +3317,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[28]:FE_Nedelec<3>(1)
+DEAL::fe_data[30]:FE_Nedelec<3>(1)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=2
 DEAL::dofs_per_quad=4
@@ -3160,7 +3352,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[29]:FE_RaviartThomasNodal<3>(0)
+DEAL::fe_data[31]:FE_RaviartThomasNodal<3>(0)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
 DEAL::dofs_per_quad=1
@@ -3195,7 +3387,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[30]:FE_RaviartThomasNodal<3>(1)
+DEAL::fe_data[32]:FE_RaviartThomasNodal<3>(1)
 DEAL::dofs_per_vertex=0
 DEAL::dofs_per_line=0
 DEAL::dofs_per_quad=4
@@ -3230,7 +3422,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[31]:FESystem<3>[FE_RaviartThomasNodal<3>(1)-FE_DGQ<3>(1)]
+DEAL::fe_data[33]: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
@@ -3265,7 +3457,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[32]:FESystem<3>[FESystem<3>[FE_Q<3>(1)^2]^2]
+DEAL::fe_data[34]: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
@@ -3300,7 +3492,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[33]:FESystem<3>[FESystem<3>[FE_Q<3>(1)^2]-FESystem<3>[FE_DGQ<3>(1)^2]]
+DEAL::fe_data[35]: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
@@ -3335,7 +3527,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[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::fe_data[36]: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
@@ -3370,7 +3562,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[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::fe_data[37]: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

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.