From: Peter Munch Date: Thu, 26 Nov 2020 09:32:03 +0000 (+0100) Subject: MatrixFree: setup SM capabilities only in non-hp mode X-Git-Tag: v9.3.0-rc1~846^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F11255%2Fhead;p=dealii.git MatrixFree: setup SM capabilities only in non-hp mode --- diff --git a/include/deal.II/matrix_free/matrix_free.templates.h b/include/deal.II/matrix_free/matrix_free.templates.h index 2a32d279da..f2e2325cf1 100644 --- a/include/deal.II/matrix_free/matrix_free.templates.h +++ b/include/deal.II/matrix_free/matrix_free.templates.h @@ -1539,13 +1539,16 @@ MatrixFree::initialize_indices( #ifdef DEAL_II_WITH_MPI { // non-buffering mode is only supported if the indices of all cells are - // contiguous for all dof_info objects. + // contiguous for all dof_info objects and hp is not enabled. bool is_non_buffering_sm_supported = true; for (const auto &di : dof_info) - for (const auto &v : di.index_storage_variants[2]) - is_non_buffering_sm_supported &= - (v == internal::MatrixFreeFunctions::DoFInfo::IndexStorageVariants:: - contiguous); + { + is_non_buffering_sm_supported &= di.dofs_per_cell.size() == 1; + for (const auto &v : di.index_storage_variants[2]) + is_non_buffering_sm_supported &= + (v == internal::MatrixFreeFunctions::DoFInfo::IndexStorageVariants:: + contiguous); + } is_non_buffering_sm_supported = Utilities::MPI::min(static_cast(