From 117cbdbd726bd02ef6e5356445ea9f98df67ff53 Mon Sep 17 00:00:00 2001 From: Peter Munch Date: Tue, 3 Aug 2021 21:10:01 +0200 Subject: [PATCH] Continue --- include/deal.II/matrix_free/hanging_nodes_internal.h | 10 +++++++--- 1 file changed, 7 insertions(+), 3 deletions(-) diff --git a/include/deal.II/matrix_free/hanging_nodes_internal.h b/include/deal.II/matrix_free/hanging_nodes_internal.h index bd81181863..88a1c7c1e9 100644 --- a/include/deal.II/matrix_free/hanging_nodes_internal.h +++ b/include/deal.II/matrix_free/hanging_nodes_internal.h @@ -894,7 +894,7 @@ namespace CUDAWrappers const unsigned int x_idx = threadIdx.x % (fe_degree + 1); const unsigned int y_idx = threadIdx.y; - const unsigned int this_type = + const auto this_type = (direction == 0) ? dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_x : dealii::internal::MatrixFreeFunctions::ConstraintKinds::type_y; @@ -1074,7 +1074,9 @@ namespace CUDAWrappers ConstraintKinds::unconstrained) && on_face1 && on_face2)); - if (constrained_face && constrained_dof) + if ((constrained_face != dealii::internal::MatrixFreeFunctions:: + ConstraintKinds::unconstrained) && + constrained_dof) { const bool type = (constraint_mask & this_type) != dealii::internal::MatrixFreeFunctions:: @@ -1120,7 +1122,9 @@ namespace CUDAWrappers // each block being assigned to one element. __syncthreads(); - if (constrained_face && constrained_dof) + if ((constrained_face != dealii::internal::MatrixFreeFunctions:: + ConstraintKinds::unconstrained) && + constrained_dof) values[index3(x_idx, y_idx, z_idx)] = t; __syncthreads(); -- 2.39.5