From 4194e0162bffbbf8800d9d3eda33eecb59742680 Mon Sep 17 00:00:00 2001
From: kronbichler <kronbichler@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Tue, 17 Sep 2013 13:29:50 +0000
Subject: [PATCH] Fix a few issues. Still not correct.

git-svn-id: https://svn.dealii.org/trunk@30763 0785d39b-7218-0410-832d-ea1e28bc413d
---
 deal.II/source/fe/fe_face.cc      | 5 ++++-
 deal.II/source/fe/fe_face.inst.in | 3 ++-
 2 files changed, 6 insertions(+), 2 deletions(-)

diff --git a/deal.II/source/fe/fe_face.cc b/deal.II/source/fe/fe_face.cc
index 0478dfbad8..a362e3cdec 100644
--- a/deal.II/source/fe/fe_face.cc
+++ b/deal.II/source/fe/fe_face.cc
@@ -441,6 +441,8 @@ get_subface_interpolation_matrix (const FiniteElement<dim,spacedim> &x_source_fe
       for (unsigned int k = 0; k < face_quadrature.size(); ++k)
         {
           const Point<dim-1> p =
+            subface == numbers::invalid_unsigned_int ?
+            face_quadrature.point(k) :
             GeometryInfo<dim-1>::child_to_cell_coordinates (face_quadrature.point(k),
                                                             subface);
 
@@ -459,7 +461,8 @@ get_subface_interpolation_matrix (const FiniteElement<dim,spacedim> &x_source_fe
         {
           for (unsigned int k = 0; k < face_quadrature.size(); ++k)
             {
-              const Point<dim-1> p =
+              const Point<dim-1> p = numbers::invalid_unsigned_int ?
+                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);
diff --git a/deal.II/source/fe/fe_face.inst.in b/deal.II/source/fe/fe_face.inst.in
index 8d99474428..513eb0030b 100644
--- a/deal.II/source/fe/fe_face.inst.in
+++ b/deal.II/source/fe/fe_face.inst.in
@@ -20,9 +20,10 @@ for (deal_II_dimension : DIMENSIONS)
   {
 #if deal_II_dimension > 1
     template class FE_PolyFace<TensorProductPolynomials<deal_II_dimension-1> >;
-    //template class FE_PolyFace<PolynomialSpace<deal_II_dimension>, deal_II_dimension>;
+    template class FE_PolyFace<PolynomialSpace<deal_II_dimension-1>, deal_II_dimension>;
     //template class FE_PolyFace<PolynomialsP<deal_II_dimension>, deal_II_dimension>;
     template class FE_FaceQ<deal_II_dimension>;
+    template class FE_FaceP<deal_II_dimension>;
 #endif
   }
 
-- 
2.39.5