From: Peter Munch <peterrmuench@gmail.com>
Date: Wed, 19 Jan 2022 16:19:18 +0000 (+0100)
Subject: Merge pull request #13259 from peterrum/fe_eval_n_vectorization_lanes_filled
X-Git-Tag: v9.4.0-rc1~587
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=0052b8588f525f25c949de13ec718d496ff76735;p=dealii.git

Merge pull request #13259 from peterrum/fe_eval_n_vectorization_lanes_filled

FEEval: replace usages of n_vectorization_actual
---

0052b8588f525f25c949de13ec718d496ff76735
diff --cc include/deal.II/matrix_free/fe_evaluation.h
index a5c4e384ba,f8964cbbb3..c8ddf1ea76
--- a/include/deal.II/matrix_free/fe_evaluation.h
+++ b/include/deal.II/matrix_free/fe_evaluation.h
@@@ -6864,14 -6893,15 +6884,18 @@@ FEEvaluation<dim
      this->mapping_data->data_index_offsets[cell_index];
    this->jacobian = &this->mapping_data->jacobians[0][offsets];
    this->J_value  = &this->mapping_data->JxW_values[offsets];
 +  this->jacobian_gradients =
 +    this->mapping_data->jacobian_gradients[0].data() + offsets;
  
-   for (unsigned int i = 0; i < VectorizedArrayType::size(); ++i)
+   unsigned int i = 0;
+   for (; i < this->matrix_free->n_active_entries_per_cell_batch(this->cell);
+        ++i)
      this->cell_ids[i] = cell_index * VectorizedArrayType::size() + i;
+   for (; i < VectorizedArrayType::size(); ++i)
+     this->cell_ids[i] = numbers::invalid_unsigned_int;
  
  #  ifdef DEBUG
 +  this->is_reinitialized           = true;
    this->dof_values_initialized     = false;
    this->values_quad_initialized    = false;
    this->gradients_quad_initialized = false;