]> https://gitweb.dealii.org/ - dealii.git/commitdiff
MatrixFree: setup SM capabilities only in non-hp mode 11255/head
authorPeter Munch <peterrmuench@gmail.com>
Thu, 26 Nov 2020 09:32:03 +0000 (10:32 +0100)
committerPeter Munch <peterrmuench@gmail.com>
Thu, 26 Nov 2020 11:19:45 +0000 (12:19 +0100)
include/deal.II/matrix_free/matrix_free.templates.h

index 2a32d279dafd6adf7ea7a025392c59276c7f9e55..f2e2325cf1941ad0603850804fac4f16032cb8b4 100644 (file)
@@ -1539,13 +1539,16 @@ MatrixFree<dim, Number, VectorizedArrayType>::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<unsigned int>(

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.