]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Fix the one last testcase, where we ran into trouble generating constraints that...
authorbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 15 Dec 2006 21:13:01 +0000 (21:13 +0000)
committerbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 15 Dec 2006 21:13:01 +0000 (21:13 +0000)
git-svn-id: https://svn.dealii.org/trunk@14258 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/source/dofs/dof_tools.cc
tests/hp/hp_constraints_q_system_06/cmp/generic

index 3aeaad1958f3d3fb41881e31551a06adbd23b7dd..9352ee66d5ca1107089facf5b0197c55e3a41a1d 100644 (file)
@@ -1742,76 +1742,42 @@ namespace internal
                                         // always false warning for the cases
                                         // at the bottom in 1d and 2d)
                                         //
-                                        // as mentioned in the paper,
-                                        // it is not always easy to
-                                        // find a set of master dofs
-                                        // that produces an
-                                        // invertible matrix. to this
-                                        // end, we check in each step
-                                        // whether the matrix is
-                                        // still invertible. we have
-                                        // never had a problem with
-                                        // vertex and line dofs, so
-                                        // we only check inside an
-                                        // Assert statement, without
-                                        // providing any code that
-                                        // can cope with the
-                                        // situation if this goes
-                                        // wrong.
+                                        // as mentioned in the paper, it is
+                                        // not always easy to find a set of
+                                        // master dofs that produces an
+                                        // invertible matrix. to this end, we
+                                        // check in each step whether the
+                                        // matrix is still invertible and
+                                        // simply discard this dof if the
+                                        // matrix is not invertible anymore.
                                         //
-                                        // on the other hand, we did
-                                        // have trouble with adding
-                                        // more quad dofs when Q3 and
-                                        // Q4 elements meet at a
-                                        // refined face in 3d (see
-                                        // the hp/crash_12 test that
-                                        // tests that we can do
-                                        // exactly this, and failed
-                                        // before we had code to
-                                        // compensate for this
-                                        // case). therefore, for quad
-                                        // dofs (which appear only in
-                                        // 3d anyway), we have a
-                                        // piece of code that can add
-                                        // dofs until we have finally
-                                        // found a suitable set
+                                        // the cases where we did have
+                                        // trouble in the past were with
+                                        // adding more quad dofs when Q3 and
+                                        // Q4 elements meet at a refined face
+                                        // in 3d (see the hp/crash_12 test
+                                        // that tests that we can do exactly
+                                        // this, and failed before we had
+                                        // code to compensate for this
+                                        // case). the other case are system
+                                        // elements: if we have say a Q1Q2 vs
+                                        // a Q2Q3 element, then we can't just
+                                        // take all master dofs on a line
+                                        // from a single base element, since
+                                        // the shape functions of that base
+                                        // element are independent of that of
+                                        // the other one. this latter case
+                                        // shows up when running
+                                        // hp/hp_constraints_q_system_06
        std::vector<unsigned int> master_dof_list;
        unsigned int index = 0;
        for (int v=0;
             v<static_cast<signed int>(GeometryInfo<dim>::vertices_per_face);
             ++v)
