]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Tests for face interpolation in FE_FaceQ
authorMartin Kronbichler <kronbichler@lnm.mw.tum.de>
Wed, 7 Aug 2013 15:49:45 +0000 (15:49 +0000)
committerMartin Kronbichler <kronbichler@lnm.mw.tum.de>
Wed, 7 Aug 2013 15:49:45 +0000 (15:49 +0000)
git-svn-id: https://svn.dealii.org/trunk@30249 0785d39b-7218-0410-832d-ea1e28bc413d

tests/bits/faceq_face_interpolation.cc [new file with mode: 0644]
tests/bits/faceq_face_interpolation/cmp/generic [new file with mode: 0644]
tests/bits/faceq_subface_interpolation.cc [new file with mode: 0644]
tests/bits/faceq_subface_interpolation/cmp/generic [new file with mode: 0644]

diff --git a/tests/bits/faceq_face_interpolation.cc b/tests/bits/faceq_face_interpolation.cc
new file mode 100644 (file)
index 0000000..fe8f50e
--- /dev/null
@@ -0,0 +1,191 @@
+// ---------------------------------------------------------------------
+// $Id$
+//
+// Copyright (C) 2005 - 2013 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 <deal.II/fe/fe_face.h>
+
+// check
+//   FE_FaceQ::face_interpolation
+
+
+std::string output_file_name = "faceq_face_interpolation/output";
+
+
+void
+output_matrix (const FullMatrix<double> &m)
+{
+  if ((m.m() == 0) || (m.n() == 0))
+    {
+      deallog << "(Empty matrix)" << std::endl;
+      return;
+    }
+
+  deallog << m.l1_norm() << ' ' << m.linfty_norm()
+          << std::endl;
+  if (m.m() == m.n())
+    deallog << m.frobenius_norm() << std::endl;
+
+  for (unsigned int i=0; i<std::min(m.m(),m.n()); ++i)
+    deallog << m(i,i) << ' ' << m(i,std::min(m.m(),m.n())-i-1) << ' ';
+  deallog << std::endl;
+}
+
+
+template <int dim>
+void
+check_this (const FiniteElement<dim> &fe1,
+            const FiniteElement<dim> &fe2)
+{
+  // check all combinations of fe1 and fe2
+  FullMatrix<double> face_constraints;
+  try
+    {
+      face_constraints.reinit (fe1.dofs_per_face,
+                               fe1.dofs_per_face);
+      fe1.get_face_interpolation_matrix (fe1, face_constraints);
+
+      deallog << fe1.get_name()
+              << "  vs.  "
+              << fe1.get_name()
+              << std::endl;
+      output_matrix (face_constraints);
+    }
+  catch (...)
+    {
+    }
+
+  try
+    {
+      face_constraints.reinit (fe2.dofs_per_face,
+                               fe2.dofs_per_face);
+      fe2.get_face_interpolation_matrix (fe2, face_constraints);
+
+      deallog << fe2.get_name()
+              << "  vs.  "
+              << fe2.get_name()
+              << std::endl;
+      output_matrix (face_constraints);
+    }
+  catch (...)
+    {
+    }
+
+  if (fe1.dofs_per_face <= fe2.dofs_per_face)
+    try
+      {
+        face_constraints.reinit (fe2.dofs_per_face,
+                                 fe1.dofs_per_face);
+        fe1.get_face_interpolation_matrix (fe2, face_constraints);
+
+        deallog << fe1.get_name()
+                << "  vs.  "
+                << fe2.get_name()
+                << std::endl;
+        output_matrix (face_constraints);
+      }
+    catch (...)
+      {
+      }
+
+  if (fe2.dofs_per_face <= fe1.dofs_per_face)
+    try
+      {
+        face_constraints.reinit (fe1.dofs_per_face,
+                                 fe2.dofs_per_face);
+        fe2.get_face_interpolation_matrix (fe1, face_constraints);
+
+        deallog << fe2.get_name()
+                << "  vs.  "
+                << fe1.get_name()
+                << std::endl;
+        output_matrix (face_constraints);
+      }
+    catch (...)
+      {
+      }
+}
+
+
+template <int dim>
+void
+check (const unsigned int degree1,
+       const unsigned int degree2)
+{
+  deallog << "Checking FE_FaceQ in " << dim << "d:"
+          << std::endl;
+
+  FE_FaceQ<dim> fe1 (degree1), fe2 (degree2);
+
+  check_this (fe1, fe2);
+}
+
+
+template <int dim>
+void
+check ()
+{
+  check<dim>(0,0);
+  check<dim>(0,1);
+  check<dim>(0,2);
+  check<dim>(0,3);
+  check<dim>(1,2);
+  check<dim>(1,3);
+  check<dim>(2,3);
+}
+
+
+
+int
+main()
+{
+  try
+    {
+      std::ofstream logfile(output_file_name.c_str());
+      deallog << std::setprecision (6);
+      deallog.attach(logfile);
+      deallog.depth_console(0);
+      deallog.threshold_double(1.e-10);
+
+      check<2>();
+      check<3>();
+
+      return 0;
+    }
+  catch (std::exception &exc)
+    {
+      deallog << std::endl << std::endl
+              << "----------------------------------------------------"
+              << std::endl;
+      deallog << "Exception on processing: " << std::endl
+              << exc.what() << std::endl
+              << "Aborting!" << std::endl
+              << "----------------------------------------------------"
+              << std::endl;
+      return 1;
+    }
+  catch (...)
+    {
+      deallog << std::endl << std::endl
+              << "----------------------------------------------------"
+              << std::endl;
+      deallog << "Unknown exception!" << std::endl
+              << "Aborting!" << std::endl
+              << "----------------------------------------------------"
+              << std::endl;
+      return 1;
+    };
+}
diff --git a/tests/bits/faceq_face_interpolation/cmp/generic b/tests/bits/faceq_face_interpolation/cmp/generic
new file mode 100644 (file)
index 0000000..5b2aa67
--- /dev/null
@@ -0,0 +1,179 @@
+
+DEAL::Checking FE_FaceQ in 2d:
+DEAL::FE_FaceQ<2>(0)  vs.  FE_FaceQ<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<2>(0)  vs.  FE_FaceQ<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<2>(0)  vs.  FE_FaceQ<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<2>(0)  vs.  FE_FaceQ<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::Checking FE_FaceQ in 2d:
+DEAL::FE_FaceQ<2>(0)  vs.  FE_FaceQ<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<2>(1)  vs.  FE_FaceQ<2>(1)
+DEAL::1.00000 1.00000
+DEAL::1.41421
+DEAL::1.00000 0 1.00000 0 
+DEAL::FE_FaceQ<2>(0)  vs.  FE_FaceQ<2>(1)
+DEAL::2.00000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::Checking FE_FaceQ in 2d:
+DEAL::FE_FaceQ<2>(0)  vs.  FE_FaceQ<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<2>(2)  vs.  FE_FaceQ<2>(2)
+DEAL::1.00000 1.00000
+DEAL::1.73205
+DEAL::1.00000 0 1.00000 1.00000 1.00000 0 
+DEAL::FE_FaceQ<2>(0)  vs.  FE_FaceQ<2>(2)
+DEAL::3.00000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::Checking FE_FaceQ in 2d:
+DEAL::FE_FaceQ<2>(0)  vs.  FE_FaceQ<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<2>(3)  vs.  FE_FaceQ<2>(3)
+DEAL::1.00000 1.00000
+DEAL::2.00000
+DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 
+DEAL::FE_FaceQ<2>(0)  vs.  FE_FaceQ<2>(3)
+DEAL::4.00000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::Checking FE_FaceQ in 2d:
+DEAL::FE_FaceQ<2>(1)  vs.  FE_FaceQ<2>(1)
+DEAL::1.00000 1.00000
+DEAL::1.41421
+DEAL::1.00000 0 1.00000 0 
+DEAL::FE_FaceQ<2>(2)  vs.  FE_FaceQ<2>(2)
+DEAL::1.00000 1.00000
+DEAL::1.73205
+DEAL::1.00000 0 1.00000 1.00000 1.00000 0 
+DEAL::FE_FaceQ<2>(1)  vs.  FE_FaceQ<2>(2)
+DEAL::1.50000 1.00000
+DEAL::1.00000 0 0.500000 0.500000 
+DEAL::Checking FE_FaceQ in 2d:
+DEAL::FE_FaceQ<2>(1)  vs.  FE_FaceQ<2>(1)
+DEAL::1.00000 1.00000
+DEAL::1.41421
+DEAL::1.00000 0 1.00000 0 
+DEAL::FE_FaceQ<2>(3)  vs.  FE_FaceQ<2>(3)
+DEAL::1.00000 1.00000
+DEAL::2.00000
+DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 
+DEAL::FE_FaceQ<2>(1)  vs.  FE_FaceQ<2>(3)
+DEAL::2.00000 1.00000
+DEAL::1.00000 0 0.276393 0.723607 
+DEAL::Checking FE_FaceQ in 2d:
+DEAL::FE_FaceQ<2>(2)  vs.  FE_FaceQ<2>(2)
+DEAL::1.00000 1.00000
+DEAL::1.73205
+DEAL::1.00000 0 1.00000 1.00000 1.00000 0 
+DEAL::FE_FaceQ<2>(3)  vs.  FE_FaceQ<2>(3)
+DEAL::1.00000 1.00000
+DEAL::2.00000
+DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 
+DEAL::FE_FaceQ<2>(2)  vs.  FE_FaceQ<2>(3)
+DEAL::1.60000 1.24721
+DEAL::1.00000 0 0.800000 0.800000 0.323607 -0.123607 
+DEAL::Checking FE_FaceQ in 3d:
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::Checking FE_FaceQ in 3d:
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(1)  vs.  FE_FaceQ<3>(1)
+DEAL::1.00000 1.00000
+DEAL::2.00000
+DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(1)
+DEAL::4.00000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::Checking FE_FaceQ in 3d:
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(2)  vs.  FE_FaceQ<3>(2)
+DEAL::1.00000 1.00000
+DEAL::3.00000
+DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 1.00000 1.00000 0 1.00000 0 1.00000 0 1.00000 0 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(2)
+DEAL::9.00000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::Checking FE_FaceQ in 3d:
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(3)  vs.  FE_FaceQ<3>(3)
+DEAL::1.00000 1.00000
+DEAL::4.00000
+DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(3)
+DEAL::16.0000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::Checking FE_FaceQ in 3d:
+DEAL::FE_FaceQ<3>(1)  vs.  FE_FaceQ<3>(1)
+DEAL::1.00000 1.00000
+DEAL::2.00000
+DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 
+DEAL::FE_FaceQ<3>(2)  vs.  FE_FaceQ<3>(2)
+DEAL::1.00000 1.00000
+DEAL::3.00000
+DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 1.00000 1.00000 0 1.00000 0 1.00000 0 1.00000 0 
+DEAL::FE_FaceQ<3>(1)  vs.  FE_FaceQ<3>(2)
+DEAL::2.25000 1.00000
+DEAL::1.00000 0 0.500000 0 0 1.00000 0 0.500000 
+DEAL::Checking FE_FaceQ in 3d:
+DEAL::FE_FaceQ<3>(1)  vs.  FE_FaceQ<3>(1)
+DEAL::1.00000 1.00000
+DEAL::2.00000
+DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 
+DEAL::FE_FaceQ<3>(3)  vs.  FE_FaceQ<3>(3)
+DEAL::1.00000 1.00000
+DEAL::4.00000
+DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 
+DEAL::FE_FaceQ<3>(1)  vs.  FE_FaceQ<3>(3)
+DEAL::4.00000 1.00000
+DEAL::1.00000 0 0.276393 0 0 0.723607 0 0 
+DEAL::Checking FE_FaceQ in 3d:
+DEAL::FE_FaceQ<3>(2)  vs.  FE_FaceQ<3>(2)
+DEAL::1.00000 1.00000
+DEAL::3.00000
+DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 1.00000 1.00000 0 1.00000 0 1.00000 0 1.00000 0 
+DEAL::FE_FaceQ<3>(3)  vs.  FE_FaceQ<3>(3)
+DEAL::1.00000 1.00000
+DEAL::4.00000
+DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 
+DEAL::FE_FaceQ<3>(2)  vs.  FE_FaceQ<3>(3)
+DEAL::2.56000 1.55554
+DEAL::1.00000 0 0.800000 0 0.323607 0 0 0 0 0 -0.0988854 0.258885 0.0152786 0.104721 0 0 0 -0.123607 
diff --git a/tests/bits/faceq_subface_interpolation.cc b/tests/bits/faceq_subface_interpolation.cc
new file mode 100644 (file)
index 0000000..c8bfa05
--- /dev/null
@@ -0,0 +1,194 @@
+// ---------------------------------------------------------------------
+// $Id$
+//
+// Copyright (C) 2005 - 2013 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 <deal.II/fe/fe_face.h>
+
+// check
+//   FE_FaceQ::subface_interpolation
+
+
+std::string output_file_name = "faceq_subface_interpolation/output";
+
+
+void
+output_matrix (const FullMatrix<double> &m)
+{
+  if ((m.m() == 0) || (m.n() == 0))
+    {
+      deallog << "(Empty matrix)" << std::endl;
+      return;
+    }
+
+  deallog << m.l1_norm() << ' ' << m.linfty_norm()
+          << std::endl;
+  if (m.m() == m.n())
+    deallog << m.frobenius_norm() << std::endl;
+
+  for (unsigned int i=0; i<std::min(m.m(),m.n()); ++i)
+    deallog << m(i,i) << ' ' << m(i,std::min(m.m(),m.n())-i-1) << ' ';
+  deallog << std::endl;
+}
+
+
+template <int dim>
+void
+check_this (const FiniteElement<dim> &fe1,
+            const FiniteElement<dim> &fe2)
+{
+  // check all combinations of fe1 and fe2
+  for (unsigned int subface=0; subface<GeometryInfo<dim>::max_children_per_face; ++subface)
+    {
+      FullMatrix<double> face_constraints;
+      try
+        {
+          face_constraints.reinit (fe1.dofs_per_face,
+                                   fe1.dofs_per_face);
+          fe1.get_subface_interpolation_matrix (fe1, subface, face_constraints);
+
+          deallog << fe1.get_name()
+                  << "  vs.  "
+                  << fe1.get_name()
+                  << std::endl;
+          output_matrix (face_constraints);
+        }
+      catch (...)
+        {
+        }
+
+      try
+        {
+          face_constraints.reinit (fe2.dofs_per_face,
+                                   fe2.dofs_per_face);
+          fe2.get_subface_interpolation_matrix (fe2, subface, face_constraints);
+
+          deallog << fe2.get_name()
+                  << "  vs.  "
+                  << fe2.get_name()
+                  << std::endl;
+          output_matrix (face_constraints);
+        }
+      catch (...)
+        {
+        }
+
+      if (fe1.dofs_per_face <= fe2.dofs_per_face)
+        try
+          {
+            face_constraints.reinit (fe2.dofs_per_face,
+                                     fe1.dofs_per_face);
+            fe1.get_subface_interpolation_matrix (fe2, subface, face_constraints);
+
+            deallog << fe1.get_name()
+                    << "  vs.  "
+                    << fe2.get_name()
+                    << std::endl;
+            output_matrix (face_constraints);
+          }
+        catch (...)
+          {
+          }
+
+      if (fe2.dofs_per_face <= fe1.dofs_per_face)
+        try
+          {
+            face_constraints.reinit (fe1.dofs_per_face,
+                                     fe2.dofs_per_face);
+            fe2.get_subface_interpolation_matrix (fe1, subface, face_constraints);
+
+            deallog << fe2.get_name()
+                    << "  vs.  "
+                    << fe1.get_name()
+                    << std::endl;
+            output_matrix (face_constraints);
+          }
+        catch (...)
+          {
+          }
+    }
+}
+
+
+template <int dim>
+void
+check (const unsigned int degree1,
+       const unsigned int degree2)
+{
+  deallog << "Checking FE_FaceQ in " << dim << "d:"
+          << std::endl;
+
+  FE_FaceQ<dim> fe1 (degree1), fe2 (degree2);
+
+  check_this (fe1, fe2);
+}
+
+
+template <int dim>
+void
+check ()
+{
+  check<dim>(0,0);
+  check<dim>(0,1);
+  check<dim>(0,2);
+  check<dim>(0,3);
+  check<dim>(1,2);
+  check<dim>(1,3);
+  check<dim>(2,3);
+}
+
+
+
+int
+main()
+{
+  try
+    {
+      std::ofstream logfile(output_file_name.c_str());
+      deallog << std::setprecision (6);
+      deallog.attach(logfile);
+      deallog.depth_console(0);
+      deallog.threshold_double(1.e-10);
+
+      check<2>();
+      check<3>();
+
+      return 0;
+    }
+  catch (std::exception &exc)
+    {
+      deallog << std::endl << std::endl
+              << "----------------------------------------------------"
+              << std::endl;
+      deallog << "Exception on processing: " << std::endl
+              << exc.what() << std::endl
+              << "Aborting!" << std::endl
+              << "----------------------------------------------------"
+              << std::endl;
+      return 1;
+    }
+  catch (...)
+    {
+      deallog << std::endl << std::endl
+              << "----------------------------------------------------"
+              << std::endl;
+      deallog << "Unknown exception!" << std::endl
+              << "Aborting!" << std::endl
+              << "----------------------------------------------------"
+              << std::endl;
+      return 1;
+    };
+}
diff --git a/tests/bits/faceq_subface_interpolation/cmp/generic b/tests/bits/faceq_subface_interpolation/cmp/generic
new file mode 100644 (file)
index 0000000..df01dec
--- /dev/null
@@ -0,0 +1,507 @@
+
+DEAL::Checking FE_FaceQ in 2d:
+DEAL::FE_FaceQ<2>(0)  vs.  FE_FaceQ<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<2>(0)  vs.  FE_FaceQ<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<2>(0)  vs.  FE_FaceQ<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<2>(0)  vs.  FE_FaceQ<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<2>(0)  vs.  FE_FaceQ<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<2>(0)  vs.  FE_FaceQ<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<2>(0)  vs.  FE_FaceQ<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<2>(0)  vs.  FE_FaceQ<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::Checking FE_FaceQ in 2d:
+DEAL::FE_FaceQ<2>(0)  vs.  FE_FaceQ<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<2>(1)  vs.  FE_FaceQ<2>(1)
+DEAL::1.50000 1.00000
+DEAL::1.22474
+DEAL::1.00000 0 0.500000 0.500000 
+DEAL::FE_FaceQ<2>(0)  vs.  FE_FaceQ<2>(1)
+DEAL::2.00000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<2>(0)  vs.  FE_FaceQ<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<2>(1)  vs.  FE_FaceQ<2>(1)
+DEAL::1.50000 1.00000
+DEAL::1.22474
+DEAL::0.500000 0.500000 1.00000 0 
+DEAL::FE_FaceQ<2>(0)  vs.  FE_FaceQ<2>(1)
+DEAL::2.00000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::Checking FE_FaceQ in 2d:
+DEAL::FE_FaceQ<2>(0)  vs.  FE_FaceQ<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<2>(2)  vs.  FE_FaceQ<2>(2)
+DEAL::1.75000 1.25000
+DEAL::1.64886
+DEAL::1.00000 0 0.750000 0.750000 0 0 
+DEAL::FE_FaceQ<2>(0)  vs.  FE_FaceQ<2>(2)
+DEAL::3.00000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<2>(0)  vs.  FE_FaceQ<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<2>(2)  vs.  FE_FaceQ<2>(2)
+DEAL::1.75000 1.25000
+DEAL::1.64886
+DEAL::0 0 0.750000 0.750000 1.00000 0 
+DEAL::FE_FaceQ<2>(0)  vs.  FE_FaceQ<2>(2)
+DEAL::3.00000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::Checking FE_FaceQ in 2d:
+DEAL::FE_FaceQ<2>(0)  vs.  FE_FaceQ<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<2>(3)  vs.  FE_FaceQ<2>(3)
+DEAL::2.33853 1.50000
+DEAL::1.87417
+DEAL::1.00000 0 0.779508 -0.184017 0.220492 0.934017 -0.125000 -0.125000 
+DEAL::FE_FaceQ<2>(0)  vs.  FE_FaceQ<2>(3)
+DEAL::4.00000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<2>(0)  vs.  FE_FaceQ<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<2>(3)  vs.  FE_FaceQ<2>(3)
+DEAL::2.33853 1.50000
+DEAL::1.87417
+DEAL::-0.125000 -0.125000 0.220492 0.934017 0.779508 -0.184017 1.00000 0 
+DEAL::FE_FaceQ<2>(0)  vs.  FE_FaceQ<2>(3)
+DEAL::4.00000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::Checking FE_FaceQ in 2d:
+DEAL::FE_FaceQ<2>(1)  vs.  FE_FaceQ<2>(1)
+DEAL::1.50000 1.00000
+DEAL::1.22474
+DEAL::1.00000 0 0.500000 0.500000 
+DEAL::FE_FaceQ<2>(2)  vs.  FE_FaceQ<2>(2)
+DEAL::1.75000 1.25000
+DEAL::1.64886
+DEAL::1.00000 0 0.750000 0.750000 0 0 
+DEAL::FE_FaceQ<2>(1)  vs.  FE_FaceQ<2>(2)
+DEAL::2.25000 1.00000
+DEAL::1.00000 0 0.250000 0.750000 
+DEAL::FE_FaceQ<2>(1)  vs.  FE_FaceQ<2>(1)
+DEAL::1.50000 1.00000
+DEAL::1.22474
+DEAL::0.500000 0.500000 1.00000 0 
+DEAL::FE_FaceQ<2>(2)  vs.  FE_FaceQ<2>(2)
+DEAL::1.75000 1.25000
+DEAL::1.64886
+DEAL::0 0 0.750000 0.750000 1.00000 0 
+DEAL::FE_FaceQ<2>(1)  vs.  FE_FaceQ<2>(2)
+DEAL::2.25000 1.00000
+DEAL::0.500000 0.500000 0.750000 0.250000 
+DEAL::Checking FE_FaceQ in 2d:
+DEAL::FE_FaceQ<2>(1)  vs.  FE_FaceQ<2>(1)
+DEAL::1.50000 1.00000
+DEAL::1.22474
+DEAL::1.00000 0 0.500000 0.500000 
+DEAL::FE_FaceQ<2>(3)  vs.  FE_FaceQ<2>(3)
+DEAL::2.33853 1.50000
+DEAL::1.87417
+DEAL::1.00000 0 0.779508 -0.184017 0.220492 0.934017 -0.125000 -0.125000 
+DEAL::FE_FaceQ<2>(1)  vs.  FE_FaceQ<2>(3)
+DEAL::3.00000 1.00000
+DEAL::1.00000 0 0.138197 0.861803 
+DEAL::FE_FaceQ<2>(1)  vs.  FE_FaceQ<2>(1)
+DEAL::1.50000 1.00000
+DEAL::1.22474
+DEAL::0.500000 0.500000 1.00000 0 
+DEAL::FE_FaceQ<2>(3)  vs.  FE_FaceQ<2>(3)
+DEAL::2.33853 1.50000
+DEAL::1.87417
+DEAL::-0.125000 -0.125000 0.220492 0.934017 0.779508 -0.184017 1.00000 0 
+DEAL::FE_FaceQ<2>(1)  vs.  FE_FaceQ<2>(3)
+DEAL::3.00000 1.00000
+DEAL::0.500000 0.500000 0.638197 0.361803 
+DEAL::Checking FE_FaceQ in 2d:
+DEAL::FE_FaceQ<2>(2)  vs.  FE_FaceQ<2>(2)
+DEAL::1.75000 1.25000
+DEAL::1.64886
+DEAL::1.00000 0 0.750000 0.750000 0 0 
+DEAL::FE_FaceQ<2>(3)  vs.  FE_FaceQ<2>(3)
+DEAL::2.33853 1.50000
+DEAL::1.87417
+DEAL::1.00000 0 0.779508 -0.184017 0.220492 0.934017 -0.125000 -0.125000 
+DEAL::FE_FaceQ<2>(2)  vs.  FE_FaceQ<2>(3)
+DEAL::2.40000 1.20000
+DEAL::1.00000 0 0.476393 0.476393 -0.100000 0.176393 
+DEAL::FE_FaceQ<2>(2)  vs.  FE_FaceQ<2>(2)
+DEAL::1.75000 1.25000
+DEAL::1.64886
+DEAL::0 0 0.750000 0.750000 1.00000 0 
+DEAL::FE_FaceQ<2>(3)  vs.  FE_FaceQ<2>(3)
+DEAL::2.33853 1.50000
+DEAL::1.87417
+DEAL::-0.125000 -0.125000 0.220492 0.934017 0.779508 -0.184017 1.00000 0 
+DEAL::FE_FaceQ<2>(2)  vs.  FE_FaceQ<2>(3)
+DEAL::2.40000 1.20000
+DEAL::0 0 0.923607 0.923607 0.623607 -0.100000 
+DEAL::Checking FE_FaceQ in 3d:
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::Checking FE_FaceQ in 3d:
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(1)  vs.  FE_FaceQ<3>(1)
+DEAL::2.25000 1.00000
+DEAL::1.50000
+DEAL::1.00000 0 0.500000 0 0.500000 0 0.250000 0.250000 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(1)
+DEAL::4.00000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(1)  vs.  FE_FaceQ<3>(1)
+DEAL::2.25000 1.00000
+DEAL::1.50000
+DEAL::0.500000 0 1.00000 0 0.250000 0.250000 0.500000 0 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(1)
+DEAL::4.00000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(1)  vs.  FE_FaceQ<3>(1)
+DEAL::2.25000 1.00000
+DEAL::1.50000
+DEAL::0.500000 0 0.250000 0.250000 1.00000 0 0.500000 0 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(1)
+DEAL::4.00000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(1)  vs.  FE_FaceQ<3>(1)
+DEAL::2.25000 1.00000
+DEAL::1.50000
+DEAL::0.250000 0.250000 0.500000 0 0.500000 0 1.00000 0 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(1)
+DEAL::4.00000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::Checking FE_FaceQ in 3d:
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(2)  vs.  FE_FaceQ<3>(2)
+DEAL::3.06250 1.56250
+DEAL::2.71875
+DEAL::1.00000 0 0.750000 0 0 0 0.750000 0 0.562500 0.562500 0 0 0 0 0 0 0 0 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(2)
+DEAL::9.00000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(2)  vs.  FE_FaceQ<3>(2)
+DEAL::3.06250 1.56250
+DEAL::2.71875
+DEAL::0 0 0.750000 0 1.00000 0 0 0 0.562500 0.562500 0.750000 0 0 0 0 0 0 0 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(2)
+DEAL::9.00000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(2)  vs.  FE_FaceQ<3>(2)
+DEAL::3.06250 1.56250
+DEAL::2.71875
+DEAL::0 0 0 0 0 0 0.750000 0 0.562500 0.562500 0 0 1.00000 0 0.750000 0 0 0 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(2)
+DEAL::9.00000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(2)  vs.  FE_FaceQ<3>(2)
+DEAL::3.06250 1.56250
+DEAL::2.71875
+DEAL::0 0 0 0 0 0 0 0 0.562500 0.562500 0.750000 0 0 0 0.750000 0 1.00000 0 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(2)
+DEAL::9.00000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::Checking FE_FaceQ in 3d:
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(3)  vs.  FE_FaceQ<3>(3)
+DEAL::5.46870 2.25000
+DEAL::3.51250
+DEAL::1.00000 0 0.779508 0 0.220492 0 -0.125000 0 0.779508 0 0.607633 0.0338623 0.171875 -0.171875 -0.0974386 0.0230021 0.220492 0 0.171875 -0.171875 0.0486165 0.872388 -0.0275614 -0.116752 -0.125000 0 -0.0974386 0.0230021 -0.0275614 -0.116752 0.0156250 0.0156250 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(3)
+DEAL::16.0000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(3)  vs.  FE_FaceQ<3>(3)
+DEAL::5.46870 2.25000
+DEAL::3.51250
+DEAL::-0.125000 0 0.220492 0 0.779508 0 1.00000 0 -0.0974386 0.0230021 0.171875 -0.171875 0.607633 0.0338623 0.779508 0 -0.0275614 -0.116752 0.0486165 0.872388 0.171875 -0.171875 0.220492 0 0.0156250 0.0156250 -0.0275614 -0.116752 -0.0974386 0.0230021 -0.125000 0 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(3)
+DEAL::16.0000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(3)  vs.  FE_FaceQ<3>(3)
+DEAL::5.46870 2.25000
+DEAL::3.51250
+DEAL::-0.125000 0 -0.0974386 0.0230021 -0.0275614 -0.116752 0.0156250 0.0156250 0.220492 0 0.171875 -0.171875 0.0486165 0.872388 -0.0275614 -0.116752 0.779508 0 0.607633 0.0338623 0.171875 -0.171875 -0.0974386 0.0230021 1.00000 0 0.779508 0 0.220492 0 -0.125000 0 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(3)
+DEAL::16.0000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceQ<3>(3)  vs.  FE_FaceQ<3>(3)
+DEAL::5.46870 2.25000
+DEAL::3.51250
+DEAL::0.0156250 0.0156250 -0.0275614 -0.116752 -0.0974386 0.0230021 -0.125000 0 -0.0275614 -0.116752 0.0486165 0.872388 0.171875 -0.171875 0.220492 0 -0.0974386 0.0230021 0.171875 -0.171875 0.607633 0.0338623 0.779508 0 -0.125000 0 0.220492 0 0.779508 0 1.00000 0 
+DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(3)
+DEAL::16.0000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::Checking FE_FaceQ in 3d:
+DEAL::FE_FaceQ<3>(1)  vs.  FE_FaceQ<3>(1)
+DEAL::2.25000 1.00000
+DEAL::1.50000
+DEAL::1.00000 0 0.500000 0 0.500000 0 0.250000 0.250000 
+DEAL::FE_FaceQ<3>(2)  vs.  FE_FaceQ<3>(2)
+DEAL::3.06250 1.56250
+DEAL::2.71875
+DEAL::1.00000 0 0.750000 0 0 0 0.750000 0 0.562500 0.562500 0 0 0 0 0 0 0 0 
+DEAL::FE_FaceQ<3>(1)  vs.  FE_FaceQ<3>(2)
+DEAL::5.06250 1.00000
+DEAL::1.00000 0 0.250000 0 0 0.500000 0 0.750000 
+DEAL::FE_FaceQ<3>(1)  vs.  FE_FaceQ<3>(1)
+DEAL::2.25000 1.00000
+DEAL::1.50000
+DEAL::0.500000 0 1.00000 0 0.250000 0.250000 0.500000 0 
+DEAL::FE_FaceQ<3>(2)  vs.  FE_FaceQ<3>(2)
+DEAL::3.06250 1.56250
+DEAL::2.71875
+DEAL::0 0 0.750000 0 1.00000 0 0 0 0.562500 0.562500 0.750000 0 0 0 0 0 0 0 
+DEAL::FE_FaceQ<3>(1)  vs.  FE_FaceQ<3>(2)
+DEAL::5.06250 1.00000
+DEAL::0.500000 0 0.750000 0 0 1.00000 0.125000 0.375000 
+DEAL::FE_FaceQ<3>(1)  vs.  FE_FaceQ<3>(1)
+DEAL::2.25000 1.00000
+DEAL::1.50000
+DEAL::0.500000 0 0.250000 0.250000 1.00000 0 0.500000 0 
+DEAL::FE_FaceQ<3>(2)  vs.  FE_FaceQ<3>(2)
+DEAL::3.06250 1.56250
+DEAL::2.71875
+DEAL::0 0 0 0 0 0 0.750000 0 0.562500 0.562500 0 0 1.00000 0 0.750000 0 0 0 
+DEAL::FE_FaceQ<3>(1)  vs.  FE_FaceQ<3>(2)
+DEAL::5.06250 1.00000
+DEAL::0.500000 0 0.125000 0.375000 0.250000 0.250000 0 0.250000 
+DEAL::FE_FaceQ<3>(1)  vs.  FE_FaceQ<3>(1)
+DEAL::2.25000 1.00000
+DEAL::1.50000
+DEAL::0.250000 0.250000 0.500000 0 0.500000 0 1.00000 0 
+DEAL::FE_FaceQ<3>(2)  vs.  FE_FaceQ<3>(2)
+DEAL::3.06250 1.56250
+DEAL::2.71875
+DEAL::0 0 0 0 0 0 0 0 0.562500 0.562500 0.750000 0 0 0 0.750000 0 1.00000 0 
+DEAL::FE_FaceQ<3>(1)  vs.  FE_FaceQ<3>(2)
+DEAL::5.06250 1.00000
+DEAL::0.250000 0.250000 0.375000 0.125000 0 0.500000 0.375000 0.125000 
+DEAL::Checking FE_FaceQ in 3d:
+DEAL::FE_FaceQ<3>(1)  vs.  FE_FaceQ<3>(1)
+DEAL::2.25000 1.00000
+DEAL::1.50000
+DEAL::1.00000 0 0.500000 0 0.500000 0 0.250000 0.250000 
+DEAL::FE_FaceQ<3>(3)  vs.  FE_FaceQ<3>(3)
+DEAL::5.46870 2.25000
+DEAL::3.51250
+DEAL::1.00000 0 0.779508 0 0.220492 0 -0.125000 0 0.779508 0 0.607633 0.0338623 0.171875 -0.171875 -0.0974386 0.0230021 0.220492 0 0.171875 -0.171875 0.0486165 0.872388 -0.0275614 -0.116752 -0.125000 0 -0.0974386 0.0230021 -0.0275614 -0.116752 0.0156250 0.0156250 
+DEAL::FE_FaceQ<3>(1)  vs.  FE_FaceQ<3>(3)
+DEAL::9.00000 1.00000
+DEAL::1.00000 0 0.138197 0 0 0.361803 0 0.500000 
+DEAL::FE_FaceQ<3>(1)  vs.  FE_FaceQ<3>(1)
+DEAL::2.25000 1.00000
+DEAL::1.50000
+DEAL::0.500000 0 1.00000 0 0.250000 0.250000 0.500000 0 
+DEAL::FE_FaceQ<3>(3)  vs.  FE_FaceQ<3>(3)
+DEAL::5.46870 2.25000
+DEAL::3.51250
+DEAL::-0.125000 0 0.220492 0 0.779508 0 1.00000 0 -0.0974386 0.0230021 0.171875 -0.171875 0.607633 0.0338623 0.779508 0 -0.0275614 -0.116752 0.0486165 0.872388 0.171875 -0.171875 0.220492 0 0.0156250 0.0156250 -0.0275614 -0.116752 -0.0974386 0.0230021 -0.125000 0 
+DEAL::FE_FaceQ<3>(1)  vs.  FE_FaceQ<3>(3)
+DEAL::9.00000 1.00000
+DEAL::0.500000 0 0.638197 0 0 0.861803 0 0 
+DEAL::FE_FaceQ<3>(1)  vs.  FE_FaceQ<3>(1)
+DEAL::2.25000 1.00000
+DEAL::1.50000
+DEAL::0.500000 0 0.250000 0.250000 1.00000 0 0.500000 0 
+DEAL::FE_FaceQ<3>(3)  vs.  FE_FaceQ<3>(3)
+DEAL::5.46870 2.25000
+DEAL::3.51250
+DEAL::-0.125000 0 -0.0974386 0.0230021 -0.0275614 -0.116752 0.0156250 0.0156250 0.220492 0 0.171875 -0.171875 0.0486165 0.872388 -0.0275614 -0.116752 0.779508 0 0.607633 0.0338623 0.171875 -0.171875 -0.0974386 0.0230021 1.00000 0 0.779508 0 0.220492 0 -0.125000 0 
+DEAL::FE_FaceQ<3>(1)  vs.  FE_FaceQ<3>(3)
+DEAL::9.00000 1.00000
+DEAL::0.500000 0 0.0690983 0.430902 0.319098 0.180902 0.250000 0.250000 
+DEAL::FE_FaceQ<3>(1)  vs.  FE_FaceQ<3>(1)
+DEAL::2.25000 1.00000
+DEAL::1.50000
+DEAL::0.250000 0.250000 0.500000 0 0.500000 0 1.00000 0 
+DEAL::FE_FaceQ<3>(3)  vs.  FE_FaceQ<3>(3)
+DEAL::5.46870 2.25000
+DEAL::3.51250
+DEAL::0.0156250 0.0156250 -0.0275614 -0.116752 -0.0974386 0.0230021 -0.125000 0 -0.0275614 -0.116752 0.0486165 0.872388 0.171875 -0.171875 0.220492 0 -0.0974386 0.0230021 0.171875 -0.171875 0.607633 0.0338623 0.779508 0 -0.125000 0 0.220492 0 0.779508 0 1.00000 0 
+DEAL::FE_FaceQ<3>(1)  vs.  FE_FaceQ<3>(3)
+DEAL::9.00000 1.00000
+DEAL::0.250000 0.250000 0.319098 0.180902 0.0690983 0.430902 0.500000 0 
+DEAL::Checking FE_FaceQ in 3d:
+DEAL::FE_FaceQ<3>(2)  vs.  FE_FaceQ<3>(2)
+DEAL::3.06250 1.56250
+DEAL::2.71875
+DEAL::1.00000 0 0.750000 0 0 0 0.750000 0 0.562500 0.562500 0 0 0 0 0 0 0 0 
+DEAL::FE_FaceQ<3>(3)  vs.  FE_FaceQ<3>(3)
+DEAL::5.46870 2.25000
+DEAL::3.51250
+DEAL::1.00000 0 0.779508 0 0.220492 0 -0.125000 0 0.779508 0 0.607633 0.0338623 0.171875 -0.171875 -0.0974386 0.0230021 0.220492 0 0.171875 -0.171875 0.0486165 0.872388 -0.0275614 -0.116752 -0.125000 0 -0.0974386 0.0230021 -0.0275614 -0.116752 0.0156250 0.0156250 
+DEAL::FE_FaceQ<3>(2)  vs.  FE_FaceQ<3>(3)
+DEAL::5.76000 1.44000
+DEAL::1.00000 0 0.476393 0 -0.100000 0 0 0 0 0 -0.0476393 0.297082 -0.0176393 -0.0623607 -0.100000 0.623607 0 0.176393 
+DEAL::FE_FaceQ<3>(2)  vs.  FE_FaceQ<3>(2)
+DEAL::3.06250 1.56250
+DEAL::2.71875
+DEAL::0 0 0.750000 0 1.00000 0 0 0 0.562500 0.562500 0.750000 0 0 0 0 0 0 0 
+DEAL::FE_FaceQ<3>(3)  vs.  FE_FaceQ<3>(3)
+DEAL::5.46870 2.25000
+DEAL::3.51250
+DEAL::-0.125000 0 0.220492 0 0.779508 0 1.00000 0 -0.0974386 0.0230021 0.171875 -0.171875 0.607633 0.0338623 0.779508 0 -0.0275614 -0.116752 0.0486165 0.872388 0.171875 -0.171875 0.220492 0 0.0156250 0.0156250 -0.0275614 -0.116752 -0.0974386 0.0230021 -0.125000 0 
+DEAL::FE_FaceQ<3>(2)  vs.  FE_FaceQ<3>(3)
+DEAL::5.76000 1.44000
+DEAL::0 0 0.923607 0 0.623607 0 0 0 0.476393 0.476393 0.0840325 -0.0476393 0.0100000 0.388885 0 0 0 0 
+DEAL::FE_FaceQ<3>(2)  vs.  FE_FaceQ<3>(2)
+DEAL::3.06250 1.56250
+DEAL::2.71875
+DEAL::0 0 0 0 0 0 0.750000 0 0.562500 0.562500 0 0 1.00000 0 0.750000 0 0 0 
+DEAL::FE_FaceQ<3>(3)  vs.  FE_FaceQ<3>(3)
+DEAL::5.46870 2.25000
+DEAL::3.51250
+DEAL::-0.125000 0 -0.0974386 0.0230021 -0.0275614 -0.116752 0.0156250 0.0156250 0.220492 0 0.171875 -0.171875 0.0486165 0.872388 -0.0275614 -0.116752 0.779508 0 0.607633 0.0338623 0.171875 -0.171875 -0.0974386 0.0230021 1.00000 0 0.779508 0 0.220492 0 -0.125000 0 
+DEAL::FE_FaceQ<3>(2)  vs.  FE_FaceQ<3>(3)
+DEAL::5.76000 1.44000
+DEAL::0 0 0 0 0 0 0 0 0 0 -0.0923607 0.575967 0.0311146 0.0100000 0.176393 -0.100000 0 -0.100000 
+DEAL::FE_FaceQ<3>(2)  vs.  FE_FaceQ<3>(2)
+DEAL::3.06250 1.56250
+DEAL::2.71875
+DEAL::0 0 0 0 0 0 0 0 0.562500 0.562500 0.750000 0 0 0 0.750000 0 1.00000 0 
+DEAL::FE_FaceQ<3>(3)  vs.  FE_FaceQ<3>(3)
+DEAL::5.46870 2.25000
+DEAL::3.51250
+DEAL::0.0156250 0.0156250 -0.0275614 -0.116752 -0.0974386 0.0230021 -0.125000 0 -0.0275614 -0.116752 0.0486165 0.872388 0.171875 -0.171875 0.220492 0 -0.0974386 0.0230021 0.171875 -0.171875 0.607633 0.0338623 0.779508 0 -0.125000 0 0.220492 0 0.779508 0 1.00000 0 
+DEAL::FE_FaceQ<3>(2)  vs.  FE_FaceQ<3>(3)
+DEAL::5.76000 1.44000
+DEAL::0 0 0 0 0 0 0 1.00000 0.923607 0.923607 0.162918 -0.0923607 -0.0176393 -0.0623607 0 0 0 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.