Merge pull request #13253 from peterrum/fe_eval_jacobian_gradients

Move jacobian_gradients into FEEvalData

diff --cc include/deal.II/matrix_free/fe_evaluation.h
index 2898250ccba0d4746502b1a65c404bd3b23ec2ae,7b1d2a1482f6d6f1c6acc92928ad720dba99563e..a5c4e384ba006f6d9bb6669060aea10ba3a5c26a
--- 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
@@@ -7618,9 -7631,11 +7611,12 @@@ FEFaceEvaluation<dim
    this->normal_x_jacobian =
      &this->mapping_data
         ->normals_times_jacobians[!this->is_interior_face][offsets];
+   this->jacobian_gradients =
+     this->mapping_data->jacobian_gradients[!this->is_interior_face].data() +
+     offsets;
  
  #  ifdef DEBUG
 +  this->is_reinitialized           = true;
    this->dof_values_initialized     = false;
    this->values_quad_initialized    = false;
    this->gradients_quad_initialized = false;
@@@ -7755,9 -7770,11 +7751,12 @@@ FEFaceEvaluation<dim
      &this->matrix_free->get_mapping_info()
         .face_data_by_cells[this->quad_no]
         .normals_times_jacobians[!this->is_interior_face][offsets];
+   this->jacobian_gradients =
+     this->mapping_data->jacobian_gradients[!this->is_interior_face].data() +
+     offsets;
  
  #  ifdef DEBUG
 +  this->is_reinitialized           = true;
    this->dof_values_initialized     = false;
    this->values_quad_initialized    = false;
    this->gradients_quad_initialized = false;

diff --cc include/deal.II/matrix_free/fe_evaluation_data.h
Simple merge