Fix bug in eulerian mapping regarding preservation of vertex locations

diff --git a/include/deal.II/fe/mapping_q_eulerian.h b/include/deal.II/fe/mapping_q_eulerian.h
index d21ca6ce4bc1679c39c43f7a70fbed298a015a47..96613488babf2cf80835ef20f267d85aabbab18a 100644 (file)
--- a/include/deal.II/fe/mapping_q_eulerian.h
+++ b/include/deal.II/fe/mapping_q_eulerian.h
@@ -129,6 +129,7 @@ public:
    * preserve vertex locations (unless the translation vector happens to
    * provide for zero displacements at vertex locations).
    */
+  virtual
   bool preserves_vertex_locations () const;
 
   /**
@@ -200,6 +201,14 @@ private:
     std::vector<Point<spacedim> >
     compute_mapping_support_points(const typename Triangulation<dim,spacedim>::cell_iterator &cell) const;
 
+    /**
+     * Always returns @p false because MappingQ1Eulerian does not in general
+     * preserve vertex locations (unless the translation vector happens to
+     * provide for zero displacements at vertex locations).
+     */
+    virtual
+    bool preserves_vertex_locations () const;
+
   private:
     /**
      * Reference to the surrounding object off of which we live.

diff --git a/source/fe/mapping_q_eulerian.cc b/source/fe/mapping_q_eulerian.cc
index bd9f65ffc23ab827e6e2fd1a0e6e941889f2e952..39250e5d7a71546276c339d85fab5305dea6b441 100644 (file)
--- a/source/fe/mapping_q_eulerian.cc
+++ b/source/fe/mapping_q_eulerian.cc
@@ -153,6 +153,15 @@ get_vertices
 
 
 
+template <int dim, class VectorType, int spacedim>
+bool
+MappingQEulerian<dim, VectorType, spacedim>::MappingQEulerianGeneric::
+preserves_vertex_locations () const
+{
+  return false;
+}
+
+
 
 template <int dim, class VectorType, int spacedim>
 std::vector<Point<spacedim> >