-         {
-           for (unsigned int i=0; i<fe2.dofs_per_vertex; ++i)
-             {
-               master_dof_list.push_back (index+i);
-               Assert (check_master_dof_list (face_interpolation_matrix,
-                                              master_dof_list),
-                       ExcInternalError());
-             }
-           index += fe1.dofs_per_vertex;
-         }
-
-       for (int l=0;
-            l<static_cast<signed int>(GeometryInfo<dim>::lines_per_face);
-            ++l)
-         {
-           for (unsigned int i=0; i<fe2.dofs_per_line; ++i)
-             {
-               master_dof_list.push_back (index+i);
-               Assert (check_master_dof_list (face_interpolation_matrix,
-                                              master_dof_list),
-                       ExcInternalError());
-             }
-           index += fe1.dofs_per_line;
-         }
-       
-       for (int q=0;
-            q<static_cast<signed int>(GeometryInfo<dim>::quads_per_face);
-            ++q)
          {
            unsigned int dofs_added = 0;
            unsigned int i          = 0;
-           while (dofs_added < fe2.dofs_per_quad)
+           while (dofs_added < fe2.dofs_per_vertex)
              {
                                                 // make sure that we
                                                 // were able to find
@@ -1820,10 +1786,10 @@ namespace internal
                                                 // code down below
                                                 // didn't just reject
                                                 // all our efforts
-               Assert (i < fe1.dofs_per_quad,
+               Assert (i < fe1.dofs_per_vertex,
                        ExcInternalError());
                                                 // tentatively push
-                                                // this quad dof
+                                                // this vertex dof
                master_dof_list.push_back (index+i);
 
                                                 // then see what
@@ -1847,6 +1813,56 @@ namespace internal
                                                 // one
                ++i;
              }
+           index += fe1.dofs_per_vertex;           
+         }
+
+       for (int l=0;
+            l<static_cast<signed int>(GeometryInfo<dim>::lines_per_face);
+            ++l)
+         {
+                                            // same algorithm as above
+           unsigned int dofs_added = 0;
+           unsigned int i          = 0;
+           while (dofs_added < fe2.dofs_per_line)
+             {
+               Assert (i < fe1.dofs_per_line,
+                       ExcInternalError());
+               
+               master_dof_list.push_back (index+i);
+               if (check_master_dof_list (face_interpolation_matrix,
+                                          master_dof_list)
+                   == true)
+                 ++dofs_added;
+               else
+                 master_dof_list.pop_back ();
+
+               ++i;
+             }
+           index += fe1.dofs_per_line;
+         }
+       
+       for (int q=0;
+            q<static_cast<signed int>(GeometryInfo<dim>::quads_per_face);
+            ++q)
+         {
+                                            // same algorithm as above
+           unsigned int dofs_added = 0;
+           unsigned int i          = 0;
+           while (dofs_added < fe2.dofs_per_quad)
+             {
+               Assert (i < fe1.dofs_per_quad,
+                       ExcInternalError());
+               
+               master_dof_list.push_back (index+i);
+               if (check_master_dof_list (face_interpolation_matrix,
+                                          master_dof_list)
+                   == true)
+                 ++dofs_added;
+               else
+                 master_dof_list.pop_back ();
+
+               ++i;
+             }
            index += fe1.dofs_per_quad;
          }
 
index 81cd13989996f79da1f86a6b1dcb38dd9ac09317..e3d65bbe7f2e8ed5e8cf16f1b8bc37968e264225 100644 (file)
 
-DEAL::n_dofs=140
-    27 0:  0.500
-    27 46:  0.500
-    48 0:  0.500
-    48 46:  0.500
-    61 47:  0.500
-    61 49:  0.500
-    82 0:  0.667
-    82 46:  0.333
-    83 0:  0.833
-    83 46:  0.167
-    84 6:  -0.111
-    84 49:  0.222
-    84 120:  0.889
-    85 6:  0.222
-    85 49:  -0.111
-    85 120:  0.889
-    86 0:  0.333
-    86 46:  0.333
-    86 49:  0.333
-    87 0:  0.167
-    87 46:  0.167
-    87 49:  0.667
-    88 0:  0.556
-    88 1:  -0.111
-    88 6:  0.556
-    89 0:  0.222
-    89 1:  -0.111
-    89 6:  0.889
-    119 6:  1.00
-    122 0:  -0.125
-    122 1:  0.375
-    122 6:  0.750
-DEAL::n_dofs=129
-    2 23:  0.500
-    2 24:  0.500
-    5 1:  0.500
-    5 24:  0.500
-    36 1:  0.500
-    36 24:  0.500
-    37 24:  0.500
-    37 26:  0.500
-    49 1:  0.250
-    49 24:  0.250
-    49 38:  0.500
-    51 1:  0.750
-    51 24:  0.250
-    70 24:  0.333
-    70 26:  0.667
-    71 24:  0.167
-    71 26:  0.833
-    72 38:  0.222
-    72 107:  -0.111
-    72 109:  0.889
-    73 38:  -0.111
-    73 107:  0.222
-    73 109:  0.889
-    74 24:  0.333
-    74 26:  0.333
-    74 38:  0.333
-    75 24:  0.167
-    75 26:  0.167
-    75 38:  0.667
-DEAL::n_dofs=146
-    3 0:  -0.333
-    3 2:  1.00
-    3 39:  0.333
-    8 39:  0.222
-    8 40:  -0.111
-    8 44:  0.889
-    9 39:  -0.111
-    9 40:  0.222
-    9 44:  0.889
-    42 39:  0.500
-    42 136:  0.500
-    65 0:  -0.0556
-    65 2:  0.250
-    65 39:  0.0278
-    65 101:  0.889
-    65 123:  -0.111
-    66 0:  0.0278
-    66 2:  -0.125
-    66 39:  -0.0139
-    66 101:  0.889
-    66 123:  0.222
-    67 0:  0.417
-    67 2:  0.625
-    67 39:  -0.0417
-    68 0:  -2.78e-17
-    68 2:  1.00
-    69 122:  0.667
-    69 123:  0.333
-    70 122:  0.333
-    70 123:  0.667
-    100 0:  -0.250
-    100 2:  1.12
-    100 39:  0.125
-    102 39:  0.750
-    102 136:  0.250
-    103 0:  -0.312
-    103 2:  0.844
-    103 39:  0.469
-    104 39:  0.250
-    104 123:  0.500
-    104 136:  0.250
-    135 39:  0.500
-    135 136:  0.500
+DEAL::Testing FESystem<1>[FE_Q<1>(1)-FE_Q<1>(2)] vs. FESystem<1>[FE_Q<1>(1)-FE_Q<1>(2)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(1)-FE_Q<1>(2)] vs. FESystem<1>[FE_Q<1>(2)-FE_Q<1>(3)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(1)-FE_Q<1>(2)] vs. FESystem<1>[FE_Q<1>(3)-FE_Q<1>(4)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(1)-FE_Q<1>(2)] vs. FESystem<1>[FE_Q<1>(4)-FE_Q<1>(5)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(1)-FE_Q<1>(2)] vs. FESystem<1>[FE_Q<1>(5)-FE_Q<1>(6)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(2)-FE_Q<1>(3)] vs. FESystem<1>[FE_Q<1>(1)-FE_Q<1>(2)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(2)-FE_Q<1>(3)] vs. FESystem<1>[FE_Q<1>(2)-FE_Q<1>(3)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(2)-FE_Q<1>(3)] vs. FESystem<1>[FE_Q<1>(3)-FE_Q<1>(4)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(2)-FE_Q<1>(3)] vs. FESystem<1>[FE_Q<1>(4)-FE_Q<1>(5)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(2)-FE_Q<1>(3)] vs. FESystem<1>[FE_Q<1>(5)-FE_Q<1>(6)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(3)-FE_Q<1>(4)] vs. FESystem<1>[FE_Q<1>(1)-FE_Q<1>(2)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(3)-FE_Q<1>(4)] vs. FESystem<1>[FE_Q<1>(2)-FE_Q<1>(3)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(3)-FE_Q<1>(4)] vs. FESystem<1>[FE_Q<1>(3)-FE_Q<1>(4)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(3)-FE_Q<1>(4)] vs. FESystem<1>[FE_Q<1>(4)-FE_Q<1>(5)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(3)-FE_Q<1>(4)] vs. FESystem<1>[FE_Q<1>(5)-FE_Q<1>(6)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(4)-FE_Q<1>(5)] vs. FESystem<1>[FE_Q<1>(1)-FE_Q<1>(2)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(4)-FE_Q<1>(5)] vs. FESystem<1>[FE_Q<1>(2)-FE_Q<1>(3)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(4)-FE_Q<1>(5)] vs. FESystem<1>[FE_Q<1>(3)-FE_Q<1>(4)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(4)-FE_Q<1>(5)] vs. FESystem<1>[FE_Q<1>(4)-FE_Q<1>(5)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(4)-FE_Q<1>(5)] vs. FESystem<1>[FE_Q<1>(5)-FE_Q<1>(6)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(5)-FE_Q<1>(6)] vs. FESystem<1>[FE_Q<1>(1)-FE_Q<1>(2)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(5)-FE_Q<1>(6)] vs. FESystem<1>[FE_Q<1>(2)-FE_Q<1>(3)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(5)-FE_Q<1>(6)] vs. FESystem<1>[FE_Q<1>(3)-FE_Q<1>(4)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(5)-FE_Q<1>(6)] vs. FESystem<1>[FE_Q<1>(4)-FE_Q<1>(5)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(5)-FE_Q<1>(6)] vs. FESystem<1>[FE_Q<1>(5)-FE_Q<1>(6)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(1)-FE_Q<1>(2)] vs. FESystem<1>[FE_Q<1>(1)-FE_Q<1>(2)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(1)-FE_Q<1>(2)] vs. FESystem<1>[FE_Q<1>(2)-FE_Q<1>(3)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(1)-FE_Q<1>(2)] vs. FESystem<1>[FE_Q<1>(3)-FE_Q<1>(4)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(1)-FE_Q<1>(2)] vs. FESystem<1>[FE_Q<1>(4)-FE_Q<1>(5)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(1)-FE_Q<1>(2)] vs. FESystem<1>[FE_Q<1>(5)-FE_Q<1>(6)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(2)-FE_Q<1>(3)] vs. FESystem<1>[FE_Q<1>(1)-FE_Q<1>(2)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(2)-FE_Q<1>(3)] vs. FESystem<1>[FE_Q<1>(2)-FE_Q<1>(3)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(2)-FE_Q<1>(3)] vs. FESystem<1>[FE_Q<1>(3)-FE_Q<1>(4)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(2)-FE_Q<1>(3)] vs. FESystem<1>[FE_Q<1>(4)-FE_Q<1>(5)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(2)-FE_Q<1>(3)] vs. FESystem<1>[FE_Q<1>(5)-FE_Q<1>(6)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(3)-FE_Q<1>(4)] vs. FESystem<1>[FE_Q<1>(1)-FE_Q<1>(2)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(3)-FE_Q<1>(4)] vs. FESystem<1>[FE_Q<1>(2)-FE_Q<1>(3)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(3)-FE_Q<1>(4)] vs. FESystem<1>[FE_Q<1>(3)-FE_Q<1>(4)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(3)-FE_Q<1>(4)] vs. FESystem<1>[FE_Q<1>(4)-FE_Q<1>(5)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(3)-FE_Q<1>(4)] vs. FESystem<1>[FE_Q<1>(5)-FE_Q<1>(6)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(4)-FE_Q<1>(5)] vs. FESystem<1>[FE_Q<1>(1)-FE_Q<1>(2)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(4)-FE_Q<1>(5)] vs. FESystem<1>[FE_Q<1>(2)-FE_Q<1>(3)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(4)-FE_Q<1>(5)] vs. FESystem<1>[FE_Q<1>(3)-FE_Q<1>(4)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(4)-FE_Q<1>(5)] vs. FESystem<1>[FE_Q<1>(4)-FE_Q<1>(5)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(4)-FE_Q<1>(5)] vs. FESystem<1>[FE_Q<1>(5)-FE_Q<1>(6)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(5)-FE_Q<1>(6)] vs. FESystem<1>[FE_Q<1>(1)-FE_Q<1>(2)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(5)-FE_Q<1>(6)] vs. FESystem<1>[FE_Q<1>(2)-FE_Q<1>(3)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(5)-FE_Q<1>(6)] vs. FESystem<1>[FE_Q<1>(3)-FE_Q<1>(4)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(5)-FE_Q<1>(6)] vs. FESystem<1>[FE_Q<1>(4)-FE_Q<1>(5)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<1>[FE_Q<1>(5)-FE_Q<1>(6)] vs. FESystem<1>[FE_Q<1>(5)-FE_Q<1>(6)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<2>[FE_Q<2>(1)-FE_Q<2>(2)] vs. FESystem<2>[FE_Q<2>(1)-FE_Q<2>(2)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<2>[FE_Q<2>(1)-FE_Q<2>(2)] vs. FESystem<2>[FE_Q<2>(2)-FE_Q<2>(3)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<2>[FE_Q<2>(1)-FE_Q<2>(2)] vs. FESystem<2>[FE_Q<2>(3)-FE_Q<2>(4)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<2>[FE_Q<2>(1)-FE_Q<2>(2)] vs. FESystem<2>[FE_Q<2>(4)-FE_Q<2>(5)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<2>[FE_Q<2>(2)-FE_Q<2>(3)] vs. FESystem<2>[FE_Q<2>(1)-FE_Q<2>(2)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<2>[FE_Q<2>(2)-FE_Q<2>(3)] vs. FESystem<2>[FE_Q<2>(2)-FE_Q<2>(3)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<2>[FE_Q<2>(2)-FE_Q<2>(3)] vs. FESystem<2>[FE_Q<2>(3)-FE_Q<2>(4)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<2>[FE_Q<2>(2)-FE_Q<2>(3)] vs. FESystem<2>[FE_Q<2>(4)-FE_Q<2>(5)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<2>[FE_Q<2>(3)-FE_Q<2>(4)] vs. FESystem<2>[FE_Q<2>(1)-FE_Q<2>(2)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<2>[FE_Q<2>(3)-FE_Q<2>(4)] vs. FESystem<2>[FE_Q<2>(2)-FE_Q<2>(3)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<2>[FE_Q<2>(3)-FE_Q<2>(4)] vs. FESystem<2>[FE_Q<2>(3)-FE_Q<2>(4)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<2>[FE_Q<2>(3)-FE_Q<2>(4)] vs. FESystem<2>[FE_Q<2>(4)-FE_Q<2>(5)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<2>[FE_Q<2>(4)-FE_Q<2>(5)] vs. FESystem<2>[FE_Q<2>(1)-FE_Q<2>(2)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<2>[FE_Q<2>(4)-FE_Q<2>(5)] vs. FESystem<2>[FE_Q<2>(2)-FE_Q<2>(3)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<2>[FE_Q<2>(4)-FE_Q<2>(5)] vs. FESystem<2>[FE_Q<2>(3)-FE_Q<2>(4)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<2>[FE_Q<2>(4)-FE_Q<2>(5)] vs. FESystem<2>[FE_Q<2>(4)-FE_Q<2>(5)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<2>[FE_Q<2>(1)-FE_Q<2>(2)] vs. FESystem<2>[FE_Q<2>(1)-FE_Q<2>(2)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<2>[FE_Q<2>(1)-FE_Q<2>(2)] vs. FESystem<2>[FE_Q<2>(2)-FE_Q<2>(3)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<2>[FE_Q<2>(1)-FE_Q<2>(2)] vs. FESystem<2>[FE_Q<2>(3)-FE_Q<2>(4)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<2>[FE_Q<2>(1)-FE_Q<2>(2)] vs. FESystem<2>[FE_Q<2>(4)-FE_Q<2>(5)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<2>[FE_Q<2>(2)-FE_Q<2>(3)] vs. FESystem<2>[FE_Q<2>(1)-FE_Q<2>(2)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<2>[FE_Q<2>(2)-FE_Q<2>(3)] vs. FESystem<2>[FE_Q<2>(2)-FE_Q<2>(3)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<2>[FE_Q<2>(2)-FE_Q<2>(3)] vs. FESystem<2>[FE_Q<2>(3)-FE_Q<2>(4)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<2>[FE_Q<2>(2)-FE_Q<2>(3)] vs. FESystem<2>[FE_Q<2>(4)-FE_Q<2>(5)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<2>[FE_Q<2>(3)-FE_Q<2>(4)] vs. FESystem<2>[FE_Q<2>(1)-FE_Q<2>(2)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<2>[FE_Q<2>(3)-FE_Q<2>(4)] vs. FESystem<2>[FE_Q<2>(2)-FE_Q<2>(3)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<2>[FE_Q<2>(3)-FE_Q<2>(4)] vs. FESystem<2>[FE_Q<2>(3)-FE_Q<2>(4)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<2>[FE_Q<2>(3)-FE_Q<2>(4)] vs. FESystem<2>[FE_Q<2>(4)-FE_Q<2>(5)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<2>[FE_Q<2>(4)-FE_Q<2>(5)] vs. FESystem<2>[FE_Q<2>(1)-FE_Q<2>(2)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<2>[FE_Q<2>(4)-FE_Q<2>(5)] vs. FESystem<2>[FE_Q<2>(2)-FE_Q<2>(3)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<2>[FE_Q<2>(4)-FE_Q<2>(5)] vs. FESystem<2>[FE_Q<2>(3)-FE_Q<2>(4)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<2>[FE_Q<2>(4)-FE_Q<2>(5)] vs. FESystem<2>[FE_Q<2>(4)-FE_Q<2>(5)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<3>[FE_Q<3>(1)-FE_Q<3>(2)] vs. FESystem<3>[FE_Q<3>(1)-FE_Q<3>(2)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<3>[FE_Q<3>(1)-FE_Q<3>(2)] vs. FESystem<3>[FE_Q<3>(2)-FE_Q<3>(3)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<3>[FE_Q<3>(1)-FE_Q<3>(2)] vs. FESystem<3>[FE_Q<3>(3)-FE_Q<3>(4)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<3>[FE_Q<3>(2)-FE_Q<3>(3)] vs. FESystem<3>[FE_Q<3>(1)-FE_Q<3>(2)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<3>[FE_Q<3>(2)-FE_Q<3>(3)] vs. FESystem<3>[FE_Q<3>(2)-FE_Q<3>(3)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<3>[FE_Q<3>(2)-FE_Q<3>(3)] vs. FESystem<3>[FE_Q<3>(3)-FE_Q<3>(4)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<3>[FE_Q<3>(3)-FE_Q<3>(4)] vs. FESystem<3>[FE_Q<3>(1)-FE_Q<3>(2)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<3>[FE_Q<3>(3)-FE_Q<3>(4)] vs. FESystem<3>[FE_Q<3>(2)-FE_Q<3>(3)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<3>[FE_Q<3>(3)-FE_Q<3>(4)] vs. FESystem<3>[FE_Q<3>(3)-FE_Q<3>(4)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<3>[FE_Q<3>(1)-FE_Q<3>(2)] vs. FESystem<3>[FE_Q<3>(1)-FE_Q<3>(2)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<3>[FE_Q<3>(1)-FE_Q<3>(2)] vs. FESystem<3>[FE_Q<3>(2)-FE_Q<3>(3)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<3>[FE_Q<3>(1)-FE_Q<3>(2)] vs. FESystem<3>[FE_Q<3>(3)-FE_Q<3>(4)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<3>[FE_Q<3>(2)-FE_Q<3>(3)] vs. FESystem<3>[FE_Q<3>(1)-FE_Q<3>(2)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<3>[FE_Q<3>(2)-FE_Q<3>(3)] vs. FESystem<3>[FE_Q<3>(2)-FE_Q<3>(3)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<3>[FE_Q<3>(2)-FE_Q<3>(3)] vs. FESystem<3>[FE_Q<3>(3)-FE_Q<3>(4)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<3>[FE_Q<3>(3)-FE_Q<3>(4)] vs. FESystem<3>[FE_Q<3>(1)-FE_Q<3>(2)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<3>[FE_Q<3>(3)-FE_Q<3>(4)] vs. FESystem<3>[FE_Q<3>(2)-FE_Q<3>(3)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::Testing FESystem<3>[FE_Q<3>(3)-FE_Q<3>(4)] vs. FESystem<3>[FE_Q<3>(3)-FE_Q<3>(4)]
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0
+DEAL::  Relative interpolation error before constraints: 0
+DEAL::  Relative difference after constraints: 0

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.