]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Rename
authorPeter Munch <peterrmuench@gmail.com>
Tue, 3 Aug 2021 08:12:26 +0000 (10:12 +0200)
committerPeter Munch <peterrmuench@gmail.com>
Tue, 3 Aug 2021 08:12:26 +0000 (10:12 +0200)
14 files changed:
include/deal.II/matrix_free/cuda_fe_evaluation.h
include/deal.II/matrix_free/cuda_matrix_free.h
include/deal.II/matrix_free/cuda_matrix_free.templates.h
include/deal.II/matrix_free/dof_info.h
include/deal.II/matrix_free/dof_info.templates.h
include/deal.II/matrix_free/evaluation_kernels.h
include/deal.II/matrix_free/evaluation_template_factory.h
include/deal.II/matrix_free/evaluation_template_factory.templates.h
include/deal.II/matrix_free/fe_evaluation.h
include/deal.II/matrix_free/hanging_nodes_internal.h
include/deal.II/matrix_free/matrix_free.templates.h
include/deal.II/matrix_free/tools.h
source/matrix_free/dof_info.cc
tests/matrix_free/hanging_node_kernels_01.cc

index 06efd8e5ba84c7038c0c1683e54902c83cc8996b..965ff1c27ef67c76285edb279fabeb52cd7ff9c7 100644 (file)
@@ -238,7 +238,7 @@ namespace CUDAWrappers
     unsigned int             padding_length;
     const unsigned int       mf_object_id;
 
-    const dealii::internal::MatrixFreeFunctions::ConstraintTypes
+    const dealii::internal::MatrixFreeFunctions::ConstraintKinds
       constraint_mask;
 
     const bool use_coloring;
index 3557fd4218d21f9a8b0bebc971a8fbe3381df0a7..2fa34f40c161f700dc3fb61927fe801b08e3d3af 100644 (file)
@@ -212,7 +212,7 @@ namespace CUDAWrappers
       /**
        * Mask deciding where constraints are set on a given cell.
        */
-      dealii::internal::MatrixFreeFunctions::ConstraintTypes *constraint_mask;
+      dealii::internal::MatrixFreeFunctions::ConstraintKinds *constraint_mask;
 
       /**
        * If true, use graph coloring has been used and we can simply add into
@@ -591,7 +591,7 @@ namespace CUDAWrappers
     /**
      * Mask deciding where constraints are set on a given cell.
      */
-    std::vector<dealii::internal::MatrixFreeFunctions::ConstraintTypes *>
+    std::vector<dealii::internal::MatrixFreeFunctions::ConstraintKinds *>
       constraint_mask;
 
     /**
@@ -818,7 +818,7 @@ namespace CUDAWrappers
     /**
      * Mask deciding where constraints are set on a given cell.
      */
-    std::vector<dealii::internal::MatrixFreeFunctions::ConstraintTypes>
+    std::vector<dealii::internal::MatrixFreeFunctions::ConstraintKinds>
       constraint_mask;
 
     /**
index 1f6750010c14032159a2cc2c70a8ee4d0e44aa49..bffcdab551b78abd26f9384f44f246849d6d109a 100644 (file)
@@ -229,7 +229,7 @@ namespace CUDAWrappers
       std::vector<Point<dim, Number>>      q_points_host;
       std::vector<Number>                  JxW_host;
       std::vector<Number>                  inv_jacobian_host;
-      std::vector<dealii::internal::MatrixFreeFunctions::ConstraintTypes>
+      std::vector<dealii::internal::MatrixFreeFunctions::ConstraintKinds>
         constraint_mask_host;
       // Local buffer
       std::vector<types::global_dof_index> local_dof_indices;
@@ -378,7 +378,7 @@ namespace CUDAWrappers
       for (unsigned int i = 0; i < dofs_per_cell; ++i)
         lexicographic_dof_indices[i] = local_dof_indices[lexicographic_inv[i]];
 
-      const ArrayView<dealii::internal::MatrixFreeFunctions::ConstraintTypes>
+      const ArrayView<dealii::internal::MatrixFreeFunctions::ConstraintKinds>
         cell_id_view(constraint_mask_host[cell_id]);
 
       hanging_nodes.setup_constraints(cell,
@@ -496,7 +496,7 @@ namespace CUDAWrappers
 
       alloc_and_copy(
         &data->constraint_mask[color],
-        ArrayView<const dealii::internal::MatrixFreeFunctions::ConstraintTypes>(
+        ArrayView<const dealii::internal::MatrixFreeFunctions::ConstraintKinds>(
           constraint_mask_host.data(), constraint_mask_host.size()),
         n_cells);
     }
index a465e8537bc149385ab2a615b6c66ecd0c21c561..d4b69bce4c4285e71cd382b4ecf2322f8e9ca160 100644 (file)
@@ -528,7 +528,7 @@ namespace internal
        * Masks indicating for each cell and component if the optimized
        * hanging-node constraint is applicable and if yes which type.
        */
