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 (file)
index 0000000..aac80a2
--- /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.
+<br>
+(Marc Fehling, 2018/07/20)

diff --git a/source/hp/dof_handler.cc b/source/hp/dof_handler.cc
index 2b8f9c6552be8d9fba432056ab7a8d37674498db..583be7eb3a0c20f862132367551c0501f6a06b0d 100644 (file)
--- 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<dim>::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<dim>::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