]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Replace level_mg_handler by mg_level 8919/head
authorPeter Munch <peterrmuench@gmail.com>
Thu, 17 Oct 2019 10:47:24 +0000 (12:47 +0200)
committerPeter Munch <peterrmuench@gmail.com>
Thu, 17 Oct 2019 15:53:40 +0000 (17:53 +0200)
32 files changed:
doc/news/changes/incompatibilities/20191017PeterMunch [new file with mode: 0644]
examples/step-37/step-37.cc
examples/step-59/step-59.cc
include/deal.II/matrix_free/face_setup_internal.h
include/deal.II/matrix_free/matrix_free.h
include/deal.II/matrix_free/matrix_free.templates.h
tests/mappings/mapping_q_eulerian_08.cc
tests/matrix_free/cell_categorization_02.cc
tests/matrix_free/cell_level_and_index.cc
tests/matrix_free/compare_faces_by_cells.cc
tests/matrix_free/dg_pbc_02.cc
tests/matrix_free/matrix_vector_mg.cc
tests/matrix_free/multigrid_dg_periodic.cc
tests/matrix_free/multigrid_dg_sip_01.cc
tests/matrix_free/multigrid_dg_sip_02.cc
tests/matrix_free/parallel_multigrid.cc
tests/matrix_free/parallel_multigrid_02.cc
tests/matrix_free/parallel_multigrid_adaptive_01.cc
tests/matrix_free/parallel_multigrid_adaptive_02.cc
tests/matrix_free/parallel_multigrid_adaptive_03.cc
tests/matrix_free/parallel_multigrid_adaptive_04.cc
tests/matrix_free/parallel_multigrid_adaptive_05.cc
tests/matrix_free/parallel_multigrid_adaptive_06.cc
tests/matrix_free/parallel_multigrid_adaptive_06ref.cc
tests/matrix_free/parallel_multigrid_adaptive_07.cc
tests/matrix_free/parallel_multigrid_adaptive_08.cc
tests/matrix_free/parallel_multigrid_fullydistributed.cc
tests/matrix_free/parallel_multigrid_mf.cc
tests/matrix_free/step-37-inhomogeneous-1.cc
tests/matrix_free/step-37.cc
tests/matrix_free/stokes_computation.cc
tests/mpi/step-37.cc

diff --git a/doc/news/changes/incompatibilities/20191017PeterMunch b/doc/news/changes/incompatibilities/20191017PeterMunch
new file mode 100644 (file)
index 0000000..0fe663a
--- /dev/null
@@ -0,0 +1,4 @@
+Changed: The variable level_mg_handler has been renamed
+to mg_level in MatrixFree and its associated AdditionalData.
+<br>
+(Peter Munch, 2019/10/17)
index 9495fcd8fe0b3f792af547ba4ccb05ee64862099..b08866d4a7937f95d0ae2b403179a6823f1ae889 100644 (file)
@@ -866,7 +866,7 @@ namespace Step37
           MatrixFree<dim, float>::AdditionalData::none;
         additional_data.mapping_update_flags =
           (update_gradients | update_JxW_values | update_quadrature_points);
