Add a test.

diff --git a/tests/bits/fe_tools_04_overwrite.cc b/tests/bits/fe_tools_04_overwrite.cc
new file mode 100644 (file)
index 0000000..379ca5d
--- /dev/null
+++ b/tests/bits/fe_tools_04_overwrite.cc
@@ -0,0 +1,52 @@
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 2003 - 2017 by the deal.II authors
+//
+// This file is part of the deal.II library.
+//
+// The deal.II library is free software; you can use it, redistribute
+// it, and/or modify it under the terms of the GNU Lesser General
+// Public License as published by the Free Software Foundation; either
+// version 2.1 of the License, or (at your option) any later version.
+// The full text of the license can be found in the file LICENSE at
+// the top level of the deal.II distribution.
+//
+// ---------------------------------------------------------------------
+
+
+#include "../tests.h"
+#include "fe_tools_common.h"
+#include <deal.II/lac/sparsity_pattern.h>
+
+// check
+//   FETools::get_interpolation_difference_matrix
+// we used to have a bug wherein we did not clear previous content of the matrix
+// passed to this function. check that this is now fixed
+
+
+
+template <int dim>
+void
+check_this (const FiniteElement<dim> &fe1,
+            const FiniteElement<dim> &fe2)
+{
+  // only check if both elements have
+  // support points. otherwise,
+  // interpolation doesn't really
+  // work
+  if ((fe1.get_unit_support_points().size() == 0) ||
+      (fe2.get_unit_support_points().size() == 0))
+    return;
+  //  likewise for non-primitive elements
+  if (!fe1.is_primitive() || !fe2.is_primitive())
+    return;
+
+  // start with a matrix that is nonzero
+  FullMatrix<double> m = IdentityMatrix (fe1.dofs_per_cell);
+  m *= numbers::PI;
+
+  FETools::get_interpolation_difference_matrix (fe1, fe2, m);
+
+  output_matrix (m);
+}
+

