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 4dfaf968dc6c651b2ca121a9f2737a895b5f81d3..d7d8e96d88f2bc2fca58bb1d548e80971a3cc98b 100644 (file)
--- 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