]> https://gitweb.dealii.org/ - dealii.git/commitdiff
FEEval: rename get_is_interior_face() to is_interior_face() 13311/head
authorPeter Munch <peterrmuench@gmail.com>
Mon, 31 Jan 2022 07:59:26 +0000 (08:59 +0100)
committerPeter Munch <peterrmuench@gmail.com>
Mon, 31 Jan 2022 07:59:26 +0000 (08:59 +0100)
include/deal.II/matrix_free/evaluation_kernels.h
include/deal.II/matrix_free/fe_evaluation.h
include/deal.II/matrix_free/fe_evaluation_data.h

index 25b2120a0cc6758aa55af42161d3ffbecd195250..d095c03ce23cd03d64d6533fe1e147aee70f8315 100644 (file)
@@ -3206,7 +3206,7 @@ namespace internal
     std::array<const unsigned int *, n_face_orientations> orientation = {};
 
     if (dim == 3 && n_face_orientations == n_lanes && !all_faces_are_same &&
-        fe_eval.get_is_interior_face() == 0)
+        fe_eval.is_interior_face() == 0)
       for (unsigned int v = 0; v < n_lanes; ++v)
         {
           // the loop breaks once an invalid_unsigned_int is hit for
@@ -3749,7 +3749,7 @@ namespace internal
 
       if (fe_eval.get_dof_access_index() ==
             MatrixFreeFunctions::DoFInfo::dof_access_cell &&
-          fe_eval.get_is_interior_face() == false)
+          fe_eval.is_interior_face() == false)
         fe_face_evaluation_process_and_io<VectorizedArrayType::size()>(
           p, n_components, evaluation_flag, src_ptr, sm_ptr, fe_eval, temp);
       else
@@ -3794,7 +3794,7 @@ namespace internal
         {
           if (fe_eval.get_dof_access_index() ==
                 MatrixFreeFunctions::DoFInfo::dof_access_cell &&
-              fe_eval.get_is_interior_face() == false)
+              fe_eval.is_interior_face() == false)
             {
               for (unsigned int v = 0; v < VectorizedArrayType::size(); ++v)
                 {
@@ -3972,7 +3972,7 @@ namespace internal
         {
           if (fe_eval.get_dof_access_index() ==
                 MatrixFreeFunctions::DoFInfo::dof_access_cell &&
-              fe_eval.get_is_interior_face() == false)
+              fe_eval.is_interior_face() == false)
             for (unsigned int v = 0; v < VectorizedArrayType::size(); ++v)
               {
                 // the loop breaks once an invalid_unsigned_int is hit for
@@ -4047,7 +4047,7 @@ namespace internal
 
       if (fe_eval.get_dof_access_index() ==
             MatrixFreeFunctions::DoFInfo::dof_access_cell &&
-          fe_eval.get_is_interior_face() == false)
+          fe_eval.is_interior_face() == false)
         fe_face_evaluation_process_and_io<VectorizedArrayType::size()>(
           p, n_components, integration_flag, dst_ptr, sm_ptr, fe_eval, temp);
       else
index 1e2c7760a275399395dfff63e1f363e4acc2ab62..0c899d8719552f40fa39fc89d51b3db650e2a3e4 100644 (file)
@@ -3874,8 +3874,8 @@ FEEvaluationBase<dim, n_components_, Number, is_face, VectorizedArrayType>::
       !(is_face &&
         this->dof_access_index ==
           internal::MatrixFreeFunctions::DoFInfo::dof_access_cell &&
-        this->is_interior_face == false) &&
-      !(!is_face && !this->is_interior_face))
+        this->is_interior_face() == false) &&
+      !(!is_face && !this->is_interior_face()))
     {
       const unsigned int dof_index =
         dof_indices_cont[this->cell * VectorizedArrayType::size()] +
@@ -3910,8 +3910,8 @@ FEEvaluationBase<dim, n_components_, Number, is_face, VectorizedArrayType>::
   const bool is_ecl =
     (this->dof_access_index ==
        internal::MatrixFreeFunctions::DoFInfo::dof_access_cell &&
-     this->is_interior_face == false) ||
-    (!is_face && !this->is_interior_face);
+     this->is_interior_face() == false) ||
+    (!is_face && !this->is_interior_face());
 
   if (vectors_sm[0] != nullptr)
     {
@@ -7563,7 +7563,7 @@ FEFaceEvaluation<dim,
 
   this->cell = face_index;
   this->dof_access_index =
-    this->is_interior_face ?
+    this->is_interior_face() ?
       internal::MatrixFreeFunctions::DoFInfo::dof_access_face_interior :
       internal::MatrixFreeFunctions::DoFInfo::dof_access_face_exterior;
   Assert(this->mapping_data != nullptr, ExcNotInitialized());
@@ -7572,7 +7572,7 @@ FEFaceEvaluation<dim,
         this->matrix_free->get_task_info().face_partition_data.back() &&
       face_index <
         this->matrix_free->get_task_info().boundary_partition_data.back())
