From b6892be846dbe7258e5253a47734bd69e3f22a0b Mon Sep 17 00:00:00 2001 From: Peter Munch Date: Sun, 5 Jan 2020 00:51:27 +0100 Subject: [PATCH] Fix check of neigboring cells --- include/deal.II/matrix_free/mapping_info.templates.h | 8 ++++++-- 1 file changed, 6 insertions(+), 2 deletions(-) diff --git a/include/deal.II/matrix_free/mapping_info.templates.h b/include/deal.II/matrix_free/mapping_info.templates.h index 39be9ccbbe..e7357c76bb 100644 --- a/include/deal.II/matrix_free/mapping_info.templates.h +++ b/include/deal.II/matrix_free/mapping_info.templates.h @@ -2048,14 +2048,18 @@ namespace internal cells[cell * vectorization_width + v].second); fe_val.reinit(cell_it, face); - const bool is_local = cell_it->is_locally_owned(); + const bool is_local = + cell_it->is_locally_owned() && + (!cell_it->at_boundary(face) || + (cell_it->at_boundary(face) && + cell_it->has_periodic_neighbor(face))); if (is_local) { auto cell_it_neigh = cell_it->neighbor_or_periodic_neighbor(face); fe_val_neigh.reinit(cell_it_neigh, - cell_it->at_boundary() ? + cell_it->at_boundary(face) ? cell_it->periodic_neighbor_face_no( face) : cell_it->neighbor_face_no(face)); -- 2.39.5