]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Merge pull request #13259 from peterrum/fe_eval_n_vectorization_lanes_filled
authorPeter Munch <peterrmuench@gmail.com>
Wed, 19 Jan 2022 16:19:18 +0000 (17:19 +0100)
committerGitHub <noreply@github.com>
Wed, 19 Jan 2022 16:19:18 +0000 (17:19 +0100)
FEEval: replace usages of n_vectorization_actual

1  2 
include/deal.II/matrix_free/fe_evaluation.h

index a5c4e384ba006f6d9bb6669060aea10ba3a5c26a,f8964cbbb38821b5c2888b9c3c534b93f241d33e..c8ddf1ea76545b089e7de650b553a19af4ebc3c4
@@@ -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;

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.