-        additional_data.level_mg_handler = level;
+        additional_data.mg_level = level;
         std::shared_ptr<MatrixFree<dim, float>> mg_mf_storage_level(
           new MatrixFree<dim, float>());
         mg_mf_storage_level->reinit(dof_handler,
index 59ec4472c3cce61c5329f213e21ef760430a723b..9378c1f8cdaced12d8904bbcebb5e82256167fa2 100644 (file)
@@ -1041,7 +1041,7 @@ namespace Step59
           (update_gradients | update_JxW_values);
         additional_data.mapping_update_flags_boundary_faces =
           (update_gradients | update_JxW_values);
-        additional_data.level_mg_handler = level;
+        additional_data.mg_level = level;
         const auto mg_mf_storage_level =
           std::make_shared<MatrixFree<dim, float>>();
         mg_mf_storage_level->reinit(dof_handler,
index b438b23a00d0258153c0d4447dd63332741365a7..76ebe7f3dd23715ef6ed7bf8692c202f23a9e1d0 100644 (file)
@@ -171,7 +171,7 @@ namespace internal
       std::vector<std::pair<unsigned int, unsigned int>> &cell_levels)
     {
       use_active_cells =
-        additional_data.level_mg_handler == numbers::invalid_unsigned_int;
+        additional_data.mg_level == numbers::invalid_unsigned_int;
 
 #  ifdef DEBUG
       // safety check
@@ -582,7 +582,7 @@ namespace internal
               // inside the domain in case of multigrid separately
               else if ((dcell->at_boundary(f) == false ||
                         dcell->has_periodic_neighbor(f)) &&
-                       additional_data.level_mg_handler !=
+                       additional_data.mg_level !=
                          numbers::invalid_unsigned_int &&
                        dcell->neighbor_or_periodic_neighbor(f)->level() <
                          dcell->level())
index 4c7ed65a9271c6a929070234eee238221c0b53a3..d2977679f997893c36ae9623be737fbdc9606cf0 100644 (file)
@@ -166,10 +166,10 @@ public:
    * local time stepping where one would like to control which elements
    * together form a batch of cells. The array `cell_vectorization_categories`
    * is accessed by the number given by cell->active_cell_index() when working
-   * on the active cells with `level_mg_handler` set to `-1` and by
-   * cell->index() for the level cells. By default, the different categories
-   * in `cell_vectorization_category` can be mixed and the algorithm is
-   * allowed to merge lower category numbers with the next higher categories
+   * on the active cells with `mg_level` set to `numbers::invalid_unsigned_int`
+   * and by cell->index() for the level cells. By default, the different
+   * categories in `cell_vectorization_category` can be mixed and the algorithm
+   * is allowed to merge lower category numbers with the next higher categories
    * if it is necessary inside the algorithm, in order to avoid partially
    * filled SIMD lanes as much as possible. This gives a better utilization of
    * the vectorization but might need special treatment, in particular for
@@ -218,7 +218,7 @@ public:
       const UpdateFlags  mapping_update_flags_boundary_faces = update_default,
       const UpdateFlags  mapping_update_flags_inner_faces    = update_default,
       const UpdateFlags  mapping_update_flags_faces_by_cells = update_default,
-      const unsigned int level_mg_handler    = numbers::invalid_unsigned_int,
+      const unsigned int mg_level            = numbers::invalid_unsigned_int,
       const bool         store_plain_indices = true,
       const bool         initialize_indices  = true,
       const bool         initialize_mapping  = true,
@@ -231,7 +231,8 @@ public:
       , mapping_update_flags_boundary_faces(mapping_update_flags_boundary_faces)
       , mapping_update_flags_inner_faces(mapping_update_flags_inner_faces)
       , mapping_update_flags_faces_by_cells(mapping_update_flags_faces_by_cells)
-      , level_mg_handler(level_mg_handler)
+      , mg_level(mg_level)
+      , level_mg_handler(this->mg_level)
       , store_plain_indices(store_plain_indices)
       , initialize_indices(initialize_indices)
       , initialize_mapping(initialize_mapping)
@@ -381,7 +382,14 @@ public:
      * to work on a level, its dofs must be distributed by using
      * <code>dof_handler.distribute_mg_dofs(fe);</code>.
      */
-    unsigned int level_mg_handler;
+    unsigned int mg_level;
+
+    /**
+     * Alias for mg_level
+     *
+     * @deprecated Use mg_level instead.
+     */
+    DEAL_II_DEPRECATED unsigned int &level_mg_handler;
 
     /**
      * Controls whether to allow reading from vectors without resolving
@@ -442,7 +450,7 @@ public:
      * control which elements together form a batch of cells.
      *
      * This array is accessed by the number given by cell->active_cell_index()
-     * when working on the active cells with @p level_mg_handler set to -1 and
+     * when working on the active cells with @p mg_level set to numbers::invalid_unsigned_int and
      * by cell->index() for the level cells.
      *
      * @note This field is empty upon construction of AdditionalData. It is
@@ -1602,10 +1610,11 @@ public:
   mapping_initialized() const;
 
   /**
-   * Return the level of the mesh to be worked on.
+   * Return the level of the mesh to be worked on. Returns
+   * numbers::invalid_unsigned_int if working on active cells.
    */
   unsigned int
-  get_level_mg_handler() const;
+  get_mg_level() const;
 
   /**
    * Return an approximation of the memory consumption of this class in
@@ -1938,7 +1947,7 @@ private:
   /**
    * Stored the level of the mesh to be worked on.
    */
-  unsigned int level_mg_handler;
+  unsigned int mg_level;
 };
 
 
@@ -2494,9 +2503,9 @@ MatrixFree<dim, Number, VectorizedArrayType>::mapping_initialized() const
 
 template <int dim, typename Number, typename VectorizedArrayType>
 inline unsigned int
-MatrixFree<dim, Number, VectorizedArrayType>::get_level_mg_handler() const
+MatrixFree<dim, Number, VectorizedArrayType>::get_mg_level() const
 {
-  return level_mg_handler;
+  return mg_level;
 }
 
 
@@ -2644,7 +2653,7 @@ MatrixFree<dim, Number, VectorizedArrayType>::reinit(
 
   std::vector<IndexSet> locally_owned_sets =
     internal::MatrixFreeImplementation::extract_locally_owned_index_sets(
-      dof_handlers, additional_data.level_mg_handler);
+      dof_handlers, additional_data.mg_level);
 
   std::vector<hp::QCollection<1>> quad_hp;
   quad_hp.emplace_back(quad);
@@ -2678,7 +2687,7 @@ MatrixFree<dim, Number, VectorizedArrayType>::reinit(
 
   std::vector<IndexSet> locally_owned_sets =
     internal::MatrixFreeImplementation::extract_locally_owned_index_sets(
-      dof_handlers, additional_data.level_mg_handler);
+      dof_handlers, additional_data.mg_level);
 
   std::vector<hp::QCollection<1>> quad_hp;
   quad_hp.emplace_back(quad);
@@ -2705,7 +2714,7 @@ MatrixFree<dim, Number, VectorizedArrayType>::reinit(
 {
   std::vector<IndexSet> locally_owned_set =
     internal::MatrixFreeImplementation::extract_locally_owned_index_sets(
-      dof_handler, additional_data.level_mg_handler);
+      dof_handler, additional_data.mg_level);
   std::vector<hp::QCollection<1>> quad_hp;
   for (unsigned int q = 0; q < quad.size(); ++q)
     quad_hp.emplace_back(quad[q]);
@@ -2731,7 +2740,7 @@ MatrixFree<dim, Number, VectorizedArrayType>::reinit(
 {
   std::vector<IndexSet> locally_owned_set =
     internal::MatrixFreeImplementation::extract_locally_owned_index_sets(
-      dof_handler, additional_data.level_mg_handler);
+      dof_handler, additional_data.mg_level);
   std::vector<hp::QCollection<1>> quad_hp;
   quad_hp.emplace_back(quad);
   internal_reinit(StaticMappingQ1<dim>::mapping,
@@ -2757,7 +2766,7 @@ MatrixFree<dim, Number, VectorizedArrayType>::reinit(
 {
   std::vector<IndexSet> locally_owned_set =
     internal::MatrixFreeImplementation::extract_locally_owned_index_sets(
-      dof_handler, additional_data.level_mg_handler);
+      dof_handler, additional_data.mg_level);
   std::vector<hp::QCollection<1>> quad_hp;
   quad_hp.emplace_back(quad);
   internal_reinit(mapping,
@@ -2783,7 +2792,7 @@ MatrixFree<dim, Number, VectorizedArrayType>::reinit(
 {
   std::vector<IndexSet> locally_owned_set =
     internal::MatrixFreeImplementation::extract_locally_owned_index_sets(
-      dof_handler, additional_data.level_mg_handler);
+      dof_handler, additional_data.mg_level);
   std::vector<hp::QCollection<1>> quad_hp;
   for (unsigned int q = 0; q < quad.size(); ++q)
     quad_hp.emplace_back(quad[q]);
index a1c2aa7ee186cd173f0ef49f6cf42bf8d550fc22..0bf653412ef3b70e2b18dc118839c2b33c0729c6 100644 (file)
@@ -58,7 +58,7 @@ MatrixFree<dim, Number, VectorizedArrayType>::MatrixFree()
   : Subscriptor()
   , indices_are_initialized(false)
   , mapping_is_initialized(false)
-  , level_mg_handler(numbers::invalid_unsigned_int)
+  , mg_level(numbers::invalid_unsigned_int)
 {}
 
 
@@ -250,7 +250,7 @@ MatrixFree<dim, Number, VectorizedArrayType>::copy_from(
   task_info                 = v.task_info;
   indices_are_initialized   = v.indices_are_initialized;
   mapping_is_initialized    = v.mapping_is_initialized;
-  level_mg_handler          = v.level_mg_handler;
+  mg_level                  = v.mg_level;
 }
 
 
@@ -268,7 +268,7 @@ MatrixFree<dim, Number, VectorizedArrayType>::internal_reinit(
     &additional_data)
 {
   // Store the level of the mesh to be worked on.
-  this->level_mg_handler = additional_data.level_mg_handler;
+  this->mg_level = additional_data.mg_level;
 
   // Reads out the FE information and stores the shape function values,
   // gradients and Hessians for quadrature points.
@@ -431,7 +431,7 @@ MatrixFree<dim, Number, VectorizedArrayType>::internal_reinit(
     &additional_data)
 {
   // Store the level of the mesh to be worked on.
-  this->level_mg_handler = additional_data.level_mg_handler;
+  this->mg_level = additional_data.mg_level;
 
   // Reads out the FE information and stores the shape function values,
   // gradients and Hessians for quadrature points.
@@ -658,7 +658,7 @@ MatrixFree<dim, Number, VectorizedArrayType>::initialize_dof_handlers(
 {
   cell_level_index.clear();
   dof_handlers.active_dof_handler = DoFHandlers::usual;
-  dof_handlers.level              = additional_data.level_mg_handler;
+  dof_handlers.level              = additional_data.mg_level;
   dof_handlers.n_dof_handlers     = dof_handler.size();
   dof_handlers.dof_handler.resize(dof_handlers.n_dof_handlers);
   for (unsigned int no = 0; no < dof_handlers.n_dof_handlers; ++no)
@@ -676,7 +676,7 @@ MatrixFree<dim, Number, VectorizedArrayType>::initialize_dof_handlers(
 
   const Triangulation<dim> &tria =
     dof_handlers.dof_handler[0]->get_triangulation();
-  const unsigned int level = additional_data.level_mg_handler;
+  const unsigned int level = additional_data.mg_level;
   if (level == numbers::invalid_unsigned_int)
     {
       if (n_mpi_procs == 1)
@@ -727,7 +727,7 @@ MatrixFree<dim, Number, VectorizedArrayType>::initialize_dof_handlers(
 {
   cell_level_index.clear();
   dof_handlers.active_dof_handler = DoFHandlers::hp;
-  dof_handlers.level              = additional_data.level_mg_handler;
+  dof_handlers.level              = additional_data.mg_level;
   Assert(dof_handlers.level == numbers::invalid_unsigned_int,
          ExcNotImplemented());
   dof_handlers.n_dof_handlers = dof_handler.size();
@@ -1039,19 +1039,18 @@ MatrixFree<dim, Number, VectorizedArrayType>::initialize_indices(
     for (unsigned int no = 0; no < n_fe; ++no)
       {
         if (do_face_integrals &&
-            additional_data.level_mg_handler != numbers::invalid_unsigned_int)
+            additional_data.mg_level != numbers::invalid_unsigned_int)
           {
             // in case of adaptivity, go through the cells on the next finer
             // level and check whether we need to get read access to some of
             // those entries for the mg flux matrices
             const DoFHandler<dim> &dof_handler = *dof_handlers.dof_handler[no];
             std::vector<types::global_dof_index> dof_indices;
-            if (additional_data.level_mg_handler + 1 <
+            if (additional_data.mg_level + 1 <
                 dof_handler.get_triangulation().n_global_levels())
               for (typename DoFHandler<dim>::cell_iterator cell =
-                     dof_handler.begin(additional_data.level_mg_handler + 1);
-                   cell !=
-                   dof_handler.end(additional_data.level_mg_handler + 1);
+                     dof_handler.begin(additional_data.mg_level + 1);
+                   cell != dof_handler.end(additional_data.mg_level + 1);
                    ++cell)
                 if (cell->level_subdomain_id() == task_info.my_pid)
                   for (unsigned int f = 0;
index 1ab5851756c3eab888ba75ef795b1b237e018b96..a6afbbda8a8dd892084a224378224ac76b3ffb8e 100644 (file)
@@ -280,7 +280,7 @@ test(const unsigned int n_ref = 0)
         mg_additional_data;
       mg_additional_data.tasks_parallel_scheme =
         MatrixFree<dim, LevelNumberType>::AdditionalData::partition_color;
-      mg_additional_data.level_mg_handler = level;
+      mg_additional_data.mg_level = level;
       mg_additional_data.mapping_update_flags =
         update_values | update_gradients | update_JxW_values |
         update_quadrature_points;
index d415779c66f8ed7a10b5b2123dd9bff6fc170902..4269eeb27bf014551c517cdc4f5d133758b729dc 100644 (file)
@@ -135,7 +135,7 @@ test()
               .cell_vectorization_category[cell->index()] = cell->material_id();
           }
 
-      mg_additional_data[level].level_mg_handler = level;
+      mg_additional_data[level].mg_level = level;
     }
 
   std::vector<std::shared_ptr<MatrixFree<dim, float>>> mg_mf_data(max_level +
index ac16715d5d27fa146ad0e831be70533bba6cba9a..bd572189021b1c54a763f60477caaaaf9b18c0bb 100644 (file)
@@ -115,8 +115,8 @@ test(const bool adaptive_ref = true)
 
   {
     std::cout << "Compare level indices." << std::endl;
-    const unsigned int level         = tria.n_global_levels() - 1;
-    additional_data.level_mg_handler = level;
+    const unsigned int level = tria.n_global_levels() - 1;
+    additional_data.mg_level = level;
     mf_data->reinit(dof, constraints, quad, additional_data);
     compare_indices(mf_data.get());
   }
index e8626fdcae4205d4792fdeff61c34ecd8eec3f3d..835e8fce103d155edb971efb60163934231c1804 100644 (file)
@@ -64,7 +64,7 @@ public:
     addit_data.tasks_parallel_scheme =
       MatrixFree<dim, number>::AdditionalData::none;
     addit_data.tasks_block_size = 3;
-    addit_data.level_mg_handler = level;
+    addit_data.mg_level         = level;
     addit_data.mapping_update_flags_inner_faces =
       update_JxW_values | update_normal_vectors | update_jacobians;
     addit_data.mapping_update_flags_boundary_faces =
index 2a67eae8cdd1a2e7ab4bd2b369be154b676906fa..7fcd4dbe399a29de8b579d1147f799f0b8147367 100644 (file)
@@ -97,7 +97,7 @@ test()
   for (unsigned int level = 0; level < tria.n_global_levels(); ++level)
     {
       MatrixFree<dim> mf_data;
-      data.level_mg_handler = level;
+      data.mg_level = level;
       mf_data.reinit(dof, constraints, quad, data);
       std::vector<unsigned int> n_inner_faces(2 * dim),
         n_inner_other_faces(2 * dim), n_boundary_faces(2 * dim);
index ed1a74a51b03536275d21eaee69bdb3df2f837d2..db3be89600741310baef62bc06bea85dd885ce15 100644 (file)
@@ -109,7 +109,7 @@ test()
         mg_constraints[level].add_line(*bc_it);
       mg_constraints[level].close();
       typename MatrixFree<dim>::AdditionalData data;
-      data.level_mg_handler = level;
+      data.mg_level = level;
       mg_matrices[level].reinit(
         mapping, dof, mg_constraints[level], quad, data);
 
index a9a6550e9b6370873659cb21b4fb5668b2820269..2009ce4c56f8d99cebd859953a915ad0746e35c0 100644 (file)
@@ -77,7 +77,7 @@ public:
     addit_data.tasks_parallel_scheme =
       MatrixFree<dim, number>::AdditionalData::none;
     addit_data.tasks_block_size = 3;
-    addit_data.level_mg_handler = level;
+    addit_data.mg_level         = level;
     addit_data.mapping_update_flags_inner_faces =
       (update_gradients | update_JxW_values);
     addit_data.mapping_update_flags_boundary_faces =
index 7dd71a2fca0a140935c297a52ae67176b7b0a1dc..28a01677a51b505a3de2b7d71426e63f35ce3107 100644 (file)
@@ -76,7 +76,7 @@ public:
     addit_data.tasks_parallel_scheme =
       MatrixFree<dim, number>::AdditionalData::none;
     addit_data.tasks_block_size = 3;
-    addit_data.level_mg_handler = level;
+    addit_data.mg_level         = level;
     addit_data.mapping_update_flags_inner_faces =
       (update_gradients | update_JxW_values);
     addit_data.mapping_update_flags_boundary_faces =
index 980bc419259e0b063c5fa0411e88d9acc021c0b7..728bae208f3bf2e40ad9ed4666b664127a345972 100644 (file)
@@ -77,7 +77,7 @@ public:
     addit_data.tasks_parallel_scheme =
       MatrixFree<dim, number>::AdditionalData::none;
     addit_data.tasks_block_size = 3;
-    addit_data.level_mg_handler = level;
+    addit_data.mg_level         = level;
     addit_data.mapping_update_flags_inner_faces =
       update_JxW_values | update_normal_vectors | update_jacobians;
     addit_data.mapping_update_flags_boundary_faces =
index 9267b213f03bbe846e33a792f25a2b760244c000..a574a58661f20fa034ceeea83301f6a6aa4795eb 100644 (file)
@@ -71,7 +71,7 @@ public:
     typename MatrixFree<dim, number>::AdditionalData addit_data;
     addit_data.tasks_parallel_scheme =
       MatrixFree<dim, number>::AdditionalData::none;
-    addit_data.level_mg_handler = level;
+    addit_data.mg_level = level;
 
     // extract the constraints due to Dirichlet boundary conditions
     AffineConstraints<double>                           constraints;
index 9185a4448d0d81f5f1e602f877dc86fb067e00b5..9c7dc363e787a606a20823160407166cc544321e 100644 (file)
@@ -170,7 +170,7 @@ do_test(const DoFHandler<dim> &dof)
       typename MatrixFree<dim, number>::AdditionalData mg_additional_data;
       mg_additional_data.tasks_parallel_scheme =
         MatrixFree<dim, number>::AdditionalData::none;
-      mg_additional_data.level_mg_handler = level;
+      mg_additional_data.mg_level = level;
 
       AffineConstraints<double> level_constraints;
       IndexSet                  relevant_dofs;
index 388fc77ef5a346c1386c8fd74552b16268c42511..bc981df7d64bd8876fc577b101cafaa26ee75630 100644 (file)
@@ -75,7 +75,7 @@ public:
     addit_data.tasks_parallel_scheme =
       MatrixFree<dim, number>::AdditionalData::none;
     addit_data.tasks_block_size = 3;
-    addit_data.level_mg_handler = level;
+    addit_data.mg_level         = level;
     AffineConstraints<double> constraints;
     if (level == numbers::invalid_unsigned_int)
       {
index 2f2abb724e9fb460653dbfb43b91fe5fb54d7c23..3f4135055d673bd9fac26341feda92ad2c25a5d2 100644 (file)
@@ -220,7 +220,7 @@ do_test(const DoFHandler<dim> &dof)
       mg_additional_data.tasks_parallel_scheme =
         MatrixFree<dim, number>::AdditionalData::none;
       mg_additional_data.tasks_block_size = 3;
-      mg_additional_data.level_mg_handler = level;
+      mg_additional_data.mg_level         = level;
 
       AffineConstraints<double> level_constraints;
       IndexSet                  relevant_dofs;
index b3d2ae6033238717802c7bbbd270b6ba37959004..1f2f3bca58fec61bad7b88ef2101bfd9ad5a7c54 100644 (file)
@@ -81,7 +81,7 @@ public:
       addit_data.tasks_parallel_scheme =
         MatrixFree<dim, number>::AdditionalData::none;
     addit_data.tasks_block_size = 3;
-    addit_data.level_mg_handler = level;
+    addit_data.mg_level         = level;
     AffineConstraints<double> constraints;
     if (level == numbers::invalid_unsigned_int)
       {
index c30c263defe3cf60628592417ce274a2d64af6ab..94b2fce21e1c195debe3470575f3c7cf9fa18770 100644 (file)
@@ -220,7 +220,7 @@ do_test(const DoFHandler<dim> &dof)
       mg_additional_data.tasks_parallel_scheme =
         MatrixFree<dim, number>::AdditionalData::none;
       mg_additional_data.tasks_block_size = 3;
-      mg_additional_data.level_mg_handler = level;
+      mg_additional_data.mg_level         = level;
 
       AffineConstraints<double> level_constraints;
       IndexSet                  relevant_dofs;
index 3d2b4aee9a43c9600dae26849497f8b0130c26ae..172f4c06636d4b8888a15e607fedf18d30dabfd4 100644 (file)
@@ -78,7 +78,7 @@ public:
       addit_data.tasks_parallel_scheme =
         MatrixFree<dim, number>::AdditionalData::none;
     addit_data.tasks_block_size = 3;
-    addit_data.level_mg_handler = level;
+    addit_data.mg_level         = level;
     AffineConstraints<double> constraints;
     if (level == numbers::invalid_unsigned_int)
       {
index 3febd20874b399a7b9e4eead7ab0bf8e02ebd365..589fbf310cfdad1b63b25754849b3e1df67c71d5 100644 (file)
@@ -311,7 +311,7 @@ do_test(const DoFHandler<dim> &dof, const unsigned int nb)
       mg_additional_data.tasks_parallel_scheme =
         MatrixFree<dim, number>::AdditionalData::none;
       mg_additional_data.tasks_block_size = 3;
-      mg_additional_data.level_mg_handler = level;
+      mg_additional_data.mg_level         = level;
 
       AffineConstraints<double> level_constraints;
       IndexSet                  relevant_dofs;
index 4333ee97461081fdede2f1b871fcfb0f329027a3..319fc1ab96061d8b7d061388c84a5d2d25f1b13f 100644 (file)
@@ -221,7 +221,7 @@ do_test(const DoFHandler<dim> &dof)
       mg_additional_data.tasks_parallel_scheme =
         MatrixFree<dim, number>::AdditionalData::none;
       mg_additional_data.tasks_block_size = 3;
-      mg_additional_data.level_mg_handler = level;
+      mg_additional_data.mg_level         = level;
 
       AffineConstraints<double> level_constraints;
       IndexSet                  relevant_dofs;
index 85b3a9670e535eec4cff7e55283e9644ab71e198..e6fc902cea12e084fbfee570b7cc44dbf6ac339a 100644 (file)
@@ -338,7 +338,7 @@ do_test(const std::vector<const DoFHandler<dim> *> &dof)
       mg_additional_data.tasks_parallel_scheme =
         MatrixFree<dim, number>::AdditionalData::none;
       mg_additional_data.tasks_block_size = 3;
-      mg_additional_data.level_mg_handler = level;
+      mg_additional_data.mg_level         = level;
 
       std::vector<AffineConstraints<double>> level_constraints(dof.size());
       std::vector<const AffineConstraints<double> *> level_constraints_ptrs(
index 080c1facd2f0ae562586ce598cdd459ce204070e..e90872503e543f06ae60fd853dc344234efbe724 100644 (file)
@@ -292,7 +292,7 @@ do_test(const DoFHandler<dim> &dof)
       mg_additional_data.tasks_parallel_scheme =
         MatrixFree<dim, number>::AdditionalData::none;
       mg_additional_data.tasks_block_size = 3;
-      mg_additional_data.level_mg_handler = level;
+      mg_additional_data.mg_level         = level;
 
       AffineConstraints<double> level_constraints;
       IndexSet                  relevant_dofs;
index c6ca26874fb63810377ea78a987b5b1d69eeea00..77c57bcdd9888764445d0e717b2ae1fffd512b3f 100644 (file)
@@ -74,7 +74,7 @@ public:
     typename MatrixFree<dim, number>::AdditionalData addit_data;
     addit_data.tasks_parallel_scheme =
       MatrixFree<dim, number>::AdditionalData::none;
-    addit_data.level_mg_handler = level;
+    addit_data.mg_level = level;
 
     // extract the constraints due to Dirichlet boundary conditions
     AffineConstraints<double>                           constraints;
index 2ce591eb191cbac2efa213084e088de74c4b7317..2bf5b1db4cab57553efba8ff874568c89ceecdb7 100644 (file)
@@ -70,7 +70,7 @@ public:
     typename MatrixFree<dim, number>::AdditionalData addit_data;
     addit_data.tasks_parallel_scheme =
       MatrixFree<dim, number>::AdditionalData::none;
-    addit_data.level_mg_handler = level;
+    addit_data.mg_level = level;
 
     // extract the constraints due to Dirichlet boundary conditions
     AffineConstraints<double>                           constraints;
index 838c1900e8eca32da540a19b17628a403006e5a2..64d429a017bdd7f41acae8db3bbf3fd3c7a19be5 100644 (file)
@@ -459,7 +459,7 @@ namespace Step37
           MatrixFree<dim, float>::AdditionalData::none;
         additional_data.mapping_update_flags =
           (update_gradients | update_JxW_values | update_quadrature_points);
-        additional_data.level_mg_handler = level;
+        additional_data.mg_level = level;
         std::shared_ptr<MatrixFree<dim, float>> mg_mf_storage_level(
           new MatrixFree<dim, float>());
         mg_mf_storage_level->reinit(dof_handler,
index a63271f3c1dbfe29b3c260a57de3a09c9ac5b1bf..a3b445ae16f740ec0bc32f3bfcacd4793589d867 100644 (file)
@@ -218,7 +218,7 @@ namespace Step37
     typename MatrixFree<dim, number>::AdditionalData additional_data;
     additional_data.tasks_parallel_scheme =
       MatrixFree<dim, number>::AdditionalData::partition_color;
-    additional_data.level_mg_handler = level;
+    additional_data.mg_level = level;
     additional_data.mapping_update_flags =
       (update_gradients | update_JxW_values | update_quadrature_points);
     data.reinit(dof_handler,
index eb77cb2620f19b67c34869a6e13f45e88af7d996..67aebc6b53bd3f35dca87a077880c78db9c67c26 100644 (file)
@@ -1133,7 +1133,7 @@ namespace StokesClass
           MatrixFree<dim, double>::AdditionalData::none;
         additional_data.mapping_update_flags =
           (update_gradients | update_JxW_values | update_quadrature_points);
-        additional_data.level_mg_handler = level;
+        additional_data.mg_level = level;
         std::shared_ptr<MatrixFree<dim, double>> mg_mf_storage_level(
           new MatrixFree<dim, double>());
         mg_mf_storage_level->reinit(dof_handler_u,
index e8f940a1ef691c8294a0f79e8deb017a6dd2c0e1..5d342b55769f1e1ea9945f36c227fbf2efb360ca 100644 (file)
@@ -269,7 +269,7 @@ namespace Step37
           MatrixFree<dim, float>::AdditionalData::none;
         additional_data.mapping_update_flags =
           (update_gradients | update_JxW_values | update_quadrature_points);
-        additional_data.level_mg_handler = level;
+        additional_data.mg_level = level;
         std::shared_ptr<MatrixFree<dim, float>> mg_mf_storage_level(
           new MatrixFree<dim, float>());
         mg_mf_storage_level->reinit(dof_handler,

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.