]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Fix bug in eulerian mapping regarding preservation of vertex locations
authorMartin Kronbichler <kronbichler@lnm.mw.tum.de>
Sat, 3 Jun 2017 14:51:00 +0000 (16:51 +0200)
committerMartin Kronbichler <kronbichler@lnm.mw.tum.de>
Mon, 5 Jun 2017 14:09:35 +0000 (16:09 +0200)
include/deal.II/fe/mapping_q_eulerian.h
source/fe/mapping_q_eulerian.cc

index d21ca6ce4bc1679c39c43f7a70fbed298a015a47..96613488babf2cf80835ef20f267d85aabbab18a 100644 (file)
@@ -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.
index bd9f65ffc23ab827e6e2fd1a0e6e941889f2e952..39250e5d7a71546276c339d85fab5305dea6b441 100644 (file)
@@ -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> >

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.