Fix a bug in FE_Face.

diff --git a/source/fe/fe_face.cc b/source/fe/fe_face.cc
index bed160f792f58244f9ae7d860326cfeb64a7efc9..79c9797f1d75599877623e0a2e0a8949d26aba7e 100644 (file)
--- a/source/fe/fe_face.cc
+++ b/source/fe/fe_face.cc
@@ -956,7 +956,7 @@ FE_FaceP<dim, spacedim>::get_subface_interpolation_matrix(
           for (unsigned int k = 0; k < face_quadrature.size(); ++k)
             {
               const Point<dim - 1> p =
-                numbers::invalid_unsigned_int != 0u ?
+                subface == numbers::invalid_unsigned_int ?
                   face_quadrature.point(k) :
                   GeometryInfo<dim - 1>::child_to_cell_coordinates(
                     face_quadrature.point(k), subface);