From ed72b29f6a18bd0b38ed863d487d30ec9bb1bcb8 Mon Sep 17 00:00:00 2001 From: Marc Fehling Date: Fri, 20 Jul 2018 12:22:48 -0600 Subject: [PATCH] FIX: hp::DoFHandler::distribute_dofs on p::d::Triangulation. --- doc/news/changes/minor/20180720MarcFehling | 5 +++++ source/hp/dof_handler.cc | 10 ++++++---- 2 files changed, 11 insertions(+), 4 deletions(-) create mode 100644 doc/news/changes/minor/20180720MarcFehling 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 -- 2.39.5