From ddbb09b90e85afb1f0e1f039500036fec4681e24 Mon Sep 17 00:00:00 2001 From: Stefan Meggendorfer Date: Mon, 4 Feb 2019 14:56:57 +0100 Subject: [PATCH] MappingQGeneric: use of MatrixFree only for higher order elements --- source/fe/mapping_q_generic.cc | 4 ++++ 1 file changed, 4 insertions(+) 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::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 -- 2.39.5