Fixes

peterrum · Jul 24, 2023 · 26dcf9e · 26dcf9e
1 parent 3a4e6fc
commit 26dcf9e
Show file tree

Hide file tree

Showing 2 changed files with 9 additions and 8 deletions.
diff --git a/include/deal.II/multigrid/mg_transfer_global_coarsening.h b/include/deal.II/multigrid/mg_transfer_global_coarsening.h
@@ -1223,17 +1223,17 @@ class MGTransferBlockGlobalCoarsening
   get_matrix_free_transfer(const unsigned int b) const override;
 
 private:
-  /**
-   * Non-block version of transfer operation.
-   */
-  std::vector<SmartPointer<const MGTransferGlobalCoarsening<dim, VectorType>>>
-    transfer_operators;
-
   /**
    * Internal non-block version of transfer operation.
    */
   std::vector<MGTransferGlobalCoarsening<dim, VectorType>>
     transfer_operators_internal;
+
+  /**
+   * Non-block version of transfer operation.
+   */
+  std::vector<SmartPointer<const MGTransferGlobalCoarsening<dim, VectorType>>>
+    transfer_operators;
 };
 
 
@@ -2051,7 +2051,7 @@ MGTransferGlobalCoarsening<dim, VectorType>::initialize_dof_vector(
 
   if (external_partitioners.empty())
     {
-      vec.reinit(vec_reference);
+      partitioner = vec_reference.get_partitioner();
     }
   else
     {

diff --git a/include/deal.II/multigrid/mg_transfer_global_coarsening.templates.h b/include/deal.II/multigrid/mg_transfer_global_coarsening.templates.h
@@ -3540,6 +3540,7 @@ MGTransferBlockGlobalCoarsening<dim, VectorType>::build(
   const std::vector<const DoFHandler<dim> *> &dof_handler)
 {
   AssertDimension(transfer_operators.size(), dof_handler.size());
+  AssertDimension(transfer_operators_internal.size(), dof_handler.size());
 
   for (unsigned int i = 0; i < dof_handler.size(); ++i)
     this->transfer_operators_internal[i].build(*dof_handler[i]);
@@ -3552,7 +3553,7 @@ const MGTransferGlobalCoarsening<dim, VectorType> &
 MGTransferBlockGlobalCoarsening<dim, VectorType>::get_matrix_free_transfer(
   const unsigned int b) const
 {
-  AssertDimension(b, transfer_operators.size());
+  AssertIndexRange(b, transfer_operators.size());
   return *transfer_operators[b];
 }