]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Add mask 14090/head
authorMaximilian Bergbauer <maximilian.bergbauer@tum.de>
Fri, 1 Jul 2022 09:56:33 +0000 (11:56 +0200)
committerMaximilian Bergbauer <maximilian.bergbauer@tum.de>
Fri, 1 Jul 2022 09:56:33 +0000 (11:56 +0200)
include/deal.II/matrix_free/fe_evaluation.h

index ed5d5692671d8a490c829d7b0beda92024237221..6375c6c1130055885641b5b99acbee97fe44173e 100644 (file)
@@ -3914,13 +3914,17 @@ FEEvaluationBase<dim, n_components_, Number, is_face, VectorizedArrayType>::
 
   for (unsigned int v = 0; v < n_filled_lanes; ++v)
     {
-      Assert(cells[v] != numbers::invalid_unsigned_int, ExcNotImplemented());
-      dof_indices[v] =
-        dof_indices_cont[cells[v]] +
-        this->dof_info
-            ->component_dof_indices_offset[this->active_fe_index]
-                                          [this->first_selected_component] *
-          this->dof_info->dof_indices_interleave_strides[ind][cells[v]];
+      Assert(mask[v] == false || cells[v] != numbers::invalid_unsigned_int,
+             ExcNotImplemented());
+      if (mask[v] == true)
+        dof_indices[v] =
+          dof_indices_cont[cells[v]] +
+          this->dof_info
+              ->component_dof_indices_offset[this->active_fe_index]
+                                            [this->first_selected_component] *
+            this->dof_info->dof_indices_interleave_strides[ind][cells[v]];
+      else
+        dof_indices[v] = numbers::invalid_unsigned_int;
     }
 
   for (unsigned int v = n_filled_lanes; v < VectorizedArrayType::size(); ++v)

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.