]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Fix a few issues. Still not correct.
authorMartin Kronbichler <kronbichler@lnm.mw.tum.de>
Tue, 17 Sep 2013 13:29:50 +0000 (13:29 +0000)
committerMartin Kronbichler <kronbichler@lnm.mw.tum.de>
Tue, 17 Sep 2013 13:29:50 +0000 (13:29 +0000)
git-svn-id: https://svn.dealii.org/trunk@30763 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/source/fe/fe_face.cc
deal.II/source/fe/fe_face.inst.in

index 0478dfbad821336cfd56b8377b879168cdc987df..a362e3cdecf075d0bc0f442f9491b702f2753a68 100644 (file)
@@ -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);
index 8d9947442819aebfb2247b10b6d006a4d496e97c..513eb0030bf9060d0f2f0028a595eff2da8f39b2 100644 (file)
@@ -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
   }
 

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.