-      std::vector<ConstraintTypes> hanging_node_constraint_masks;
+      std::vector<ConstraintKinds> hanging_node_constraint_masks;
 
       /**
        * This variable describes the position of constraints in terms of the
index 6e019507dbcb941c57e3788bbf48464f77b33dc5..ffb59ed4c3770154e16301dbe2d45e2ed9cd6590 100644 (file)
@@ -283,7 +283,7 @@ namespace internal
       const TriaIterator<DoFCellAccessor<dim, dim, false>> &cell,
       std::vector<types::global_dof_index> &                dof_indices)
     {
-      const ArrayView<ConstraintTypes> mask_view(
+      const ArrayView<ConstraintKinds> mask_view(
         hanging_node_constraint_masks.data() +
           cell_number * cell->get_fe().n_components(),
         cell->get_fe().n_components());
index 804fbb9d6c2bb4a4bdb3201f8e253f5f32f6f24e..d713a7bc1754d07bca0e64078206c1b6f39beb08 100644 (file)
@@ -4434,7 +4434,7 @@ namespace internal
                                    is_face,
                                    Number> &fe_eval,
         const bool                          transpose,
-        const std::array<MatrixFreeFunctions::ConstraintTypes, Number::size()>
+        const std::array<MatrixFreeFunctions::ConstraintKinds, Number::size()>
           &     c_mask,
         Number *values)
     {
@@ -4655,7 +4655,7 @@ namespace internal
                                             typename Number::value_type,
                                             is_face,
                                             Number> &fe_eval,
-                 const std::array<MatrixFreeFunctions::ConstraintTypes,
+                 const std::array<MatrixFreeFunctions::ConstraintKinds,
                                   Number::size()> &  constraint_mask,
                  Number *                            values)
     {
@@ -4672,7 +4672,7 @@ namespace internal
       };
 
       const auto not_set = [](const auto a, const auto b) -> bool {
-        return (a & b) == MatrixFreeFunctions::ConstraintTypes::unconstrained;
+        return (a & b) == MatrixFreeFunctions::ConstraintKinds::unconstrained;
       };
 
       const unsigned int points = given_degree + 1;
@@ -4683,20 +4683,20 @@ namespace internal
             {
               const auto mask = constraint_mask[v];
 
-              if (mask == MatrixFreeFunctions::ConstraintTypes::unconstrained)
+              if (mask == MatrixFreeFunctions::ConstraintKinds::unconstrained)
                 continue;
 
               if (dim == 2) // 2D: only faces
                 {
                   // direction 0:
-                  if ((mask & MatrixFreeFunctions::ConstraintTypes::face_y) !=
-                      MatrixFreeFunctions::ConstraintTypes::unconstrained)
+                  if ((mask & MatrixFreeFunctions::ConstraintKinds::face_y) !=
+                      MatrixFreeFunctions::ConstraintKinds::unconstrained)
                     {
                       const bool is_subface_0 =
-                        (mask & MatrixFreeFunctions::ConstraintTypes::type_x) !=
-                        MatrixFreeFunctions::ConstraintTypes::unconstrained;
+                        (mask & MatrixFreeFunctions::ConstraintKinds::type_x) !=
+                        MatrixFreeFunctions::ConstraintKinds::unconstrained;
                       if (is_set(mask,
-                                 MatrixFreeFunctions::ConstraintTypes::type_y))
+                                 MatrixFreeFunctions::ConstraintKinds::type_y))
                         interpolate_2D<fe_degree, 2, transpose>(
                           given_degree,
                           is_subface_0,
@@ -4713,14 +4713,14 @@ namespace internal
                     }
 
                   // direction 1:
-                  if ((mask & MatrixFreeFunctions::ConstraintTypes::face_x) !=
-                      MatrixFreeFunctions::ConstraintTypes::unconstrained)
+                  if ((mask & MatrixFreeFunctions::ConstraintKinds::face_x) !=
+                      MatrixFreeFunctions::ConstraintKinds::unconstrained)
                     {
                       const bool is_subface_0 =
-                        (mask & MatrixFreeFunctions::ConstraintTypes::type_y) !=
-                        MatrixFreeFunctions::ConstraintTypes::unconstrained;
+                        (mask & MatrixFreeFunctions::ConstraintKinds::type_y) !=
+                        MatrixFreeFunctions::ConstraintKinds::unconstrained;
                       if (is_set(mask,
-                                 MatrixFreeFunctions::ConstraintTypes::type_x))
+                                 MatrixFreeFunctions::ConstraintKinds::type_x))
                         interpolate_2D<fe_degree, 0, transpose>(
                           given_degree,
                           is_subface_0,
@@ -4750,109 +4750,109 @@ namespace internal
 
                   const bool is_face_0 =
                     is_set(mask,
-                           MatrixFreeFunctions::ConstraintTypes::face_x) &&
-                    is_set(mask, MatrixFreeFunctions::ConstraintTypes::type_x);
+                           MatrixFreeFunctions::ConstraintKinds::face_x) &&
+                    is_set(mask, MatrixFreeFunctions::ConstraintKinds::type_x);
                   const bool is_face_1 =
                     is_set(mask,
-                           MatrixFreeFunctions::ConstraintTypes::face_x) &&
-                    not_set(mask, MatrixFreeFunctions::ConstraintTypes::type_x);
+                           MatrixFreeFunctions::ConstraintKinds::face_x) &&
+                    not_set(mask, MatrixFreeFunctions::ConstraintKinds::type_x);
                   const bool is_face_2 =
                     is_set(mask,
-                           MatrixFreeFunctions::ConstraintTypes::face_y) &&
-                    is_set(mask, MatrixFreeFunctions::ConstraintTypes::type_y);
+                           MatrixFreeFunctions::ConstraintKinds::face_y) &&
+                    is_set(mask, MatrixFreeFunctions::ConstraintKinds::type_y);
                   const bool is_face_3 =
                     is_set(mask,
-                           MatrixFreeFunctions::ConstraintTypes::face_y) &&
-                    not_set(mask, MatrixFreeFunctions::ConstraintTypes::type_y);
+                           MatrixFreeFunctions::ConstraintKinds::face_y) &&
+                    not_set(mask, MatrixFreeFunctions::ConstraintKinds::type_y);
                   const bool is_face_4 =
                     is_set(mask,
-                           MatrixFreeFunctions::ConstraintTypes::face_z) &&
-                    is_set(mask, MatrixFreeFunctions::ConstraintTypes::type_z);
+                           MatrixFreeFunctions::ConstraintKinds::face_z) &&
+                    is_set(mask, MatrixFreeFunctions::ConstraintKinds::type_z);
                   const bool is_face_5 =
                     is_set(mask,
-                           MatrixFreeFunctions::ConstraintTypes::face_z) &&
-                    not_set(mask, MatrixFreeFunctions::ConstraintTypes::type_z);
+                           MatrixFreeFunctions::ConstraintKinds::face_z) &&
+                    not_set(mask, MatrixFreeFunctions::ConstraintKinds::type_z);
 
                   const bool is_edge_2 =
                     is_set(mask,
-                           MatrixFreeFunctions::ConstraintTypes::edge_yz) &&
+                           MatrixFreeFunctions::ConstraintKinds::edge_yz) &&
                     is_set(mask,
-                           MatrixFreeFunctions::ConstraintTypes::type_y) &&
-                    is_set(mask, MatrixFreeFunctions::ConstraintTypes::type_z);
+                           MatrixFreeFunctions::ConstraintKinds::type_y) &&
+                    is_set(mask, MatrixFreeFunctions::ConstraintKinds::type_z);
                   const bool is_edge_3 =
                     is_set(mask,
-                           MatrixFreeFunctions::ConstraintTypes::edge_yz) &&
+                           MatrixFreeFunctions::ConstraintKinds::edge_yz) &&
                     not_set(mask,
-                            MatrixFreeFunctions::ConstraintTypes::type_y) &&
-                    is_set(mask, MatrixFreeFunctions::ConstraintTypes::type_z);
+                            MatrixFreeFunctions::ConstraintKinds::type_y) &&
+                    is_set(mask, MatrixFreeFunctions::ConstraintKinds::type_z);
                   const bool is_edge_6 =
                     is_set(mask,
-                           MatrixFreeFunctions::ConstraintTypes::edge_yz) &&
+                           MatrixFreeFunctions::ConstraintKinds::edge_yz) &&
                     is_set(mask,
-                           MatrixFreeFunctions::ConstraintTypes::type_y) &&
-                    not_set(mask, MatrixFreeFunctions::ConstraintTypes::type_z);
+                           MatrixFreeFunctions::ConstraintKinds::type_y) &&
+                    not_set(mask, MatrixFreeFunctions::ConstraintKinds::type_z);
                   const bool is_edge_7 =
                     is_set(mask,
-                           MatrixFreeFunctions::ConstraintTypes::edge_yz) &&
+                           MatrixFreeFunctions::ConstraintKinds::edge_yz) &&
                     not_set(mask,
-                            MatrixFreeFunctions::ConstraintTypes::type_y) &&
-                    not_set(mask, MatrixFreeFunctions::ConstraintTypes::type_z);
+                            MatrixFreeFunctions::ConstraintKinds::type_y) &&
+                    not_set(mask, MatrixFreeFunctions::ConstraintKinds::type_z);
 
                   const bool is_edge_0 =
                     is_set(mask,
-                           MatrixFreeFunctions::ConstraintTypes::edge_zx) &&
+                           MatrixFreeFunctions::ConstraintKinds::edge_zx) &&
                     is_set(mask,
-                           MatrixFreeFunctions::ConstraintTypes::type_x) &&
-                    is_set(mask, MatrixFreeFunctions::ConstraintTypes::type_z);
+                           MatrixFreeFunctions::ConstraintKinds::type_x) &&
+                    is_set(mask, MatrixFreeFunctions::ConstraintKinds::type_z);
                   const bool is_edge_1 =
                     is_set(mask,
-                           MatrixFreeFunctions::ConstraintTypes::edge_zx) &&
+                           MatrixFreeFunctions::ConstraintKinds::edge_zx) &&
                     not_set(mask,
-                            MatrixFreeFunctions::ConstraintTypes::type_x) &&
-                    is_set(mask, MatrixFreeFunctions::ConstraintTypes::type_z);
+                            MatrixFreeFunctions::ConstraintKinds::type_x) &&
+                    is_set(mask, MatrixFreeFunctions::ConstraintKinds::type_z);
                   const bool is_edge_4 =
                     is_set(mask,
-                           MatrixFreeFunctions::ConstraintTypes::edge_zx) &&
+                           MatrixFreeFunctions::ConstraintKinds::edge_zx) &&
                     is_set(mask,
-                           MatrixFreeFunctions::ConstraintTypes::type_x) &&
-                    not_set(mask, MatrixFreeFunctions::ConstraintTypes::type_z);
+                           MatrixFreeFunctions::ConstraintKinds::type_x) &&
+                    not_set(mask, MatrixFreeFunctions::ConstraintKinds::type_z);
                   const bool is_edge_5 =
                     is_set(mask,
-                           MatrixFreeFunctions::ConstraintTypes::edge_zx) &&
+                           MatrixFreeFunctions::ConstraintKinds::edge_zx) &&
                     not_set(mask,
-                            MatrixFreeFunctions::ConstraintTypes::type_x) &&
-                    not_set(mask, MatrixFreeFunctions::ConstraintTypes::type_z);
+                            MatrixFreeFunctions::ConstraintKinds::type_x) &&
+                    not_set(mask, MatrixFreeFunctions::ConstraintKinds::type_z);
 
                   const bool is_edge_8 =
                     is_set(mask,
-                           MatrixFreeFunctions::ConstraintTypes::edge_xy) &&
+                           MatrixFreeFunctions::ConstraintKinds::edge_xy) &&
                     is_set(mask,
-                           MatrixFreeFunctions::ConstraintTypes::type_x) &&
-                    is_set(mask, MatrixFreeFunctions::ConstraintTypes::type_y);
+                           MatrixFreeFunctions::ConstraintKinds::type_x) &&
+                    is_set(mask, MatrixFreeFunctions::ConstraintKinds::type_y);
                   const bool is_edge_9 =
                     is_set(mask,
-                           MatrixFreeFunctions::ConstraintTypes::edge_xy) &&
+                           MatrixFreeFunctions::ConstraintKinds::edge_xy) &&
                     not_set(mask,
-                            MatrixFreeFunctions::ConstraintTypes::type_x) &&
-                    is_set(mask, MatrixFreeFunctions::ConstraintTypes::type_y);
+                            MatrixFreeFunctions::ConstraintKinds::type_x) &&
+                    is_set(mask, MatrixFreeFunctions::ConstraintKinds::type_y);
                   const bool is_edge_10 =
                     is_set(mask,
-                           MatrixFreeFunctions::ConstraintTypes::edge_xy) &&
+                           MatrixFreeFunctions::ConstraintKinds::edge_xy) &&
                     is_set(mask,
-                           MatrixFreeFunctions::ConstraintTypes::type_x) &&
-                    not_set(mask, MatrixFreeFunctions::ConstraintTypes::type_y);
+                           MatrixFreeFunctions::ConstraintKinds::type_x) &&
+                    not_set(mask, MatrixFreeFunctions::ConstraintKinds::type_y);
                   const bool is_edge_11 =
                     is_set(mask,
-                           MatrixFreeFunctions::ConstraintTypes::edge_xy) &&
+                           MatrixFreeFunctions::ConstraintKinds::edge_xy) &&
                     not_set(mask,
-                            MatrixFreeFunctions::ConstraintTypes::type_x) &&
-                    not_set(mask, MatrixFreeFunctions::ConstraintTypes::type_y);
+                            MatrixFreeFunctions::ConstraintKinds::type_x) &&
+                    not_set(mask, MatrixFreeFunctions::ConstraintKinds::type_y);
 
                   // direction 0:
                   {
                     const bool is_subface_0 =
-                      (mask & MatrixFreeFunctions::ConstraintTypes::type_x) !=
-                      MatrixFreeFunctions::ConstraintTypes::unconstrained;
+                      (mask & MatrixFreeFunctions::ConstraintKinds::type_x) !=
+                      MatrixFreeFunctions::ConstraintKinds::unconstrained;
 
                     // ... faces
                     if (is_face_2)
@@ -4926,8 +4926,8 @@ namespace internal
                   // direction 1:
                   {
                     const bool is_subface_0 =
-                      (mask & MatrixFreeFunctions::ConstraintTypes::type_y) !=
-                      MatrixFreeFunctions::ConstraintTypes::unconstrained;
+                      (mask & MatrixFreeFunctions::ConstraintKinds::type_y) !=
+                      MatrixFreeFunctions::ConstraintKinds::unconstrained;
 
                     // ... faces
                     if (is_face_0)
@@ -5001,8 +5001,8 @@ namespace internal
                   // direction 2:
                   {
                     const bool is_subface_0 =
-                      (mask & MatrixFreeFunctions::ConstraintTypes::type_z) !=
-                      MatrixFreeFunctions::ConstraintTypes::unconstrained;
+                      (mask & MatrixFreeFunctions::ConstraintKinds::type_z) !=
+                      MatrixFreeFunctions::ConstraintKinds::unconstrained;
 
                     // ... faces
                     if (is_face_0)
index ca38ec24c531b6743d8134399895f0bebdb7eacc..59204dfa0b83da13bdff79dbff7fccd6f57d5cd0 100644 (file)
@@ -202,7 +202,7 @@ namespace internal
           const FEEvaluationBaseData<dim, Number, false, VectorizedArrayType>
             &                                            fe_eval,
           const bool                                     transpose,
-          const std::array<MatrixFreeFunctions::ConstraintTypes,
+          const std::array<MatrixFreeFunctions::ConstraintKinds,
                            VectorizedArrayType::size()> &c_mask,
           VectorizedArrayType *                          values);
 
@@ -217,7 +217,7 @@ namespace internal
           const FEEvaluationBaseData<dim, Number, true, VectorizedArrayType>
             &                                            fe_eval,
           const bool                                     transpose,
-          const std::array<MatrixFreeFunctions::ConstraintTypes,
+          const std::array<MatrixFreeFunctions::ConstraintKinds,
                            VectorizedArrayType::size()> &c_mask,
           VectorizedArrayType *                          values);
   };
index ce74b34fb8ae62895b457b16f90593e9b8de7874..74e5086c88b9b8e4e046114b953b6b444d93d39b 100644 (file)
@@ -412,7 +412,7 @@ namespace internal
     const FEEvaluationBaseData<dim, Number, false, VectorizedArrayType>
       &                                            fe_eval,
     const bool                                     transpose,
-    const std::array<MatrixFreeFunctions::ConstraintTypes,
+    const std::array<MatrixFreeFunctions::ConstraintKinds,
                      VectorizedArrayType::size()> &c_mask,
     VectorizedArrayType *                          values)
   {
@@ -437,7 +437,7 @@ namespace internal
     const unsigned int fe_degree,
     const FEEvaluationBaseData<dim, Number, true, VectorizedArrayType> &fe_eval,
     const bool                                     transpose,
-    const std::array<MatrixFreeFunctions::ConstraintTypes,
+    const std::array<MatrixFreeFunctions::ConstraintKinds,
                      VectorizedArrayType::size()> &c_mask,
     VectorizedArrayType *                          values)
   {
index 4f7654d2daaed086c8b2677b40bf904aa0cda781..6ef074383c6912270756b014a503a6f334040cb0 100644 (file)
@@ -4475,7 +4475,7 @@ FEEvaluationBase<dim, n_components_, Number, is_face, VectorizedArrayType>::
                 ->hanging_node_constraint_masks[(this->cell * n_lanes + v) *
                                                   n_fe_components +
                                                 first_selected_component] !=
-              internal::MatrixFreeFunctions::ConstraintTypes::unconstrained)
+              internal::MatrixFreeFunctions::ConstraintKinds::unconstrained)
           has_hn_constraints = true;
     }
 
@@ -4582,7 +4582,7 @@ FEEvaluationBase<dim, n_components_, Number, is_face, VectorizedArrayType>::
                   ->hanging_node_constraint_masks[(this->cell * n_lanes + v) *
                                                     n_fe_components +
                                                   first_selected_component] !=
-                internal::MatrixFreeFunctions::ConstraintTypes::unconstrained)
+                internal::MatrixFreeFunctions::ConstraintKinds::unconstrained)
             has_hn_constraints = true;
 
           Assert(my_index_start[n_components_read].first ==
@@ -4658,7 +4658,7 @@ FEEvaluationBase<dim, n_components_, Number, is_face, VectorizedArrayType>::
                .second ||
            (this->dof_info->hanging_node_constraint_masks.size() > 0 &&
             this->dof_info->hanging_node_constraint_masks[cell_dof_index] !=
-              internal::MatrixFreeFunctions::ConstraintTypes::unconstrained)))
+              internal::MatrixFreeFunctions::ConstraintKinds::unconstrained)))
         {
           Assert(this->dof_info->row_starts_plain_indices[cell_index] !=
                    numbers::invalid_unsigned_int,
@@ -5315,7 +5315,7 @@ FEEvaluationBase<dim, n_components_, Number, is_face, VectorizedArrayType>::
       ->n_vectorization_lanes_filled[this->dof_access_index][this->cell];
 
   constexpr unsigned int n_lanes = VectorizedArrayType::size();
-  std::array<internal::MatrixFreeFunctions::ConstraintTypes, n_lanes>
+  std::array<internal::MatrixFreeFunctions::ConstraintKinds, n_lanes>
     constraint_mask;
 
   bool hn_available = false;
@@ -5350,7 +5350,7 @@ FEEvaluationBase<dim, n_components_, Number, is_face, VectorizedArrayType>::
       constraint_mask[v] = mask;
 
       hn_available |=
-        (mask != internal::MatrixFreeFunctions::ConstraintTypes::unconstrained);
+        (mask != internal::MatrixFreeFunctions::ConstraintKinds::unconstrained);
     }
 
   if (hn_available == false)
@@ -5358,7 +5358,7 @@ FEEvaluationBase<dim, n_components_, Number, is_face, VectorizedArrayType>::
 
   for (unsigned int v = n_vectorization_actual; v < n_lanes; ++v)
     constraint_mask[v] =
-      internal::MatrixFreeFunctions::ConstraintTypes::unconstrained;
+      internal::MatrixFreeFunctions::ConstraintKinds::unconstrained;
 
   internal::FEEvaluationHangingNodesFactory<dim, Number, VectorizedArrayType>::
     apply(n_components,
index 156f7e9c2054b20a19c7c1c3e4f45aeaf10faaf3..1aed9095fe0d10541010a3c7ddcc30c94d617123 100644 (file)
@@ -50,7 +50,7 @@ namespace internal
      * determine which face is constrained. For example, in 2D, if
      * face_x and type are set, then x = 0 is constrained.
      */
