From f578ce96510b36c7d6feefd54ef6b36c3854e80d Mon Sep 17 00:00:00 2001
From: kronbichler <kronbichler@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Tue, 17 Sep 2013 14:24:00 +0000
Subject: [PATCH] Fix subface interpolation.

git-svn-id: https://svn.dealii.org/trunk@30765 0785d39b-7218-0410-832d-ea1e28bc413d
---
 deal.II/source/fe/fe_face.cc                  |  28 +-
 tests/bits/faceq_face_interpolation.cc        |  20 +-
 .../bits/faceq_face_interpolation/cmp/generic | 178 ++++++
 tests/bits/faceq_subface_interpolation.cc     |  20 +-
 .../faceq_subface_interpolation/cmp/generic   | 506 ++++++++++++++++++
 5 files changed, 722 insertions(+), 30 deletions(-)

diff --git a/deal.II/source/fe/fe_face.cc b/deal.II/source/fe/fe_face.cc
index a362e3cdec..641e20a684 100644
--- a/deal.II/source/fe/fe_face.cc
+++ b/deal.II/source/fe/fe_face.cc
@@ -436,7 +436,7 @@ get_subface_interpolation_matrix (const FiniteElement<dim,spacedim> &x_source_fe
       // zero.
       const double eps = 2e-13*(this->degree+1)*(dim-1);
 
-      FullMatrix<double> mass (face_quadrature.size(), this->dofs_per_face);
+      FullMatrix<double> mass (face_quadrature.size(), source_fe->dofs_per_face);
  
       for (unsigned int k = 0; k < face_quadrature.size(); ++k)
         {
@@ -446,18 +446,18 @@ get_subface_interpolation_matrix (const FiniteElement<dim,spacedim> &x_source_fe
             GeometryInfo<dim-1>::child_to_cell_coordinates (face_quadrature.point(k),
                                                             subface);
 
-          for (unsigned int j = 0; j < this->dofs_per_face; ++j)
-            mass (k , j) = this->poly_space.compute_value(j, p);
+          for (unsigned int j = 0; j < source_fe->dofs_per_face; ++j)
+            mass (k , j) = source_fe->poly_space.compute_value(j, p);
         }
 
       Householder<double> H(mass);
       Vector<double> v_in(face_quadrature.size());
-      Vector<double> v_out(this->dofs_per_face);
+      Vector<double> v_out(source_fe->dofs_per_face);
 
 
       // compute the interpolation matrix by evaluating on the fine side and
       // then solving the least squares problem
-      for (unsigned int i=0; i<source_fe->dofs_per_face; ++i)
+      for (unsigned int i=0; i<this->dofs_per_face; ++i)
         {
           for (unsigned int k = 0; k < face_quadrature.size(); ++k)
             {
@@ -465,12 +465,12 @@ get_subface_interpolation_matrix (const FiniteElement<dim,spacedim> &x_source_fe
                 face_quadrature.point(k) :
                 GeometryInfo<dim-1>::child_to_cell_coordinates (face_quadrature.point(k),
                                                                 subface);
-              v_in(k) = source_fe->poly_space.compute_value(i, p);
+              v_in(k) = this->poly_space.compute_value(i, p);
             }
           const double result = H.least_squares(v_out, v_in);
           Assert(result < 1e-12, FETools::ExcLeastSquaresError (result));
 
-          for (unsigned int j = 0; j < this->dofs_per_face; ++j)
+          for (unsigned int j = 0; j < source_fe->dofs_per_face; ++j)
             {
               double matrix_entry = v_out(j);
 
@@ -482,21 +482,9 @@ get_subface_interpolation_matrix (const FiniteElement<dim,spacedim> &x_source_fe
               if (std::fabs (matrix_entry) < eps)
                 matrix_entry = 0.0;
 
-              interpolation_matrix(i,j) = matrix_entry;
+              interpolation_matrix(j,i) = matrix_entry;
             }
         }
-
-      // make sure that the row sum of each of the matrices is 1 at this
-      // point. this must be so since the shape functions sum up to 1
-      for (unsigned int j=0; j<source_fe->dofs_per_face; ++j)
-        {
-          double sum = 0.;
-
-          for (unsigned int i=0; i<this->dofs_per_face; ++i)
-            sum += interpolation_matrix(j,i);
-
-          Assert (std::fabs(sum-1) < eps, ExcInternalError());
-        }
     }
   else if (dynamic_cast<const FE_Nothing<dim> *>(&x_source_fe) != 0)
     {
diff --git a/tests/bits/faceq_face_interpolation.cc b/tests/bits/faceq_face_interpolation.cc
index fe8f50e6dd..ca0d906738 100644
--- a/tests/bits/faceq_face_interpolation.cc
+++ b/tests/bits/faceq_face_interpolation.cc
@@ -19,7 +19,7 @@
 #include <deal.II/fe/fe_face.h>
 
 // check
-//   FE_FaceQ::face_interpolation
+//   FE_FaceQ/P::face_interpolation
 
 
 std::string output_file_name = "faceq_face_interpolation/output";
@@ -125,12 +125,22 @@ void
 check (const unsigned int degree1,
        const unsigned int degree2)
 {
-  deallog << "Checking FE_FaceQ in " << dim << "d:"
-          << std::endl;
+  {
+    deallog << "Checking FE_FaceQ in " << dim << "d:"
+            << std::endl;
+
+    FE_FaceQ<dim> fe1 (degree1), fe2 (degree2);
+
+    check_this (fe1, fe2);
+  }
+  {
+    deallog << "Checking FE_FaceP in " << dim << "d:"
+            << std::endl;
 
-  FE_FaceQ<dim> fe1 (degree1), fe2 (degree2);
+    FE_FaceP<dim> fe1 (degree1), fe2 (degree2);
 
-  check_this (fe1, fe2);
+    check_this (fe1, fe2);
+  }
 }
 
 
diff --git a/tests/bits/faceq_face_interpolation/cmp/generic b/tests/bits/faceq_face_interpolation/cmp/generic
index 5b2aa67ac2..f477f4dd29 100644
--- a/tests/bits/faceq_face_interpolation/cmp/generic
+++ b/tests/bits/faceq_face_interpolation/cmp/generic
@@ -16,6 +16,23 @@ 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_FaceP in 2d:
+DEAL::FE_FaceP<2>(0)  vs.  FE_FaceP<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<2>(0)  vs.  FE_FaceP<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<2>(0)  vs.  FE_FaceP<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<2>(0)  vs.  FE_FaceP<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
@@ -28,6 +45,18 @@ 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_FaceP in 2d:
+DEAL::FE_FaceP<2>(0)  vs.  FE_FaceP<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<2>(1)  vs.  FE_FaceP<2>(1)
+DEAL::1.00000 1.00000
+DEAL::1.41421
+DEAL::1.00000 0 1.00000 0 
+DEAL::FE_FaceP<2>(0)  vs.  FE_FaceP<2>(1)
+DEAL::1.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
@@ -40,6 +69,18 @@ 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_FaceP in 2d:
+DEAL::FE_FaceP<2>(0)  vs.  FE_FaceP<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<2>(2)  vs.  FE_FaceP<2>(2)
+DEAL::1.00000 1.00000
+DEAL::1.73205
+DEAL::1.00000 0 1.00000 1.00000 1.00000 0 
+DEAL::FE_FaceP<2>(0)  vs.  FE_FaceP<2>(2)
+DEAL::1.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
@@ -52,6 +93,18 @@ 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_FaceP in 2d:
+DEAL::FE_FaceP<2>(0)  vs.  FE_FaceP<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<2>(3)  vs.  FE_FaceP<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_FaceP<2>(0)  vs.  FE_FaceP<2>(3)
+DEAL::1.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
@@ -64,6 +117,18 @@ 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_FaceP in 2d:
+DEAL::FE_FaceP<2>(1)  vs.  FE_FaceP<2>(1)
+DEAL::1.00000 1.00000
+DEAL::1.41421
+DEAL::1.00000 0 1.00000 0 
+DEAL::FE_FaceP<2>(2)  vs.  FE_FaceP<2>(2)
+DEAL::1.00000 1.00000
+DEAL::1.73205
+DEAL::1.00000 0 1.00000 1.00000 1.00000 0 
+DEAL::FE_FaceP<2>(1)  vs.  FE_FaceP<2>(2)
+DEAL::1.00000 1.00000
+DEAL::1.00000 0 1.00000 0 
 DEAL::Checking FE_FaceQ in 2d:
 DEAL::FE_FaceQ<2>(1)  vs.  FE_FaceQ<2>(1)
 DEAL::1.00000 1.00000
@@ -76,6 +141,18 @@ 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_FaceP in 2d:
+DEAL::FE_FaceP<2>(1)  vs.  FE_FaceP<2>(1)
+DEAL::1.00000 1.00000
+DEAL::1.41421
+DEAL::1.00000 0 1.00000 0 
+DEAL::FE_FaceP<2>(3)  vs.  FE_FaceP<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_FaceP<2>(1)  vs.  FE_FaceP<2>(3)
+DEAL::1.00000 1.00000
+DEAL::1.00000 0 1.00000 0 
 DEAL::Checking FE_FaceQ in 2d:
 DEAL::FE_FaceQ<2>(2)  vs.  FE_FaceQ<2>(2)
 DEAL::1.00000 1.00000
@@ -88,6 +165,18 @@ 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_FaceP in 2d:
+DEAL::FE_FaceP<2>(2)  vs.  FE_FaceP<2>(2)
+DEAL::1.00000 1.00000
+DEAL::1.73205
+DEAL::1.00000 0 1.00000 1.00000 1.00000 0 
+DEAL::FE_FaceP<2>(3)  vs.  FE_FaceP<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_FaceP<2>(2)  vs.  FE_FaceP<2>(3)
+DEAL::1.00000 1.00000
+DEAL::1.00000 0 1.00000 1.00000 1.00000 0 
 DEAL::Checking FE_FaceQ in 3d:
 DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
 DEAL::1.00000 1.00000
@@ -105,6 +194,23 @@ 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_FaceP in 3d:
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<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
@@ -117,6 +223,18 @@ 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_FaceP in 3d:
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(1)  vs.  FE_FaceP<3>(1)
+DEAL::1.00000 1.00000
+DEAL::1.73205
+DEAL::1.00000 0 1.00000 1.00000 1.00000 0 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(1)
+DEAL::1.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
@@ -129,6 +247,18 @@ DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 1.00000 1.00000 0 1.00000
 DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(2)
 DEAL::9.00000 1.00000
 DEAL::1.00000 1.00000 
+DEAL::Checking FE_FaceP in 3d:
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(2)  vs.  FE_FaceP<3>(2)
+DEAL::1.00000 1.00000
+DEAL::2.44949
+DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(2)
+DEAL::1.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
@@ -141,6 +271,18 @@ DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00
 DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(3)
 DEAL::16.0000 1.00000
 DEAL::1.00000 1.00000 
+DEAL::Checking FE_FaceP in 3d:
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(3)  vs.  FE_FaceP<3>(3)
+DEAL::1.00000 1.00000
+DEAL::3.16228
+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 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(3)
+DEAL::1.00000 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
@@ -153,6 +295,18 @@ DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 1.00000 1.00000 0 1.00000
 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_FaceP in 3d:
+DEAL::FE_FaceP<3>(1)  vs.  FE_FaceP<3>(1)
+DEAL::1.00000 1.00000
+DEAL::1.73205
+DEAL::1.00000 0 1.00000 1.00000 1.00000 0 
+DEAL::FE_FaceP<3>(2)  vs.  FE_FaceP<3>(2)
+DEAL::1.00000 1.00000
+DEAL::2.44949
+DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 
+DEAL::FE_FaceP<3>(1)  vs.  FE_FaceP<3>(2)
+DEAL::1.00000 1.00000
+DEAL::1.00000 0 1.00000 1.00000 0 0 
 DEAL::Checking FE_FaceQ in 3d:
 DEAL::FE_FaceQ<3>(1)  vs.  FE_FaceQ<3>(1)
 DEAL::1.00000 1.00000
@@ -165,6 +319,18 @@ DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00
 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_FaceP in 3d:
+DEAL::FE_FaceP<3>(1)  vs.  FE_FaceP<3>(1)
+DEAL::1.00000 1.00000
+DEAL::1.73205
+DEAL::1.00000 0 1.00000 1.00000 1.00000 0 
+DEAL::FE_FaceP<3>(3)  vs.  FE_FaceP<3>(3)
+DEAL::1.00000 1.00000
+DEAL::3.16228
+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 
+DEAL::FE_FaceP<3>(1)  vs.  FE_FaceP<3>(3)
+DEAL::1.00000 1.00000
+DEAL::1.00000 0 1.00000 1.00000 0 0 
 DEAL::Checking FE_FaceQ in 3d:
 DEAL::FE_FaceQ<3>(2)  vs.  FE_FaceQ<3>(2)
 DEAL::1.00000 1.00000
@@ -177,3 +343,15 @@ DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00
 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 
+DEAL::Checking FE_FaceP in 3d:
+DEAL::FE_FaceP<3>(2)  vs.  FE_FaceP<3>(2)
+DEAL::1.00000 1.00000
+DEAL::2.44949
+DEAL::1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 1.00000 0 
+DEAL::FE_FaceP<3>(3)  vs.  FE_FaceP<3>(3)
+DEAL::1.00000 1.00000
+DEAL::3.16228
+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 
+DEAL::FE_FaceP<3>(2)  vs.  FE_FaceP<3>(3)
+DEAL::1.00000 1.00000
+DEAL::1.00000 0 1.00000 0 1.00000 0 0 0 0 0 0 0 
diff --git a/tests/bits/faceq_subface_interpolation.cc b/tests/bits/faceq_subface_interpolation.cc
index c8bfa05d5e..ecf6cb3196 100644
--- a/tests/bits/faceq_subface_interpolation.cc
+++ b/tests/bits/faceq_subface_interpolation.cc
@@ -19,7 +19,7 @@
 #include <deal.II/fe/fe_face.h>
 
 // check
-//   FE_FaceQ::subface_interpolation
+//   FE_FaceQ/P::subface_interpolation
 
 
 std::string output_file_name = "faceq_subface_interpolation/output";
@@ -128,12 +128,22 @@ void
 check (const unsigned int degree1,
        const unsigned int degree2)
 {
-  deallog << "Checking FE_FaceQ in " << dim << "d:"
-          << std::endl;
+  {
+    deallog << "Checking FE_FaceQ in " << dim << "d:"
+            << std::endl;
+
+    FE_FaceQ<dim> fe1 (degree1), fe2 (degree2);
+
+    check_this (fe1, fe2);
+  }
+  {
+    deallog << "Checking FE_FaceP in " << dim << "d:"
+            << std::endl;
 
-  FE_FaceQ<dim> fe1 (degree1), fe2 (degree2);
+    FE_FaceP<dim> fe1 (degree1), fe2 (degree2);
 
-  check_this (fe1, fe2);
+    check_this (fe1, fe2);
+  }
 }
 
 
diff --git a/tests/bits/faceq_subface_interpolation/cmp/generic b/tests/bits/faceq_subface_interpolation/cmp/generic
index df01dec611..5c0103af34 100644
--- a/tests/bits/faceq_subface_interpolation/cmp/generic
+++ b/tests/bits/faceq_subface_interpolation/cmp/generic
@@ -32,6 +32,39 @@ 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_FaceP in 2d:
+DEAL::FE_FaceP<2>(0)  vs.  FE_FaceP<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<2>(0)  vs.  FE_FaceP<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<2>(0)  vs.  FE_FaceP<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<2>(0)  vs.  FE_FaceP<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<2>(0)  vs.  FE_FaceP<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<2>(0)  vs.  FE_FaceP<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<2>(0)  vs.  FE_FaceP<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<2>(0)  vs.  FE_FaceP<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
@@ -55,6 +88,29 @@ 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_FaceP in 2d:
+DEAL::FE_FaceP<2>(0)  vs.  FE_FaceP<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<2>(1)  vs.  FE_FaceP<2>(1)
+DEAL::3.73205 2.73205
+DEAL::2.82843
+DEAL::1.00000 1.73205 2.00000 0 
+DEAL::FE_FaceP<2>(0)  vs.  FE_FaceP<2>(1)
+DEAL::1.00000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<2>(0)  vs.  FE_FaceP<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<2>(1)  vs.  FE_FaceP<2>(1)
+DEAL::3.73205 2.73205
+DEAL::2.82843
+DEAL::1.00000 -1.73205 2.00000 0 
+DEAL::FE_FaceP<2>(0)  vs.  FE_FaceP<2>(1)
+DEAL::1.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
@@ -78,6 +134,29 @@ 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_FaceP in 2d:
+DEAL::FE_FaceP<2>(0)  vs.  FE_FaceP<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<2>(2)  vs.  FE_FaceP<2>(2)
+DEAL::18.4542 9.74597
+DEAL::11.3578
+DEAL::1.00000 6.70820 2.00000 2.00000 4.00000 0 
+DEAL::FE_FaceP<2>(0)  vs.  FE_FaceP<2>(2)
+DEAL::1.00000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<2>(0)  vs.  FE_FaceP<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<2>(2)  vs.  FE_FaceP<2>(2)
+DEAL::18.4542 9.74597
+DEAL::11.3578
+DEAL::1.00000 6.70820 2.00000 2.00000 4.00000 0 
+DEAL::FE_FaceP<2>(0)  vs.  FE_FaceP<2>(2)
+DEAL::1.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
@@ -101,6 +180,29 @@ 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_FaceP in 2d:
+DEAL::FE_FaceP<2>(0)  vs.  FE_FaceP<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<2>(3)  vs.  FE_FaceP<2>(3)
+DEAL::97.4282 46.4066
+DEAL::54.2586
+DEAL::1.00000 29.1033 2.00000 7.74597 4.00000 0 8.00000 0 
+DEAL::FE_FaceP<2>(0)  vs.  FE_FaceP<2>(3)
+DEAL::1.00000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<2>(0)  vs.  FE_FaceP<2>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<2>(3)  vs.  FE_FaceP<2>(3)
+DEAL::97.4282 46.4066
+DEAL::54.2586
+DEAL::1.00000 -29.1033 2.00000 -7.74597 4.00000 0 8.00000 0 
+DEAL::FE_FaceP<2>(0)  vs.  FE_FaceP<2>(3)
+DEAL::1.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
@@ -124,6 +226,29 @@ 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_FaceP in 2d:
+DEAL::FE_FaceP<2>(1)  vs.  FE_FaceP<2>(1)
+DEAL::3.73205 2.73205
+DEAL::2.82843
+DEAL::1.00000 1.73205 2.00000 0 
+DEAL::FE_FaceP<2>(2)  vs.  FE_FaceP<2>(2)
+DEAL::18.4542 9.74597
+DEAL::11.3578
+DEAL::1.00000 6.70820 2.00000 2.00000 4.00000 0 
+DEAL::FE_FaceP<2>(1)  vs.  FE_FaceP<2>(2)
+DEAL::3.73205 2.73205
+DEAL::1.00000 1.73205 2.00000 0 
+DEAL::FE_FaceP<2>(1)  vs.  FE_FaceP<2>(1)
+DEAL::3.73205 2.73205
+DEAL::2.82843
+DEAL::1.00000 -1.73205 2.00000 0 
+DEAL::FE_FaceP<2>(2)  vs.  FE_FaceP<2>(2)
+DEAL::18.4542 9.74597
+DEAL::11.3578
+DEAL::1.00000 6.70820 2.00000 2.00000 4.00000 0 
+DEAL::FE_FaceP<2>(1)  vs.  FE_FaceP<2>(2)
+DEAL::3.73205 2.73205
+DEAL::1.00000 -1.73205 2.00000 0 
 DEAL::Checking FE_FaceQ in 2d:
 DEAL::FE_FaceQ<2>(1)  vs.  FE_FaceQ<2>(1)
 DEAL::1.50000 1.00000
@@ -147,6 +272,29 @@ 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_FaceP in 2d:
+DEAL::FE_FaceP<2>(1)  vs.  FE_FaceP<2>(1)
+DEAL::3.73205 2.73205
+DEAL::2.82843
+DEAL::1.00000 1.73205 2.00000 0 
+DEAL::FE_FaceP<2>(3)  vs.  FE_FaceP<2>(3)
+DEAL::97.4282 46.4066
+DEAL::54.2586
+DEAL::1.00000 29.1033 2.00000 7.74597 4.00000 0 8.00000 0 
+DEAL::FE_FaceP<2>(1)  vs.  FE_FaceP<2>(3)
+DEAL::3.73205 2.73205
+DEAL::1.00000 1.73205 2.00000 0 
+DEAL::FE_FaceP<2>(1)  vs.  FE_FaceP<2>(1)
+DEAL::3.73205 2.73205
+DEAL::2.82843
+DEAL::1.00000 -1.73205 2.00000 0 
+DEAL::FE_FaceP<2>(3)  vs.  FE_FaceP<2>(3)
+DEAL::97.4282 46.4066
+DEAL::54.2586
+DEAL::1.00000 -29.1033 2.00000 -7.74597 4.00000 0 8.00000 0 
+DEAL::FE_FaceP<2>(1)  vs.  FE_FaceP<2>(3)
+DEAL::3.73205 2.73205
+DEAL::1.00000 -1.73205 2.00000 0 
 DEAL::Checking FE_FaceQ in 2d:
 DEAL::FE_FaceQ<2>(2)  vs.  FE_FaceQ<2>(2)
 DEAL::1.75000 1.25000
@@ -170,6 +318,29 @@ 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_FaceP in 2d:
+DEAL::FE_FaceP<2>(2)  vs.  FE_FaceP<2>(2)
+DEAL::18.4542 9.74597
+DEAL::11.3578
+DEAL::1.00000 6.70820 2.00000 2.00000 4.00000 0 
+DEAL::FE_FaceP<2>(3)  vs.  FE_FaceP<2>(3)
+DEAL::97.4282 46.4066
+DEAL::54.2586
+DEAL::1.00000 29.1033 2.00000 7.74597 4.00000 0 8.00000 0 
+DEAL::FE_FaceP<2>(2)  vs.  FE_FaceP<2>(3)
+DEAL::18.4542 9.74597
+DEAL::1.00000 6.70820 2.00000 2.00000 4.00000 0 
+DEAL::FE_FaceP<2>(2)  vs.  FE_FaceP<2>(2)
+DEAL::18.4542 9.74597
+DEAL::11.3578
+DEAL::1.00000 6.70820 2.00000 2.00000 4.00000 0 
+DEAL::FE_FaceP<2>(3)  vs.  FE_FaceP<2>(3)
+DEAL::97.4282 46.4066
+DEAL::54.2586
+DEAL::1.00000 -29.1033 2.00000 -7.74597 4.00000 0 8.00000 0 
+DEAL::FE_FaceP<2>(2)  vs.  FE_FaceP<2>(3)
+DEAL::18.4542 9.74597
+DEAL::1.00000 6.70820 2.00000 2.00000 4.00000 0 
 DEAL::Checking FE_FaceQ in 3d:
 DEAL::FE_FaceQ<3>(0)  vs.  FE_FaceQ<3>(0)
 DEAL::1.00000 1.00000
@@ -235,6 +406,71 @@ 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_FaceP in 3d:
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<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
@@ -280,6 +516,51 @@ 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_FaceP in 3d:
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(1)  vs.  FE_FaceP<3>(1)
+DEAL::3.73205 4.46410
+DEAL::3.87298
+DEAL::1.00000 1.73205 2.00000 2.00000 2.00000 0 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(1)
+DEAL::1.00000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(1)  vs.  FE_FaceP<3>(1)
+DEAL::3.73205 4.46410
+DEAL::3.87298
+DEAL::1.00000 1.73205 2.00000 2.00000 2.00000 0 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(1)
+DEAL::1.00000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(1)  vs.  FE_FaceP<3>(1)
+DEAL::3.73205 4.46410
+DEAL::3.87298
+DEAL::1.00000 -1.73205 2.00000 2.00000 2.00000 0 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(1)
+DEAL::1.00000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(1)  vs.  FE_FaceP<3>(1)
+DEAL::3.73205 4.46410
+DEAL::3.87298
+DEAL::1.00000 -1.73205 2.00000 2.00000 2.00000 0 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(1)
+DEAL::1.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
@@ -325,6 +606,51 @@ 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_FaceP in 3d:
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(2)  vs.  FE_FaceP<3>(2)
+DEAL::18.4542 20.8805
+DEAL::17.4929
+DEAL::1.00000 6.70820 2.00000 3.46410 4.00000 0 2.00000 0 4.00000 0 4.00000 0 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(2)
+DEAL::1.00000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(2)  vs.  FE_FaceP<3>(2)
+DEAL::18.4542 20.8805
+DEAL::17.4929
+DEAL::1.00000 6.70820 2.00000 3.46410 4.00000 0 2.00000 0 4.00000 0 4.00000 0 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(2)
+DEAL::1.00000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(2)  vs.  FE_FaceP<3>(2)
+DEAL::18.4542 20.8805
+DEAL::17.4929
+DEAL::1.00000 6.70820 2.00000 -3.46410 4.00000 0 2.00000 0 4.00000 0 4.00000 0 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(2)
+DEAL::1.00000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(2)  vs.  FE_FaceP<3>(2)
+DEAL::18.4542 20.8805
+DEAL::17.4929
+DEAL::1.00000 6.70820 2.00000 -3.46410 4.00000 0 2.00000 0 4.00000 0 4.00000 0 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(2)
+DEAL::1.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
@@ -370,6 +696,51 @@ DEAL::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::Checking FE_FaceP in 3d:
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(3)  vs.  FE_FaceP<3>(3)
+DEAL::97.4282 102.325
+DEAL::87.3499
+DEAL::1.00000 29.1033 2.00000 13.4164 4.00000 0 8.00000 0 2.00000 3.46410 4.00000 0 8.00000 0 4.00000 0 8.00000 0 8.00000 0 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(3)
+DEAL::1.00000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(3)  vs.  FE_FaceP<3>(3)
+DEAL::97.4282 102.325
+DEAL::87.3499
+DEAL::1.00000 29.1033 2.00000 13.4164 4.00000 0 8.00000 0 2.00000 -3.46410 4.00000 0 8.00000 0 4.00000 0 8.00000 0 8.00000 0 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(3)
+DEAL::1.00000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(3)  vs.  FE_FaceP<3>(3)
+DEAL::97.4282 102.325
+DEAL::87.3499
+DEAL::1.00000 -29.1033 2.00000 13.4164 4.00000 0 8.00000 0 2.00000 3.46410 4.00000 0 8.00000 0 4.00000 0 8.00000 0 8.00000 0 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(3)
+DEAL::1.00000 1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(0)
+DEAL::1.00000 1.00000
+DEAL::1.00000
+DEAL::1.00000 1.00000 
+DEAL::FE_FaceP<3>(3)  vs.  FE_FaceP<3>(3)
+DEAL::97.4282 102.325
+DEAL::87.3499
+DEAL::1.00000 -29.1033 2.00000 13.4164 4.00000 0 8.00000 0 2.00000 -3.46410 4.00000 0 8.00000 0 4.00000 0 8.00000 0 8.00000 0 
+DEAL::FE_FaceP<3>(0)  vs.  FE_FaceP<3>(3)
+DEAL::1.00000 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
@@ -415,6 +786,51 @@ 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_FaceP in 3d:
+DEAL::FE_FaceP<3>(1)  vs.  FE_FaceP<3>(1)
+DEAL::3.73205 4.46410
+DEAL::3.87298
+DEAL::1.00000 1.73205 2.00000 2.00000 2.00000 0 
+DEAL::FE_FaceP<3>(2)  vs.  FE_FaceP<3>(2)
+DEAL::18.4542 20.8805
+DEAL::17.4929
+DEAL::1.00000 6.70820 2.00000 3.46410 4.00000 0 2.00000 0 4.00000 0 4.00000 0 
+DEAL::FE_FaceP<3>(1)  vs.  FE_FaceP<3>(2)
+DEAL::3.73205 4.46410
+DEAL::1.00000 1.73205 2.00000 2.00000 0 0 
+DEAL::FE_FaceP<3>(1)  vs.  FE_FaceP<3>(1)
+DEAL::3.73205 4.46410
+DEAL::3.87298
+DEAL::1.00000 1.73205 2.00000 2.00000 2.00000 0 
+DEAL::FE_FaceP<3>(2)  vs.  FE_FaceP<3>(2)
+DEAL::18.4542 20.8805
+DEAL::17.4929
+DEAL::1.00000 6.70820 2.00000 3.46410 4.00000 0 2.00000 0 4.00000 0 4.00000 0 
+DEAL::FE_FaceP<3>(1)  vs.  FE_FaceP<3>(2)
+DEAL::3.73205 4.46410
+DEAL::1.00000 1.73205 2.00000 2.00000 0 0 
+DEAL::FE_FaceP<3>(1)  vs.  FE_FaceP<3>(1)
+DEAL::3.73205 4.46410
+DEAL::3.87298
+DEAL::1.00000 -1.73205 2.00000 2.00000 2.00000 0 
+DEAL::FE_FaceP<3>(2)  vs.  FE_FaceP<3>(2)
+DEAL::18.4542 20.8805
+DEAL::17.4929
+DEAL::1.00000 6.70820 2.00000 -3.46410 4.00000 0 2.00000 0 4.00000 0 4.00000 0 
+DEAL::FE_FaceP<3>(1)  vs.  FE_FaceP<3>(2)
+DEAL::3.73205 4.46410
+DEAL::1.00000 -1.73205 2.00000 2.00000 0 0 
+DEAL::FE_FaceP<3>(1)  vs.  FE_FaceP<3>(1)
+DEAL::3.73205 4.46410
+DEAL::3.87298
+DEAL::1.00000 -1.73205 2.00000 2.00000 2.00000 0 
+DEAL::FE_FaceP<3>(2)  vs.  FE_FaceP<3>(2)
+DEAL::18.4542 20.8805
+DEAL::17.4929
+DEAL::1.00000 6.70820 2.00000 -3.46410 4.00000 0 2.00000 0 4.00000 0 4.00000 0 
+DEAL::FE_FaceP<3>(1)  vs.  FE_FaceP<3>(2)
+DEAL::3.73205 4.46410
+DEAL::1.00000 -1.73205 2.00000 2.00000 0 0 
 DEAL::Checking FE_FaceQ in 3d:
 DEAL::FE_FaceQ<3>(1)  vs.  FE_FaceQ<3>(1)
 DEAL::2.25000 1.00000
@@ -460,6 +876,51 @@ DEAL::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.250000 0.250000 0.319098 0.180902 0.0690983 0.430902 0.500000 0 
+DEAL::Checking FE_FaceP in 3d:
+DEAL::FE_FaceP<3>(1)  vs.  FE_FaceP<3>(1)
+DEAL::3.73205 4.46410
+DEAL::3.87298
+DEAL::1.00000 1.73205 2.00000 2.00000 2.00000 0 
+DEAL::FE_FaceP<3>(3)  vs.  FE_FaceP<3>(3)
+DEAL::97.4282 102.325
+DEAL::87.3499
+DEAL::1.00000 29.1033 2.00000 13.4164 4.00000 0 8.00000 0 2.00000 3.46410 4.00000 0 8.00000 0 4.00000 0 8.00000 0 8.00000 0 
+DEAL::FE_FaceP<3>(1)  vs.  FE_FaceP<3>(3)
+DEAL::3.73205 4.46410
+DEAL::1.00000 1.73205 2.00000 2.00000 0 0 
+DEAL::FE_FaceP<3>(1)  vs.  FE_FaceP<3>(1)
+DEAL::3.73205 4.46410
+DEAL::3.87298
+DEAL::1.00000 1.73205 2.00000 2.00000 2.00000 0 
+DEAL::FE_FaceP<3>(3)  vs.  FE_FaceP<3>(3)
+DEAL::97.4282 102.325
+DEAL::87.3499
+DEAL::1.00000 29.1033 2.00000 13.4164 4.00000 0 8.00000 0 2.00000 -3.46410 4.00000 0 8.00000 0 4.00000 0 8.00000 0 8.00000 0 
+DEAL::FE_FaceP<3>(1)  vs.  FE_FaceP<3>(3)
+DEAL::3.73205 4.46410
+DEAL::1.00000 1.73205 2.00000 2.00000 0 0 
+DEAL::FE_FaceP<3>(1)  vs.  FE_FaceP<3>(1)
+DEAL::3.73205 4.46410
+DEAL::3.87298
+DEAL::1.00000 -1.73205 2.00000 2.00000 2.00000 0 
+DEAL::FE_FaceP<3>(3)  vs.  FE_FaceP<3>(3)
+DEAL::97.4282 102.325
+DEAL::87.3499
+DEAL::1.00000 -29.1033 2.00000 13.4164 4.00000 0 8.00000 0 2.00000 3.46410 4.00000 0 8.00000 0 4.00000 0 8.00000 0 8.00000 0 
+DEAL::FE_FaceP<3>(1)  vs.  FE_FaceP<3>(3)
+DEAL::3.73205 4.46410
+DEAL::1.00000 -1.73205 2.00000 2.00000 0 0 
+DEAL::FE_FaceP<3>(1)  vs.  FE_FaceP<3>(1)
+DEAL::3.73205 4.46410
+DEAL::3.87298
+DEAL::1.00000 -1.73205 2.00000 2.00000 2.00000 0 
+DEAL::FE_FaceP<3>(3)  vs.  FE_FaceP<3>(3)
+DEAL::97.4282 102.325
+DEAL::87.3499
+DEAL::1.00000 -29.1033 2.00000 13.4164 4.00000 0 8.00000 0 2.00000 -3.46410 4.00000 0 8.00000 0 4.00000 0 8.00000 0 8.00000 0 
+DEAL::FE_FaceP<3>(1)  vs.  FE_FaceP<3>(3)
+DEAL::3.73205 4.46410
+DEAL::1.00000 -1.73205 2.00000 2.00000 0 0 
 DEAL::Checking FE_FaceQ in 3d:
 DEAL::FE_FaceQ<3>(2)  vs.  FE_FaceQ<3>(2)
 DEAL::3.06250 1.56250
@@ -505,3 +966,48 @@ DEAL::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 0 0 0 0 1.00000 0.923607 0.923607 0.162918 -0.0923607 -0.0176393 -0.0623607 0 0 0 0 
+DEAL::Checking FE_FaceP in 3d:
+DEAL::FE_FaceP<3>(2)  vs.  FE_FaceP<3>(2)
+DEAL::18.4542 20.8805
+DEAL::17.4929
+DEAL::1.00000 6.70820 2.00000 3.46410 4.00000 0 2.00000 0 4.00000 0 4.00000 0 
+DEAL::FE_FaceP<3>(3)  vs.  FE_FaceP<3>(3)
+DEAL::97.4282 102.325
+DEAL::87.3499
+DEAL::1.00000 29.1033 2.00000 13.4164 4.00000 0 8.00000 0 2.00000 3.46410 4.00000 0 8.00000 0 4.00000 0 8.00000 0 8.00000 0 
+DEAL::FE_FaceP<3>(2)  vs.  FE_FaceP<3>(3)
+DEAL::18.4542 20.8805
+DEAL::1.00000 6.70820 2.00000 3.46410 4.00000 0 0 0 3.46410 0 0 0 
+DEAL::FE_FaceP<3>(2)  vs.  FE_FaceP<3>(2)
+DEAL::18.4542 20.8805
+DEAL::17.4929
+DEAL::1.00000 6.70820 2.00000 3.46410 4.00000 0 2.00000 0 4.00000 0 4.00000 0 
+DEAL::FE_FaceP<3>(3)  vs.  FE_FaceP<3>(3)
+DEAL::97.4282 102.325
+DEAL::87.3499
+DEAL::1.00000 29.1033 2.00000 13.4164 4.00000 0 8.00000 0 2.00000 -3.46410 4.00000 0 8.00000 0 4.00000 0 8.00000 0 8.00000 0 
+DEAL::FE_FaceP<3>(2)  vs.  FE_FaceP<3>(3)
+DEAL::18.4542 20.8805
+DEAL::1.00000 6.70820 2.00000 3.46410 4.00000 0 0 0 -3.46410 0 0 0 
+DEAL::FE_FaceP<3>(2)  vs.  FE_FaceP<3>(2)
+DEAL::18.4542 20.8805
+DEAL::17.4929
+DEAL::1.00000 6.70820 2.00000 -3.46410 4.00000 0 2.00000 0 4.00000 0 4.00000 0 
+DEAL::FE_FaceP<3>(3)  vs.  FE_FaceP<3>(3)
+DEAL::97.4282 102.325
+DEAL::87.3499
+DEAL::1.00000 -29.1033 2.00000 13.4164 4.00000 0 8.00000 0 2.00000 3.46410 4.00000 0 8.00000 0 4.00000 0 8.00000 0 8.00000 0 
+DEAL::FE_FaceP<3>(2)  vs.  FE_FaceP<3>(3)
+DEAL::18.4542 20.8805
+DEAL::1.00000 6.70820 2.00000 -3.46410 4.00000 0 0 0 3.46410 0 0 0 
+DEAL::FE_FaceP<3>(2)  vs.  FE_FaceP<3>(2)
+DEAL::18.4542 20.8805
+DEAL::17.4929
+DEAL::1.00000 6.70820 2.00000 -3.46410 4.00000 0 2.00000 0 4.00000 0 4.00000 0 
+DEAL::FE_FaceP<3>(3)  vs.  FE_FaceP<3>(3)
+DEAL::97.4282 102.325
+DEAL::87.3499
+DEAL::1.00000 -29.1033 2.00000 13.4164 4.00000 0 8.00000 0 2.00000 -3.46410 4.00000 0 8.00000 0 4.00000 0 8.00000 0 8.00000 0 
+DEAL::FE_FaceP<3>(2)  vs.  FE_FaceP<3>(3)
+DEAL::18.4542 20.8805
+DEAL::1.00000 6.70820 2.00000 -3.46410 4.00000 0 0 0 -3.46410 0 0 0 
-- 
2.39.5