From f78d4ee6f82287b13c2a201ab829e2bf6cb721d1 Mon Sep 17 00:00:00 2001 From: Timo Heister Date: Wed, 21 Oct 2015 18:04:49 -0400 Subject: [PATCH] smaller updates --- source/dofs/dof_tools.cc | 2 +- source/multigrid/mg_tools.cc | 11 ++++------- 2 files changed, 5 insertions(+), 8 deletions(-) diff --git a/source/dofs/dof_tools.cc b/source/dofs/dof_tools.cc index ddef986596..55ea3f8664 100644 --- a/source/dofs/dof_tools.cc +++ b/source/dofs/dof_tools.cc @@ -992,7 +992,7 @@ namespace DoFTools endc = dof_handler.end(level); for (; cell!=endc; ++cell) { - types::subdomain_id id = cell->level_subdomain_id(); + const types::subdomain_id id = cell->level_subdomain_id(); // skip artificial and own cells (only look at ghost cells) if (id == dof_handler.get_tria().locally_owned_subdomain() diff --git a/source/multigrid/mg_tools.cc b/source/multigrid/mg_tools.cc index 81c8df6ac7..583410e078 100644 --- a/source/multigrid/mg_tools.cc +++ b/source/multigrid/mg_tools.cc @@ -1511,7 +1511,8 @@ namespace MGTools for (; cell!=endc; ++cell) { - // do not look at artificial cells + // Do not look at artificial level cells (in a serial computation we + // need to ignore the level_subdomain_id() because it is never set). if (mg_dof_handler.get_tria().locally_owned_subdomain()!=numbers::invalid_subdomain_id && cell->level_subdomain_id()==numbers::artificial_subdomain_id) continue; @@ -1529,17 +1530,13 @@ namespace MGTools const typename DoFHandler::cell_iterator neighbor = cell->neighbor(face_nr); - // only process cell pairs if one of them is mine + // only process cell pairs if one or both of them are owned by me (ignore if running in serial) if (mg_dof_handler.get_tria().locally_owned_subdomain()!=numbers::invalid_subdomain_id && neighbor->level_subdomain_id()==numbers::artificial_subdomain_id) -// neighbor->level_subdomain_id()!=mg_dof_handler.get_tria().locally_owned_subdomain() -// && -// cell->level_subdomain_id()!=mg_dof_handler.get_tria().locally_owned_subdomain()) continue; - // Do refinement face - // from the coarse side + // Do refinement face from the coarse side if (neighbor->level() < cell->level()) { for (unsigned int j=0; j