From: Stefan Meggendorfer <stefan.meggendorfer@iwr.uni-heidelberg.de>
Date: Mon, 4 Feb 2019 13:56:57 +0000 (+0100)
Subject: MappingQGeneric: use of MatrixFree only for higher order elements
X-Git-Tag: v9.1.0-rc1~369^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F7687%2Fhead;p=dealii.git

MappingQGeneric: use of MatrixFree only for higher order elements
---

diff --git a/source/fe/mapping_q_generic.cc b/source/fe/mapping_q_generic.cc
index 4dfaf968dc..d7d8e96d88 100644
--- a/source/fe/mapping_q_generic.cc
+++ b/source/fe/mapping_q_generic.cc
@@ -730,6 +730,10 @@ MappingQGeneric<dim, spacedim>::InternalData::initialize(
 
   tensor_product_quadrature = q.is_tensor_product();
 
+  // use of MatrixFree only for higher order elements
+  if (polynomial_degree < 2)
+    tensor_product_quadrature = false;
+
   if (dim > 1)
     {
       // find out if the one-dimensional formula is the same