MGTwoLevelTransfer: Make sure ot clean old vector content in prolonga…

…te_add
dealii · May 31, 2023 · bd91cb6 · bd91cb6
1 parent 2896a7e
commit bd91cb6
Show file tree

Hide file tree

Showing 3 changed files with 213 additions and 1 deletion.
diff --git a/include/deal.II/multigrid/mg_transfer_global_coarsening.templates.h b/include/deal.II/multigrid/mg_transfer_global_coarsening.templates.h
@@ -2630,7 +2630,7 @@ MGTwoLevelTransferBase<LinearAlgebra::distributed::Vector<Number>>::
 
   this->update_ghost_values(*vec_coarse_ptr);
 
-  if (fine_element_is_continuous && (use_dst_inplace == false))
+  if (use_dst_inplace == false)
     *vec_fine_ptr = Number(0.);
 
   this->prolongate_and_add_internal(*vec_fine_ptr, *vec_coarse_ptr);

diff --git a/tests/multigrid-global-coarsening/mg_transfer_a_01.cc b/tests/multigrid-global-coarsening/mg_transfer_a_01.cc
@@ -112,6 +112,10 @@ do_test(const FiniteElement<dim> &fe_fine, const FiniteElement<dim> &fe_coarse)
                   constraint_fine,
                   constraint_coarse);
 
+  deallog << "test first time" << std::endl;
+  test_transfer_operator(transfer, dof_handler_fine, dof_handler_coarse);
+
+  deallog << "test second time" << std::endl;
   test_transfer_operator(transfer, dof_handler_fine, dof_handler_coarse);
 }