]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Fix bug in MappingQGeneric::transform_real_to_unit_cell
authorMartin Kronbichler <kronbichler@lnm.mw.tum.de>
Thu, 14 Nov 2019 09:18:31 +0000 (10:18 +0100)
committerMartin Kronbichler <kronbichler@lnm.mw.tum.de>
Thu, 14 Nov 2019 09:18:31 +0000 (10:18 +0100)
source/fe/mapping_q_generic.cc

index b4e5d22f8b8307918251289332c4d916e4be7317..959445efecbe6e0b8466b13a0a3dcd14fb8f196f 100644 (file)
@@ -2534,7 +2534,7 @@ MappingQGeneric<dim, spacedim>::transform_real_to_unit_cell(
 {
   // Use an exact formula if one is available. this is only the case
   // for Q1 mappings in 1d, and in 2d if dim==spacedim
-  if ((polynomial_degree == 1) &&
+  if (this->preserves_vertex_locations() && (polynomial_degree == 1) &&
       ((dim == 1) || ((dim == 2) && (dim == spacedim))))
     {
       // The dimension-dependent algorithms are much faster (about 25-45x in

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.