]> https://gitweb.dealii.org/ - dealii.git/commitdiff
MappingQGeneric: use of MatrixFree only for higher order elements 7687/head
authorStefan Meggendorfer <stefan.meggendorfer@iwr.uni-heidelberg.de>
Mon, 4 Feb 2019 13:56:57 +0000 (14:56 +0100)
committerStefan Meggendorfer <stefan.meggendorfer@iwr.uni-heidelberg.de>
Mon, 4 Feb 2019 13:56:57 +0000 (14:56 +0100)
source/fe/mapping_q_generic.cc

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

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.