From: Timo Heister Date: Wed, 21 Oct 2015 22:04:49 +0000 (-0400) Subject: smaller updates X-Git-Tag: v8.4.0-rc2~277^2~6 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=f78d4ee6f82287b13c2a201ab829e2bf6cb721d1;p=dealii.git smaller updates --- 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