]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Merge pull request #13253 from peterrum/fe_eval_jacobian_gradients
authorMartin Kronbichler <martin.kronbichler@it.uu.se>
Tue, 18 Jan 2022 13:45:12 +0000 (14:45 +0100)
committerGitHub <noreply@github.com>
Tue, 18 Jan 2022 13:45:12 +0000 (14:45 +0100)
Move jacobian_gradients into FEEvalData

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

index 2898250ccba0d4746502b1a65c404bd3b23ec2ae,7b1d2a1482f6d6f1c6acc92928ad720dba99563e..a5c4e384ba006f6d9bb6669060aea10ba3a5c26a
@@@ -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;

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.