From: Marc Fehling Date: Fri, 20 Jul 2018 18:22:48 +0000 (-0600) Subject: FIX: hp::DoFHandler::distribute_dofs on p::d::Triangulation. X-Git-Tag: v9.1.0-rc1~905^2~1 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=ed72b29f6a18bd0b38ed863d487d30ec9bb1bcb8;p=dealii.git FIX: hp::DoFHandler::distribute_dofs on p::d::Triangulation. --- diff --git a/doc/news/changes/minor/20180720MarcFehling b/doc/news/changes/minor/20180720MarcFehling new file mode 100644 index 0000000000..aac80a269b --- /dev/null +++ b/doc/news/changes/minor/20180720MarcFehling @@ -0,0 +1,5 @@ +Fixed: hp::DoFHandler::distribute_dofs() now works on +parallel::distributed::Triangulation<3> objects that +contain artificial cells. +
+(Marc Fehling, 2018/07/20) diff --git a/source/hp/dof_handler.cc b/source/hp/dof_handler.cc index 2b8f9c6552..583be7eb3a 100644 --- a/source/hp/dof_handler.cc +++ b/source/hp/dof_handler.cc @@ -705,10 +705,12 @@ namespace internal cell = dof_handler.begin_active(); cell != dof_handler.end(); ++cell) - for (unsigned int l = 0; l < GeometryInfo::lines_per_cell; - ++l) - line_fe_association[cell->active_fe_index()] - [cell->line_index(l)] = true; + if (!cell->is_artificial()) + for (unsigned int l = 0; + l < GeometryInfo::lines_per_cell; + ++l) + line_fe_association[cell->active_fe_index()] + [cell->line_index(l)] = true; // first check which of the lines is used at all, // i.e. is associated with a finite element. we do this