Add another barrier in MatrixFreeTools::compute_diagonal()

diff --git a/include/deal.II/matrix_free/tools.h b/include/deal.II/matrix_free/tools.h
index 01c9b19ecd0476e625d904d13f1f6b5283bd0925..b634cd0144c2a5c4797859a3158b557746a670b4 100644 (file)
--- a/include/deal.II/matrix_free/tools.h
+++ b/include/deal.II/matrix_free/tools.h
@@ -1409,6 +1409,8 @@ namespace MatrixFreeTools
             Kokkos::TeamThreadRange(shared_data->team_member, dofs_per_cell),
             [&](int j) { fe_eval.submit_dof_value(i == j ? 1 : 0, j); });
 
+          shared_data->team_member.team_barrier();
+
           Portable::internal::
             resolve_hanging_nodes<dim, fe_degree, false, Number>(
               shared_data->team_member,