-    enum class ConstraintTypes : std::uint16_t
+    enum class ConstraintKinds : std::uint16_t
     {
       unconstrained = 0,
 
@@ -78,10 +78,10 @@ namespace internal
      * would be an integer which would in turn trigger a compiler warning when
      * we tried to assign it to an object of type UpdateFlags.
      */
-    DEAL_II_CUDA_HOST_DEV inline ConstraintTypes
-    operator|(const ConstraintTypes f1, const ConstraintTypes f2)
+    DEAL_II_CUDA_HOST_DEV inline ConstraintKinds
+    operator|(const ConstraintKinds f1, const ConstraintKinds f2)
     {
-      return static_cast<ConstraintTypes>(static_cast<std::uint16_t>(f1) |
+      return static_cast<ConstraintKinds>(static_cast<std::uint16_t>(f1) |
                                           static_cast<std::uint16_t>(f2));
     }
 
@@ -91,8 +91,8 @@ namespace internal
      * Global operator which sets the bits from the second argument also in the
      * first one.
      */
-    DEAL_II_CUDA_HOST_DEV inline ConstraintTypes &
-    operator|=(ConstraintTypes &f1, const ConstraintTypes f2)
+    DEAL_II_CUDA_HOST_DEV inline ConstraintKinds &
+    operator|=(ConstraintKinds &f1, const ConstraintKinds f2)
     {
       f1 = f1 | f2;
       return f1;
@@ -104,7 +104,7 @@ namespace internal
      * Global operator which checks inequality.
      */
     DEAL_II_CUDA_HOST_DEV inline bool
-    operator!=(const ConstraintTypes f1, const ConstraintTypes f2)
+    operator!=(const ConstraintKinds f1, const ConstraintKinds f2)
     {
       return static_cast<std::uint16_t>(f1) != static_cast<std::uint16_t>(f2);
     }
@@ -116,7 +116,7 @@ namespace internal
      * second.
      */
     DEAL_II_CUDA_HOST_DEV inline bool
-    operator<(const ConstraintTypes f1, const ConstraintTypes f2)
+    operator<(const ConstraintKinds f1, const ConstraintKinds f2)
     {
       return static_cast<std::uint16_t>(f1) < static_cast<std::uint16_t>(f2);
     }
@@ -126,10 +126,10 @@ namespace internal
     /**
      * Global operator which performs a binary and for the provided arguments.
      */
-    DEAL_II_CUDA_HOST_DEV inline ConstraintTypes
-    operator&(const ConstraintTypes f1, const ConstraintTypes f2)
+    DEAL_II_CUDA_HOST_DEV inline ConstraintKinds
+    operator&(const ConstraintKinds f1, const ConstraintKinds f2)
     {
-      return static_cast<ConstraintTypes>(static_cast<std::uint16_t>(f1) &
+      return static_cast<ConstraintKinds>(static_cast<std::uint16_t>(f1) &
                                           static_cast<std::uint16_t>(f2));
     }
 
@@ -162,7 +162,7 @@ namespace internal
         const std::shared_ptr<const Utilities::MPI::Partitioner> &partitioner,
         const std::vector<unsigned int> &     lexicographic_mapping,
         std::vector<types::global_dof_index> &dof_indices,
-        const ArrayView<ConstraintTypes> &    mask) const;
+        const ArrayView<ConstraintKinds> &    mask) const;
 
     private:
       /**
@@ -305,7 +305,7 @@ namespace internal
       const std::shared_ptr<const Utilities::MPI::Partitioner> &partitioner,
       const std::vector<unsigned int> &     lexicographic_mapping,
       std::vector<types::global_dof_index> &dof_indices,
-      const ArrayView<ConstraintTypes> &    masks) const
+      const ArrayView<ConstraintKinds> &    masks) const
     {
       bool cell_has_hanging_node_constraints = false;
 
@@ -344,7 +344,7 @@ namespace internal
              ++c, ++comp)
           {
             auto &mask = masks[comp];
-            mask       = ConstraintTypes::unconstrained;
+            mask       = ConstraintKinds::unconstrained;
 
             const auto &fe = cell->get_fe().base_element(base_element_index);
 
@@ -410,19 +410,19 @@ namespace internal
                           {
                             if (face < 2)
                               {
-                                mask |= ConstraintTypes::face_x;
+                                mask |= ConstraintKinds::face_x;
                                 if (face == 0)
-                                  mask |= ConstraintTypes::type_x;
+                                  mask |= ConstraintKinds::type_x;
                                 if (subface == 0)
-                                  mask |= ConstraintTypes::type_y;
+                                  mask |= ConstraintKinds::type_y;
                               }
                             else
                               {
-                                mask |= ConstraintTypes::face_y;
+                                mask |= ConstraintKinds::face_y;
                                 if (face == 2)
-                                  mask |= ConstraintTypes::type_y;
+                                  mask |= ConstraintKinds::type_y;
                                 if (subface == 0)
-                                  mask |= ConstraintTypes::type_x;
+                                  mask |= ConstraintKinds::type_x;
                               }
 
                             // Reorder neighbor_dofs and copy into faceth face
@@ -470,35 +470,35 @@ namespace internal
                             // YZ-plane
                             if (face < 2)
                               {
-                                mask |= ConstraintTypes::face_x;
+                                mask |= ConstraintKinds::face_x;
                                 if (face == 0)
-                                  mask |= ConstraintTypes::type_x;
+                                  mask |= ConstraintKinds::type_x;
                                 if (subface % 2 == 0)
-                                  mask |= ConstraintTypes::type_y;
+                                  mask |= ConstraintKinds::type_y;
                                 if (subface / 2 == 0)
-                                  mask |= ConstraintTypes::type_z;
+                                  mask |= ConstraintKinds::type_z;
                               }
                             // XZ-plane
                             else if (face < 4)
                               {
-                                mask |= ConstraintTypes::face_y;
+                                mask |= ConstraintKinds::face_y;
                                 if (face == 2)
-                                  mask |= ConstraintTypes::type_y;
+                                  mask |= ConstraintKinds::type_y;
                                 if (subface % 2 == 0)
-                                  mask |= ConstraintTypes::type_z;
+                                  mask |= ConstraintKinds::type_z;
                                 if (subface / 2 == 0)
-                                  mask |= ConstraintTypes::type_x;
+                                  mask |= ConstraintKinds::type_x;
                               }
                             // XY-plane
                             else
                               {
-                                mask |= ConstraintTypes::face_z;
+                                mask |= ConstraintKinds::face_z;
                                 if (face == 4)
-                                  mask |= ConstraintTypes::type_z;
+                                  mask |= ConstraintKinds::type_z;
                                 if (subface % 2 == 0)
-                                  mask |= ConstraintTypes::type_x;
+                                  mask |= ConstraintKinds::type_x;
                                 if (subface / 2 == 0)
-                                  mask |= ConstraintTypes::type_y;
+                                  mask |= ConstraintKinds::type_y;
                               }
 
                             // Offset if upper/right/back face
@@ -545,55 +545,55 @@ namespace internal
                 // For each line on cell, which faces does it belong to, what is
                 // the edge mask, what is the types of the faces it belong to,
                 // and what is the type along the edge.
-                const ConstraintTypes line_to_edge[12][4] = {
-                  {ConstraintTypes::face_x | ConstraintTypes::face_z,
-                   ConstraintTypes::edge_zx,
-                   ConstraintTypes::type_x | ConstraintTypes::type_z,
-                   ConstraintTypes::type_y},
-                  {ConstraintTypes::face_x | ConstraintTypes::face_z,
-                   ConstraintTypes::edge_zx,
-                   ConstraintTypes::type_z,
-                   ConstraintTypes::type_y},
-                  {ConstraintTypes::face_y | ConstraintTypes::face_z,
-                   ConstraintTypes::edge_yz,
-                   ConstraintTypes::type_y | ConstraintTypes::type_z,
-                   ConstraintTypes::type_x},
-                  {ConstraintTypes::face_y | ConstraintTypes::face_z,
-                   ConstraintTypes::edge_yz,
-                   ConstraintTypes::type_z,
-                   ConstraintTypes::type_x},
-                  {ConstraintTypes::face_x | ConstraintTypes::face_z,
-                   ConstraintTypes::edge_zx,
-                   ConstraintTypes::type_x,
-                   ConstraintTypes::type_y},
-                  {ConstraintTypes::face_x | ConstraintTypes::face_z,
-                   ConstraintTypes::edge_zx,
-                   ConstraintTypes::unconstrained,
-                   ConstraintTypes::type_y},
-                  {ConstraintTypes::face_y | ConstraintTypes::face_z,
-                   ConstraintTypes::edge_yz,
-                   ConstraintTypes::type_y,
-                   ConstraintTypes::type_x},
-                  {ConstraintTypes::face_y | ConstraintTypes::face_z,
-                   ConstraintTypes::edge_yz,
-                   ConstraintTypes::unconstrained,
-                   ConstraintTypes::type_x},
-                  {ConstraintTypes::face_x | ConstraintTypes::face_y,
-                   ConstraintTypes::edge_xy,
-                   ConstraintTypes::type_x | ConstraintTypes::type_y,
-                   ConstraintTypes::type_z},
-                  {ConstraintTypes::face_x | ConstraintTypes::face_y,
-                   ConstraintTypes::edge_xy,
-                   ConstraintTypes::type_y,
-                   ConstraintTypes::type_z},
-                  {ConstraintTypes::face_x | ConstraintTypes::face_y,
-                   ConstraintTypes::edge_xy,
-                   ConstraintTypes::type_x,
-                   ConstraintTypes::type_z},
-                  {ConstraintTypes::face_x | ConstraintTypes::face_y,
-                   ConstraintTypes::edge_xy,
-                   ConstraintTypes::unconstrained,
-                   ConstraintTypes::type_z}};
+                const ConstraintKinds line_to_edge[12][4] = {
+                  {ConstraintKinds::face_x | ConstraintKinds::face_z,
+                   ConstraintKinds::edge_zx,
+                   ConstraintKinds::type_x | ConstraintKinds::type_z,
+                   ConstraintKinds::type_y},
+                  {ConstraintKinds::face_x | ConstraintKinds::face_z,
+                   ConstraintKinds::edge_zx,
+                   ConstraintKinds::type_z,
+                   ConstraintKinds::type_y},
+                  {ConstraintKinds::face_y | ConstraintKinds::face_z,
+                   ConstraintKinds::edge_yz,
+                   ConstraintKinds::type_y | ConstraintKinds::type_z,
+                   ConstraintKinds::type_x},
+                  {ConstraintKinds::face_y | ConstraintKinds::face_z,
+                   ConstraintKinds::edge_yz,
+                   ConstraintKinds::type_z,
+                   ConstraintKinds::type_x},
+                  {ConstraintKinds::face_x | ConstraintKinds::face_z,
+                   ConstraintKinds::edge_zx,
+                   ConstraintKinds::type_x,
+                   ConstraintKinds::type_y},
+                  {ConstraintKinds::face_x | ConstraintKinds::face_z,
+                   ConstraintKinds::edge_zx,
+                   ConstraintKinds::unconstrained,
+                   ConstraintKinds::type_y},
+                  {ConstraintKinds::face_y | ConstraintKinds::face_z,
+                   ConstraintKinds::edge_yz,
+                   ConstraintKinds::type_y,
+                   ConstraintKinds::type_x},
+                  {ConstraintKinds::face_y | ConstraintKinds::face_z,
+                   ConstraintKinds::edge_yz,
+                   ConstraintKinds::unconstrained,
+                   ConstraintKinds::type_x},
+                  {ConstraintKinds::face_x | ConstraintKinds::face_y,
+                   ConstraintKinds::edge_xy,
+                   ConstraintKinds::type_x | ConstraintKinds::type_y,
+                   ConstraintKinds::type_z},
+                  {ConstraintKinds::face_x | ConstraintKinds::face_y,
+                   ConstraintKinds::edge_xy,
+                   ConstraintKinds::type_y,
+                   ConstraintKinds::type_z},
+                  {ConstraintKinds::face_x | ConstraintKinds::face_y,
+                   ConstraintKinds::edge_xy,
+                   ConstraintKinds::type_x,
+                   ConstraintKinds::type_z},
+                  {ConstraintKinds::face_x | ConstraintKinds::face_y,
+                   ConstraintKinds::edge_xy,
+                   ConstraintKinds::unconstrained,
+                   ConstraintKinds::type_z}};
 
                 for (unsigned int local_line = 0;
                      local_line < GeometryInfo<dim>::lines_per_cell;
@@ -602,7 +602,7 @@ namespace internal
                     // If we don't already have a constraint for as part of a
                     // face
                     if ((mask & line_to_edge[local_line][0]) ==
-                        ConstraintTypes::unconstrained)
+                        ConstraintKinds::unconstrained)
                       {
                         // For each cell which share that edge
                         const unsigned int line =
@@ -698,7 +698,7 @@ namespace internal
                   }
               }
             cell_has_hanging_node_constraints |=
-              mask != ConstraintTypes::unconstrained;
+              mask != ConstraintKinds::unconstrained;
           }
       return cell_has_hanging_node_constraints;
     }
@@ -887,7 +887,7 @@ namespace CUDAWrappers
               typename Number>
     __device__ inline void
     interpolate_boundary_2d(
-      const dealii::internal::MatrixFreeFunctions::ConstraintTypes
+      const dealii::internal::MatrixFreeFunctions::ConstraintKinds
               constraint_mask,
       Number *values)
     {
@@ -896,8 +896,8 @@ namespace CUDAWrappers
 
       const unsigned int this_type =
         (direction == 0) ?
-          dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_x :
-          dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_y;
+          dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_x :
+          dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_y;
 
       const unsigned int interp_idx = (direction == 0) ? x_idx : y_idx;
 
@@ -907,12 +907,12 @@ namespace CUDAWrappers
       const bool constrained_face =
         (constraint_mask &
          (((direction == 0) ?
-             dealii::internal::MatrixFreeFunctions::ConstraintTypes::face_y :
-             dealii::internal::MatrixFreeFunctions::ConstraintTypes::
+             dealii::internal::MatrixFreeFunctions::ConstraintKinds::face_y :
+             dealii::internal::MatrixFreeFunctions::ConstraintKinds::
                unconstrained) |
           ((direction == 1) ?
-             dealii::internal::MatrixFreeFunctions::ConstraintTypes::face_x :
-             dealii::internal::MatrixFreeFunctions::ConstraintTypes::
+             dealii::internal::MatrixFreeFunctions::ConstraintKinds::face_x :
+             dealii::internal::MatrixFreeFunctions::ConstraintKinds::
                unconstrained)));
 
       // Flag is true if for the given direction, the dof is constrained with
@@ -921,12 +921,12 @@ namespace CUDAWrappers
       const bool constrained_dof =
         ((direction == 0) &&
          ((constraint_mask &
-           dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_y) ?
+           dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_y) ?
             (y_idx == 0) :
             (y_idx == fe_degree))) ||
         ((direction == 1) &&
          ((constraint_mask &
-           dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_x) ?
+           dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_x) ?
             (x_idx == 0) :
             (x_idx == fe_degree)));
 
@@ -986,7 +986,7 @@ namespace CUDAWrappers
               typename Number>
     __device__ inline void
     interpolate_boundary_3d(
-      const dealii::internal::MatrixFreeFunctions::ConstraintTypes
+      const dealii::internal::MatrixFreeFunctions::ConstraintKinds
               constraint_mask,
       Number *values)
     {
@@ -996,43 +996,43 @@ namespace CUDAWrappers
 
       const unsigned int this_type =
         (direction == 0) ?
-          dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_x :
+          dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_x :
         (direction == 1) ?
-          dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_y :
-          dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_z;
+          dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_y :
+          dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_z;
       const unsigned int face1_type =
         (direction == 0) ?
-          dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_y :
+          dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_y :
         (direction == 1) ?
-          dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_z :
-          dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_x;
+          dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_z :
+          dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_x;
       const unsigned int face2_type =
         (direction == 0) ?
-          dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_z :
+          dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_z :
         (direction == 1) ?
-          dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_x :
-          dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_y;
+          dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_x :
+          dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_y;
 
       // If computing in x-direction, need to match against face_y or
       // face_z
       const unsigned int face1 =
         (direction == 0) ?
-          dealii::internal::MatrixFreeFunctions::ConstraintTypes::face_y :
+          dealii::internal::MatrixFreeFunctions::ConstraintKinds::face_y :
         (direction == 1) ?
-          dealii::internal::MatrixFreeFunctions::ConstraintTypes::face_z :
-          dealii::internal::MatrixFreeFunctions::ConstraintTypes::face_x;
+          dealii::internal::MatrixFreeFunctions::ConstraintKinds::face_z :
+          dealii::internal::MatrixFreeFunctions::ConstraintKinds::face_x;
       const unsigned int face2 =
         (direction == 0) ?
-          dealii::internal::MatrixFreeFunctions::ConstraintTypes::face_z :
+          dealii::internal::MatrixFreeFunctions::ConstraintKinds::face_z :
         (direction == 1) ?
-          dealii::internal::MatrixFreeFunctions::ConstraintTypes::face_x :
-          dealii::internal::MatrixFreeFunctions::ConstraintTypes::face_y;
+          dealii::internal::MatrixFreeFunctions::ConstraintKinds::face_x :
+          dealii::internal::MatrixFreeFunctions::ConstraintKinds::face_y;
       const unsigned int edge =
         (direction == 0) ?
-          dealii::internal::MatrixFreeFunctions::ConstraintTypes::edge_yz :
+          dealii::internal::MatrixFreeFunctions::ConstraintKinds::edge_yz :
         (direction == 1) ?
-          dealii::internal::MatrixFreeFunctions::ConstraintTypes::edge_zx :
-          dealii::internal::MatrixFreeFunctions::ConstraintTypes::edge_xy;
+          dealii::internal::MatrixFreeFunctions::ConstraintKinds::edge_zx :
+          dealii::internal::MatrixFreeFunctions::ConstraintKinds::edge_xy;
       const unsigned int constrained_face =
         constraint_mask & (face1 | face2 | edge);
 
@@ -1121,7 +1121,7 @@ namespace CUDAWrappers
     template <int dim, int fe_degree, bool transpose, typename Number>
     __device__ void
     resolve_hanging_nodes(
-      const dealii::internal::MatrixFreeFunctions::ConstraintTypes
+      const dealii::internal::MatrixFreeFunctions::ConstraintKinds
               constraint_mask,
       Number *values)
     {
index d8160a3e816a1b19f74799428c9be24f6078373d..9b0438583ecd6e9d4d56195d9d6ed87a7d75dee1 100644 (file)
@@ -1306,7 +1306,7 @@ namespace internal
         {
           auto &vec = dof_info[no].hanging_node_constraint_masks;
           if (std::all_of(vec.begin(), vec.end(), [](const auto i) {
-                return i == internal::MatrixFreeFunctions::ConstraintTypes::
+                return i == internal::MatrixFreeFunctions::ConstraintKinds::
                               unconstrained;
               }))
             vec.clear();
index 0d8d701425f9cd186af247b5ed87e3ad60da15ed..9e7461d4e31f84c24d29887c1ff9e31a0d3d70c0 100644 (file)
@@ -454,7 +454,7 @@ namespace MatrixFreeTools
 
                 // cell has hanging nodes
                 if (mask != dealii::internal::MatrixFreeFunctions::
-                              ConstraintTypes::unconstrained)
+                              ConstraintKinds::unconstrained)
                   {
                     // check if hanging node internpolation matrix has been set
                     // up
@@ -467,7 +467,7 @@ namespace MatrixFreeTools
                           dofs_per_component);
 
                         std::array<dealii::internal::MatrixFreeFunctions::
-                                     ConstraintTypes,
+                                     ConstraintKinds,
                                    VectorizedArrayType::size()>
                           constraint_mask;
                         constraint_mask[0] = mask;
@@ -768,7 +768,7 @@ namespace MatrixFreeTools
       // constraints!
       std::array<std::vector<Number>, n_lanes> diagonals_local_constrained;
 
-      std::map<dealii::internal::MatrixFreeFunctions::ConstraintTypes,
+      std::map<dealii::internal::MatrixFreeFunctions::ConstraintKinds,
                std::vector<std::tuple<unsigned int, unsigned int, Number>>>
         locally_relevant_constrains_hn_map;
     };
index 31d17618c9f1904b0d74de0c6f5777a359af330b..f445e34f1ecb2979a49c775cdb9f934162e74260 100644 (file)
@@ -105,7 +105,7 @@ namespace internal
           const bool has_constraints =
             (hanging_node_constraint_masks.size() != 0 &&
              hanging_node_constraint_masks[ib] !=
-               ConstraintTypes::unconstrained) ||
+               ConstraintKinds::unconstrained) ||
             (row_starts[ib].second != row_starts[ib + n_fe_components].second);
 
           auto do_copy = [&](const unsigned int *begin,
@@ -232,7 +232,7 @@ namespace internal
                         hanging_node_constraint_masks[boundary_cells[i] *
                                                         n_components +
                                                       comp] !=
-                        ConstraintTypes::unconstrained;
+                        ConstraintKinds::unconstrained;
 
                   if (has_hanging_nodes ||
                       row_starts[boundary_cells[i] * n_components].second !=
@@ -338,7 +338,7 @@ namespace internal
       new_dof_indices.reserve(dof_indices.size());
       new_constraint_indicator.reserve(constraint_indicator.size());
 
-      std::vector<ConstraintTypes> new_hanging_node_constraint_masks;
+      std::vector<ConstraintKinds> new_hanging_node_constraint_masks;
       new_hanging_node_constraint_masks.reserve(
         new_hanging_node_constraint_masks.size());
 
@@ -388,7 +388,7 @@ namespace internal
                         hanging_node_constraint_masks[cell_no + comp];
                       new_hanging_node_constraint_masks.push_back(mask);
                       has_hanging_nodes |=
-                        mask != ConstraintTypes::unconstrained;
+                        mask != ConstraintKinds::unconstrained;
                     }
 
                   new_dof_indices.insert(
@@ -429,7 +429,7 @@ namespace internal
 
                 if (hanging_node_constraint_masks.size() > 0)
                   new_hanging_node_constraint_masks.push_back(
-                    ConstraintTypes::unconstrained);
+                    ConstraintKinds::unconstrained);
               }
           position_cell += n_vect;
         }
index 1a5a7e3c14989239157af2b58aa57b9a37f5d4bc..ff0999b81b6bfe68ab60844a42b4f4374b06e4db 100644 (file)
@@ -51,7 +51,7 @@ namespace dealii
                                    is_face,
                                    Number> &fe_eval,
         const bool                          transpose,
-        const std::array<dealii::internal::MatrixFreeFunctions::ConstraintTypes,
+        const std::array<dealii::internal::MatrixFreeFunctions::ConstraintKinds,
                          Number::size()> &  c_mask,
         Number *                            values)
       {
@@ -109,7 +109,7 @@ namespace dealii
                                    typename Number::value_type,
                                    is_face,
                                    Number> &fe_eval,
-        const std::array<dealii::internal::MatrixFreeFunctions::ConstraintTypes,
+        const std::array<dealii::internal::MatrixFreeFunctions::ConstraintKinds,
                          Number::size()> &  constraint_mask,
         Number *                            values)
       {
@@ -131,30 +131,30 @@ namespace dealii
         for (unsigned int v = 0; v < Number::size(); ++v)
           {
             if (constraint_mask[v] == dealii::internal::MatrixFreeFunctions::
-                                        ConstraintTypes::unconstrained)
+                                        ConstraintKinds::unconstrained)
               continue;
 
             const auto this_type =
               (direction == 0) ?
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_x :
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_y;
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_x :
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_y;
 
             const bool constrained_face =
               (constraint_mask[v] &
                (((direction == 0) ? dealii::internal::MatrixFreeFunctions::
-                                      ConstraintTypes::face_y :
+                                      ConstraintKinds::face_y :
                                     dealii::internal::MatrixFreeFunctions::
-                                      ConstraintTypes::unconstrained) |
+                                      ConstraintKinds::unconstrained) |
                 ((direction == 1) ? dealii::internal::MatrixFreeFunctions::
-                                      ConstraintTypes::face_x :
+                                      ConstraintKinds::face_x :
                                     dealii::internal::MatrixFreeFunctions::
-                                      ConstraintTypes::unconstrained))) !=
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::
+                                      ConstraintKinds::unconstrained))) !=
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::
                 unconstrained;
 
             const bool type = (constraint_mask[v] & this_type) !=
                               dealii::internal::MatrixFreeFunctions::
-                                ConstraintTypes::unconstrained;
+                                ConstraintKinds::unconstrained;
 
             for (unsigned int x_idx = 0; x_idx < fe_degree + 1; ++x_idx)
               for (unsigned int y_idx = 0; y_idx < fe_degree + 1; ++y_idx)
@@ -172,16 +172,16 @@ namespace dealii
                   const bool constrained_dof =
                     ((direction == 0) &&
                      (((constraint_mask[v] &
-                        dealii::internal::MatrixFreeFunctions::ConstraintTypes::
+                        dealii::internal::MatrixFreeFunctions::ConstraintKinds::
                           type_y) != dealii::internal::MatrixFreeFunctions::
-                                       ConstraintTypes::unconstrained) ?
+                                       ConstraintKinds::unconstrained) ?
                         (y_idx == 0) :
                         (y_idx == fe_degree))) ||
                     ((direction == 1) &&
                      (((constraint_mask[v] &
-                        dealii::internal::MatrixFreeFunctions::ConstraintTypes::
+                        dealii::internal::MatrixFreeFunctions::ConstraintKinds::
                           type_x) != dealii::internal::MatrixFreeFunctions::
-                                       ConstraintTypes::unconstrained) ?
+                                       ConstraintKinds::unconstrained) ?
                         (x_idx == 0) :
                         (x_idx == fe_degree)));
 
@@ -248,7 +248,7 @@ namespace dealii
                                    typename Number::value_type,
                                    is_face,
                                    Number> &fe_eval,
-        const std::array<dealii::internal::MatrixFreeFunctions::ConstraintTypes,
+        const std::array<dealii::internal::MatrixFreeFunctions::ConstraintKinds,
                          Number::size()> &  constraint_mask,
         Number *                            values)
       {
@@ -269,49 +269,49 @@ namespace dealii
 
         const auto this_type =
           (direction == 0) ?
-            dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_x :
+            dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_x :
           (direction == 1) ?
-            dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_y :
-            dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_z;
+            dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_y :
+            dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_z;
         const auto face1_type =
           (direction == 0) ?
-            dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_y :
+            dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_y :
           (direction == 1) ?
-            dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_z :
-            dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_x;
+            dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_z :
+            dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_x;
         const auto face2_type =
           (direction == 0) ?
-            dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_z :
+            dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_z :
           (direction == 1) ?
-            dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_x :
-            dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_y;
+            dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_x :
+            dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_y;
 
         // If computing in x-direction, need to match against
-        // dealii::internal::MatrixFreeFunctions::ConstraintTypes::face_y or
-        // dealii::internal::MatrixFreeFunctions::ConstraintTypes::face_z
+        // dealii::internal::MatrixFreeFunctions::ConstraintKinds::face_y or
+        // dealii::internal::MatrixFreeFunctions::ConstraintKinds::face_z
         const auto face1 =
           (direction == 0) ?
-            dealii::internal::MatrixFreeFunctions::ConstraintTypes::face_y :
+            dealii::internal::MatrixFreeFunctions::ConstraintKinds::face_y :
           (direction == 1) ?
-            dealii::internal::MatrixFreeFunctions::ConstraintTypes::face_z :
-            dealii::internal::MatrixFreeFunctions::ConstraintTypes::face_x;
+            dealii::internal::MatrixFreeFunctions::ConstraintKinds::face_z :
+            dealii::internal::MatrixFreeFunctions::ConstraintKinds::face_x;
         const auto face2 =
           (direction == 0) ?
-            dealii::internal::MatrixFreeFunctions::ConstraintTypes::face_z :
+            dealii::internal::MatrixFreeFunctions::ConstraintKinds::face_z :
           (direction == 1) ?
-            dealii::internal::MatrixFreeFunctions::ConstraintTypes::face_x :
-            dealii::internal::MatrixFreeFunctions::ConstraintTypes::face_y;
+            dealii::internal::MatrixFreeFunctions::ConstraintKinds::face_x :
+            dealii::internal::MatrixFreeFunctions::ConstraintKinds::face_y;
         const auto edge =
           (direction == 0) ?
-            dealii::internal::MatrixFreeFunctions::ConstraintTypes::edge_yz :
+            dealii::internal::MatrixFreeFunctions::ConstraintKinds::edge_yz :
           (direction == 1) ?
-            dealii::internal::MatrixFreeFunctions::ConstraintTypes::edge_zx :
-            dealii::internal::MatrixFreeFunctions::ConstraintTypes::edge_xy;
+            dealii::internal::MatrixFreeFunctions::ConstraintKinds::edge_zx :
+            dealii::internal::MatrixFreeFunctions::ConstraintKinds::edge_xy;
 
         for (unsigned int v = 0; v < Number::size(); ++v)
           {
             if (constraint_mask[v] == dealii::internal::MatrixFreeFunctions::
-                                        ConstraintTypes::unconstrained)
+                                        ConstraintKinds::unconstrained)
               continue;
 
             const auto constrained_face =
@@ -319,7 +319,7 @@ namespace dealii
 
             const bool type = (constraint_mask[v] & this_type) !=
                               dealii::internal::MatrixFreeFunctions::
-                                ConstraintTypes::unconstrained;
+                                ConstraintKinds::unconstrained;
 
             for (unsigned int x_idx = 0; x_idx < fe_degree + 1; ++x_idx)
               for (unsigned int y_idx = 0; y_idx < fe_degree + 1; ++y_idx)
@@ -338,33 +338,33 @@ namespace dealii
                     typename Number::value_type t = 0;
                     const bool                  on_face1 =
                       ((constraint_mask[v] & face1_type) !=
-                       dealii::internal::MatrixFreeFunctions::ConstraintTypes::
+                       dealii::internal::MatrixFreeFunctions::ConstraintKinds::
                          unconstrained) ?
                                          (face1_idx == 0) :
                                          (face1_idx == fe_degree);
                     const bool on_face2 =
                       ((constraint_mask[v] & face2_type) !=
-                       dealii::internal::MatrixFreeFunctions::ConstraintTypes::
+                       dealii::internal::MatrixFreeFunctions::ConstraintKinds::
                          unconstrained) ?
                         (face2_idx == 0) :
                         (face2_idx == fe_degree);
                     const bool constrained_dof =
                       ((((constraint_mask[v] & face1) !=
                          dealii::internal::MatrixFreeFunctions::
-                           ConstraintTypes::unconstrained) &&
+                           ConstraintKinds::unconstrained) &&
                         on_face1) ||
                        (((constraint_mask[v] & face2) !=
                          dealii::internal::MatrixFreeFunctions::
-                           ConstraintTypes::unconstrained) &&
+                           ConstraintKinds::unconstrained) &&
                         on_face2) ||
                        (((constraint_mask[v] & edge) !=
                          dealii::internal::MatrixFreeFunctions::
-                           ConstraintTypes::unconstrained) &&
+                           ConstraintKinds::unconstrained) &&
                         on_face1 && on_face2));
 
                     if ((constrained_face !=
                          dealii::internal::MatrixFreeFunctions::
-                           ConstraintTypes::unconstrained) &&
+                           ConstraintKinds::unconstrained) &&
                         constrained_dof)
                       {
                         if (type)
@@ -426,7 +426,7 @@ namespace dealii
 template <int dim>
 void
 test(const unsigned int                                           degree,
-     const dealii::internal::MatrixFreeFunctions::ConstraintTypes mask_value)
+     const dealii::internal::MatrixFreeFunctions::ConstraintKinds mask_value)
 {
   Triangulation<dim> tria;
   GridGenerator::subdivided_hyper_cube(tria, 2);
@@ -455,13 +455,13 @@ test(const unsigned int                                           degree,
   FEEvaluation<dim, -1, 0, 1, double> eval(matrix_free);
   eval.reinit(0);
 
-  std::array<dealii::internal::MatrixFreeFunctions::ConstraintTypes,
+  std::array<dealii::internal::MatrixFreeFunctions::ConstraintKinds,
              VectorizedArray<double>::size()>
     cmask;
   std::fill(
     cmask.begin(),
     cmask.end(),
-    dealii::internal::MatrixFreeFunctions::ConstraintTypes::unconstrained);
+    dealii::internal::MatrixFreeFunctions::ConstraintKinds::unconstrained);
   cmask[0] = mask_value;
 
   for (unsigned int b = 0; b < 2; ++b)
@@ -515,44 +515,44 @@ main(int argc, char **argv)
     {
       test<2>(
         degree,
-        dealii::internal::MatrixFreeFunctions::ConstraintTypes::unconstrained);
+        dealii::internal::MatrixFreeFunctions::ConstraintKinds::unconstrained);
       deallog << std::endl;
 
       test<2>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::face_x |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_x |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::face_x |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_x |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::
                   type_y); // face 0/0
       test<2>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::face_x |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::face_x |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::
                   type_x); // face 0/1
       test<2>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::face_x |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::face_x |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::
                   type_y); // face 1/0
       test<2>(
         degree,
-        dealii::internal::MatrixFreeFunctions::ConstraintTypes::face_x); // face
+        dealii::internal::MatrixFreeFunctions::ConstraintKinds::face_x); // face
                                                                          // 1/1
       deallog << std::endl;
 
       test<2>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::face_y |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_y |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::face_y |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_y |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::
                   type_x); // face 2/0
       test<2>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::face_y |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::face_y |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::
                   type_y); // face 2/1
       test<2>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::face_y |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::face_y |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::
                   type_x); // face 3/0
       test<2>(
         degree,
-        dealii::internal::MatrixFreeFunctions::ConstraintTypes::face_y); // face
+        dealii::internal::MatrixFreeFunctions::ConstraintKinds::face_y); // face
                                                                          // 3/1
       deallog << std::endl;
     }
@@ -561,140 +561,140 @@ main(int argc, char **argv)
     {
       // edge 2
       test<3>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::edge_yz |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_y |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_z);
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::edge_yz |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_y |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_z);
       test<3>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::edge_yz |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_y |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_z |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_x);
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::edge_yz |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_y |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_z |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_x);
 
       // edge 3
       test<3>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::edge_yz |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_z);
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::edge_yz |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_z);
       test<3>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::edge_yz |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_z |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_x);
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::edge_yz |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_z |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_x);
 
       // edge 6
       test<3>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::edge_yz |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_y);
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::edge_yz |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_y);
       test<3>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::edge_yz |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_y |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_x);
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::edge_yz |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_y |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_x);
 
       // edge 7
       test<3>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::edge_yz);
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::edge_yz);
       test<3>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::edge_yz |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_x);
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::edge_yz |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_x);
 
 
       // edge 0
       test<3>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::edge_zx |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_x |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_z);
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::edge_zx |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_x |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_z);
       test<3>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::edge_zx |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_x |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_z |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_y);
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::edge_zx |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_x |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_z |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_y);
 
       // edge 1
       test<3>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::edge_zx |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_z);
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::edge_zx |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_z);
       test<3>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::edge_zx |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_z |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_y);
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::edge_zx |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_z |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_y);
 
       // edge 4
       test<3>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::edge_zx |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_x);
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::edge_zx |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_x);
       test<3>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::edge_zx |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_x |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_y);
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::edge_zx |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_x |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_y);
 
       // edge 5
       test<3>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::edge_zx);
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::edge_zx);
       test<3>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::edge_zx |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_y);
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::edge_zx |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_y);
 
 
       // edge 8
       test<3>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::edge_xy |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_x |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_y);
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::edge_xy |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_x |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_y);
       test<3>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::edge_xy |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_x |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_y |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_z);
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::edge_xy |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_x |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_y |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_z);
 
       // edge 9
       test<3>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::edge_xy |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_y);
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::edge_xy |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_y);
       test<3>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::edge_xy |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_y |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_z);
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::edge_xy |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_y |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_z);
 
       // edge 10
       test<3>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::edge_xy |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_x);
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::edge_xy |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_x);
       test<3>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::edge_xy |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_x |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_z);
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::edge_xy |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_x |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_z);
 
       // edge 11
       test<3>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::edge_xy);
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::edge_xy);
       test<3>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::edge_xy |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_z);
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::edge_xy |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_z);
 
 
       // face 0
       test<3>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::face_x |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_x);
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::face_x |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_x);
 
       // face 1
       test<3>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::face_x);
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::face_x);
 
       // face 2
       test<3>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::face_y |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_y);
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::face_y |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_y);
 
       // face 3
       test<3>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::face_y);
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::face_y);
 
       // face 4
       test<3>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::face_z |
-                dealii::internal::MatrixFreeFunctions::ConstraintTypes::type_z);
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::face_z |
+                dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_z);
 
       // face 5
       test<3>(degree,
-              dealii::internal::MatrixFreeFunctions::ConstraintTypes::face_z);
+              dealii::internal::MatrixFreeFunctions::ConstraintKinds::face_z);
     }
 }

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.