Merge pull request #13252 from peterrum/fe_eval_quadrature_points

Move quadrature_points into FEEvalData

diff --cc include/deal.II/matrix_free/fe_evaluation.h
index c8ddf1ea76545b089e7de650b553a19af4ebc3c4,7e4bcae909b6ade95083de39c48c783a71572de6..6d9b4290a7e4c4c4e763ca388fddf4e17d1f1222
--- 1/include/deal.II/matrix_free/fe_evaluation.h
--- 2/include/deal.II/matrix_free/fe_evaluation.h
+++ b/include/deal.II/matrix_free/fe_evaluation.h
@@@ -6887,13 -6861,14 +6881,18 @@@ FEEvaluation<dim
    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;
  
+   if (this->mapping_data->quadrature_points.empty() == false)
+     this->quadrature_points =
+       &this->mapping_data->quadrature_points
+          [this->mapping_data->quadrature_point_offsets[this->cell]];
+ 
  #  ifdef DEBUG
    this->is_reinitialized           = true;
    this->dof_values_initialized     = false;