Fix a few issues. Still not correct.

git-svn-id: https://svn.dealii.org/trunk@30763 0785d39b-7218-0410-832d-ea1e28bc413d

diff --git a/deal.II/source/fe/fe_face.cc b/deal.II/source/fe/fe_face.cc
index 0478dfbad821336cfd56b8377b879168cdc987df..a362e3cdecf075d0bc0f442f9491b702f2753a68 100644 (file)
--- 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 8d9947442819aebfb2247b10b6d006a4d496e97c..513eb0030bf9060d0f2f0028a595eff2da8f39b2 100644 (file)
--- 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
   }