From: Marc Fehling <marc.fehling@gmx.net>
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.
+<br>
+(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<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