diff --git a/tests/bits/fe_tools_04_overwrite.output b/tests/bits/fe_tools_04_overwrite.output
new file mode 100644 (file)
index 0000000..9a63acc
--- /dev/null
+++ b/tests/bits/fe_tools_04_overwrite.output
@@ -0,0 +1,503 @@
+
+DEAL::Checking Q1 against Q1 in 1d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 
+DEAL::Checking Q1 against Q1 in 1d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 
+DEAL::Checking Q1 against Q1 in 2d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 
+DEAL::Checking Q1 against Q1 in 2d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 
+DEAL::Checking Q1 against Q1 in 3d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking Q1 against Q1 in 3d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking Q1 against Q2 in 1d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 
+DEAL::Checking Q2 against Q1 in 1d:
+DEAL::1.0000000 2.0000000
+DEAL::1.2247449
+DEAL::0 0 0 0 1.0000000 -0.50000000 
+DEAL::Checking Q1 against Q2 in 2d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 
+DEAL::Checking Q2 against Q1 in 2d:
+DEAL::1.2500000 2.0000000
+DEAL::2.6925824
+DEAL::0 0 0 0 0 0 0 0 1.0000000 1.0000000 1.0000000 -0.50000000 1.0000000 0 1.0000000 0 1.0000000 -0.25000000 
+DEAL::Checking Q1 against Q2 in 3d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking Q2 against Q1 in 3d:
+DEAL::2.3750000 2.0000000
+DEAL::5.1599419
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 1.0000000 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 -0.50000000 1.0000000 -0.25000000 1.0000000 -0.25000000 1.0000000 -0.25000000 1.0000000 -0.25000000 1.0000000 -0.25000000 1.0000000 0 1.0000000 -0.12500000 
+DEAL::Checking Q1 against Q3 in 1d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 
+DEAL::Checking Q3 against Q1 in 1d:
+DEAL::1.0000000 2.0000000
+DEAL::1.7888544
+DEAL::0 0 0 0 1.0000000 -0.27639320 1.0000000 -0.27639320 
+DEAL::Checking Q1 against Q3 in 2d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 
+DEAL::Checking Q3 against Q1 in 2d:
+DEAL::3.0000000 2.0000000
+DEAL::4.2708313
+DEAL::0 0 0 0 0 0 0 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 -0.076393202 1.0000000 -0.076393202 1.0000000 -0.076393202 1.0000000 -0.076393202 
+DEAL::Checking Q1 against Q3 in 3d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking Q3 against Q1 in 3d:
+DEAL::7.0000000 2.0000000
+DEAL::8.9871019
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 -0.021114562 1.0000000 -0.021114562 1.0000000 -0.021114562 1.0000000 -0.021114562 1.0000000 -0.021114562 1.0000000 -0.021114562 1.0000000 -0.021114562 1.0000000 -0.021114562 
+DEAL::Checking Q2 against Q2 in 1d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 
+DEAL::Checking Q2 against Q2 in 1d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 
+DEAL::Checking Q2 against Q2 in 2d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking Q2 against Q2 in 2d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking Q2 against Q2 in 3d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking Q2 against Q2 in 3d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking Q2 against Q3 in 1d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 
+DEAL::Checking Q3 against Q2 in 1d:
+DEAL::1.0000000 1.4472136
+DEAL::1.0954451
+DEAL::0 0 0 0 0.50000000 0.22360680 0.50000000 0.22360680 
+DEAL::Checking Q2 against Q3 in 2d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking Q3 against Q2 in 2d:
+DEAL::1.5000000 2.5944272
+DEAL::2.8705400
+DEAL::0 0 0 0 0 0 0 0 0.50000000 0 0.50000000 0 0.50000000 0 0.50000000 0 0.50000000 0 0.50000000 0 0.50000000 0 0.50000000 0 0.75000000 -0.050000000 0.75000000 -0.050000000 0.75000000 -0.050000000 0.75000000 -0.050000000 
+DEAL::Checking Q2 against Q3 in 3d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking Q3 against Q2 in 3d:
+DEAL::2.0944272 3.7810835
+DEAL::6.5397248
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.50000000 0 0.50000000 0 0.50000000 0 0.50000000 0 0.50000000 0 0.50000000 0 0.50000000 0 0.50000000 0 0.50000000 0 0.50000000 0 0.50000000 0 0.50000000 0 0.50000000 0 0.50000000 0 0.50000000 0 0.50000000 0 0.50000000 0 0.50000000 0 0.50000000 0 0.50000000 0 0.50000000 0 0.50000000 0 0.50000000 0 0.50000000 0 0.75000000 0 0.75000000 0 0.75000000 0.11180340 0.75000000 -0.11180340 0.75000000 -0.11180340 0.75000000 0.11180340 0.75000000 0 0.75000000 0 0.75000000 0 0.75000000 0 0.75000000 0.11180340 0.75000000 -0.11180340 0.75000000 0 0.75000000 0 0.75000000 0 0.75000000 0 0.75000000 0.11180340 0.75000000 0.11180340 0.75000000 0.11180340 0.75000000 0.11180340 0.75000000 0 0.75000000 0 0.75000000 0 0.75000000 0 0.87500000 0.011180340 0.87500000 0.011180340 0.87500000 0.011180340 0.87500000 0.011180340 0.87500000 0.011180340 0.87500000 0.011180340 0.87500000 0.011180340 0.87500000 0.011180340 
+DEAL::Checking Q3 against Q3 in 1d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 
+DEAL::Checking Q3 against Q3 in 1d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 
+DEAL::Checking Q3 against Q3 in 2d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking Q3 against Q3 in 2d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking Q3 against Q3 in 3d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking Q3 against Q3 in 3d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking DGQ0 against DGQ0 in 1d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 
+DEAL::Checking DGQ0 against DGQ0 in 1d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 
+DEAL::Checking DGQ0 against DGQ0 in 2d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 
+DEAL::Checking DGQ0 against DGQ0 in 2d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 
+DEAL::Checking DGQ0 against DGQ0 in 3d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 
+DEAL::Checking DGQ0 against DGQ0 in 3d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 
+DEAL::Checking DGQ0 against DGQ1 in 1d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 
+DEAL::Checking DGQ1 against DGQ0 in 1d:
+DEAL::1.0000000 1.0000000
+DEAL::1.0000000
+DEAL::0.50000000 -0.50000000 0.50000000 -0.50000000 
+DEAL::Checking DGQ0 against DGQ1 in 2d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 
+DEAL::Checking DGQ1 against DGQ0 in 2d:
+DEAL::1.5000000 1.5000000
+DEAL::1.7320508
+DEAL::0.75000000 -0.25000000 0.75000000 -0.25000000 0.75000000 -0.25000000 0.75000000 -0.25000000 
+DEAL::Checking DGQ0 against DGQ1 in 3d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 
+DEAL::Checking DGQ1 against DGQ0 in 3d:
+DEAL::1.7500000 1.7500000
+DEAL::2.6457513
+DEAL::0.87500000 -0.12500000 0.87500000 -0.12500000 0.87500000 -0.12500000 0.87500000 -0.12500000 0.87500000 -0.12500000 0.87500000 -0.12500000 0.87500000 -0.12500000 0.87500000 -0.12500000 
+DEAL::Checking DGQ0 against DGQ2 in 1d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 
+DEAL::Checking DGQ2 against DGQ0 in 1d:
+DEAL::2.0000000 2.0000000
+DEAL::2.0000000
+DEAL::1.0000000 0 0 0 1.0000000 0 
+DEAL::Checking DGQ0 against DGQ2 in 2d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 
+DEAL::Checking DGQ2 against DGQ0 in 2d:
+DEAL::8.0000000 2.0000000
+DEAL::4.0000000
+DEAL::1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 0 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 
+DEAL::Checking DGQ0 against DGQ2 in 3d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 
+DEAL::Checking DGQ2 against DGQ0 in 3d:
+DEAL::26.000000 2.0000000
+DEAL::7.2111026
+DEAL::1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 0 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 
+DEAL::Checking DGQ0 against DGQ4 in 1d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 
+DEAL::Checking DGQ4 against DGQ0 in 1d:
+DEAL::4.0000000 2.0000000
+DEAL::2.8284271
+DEAL::1.0000000 0 1.0000000 0 0 0 1.0000000 0 1.0000000 0 
+DEAL::Checking DGQ0 against DGQ4 in 2d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 
+DEAL::Checking DGQ4 against DGQ0 in 2d:
+DEAL::24.000000 2.0000000
+DEAL::6.9282032
+DEAL::1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 0 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 
+DEAL::Checking DGQ0 against DGQ4 in 3d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 
+DEAL::Checking DGQ4 against DGQ0 in 3d:
+DEAL::124.00000 2.0000000
+DEAL::15.748016
+DEAL::1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 0 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 
+DEAL::Checking DGQ1 against DGQ1 in 1d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 
+DEAL::Checking DGQ1 against DGQ1 in 1d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 
+DEAL::Checking DGQ1 against DGQ1 in 2d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 
+DEAL::Checking DGQ1 against DGQ1 in 2d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 
+DEAL::Checking DGQ1 against DGQ1 in 3d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking DGQ1 against DGQ1 in 3d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking DGQ1 against DGQ3 in 1d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 
+DEAL::Checking DGQ3 against DGQ1 in 1d:
+DEAL::1.0000000 2.0000000
+DEAL::1.7888544
+DEAL::0 0 1.0000000 0 1.0000000 0 0 0 
+DEAL::Checking DGQ1 against DGQ3 in 2d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 
+DEAL::Checking DGQ3 against DGQ1 in 2d:
+DEAL::3.0000000 2.0000000
+DEAL::4.2708313
+DEAL::0 0 1.0000000 0 1.0000000 0 0 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 0 0 1.0000000 0 1.0000000 0 0 0 
+DEAL::Checking DGQ1 against DGQ3 in 3d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking DGQ3 against DGQ1 in 3d:
+DEAL::7.0000000 2.0000000
+DEAL::8.9871019
+DEAL::0 0 1.0000000 0 1.0000000 0 0 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 0 0 1.0000000 0 1.0000000 0 0 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 0 0 1.0000000 0 1.0000000 0 0 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 0 0 1.0000000 0 1.0000000 0 0 0 
+DEAL::Checking DGQ2 against DGQ2 in 1d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 
+DEAL::Checking DGQ2 against DGQ2 in 1d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 
+DEAL::Checking DGQ2 against DGQ2 in 2d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking DGQ2 against DGQ2 in 2d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking DGQ2 against DGQ2 in 3d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking DGQ2 against DGQ2 in 3d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking DGQ2 against DGQ2 in 1d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 
+DEAL::Checking DGQ2 against DGQ2 in 1d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 
+DEAL::Checking DGQ2 against DGQ2 in 2d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking DGQ2 against DGQ2 in 2d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking DGQ2 against DGQ2 in 3d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking DGQ2 against DGQ2 in 3d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking DGQ2 against DGQ4 in 1d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 
+DEAL::Checking DGQ4 against DGQ2 in 1d:
+DEAL::1.1428571 2.2260822
+DEAL::1.8070158
+DEAL::0 0 1.0000000 0 0 0 1.0000000 0 0 0 
+DEAL::Checking DGQ2 against DGQ4 in 2d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking DGQ4 against DGQ2 in 2d:
+DEAL::3.5918367 2.5032777
+DEAL::5.0192466
+DEAL::0 0 1.0000000 0 0 0 1.0000000 0 0 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 0 0 1.0000000 0 0 0 1.0000000 0 0 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 0 0 1.0000000 0 0 0 1.0000000 0 0 0 
+DEAL::Checking DGQ2 against DGQ4 in 3d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking DGQ4 against DGQ2 in 3d:
+DEAL::8.8396501 2.8431420
+DEAL::12.190079
+DEAL::0 0 1.0000000 0 0 0 1.0000000 0 0 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 0 0 1.0000000 0 0 0 1.0000000 0 0 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 0 0 1.0000000 0 0 0 1.0000000 0 0 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 0 0 1.0000000 0 0 0 1.0000000 0 0 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 0 0 1.0000000 0 0 0 1.0000000 0 0 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 0 0 1.0000000 0 0 0 1.0000000 0 0 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 0 0 1.0000000 0 0 0 1.0000000 0 0 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 0 0 1.0000000 0 0 0 1.0000000 0 0 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 0 0 1.0000000 0 0 0 1.0000000 0 0 0 
+DEAL::Checking DGQ3 against DGQ3 in 1d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 
+DEAL::Checking DGQ3 against DGQ3 in 1d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 
+DEAL::Checking DGQ3 against DGQ3 in 2d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking DGQ3 against DGQ3 in 2d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking DGQ3 against DGQ3 in 3d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking DGQ3 against DGQ3 in 3d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking DGP0 against DGP0 in 1d:
+DEAL::Checking DGP0 against DGP0 in 1d:
+DEAL::Checking DGP0 against DGP0 in 2d:
+DEAL::Checking DGP0 against DGP0 in 2d:
+DEAL::Checking DGP0 against DGP0 in 3d:
+DEAL::Checking DGP0 against DGP0 in 3d:
+DEAL::Checking DGP0 against DGP1 in 1d:
+DEAL::Checking DGP1 against DGP0 in 1d:
+DEAL::Checking DGP0 against DGP1 in 2d:
+DEAL::Checking DGP1 against DGP0 in 2d:
+DEAL::Checking DGP0 against DGP1 in 3d:
+DEAL::Checking DGP1 against DGP0 in 3d:
+DEAL::Checking DGP0 against DGP2 in 1d:
+DEAL::Checking DGP2 against DGP0 in 1d:
+DEAL::Checking DGP0 against DGP2 in 2d:
+DEAL::Checking DGP2 against DGP0 in 2d:
+DEAL::Checking DGP0 against DGP2 in 3d:
+DEAL::Checking DGP2 against DGP0 in 3d:
+DEAL::Checking DGP0 against DGP4 in 1d:
+DEAL::Checking DGP4 against DGP0 in 1d:
+DEAL::Checking DGP0 against DGP4 in 2d:
+DEAL::Checking DGP4 against DGP0 in 2d:
+DEAL::Checking DGP0 against DGP4 in 3d:
+DEAL::Checking DGP4 against DGP0 in 3d:
+DEAL::Checking DGP1 against DGP1 in 1d:
+DEAL::Checking DGP1 against DGP1 in 1d:
+DEAL::Checking DGP1 against DGP1 in 2d:
+DEAL::Checking DGP1 against DGP1 in 2d:
+DEAL::Checking DGP1 against DGP1 in 3d:
+DEAL::Checking DGP1 against DGP1 in 3d:
+DEAL::Checking DGP1 against DGP3 in 1d:
+DEAL::Checking DGP3 against DGP1 in 1d:
+DEAL::Checking DGP1 against DGP3 in 2d:
+DEAL::Checking DGP3 against DGP1 in 2d:
+DEAL::Checking DGP1 against DGP3 in 3d:
+DEAL::Checking DGP3 against DGP1 in 3d:
+DEAL::Checking DGP2 against DGP2 in 1d:
+DEAL::Checking DGP2 against DGP2 in 1d:
+DEAL::Checking DGP2 against DGP2 in 2d:
+DEAL::Checking DGP2 against DGP2 in 2d:
+DEAL::Checking DGP2 against DGP2 in 3d:
+DEAL::Checking DGP2 against DGP2 in 3d:
+DEAL::Checking DGP2 against DGP2 in 1d:
+DEAL::Checking DGP2 against DGP2 in 1d:
+DEAL::Checking DGP2 against DGP2 in 2d:
+DEAL::Checking DGP2 against DGP2 in 2d:
+DEAL::Checking DGP2 against DGP2 in 3d:
+DEAL::Checking DGP2 against DGP2 in 3d:
+DEAL::Checking DGP2 against DGP4 in 1d:
+DEAL::Checking DGP4 against DGP2 in 1d:
+DEAL::Checking DGP2 against DGP4 in 2d:
+DEAL::Checking DGP4 against DGP2 in 2d:
+DEAL::Checking DGP2 against DGP4 in 3d:
+DEAL::Checking DGP4 against DGP2 in 3d:
+DEAL::Checking DGP3 against DGP3 in 1d:
+DEAL::Checking DGP3 against DGP3 in 1d:
+DEAL::Checking DGP3 against DGP3 in 2d:
+DEAL::Checking DGP3 against DGP3 in 2d:
+DEAL::Checking DGP3 against DGP3 in 3d:
+DEAL::Checking DGP3 against DGP3 in 3d:
+DEAL::Checking Nedelec0 against Nedelec0 in 2d:
+DEAL::Checking Nedelec0 against Nedelec0 in 2d:
+DEAL::Checking Nedelec0 against Nedelec0 in 3d:
+DEAL::Checking Nedelec0 against Nedelec0 in 3d:
+DEAL::Checking FE_Q<1>(1)3 against FE_Q<1>(2)3 in 1d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking FE_Q<1>(2)3 against FE_Q<1>(1)3 in 1d:
+DEAL::1.0000000 2.0000000
+DEAL::2.1213203
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 1.0000000 0 1.0000000 -0.50000000 1.0000000 0 
+DEAL::Checking FE_DGQ<1>(2)2 against FE_DGQ<1>(3)2 in 1d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking FE_DGQ<1>(3)2 against FE_DGQ<1>(2)2 in 1d:
+DEAL::1.0000000 1.4472136
+DEAL::1.5491933
+DEAL::0 0 0.50000000 0 0.50000000 0 0 0 0 0 0.50000000 0 0.50000000 0 0 0 
+DEAL::Checking FE_DGP<1>(3)1 against FE_DGP<1>(1)1 in 1d:
+DEAL::Checking FE_DGP<1>(1)1 against FE_DGP<1>(3)1 in 1d:
+DEAL::Checking FE_Q<2>(1)3 against FE_Q<2>(2)3 in 2d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking FE_Q<2>(2)3 against FE_Q<2>(1)3 in 2d:
+DEAL::1.2500000 2.0000000
+DEAL::4.6636895
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1.0000000 0 1.0000000 1.0000000 1.0000000 0 1.0000000 0 1.0000000 -0.50000000 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 -0.25000000 1.0000000 0 
+DEAL::Checking FE_DGQ<2>(2)2 against FE_DGQ<2>(3)2 in 2d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking FE_DGQ<2>(3)2 against FE_DGQ<2>(2)2 in 2d:
+DEAL::1.5000000 2.5944272
+DEAL::4.0595566
+DEAL::0 0 0.50000000 0 0.50000000 0 0 0 0.50000000 0 0.75000000 0 0.75000000 0 0.50000000 0 0.50000000 0 0.75000000 0 0.75000000 0 0.50000000 0 0 0 0.50000000 0 0.50000000 0 0 0 0 0 0.50000000 0 0.50000000 0 0 0 0.50000000 0 0.75000000 0 0.75000000 0 0.50000000 0 0.50000000 0 0.75000000 0 0.75000000 0 0.50000000 0 0 0 0.50000000 0 0.50000000 0 0 0 
+DEAL::Checking FE_DGP<2>(3)1 against FE_DGP<2>(1)1 in 2d:
+DEAL::Checking FE_DGP<2>(1)1 against FE_DGP<2>(3)1 in 2d:
+DEAL::Checking FE_Q<3>(1)3 against FE_Q<3>(2)3 in 3d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking FE_Q<3>(2)3 against FE_Q<3>(1)3 in 3d:
+DEAL::2.3750000 2.0000000
+DEAL::8.9372815
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 1.0000000 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 -0.50000000 1.0000000 0 1.0000000 0 1.0000000 -0.25000000 1.0000000 0 1.0000000 0 1.0000000 -0.25000000 1.0000000 0 1.0000000 0 1.0000000 -0.25000000 1.0000000 0 1.0000000 0 1.0000000 -0.25000000 1.0000000 0 1.0000000 0 1.0000000 -0.25000000 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 0 1.0000000 -0.12500000 1.0000000 0 
+DEAL::Checking FE_DGQ<3>(2)2 against FE_DGQ<3>(3)2 in 3d:
+DEAL::0 0
+DEAL::0
+DEAL::0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+DEAL::Checking FE_DGQ<3>(3)2 against FE_DGQ<3>(2)2 in 3d:
+DEAL::2.0944272 3.7810835
+DEAL::9.2485675
+DEAL::0 0 0.50000000 0 0.50000000 0 0 0 0.50000000 0 0.75000000 0 0.75000000 0 0.50000000 0 0.50000000 0 0.75000000 0 0.75000000 0 0.50000000 0 0 0 0.50000000 0 0.50000000 0 0 0 0.50000000 0 0.75000000 0 0.75000000 0 0.50000000 0 0.75000000 0 0.87500000 0 0.87500000 0 0.75000000 0 0.75000000 0 0.87500000 0 0.87500000 0 0.75000000 0 0.50000000 0 0.75000000 0 0.75000000 0 0.50000000 0 0.50000000 0 0.75000000 0 0.75000000 0 0.50000000 0 0.75000000 0 0.87500000 0 0.87500000 0 0.75000000 0 0.75000000 0 0.87500000 0 0.87500000 0 0.75000000 0 0.50000000 0 0.75000000 0 0.75000000 0 0.50000000 0 0 0 0.50000000 0 0.50000000 0 0 0 0.50000000 0 0.75000000 0 0.75000000 0 0.50000000 0 0.50000000 0 0.75000000 0 0.75000000 0 0.50000000 0 0 0 0.50000000 0 0.50000000 0 0 0 0 0 0.50000000 0 0.50000000 0 0 0 0.50000000 0 0.75000000 0 0.75000000 0 0.50000000 0 0.50000000 0 0.75000000 0 0.75000000 0 0.50000000 0 0 0 0.50000000 0 0.50000000 0 0 0 0.50000000 0 0.75000000 0 0.75000000 0 0.50000000 0 0.75000000 0 0.87500000 0 0.87500000 0 0.75000000 0 0.75000000 0 0.87500000 0 0.87500000 0 0.75000000 0 0.50000000 0 0.75000000 0 0.75000000 0 0.50000000 0 0.50000000 0 0.75000000 0 0.75000000 0 0.50000000 0 0.75000000 0 0.87500000 0 0.87500000 0 0.75000000 0 0.75000000 0 0.87500000 0 0.87500000 0 0.75000000 0 0.50000000 0 0.75000000 0 0.75000000 0 0.50000000 0 0 0 0.50000000 0 0.50000000 0 0 0 0.50000000 0 0.75000000 0 0.75000000 0 0.50000000 0 0.50000000 0 0.75000000 0 0.75000000 0 0.50000000 0 0 0 0.50000000 0 0.50000000 0 0 0 
+DEAL::Checking FE_DGP<3>(3)1 against FE_DGP<3>(1)1 in 3d:
+DEAL::Checking FE_DGP<3>(1)1 against FE_DGP<3>(3)1 in 3d: