Change MappingManifold<1, spacedim> to use face manifolds.

As of db5ea0f52db we support get_manifold on 1D manifolds, so we can look up the
manifold on a 1D face in a dimension-independent way. This is a slight change
from the current behavior, but using Manifolds in 1D is a very obscure use case.

diff --git a/doc/news/changes/incompatibilities/20170708DavidWells b/doc/news/changes/incompatibilities/20170708DavidWells
new file mode 100644 (file)
index 0000000..818dece
--- /dev/null
+++ b/doc/news/changes/incompatibilities/20170708DavidWells
@@ -0,0 +1,3 @@
+Changed: In 1D, the face values calculated by <code>MappingManifold</code> now use the manifold attached to the relevant vertex instead of the manifold on the current face.
+<br>
+(David Wells, 2017/07/08)

diff --git a/source/fe/mapping_manifold.cc b/source/fe/mapping_manifold.cc
index 2c8297ca8c440da734ed1c0b6a3284eac7100364..2ac37d2cb988651377f8be7f622c61064c8ecc30 100644 (file)
--- a/source/fe/mapping_manifold.cc
+++ b/source/fe/mapping_manifold.cc
@@ -351,37 +351,6 @@ namespace internal
 {
   namespace
   {
-    /**
-     * Some specialization for face Manifolds. In one dimension, there
-     * are no Manifolds associated to faces. The mapping argument is
-     * only used to help the compiler infer dim and spacedim.
-     */
-    template <int spacedim>
-    const dealii::Manifold<1, spacedim> &
-    get_face_manifold(const MappingManifold<1,spacedim> &,
-                      const typename dealii::Triangulation<1,spacedim>::cell_iterator &cell,
-                      const unsigned int &)
-    {
-      return cell->get_manifold();
-    }
-
-
-
-    /**
-     * Some specialization for face Manifolds. The mapping argument is
-     * only used to help the compiler infer dim and spacedim.
-     */
-    template <int dim, int spacedim>
-    const dealii::Manifold<dim,spacedim> &
-    get_face_manifold(const MappingManifold<dim,spacedim> &,
-                      const typename dealii::Triangulation<dim,spacedim>::cell_iterator &cell,
-                      const unsigned int face_no)
-    {
-      return cell->face(face_no)->get_manifold();
-    }
-
-
-
     /**
      * Compute the locations of quadrature points on the object described by
      * the first argument (and the cell for which the mapping support points
@@ -822,7 +791,7 @@ namespace internal
     {
       data.store_vertices(cell);
 
-      data.manifold = &get_face_manifold(mapping, cell, face_no);
+      data.manifold = &cell->face(face_no)->get_manifold();
 
       maybe_compute_q_points<dim,spacedim> (data_set,
                                             data,