-    Assert(this->is_interior_face,
+    Assert(this->is_interior_face(),
            ExcMessage(
              "Boundary faces do not have a neighbor. When looping over "
              "boundary faces use FEFaceEvaluation with the parameter "
@@ -7593,12 +7593,12 @@ FEFaceEvaluation<dim,
   this->J_value        = &this->mapping_data->JxW_values[offsets];
   this->normal_vectors = &this->mapping_data->normal_vectors[offsets];
   this->jacobian =
-    &this->mapping_data->jacobians[!this->is_interior_face][offsets];
+    &this->mapping_data->jacobians[!this->is_interior_face()][offsets];
   this->normal_x_jacobian =
     &this->mapping_data
-       ->normals_times_jacobians[!this->is_interior_face][offsets];
+       ->normals_times_jacobians[!this->is_interior_face()][offsets];
   this->jacobian_gradients =
-    this->mapping_data->jacobian_gradients[!this->is_interior_face].data() +
+    this->mapping_data->jacobian_gradients[!this->is_interior_face()].data() +
     offsets;
 
   if (this->mapping_data->quadrature_point_offsets.empty() == false)
@@ -7659,7 +7659,7 @@ FEFaceEvaluation<dim,
 
   constexpr unsigned int n_lanes = VectorizedArrayType::size();
 
-  if (this->is_interior_face == false)
+  if (this->is_interior_face() == false)
     {
       // for this case, we need to look into the FaceInfo field that collects
       // information from both sides of a face once for the global mesh, and
@@ -7741,13 +7741,13 @@ FEFaceEvaluation<dim,
                             .normal_vectors[offsets];
   this->jacobian = &this->matrix_free->get_mapping_info()
                       .face_data_by_cells[this->quad_no]
-                      .jacobians[!this->is_interior_face][offsets];
+                      .jacobians[!this->is_interior_face()][offsets];
   this->normal_x_jacobian =
     &this->matrix_free->get_mapping_info()
        .face_data_by_cells[this->quad_no]
-       .normals_times_jacobians[!this->is_interior_face][offsets];
+       .normals_times_jacobians[!this->is_interior_face()][offsets];
   this->jacobian_gradients =
-    this->mapping_data->jacobian_gradients[!this->is_interior_face].data() +
+    this->mapping_data->jacobian_gradients[!this->is_interior_face()].data() +
     offsets;
 
   if (this->matrix_free->get_mapping_info()
@@ -8198,7 +8198,7 @@ FEFaceEvaluation<dim,
 {
   Assert((this->dof_access_index ==
             internal::MatrixFreeFunctions::DoFInfo::dof_access_cell &&
-          this->is_interior_face == false) == false,
+          this->is_interior_face() == false) == false,
          ExcNotImplemented());
 
   const auto shared_vector_data = internal::get_shared_vector_data(
index 7de101a44ff575dc47360e85162b83cb9ed6124f..7ff407e8e1315d9fec4f54f9c95ec58dda116054 100644 (file)
@@ -469,7 +469,7 @@ public:
    * internal use.
    */
   bool
-  get_is_interior_face() const;
+  is_interior_face() const;
 
   /**
    * Return the id of the cells this FEEvaluation or FEFaceEvaluation is
@@ -800,7 +800,7 @@ protected:
    * the dof values should be read from the actual cell corresponding to the
    * interior face or the neighboring cell corresponding to the exterior face.
    */
-  bool is_interior_face;
+  bool interior_face;
 
   /**
    * Stores the index an FEFaceEvaluation object is currently pointing into
@@ -940,7 +940,7 @@ inline FEEvaluationData<dim, Number, is_face>::FEEvaluationData(
   , gradients_quad_submitted(false)
 #  endif
   , cell(numbers::invalid_unsigned_int)
-  , is_interior_face(is_interior_face)
+  , interior_face(is_interior_face)
   , dof_access_index(
       is_face ?
         (is_interior_face ?
@@ -980,7 +980,7 @@ inline FEEvaluationData<dim, Number, is_face>::FEEvaluationData(
   , quadrature_weights(nullptr)
   , cell(0)
   , cell_type(internal::MatrixFreeFunctions::general)
-  , is_interior_face(true)
+  , interior_face(true)
   , dof_access_index(internal::MatrixFreeFunctions::DoFInfo::dof_access_cell)
   , mapped_geometry(mapped_geometry)
   , is_reinitialized(false)
@@ -1029,11 +1029,11 @@ FEEvaluationData<dim, Number, is_face>::operator=(const FEEvaluationData &other)
   gradients_quad_submitted   = false;
 #  endif
 
-  cell             = numbers::invalid_unsigned_int;
-  is_interior_face = other.is_interior_face;
+  cell          = numbers::invalid_unsigned_int;
+  interior_face = other.is_interior_face();
   dof_access_index =
     is_face ?
-      (is_interior_face ?
+      (is_interior_face() ?
          internal::MatrixFreeFunctions::DoFInfo::dof_access_face_interior :
          internal::MatrixFreeFunctions::DoFInfo::dof_access_face_exterior) :
       internal::MatrixFreeFunctions::DoFInfo::dof_access_cell;
@@ -1097,8 +1097,8 @@ FEEvaluationData<dim, Number, is_face>::reinit_face(
          ExcMessage("Faces can only be set if the is_face template parameter "
                     "is true"));
   face_numbers[0] =
-    (is_interior_face ? face.interior_face_no : face.exterior_face_no);
-  subface_index = is_interior_face == true ?
+    (is_interior_face() ? face.interior_face_no : face.exterior_face_no);
+  subface_index = is_interior_face() == true ?
                     GeometryInfo<dim>::max_children_per_cell :
                     face.subface_index;
 
@@ -1107,11 +1107,11 @@ FEEvaluationData<dim, Number, is_face>::reinit_face(
   // standard-orientation else copy the first three bits
   // (which is equivalent to modulo 8). See also the documentation of
   // internal::MatrixFreeFunctions::FaceToCellTopology::face_orientation.
-  face_orientations[0] = (is_interior_face == (face.face_orientation >= 8)) ?
+  face_orientations[0] = (is_interior_face() == (face.face_orientation >= 8)) ?
                            (face.face_orientation % 8) :
                            0;
 
-  if (is_interior_face)
+  if (is_interior_face())
     cell_ids = face.cells_interior;
   else
     cell_ids = face.cells_exterior;
@@ -1431,7 +1431,7 @@ FEEvaluationData<dim, Number, is_face>::get_face_no(const unsigned int v) const
   Assert(v == 0 || (cell != numbers::invalid_unsigned_int &&
                     dof_access_index ==
                       internal::MatrixFreeFunctions::DoFInfo::dof_access_cell &&
-                    is_interior_face == false),
+                    is_interior_face() == false),
          ExcMessage("All face numbers can only be queried for ECL at exterior "
                     "faces. Use get_face_no() in other cases."));
 
@@ -1458,7 +1458,7 @@ FEEvaluationData<dim, Number, is_face>::get_face_orientation(
   Assert(v == 0 || (cell != numbers::invalid_unsigned_int &&
                     dof_access_index ==
                       internal::MatrixFreeFunctions::DoFInfo::dof_access_cell &&
-                    is_interior_face == false),
+                    is_interior_face() == false),
          ExcMessage("All face numbers can only be queried for ECL at exterior "
                     "faces. Use get_face_no() in other cases."));
 
@@ -1478,9 +1478,9 @@ FEEvaluationData<dim, Number, is_face>::get_dof_access_index() const
 
 template <int dim, typename Number, bool is_face>
 inline bool
-FEEvaluationData<dim, Number, is_face>::get_is_interior_face() const
+FEEvaluationData<dim, Number, is_face>::is_interior_face() const
 {
-  return is_interior_face;
+  return interior_face;
 }
 
 

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.