From: Wolfgang Bangerth <bangerth@math.tamu.edu>
Date: Sun, 6 Sep 2015 23:26:09 +0000 (-0500)
Subject: Relegate getting an initial guess to StaticMappingQ1.
X-Git-Tag: v8.4.0-rc2~439^2~4
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=27e4dabe6a38383b3d671c41ea8417346e6236af;p=dealii.git

Relegate getting an initial guess to StaticMappingQ1.
---

diff --git a/source/fe/mapping_q.cc b/source/fe/mapping_q.cc
index d1aeb6b850..4447272c2b 100644
--- a/source/fe/mapping_q.cc
+++ b/source/fe/mapping_q.cc
@@ -63,8 +63,17 @@ MappingQ<dim,spacedim>::MappingQ (const unsigned int degree,
           ((dim==2) ?
            4+4*(degree-1) :
            8+12*(degree-1)+6*(degree-1)*(degree-1))),
-  use_mapping_q_on_all_cells (use_mapping_q_on_all_cells
-                              || (dim != spacedim)),
+
+  // see whether we want to use *this* mapping objects on *all* cells,
+  // or defer to an explicit Q1 mapping on interior cells. if
+  // degree==1, then we are already that Q1 mapping, so we don't need
+  // it; if dim!=spacedim, there is also no need for anything because
+  // we're most likely on a curved manifold
+  use_mapping_q_on_all_cells (degree==1
+                              ||
+                              use_mapping_q_on_all_cells
+                              ||
+                              (dim != spacedim)),
   feq(degree),
   q1_mapping (this->use_mapping_q_on_all_cells
               ?
@@ -1028,7 +1037,7 @@ transform_real_to_unit_cell (const typename Triangulation<dim,spacedim>::cell_it
   try
     {
       initial_p_unit
-        = MappingQ1<dim,spacedim>::transform_real_to_unit_cell(cell, p);
+        = StaticMappingQ1<dim,spacedim>::mapping.transform_real_to_unit_cell(cell, p);
     }
   catch (const typename Mapping<dim,spacedim>::ExcTransformationFailed &)
     {