Fix FEInterface for MG support

Indetation

diff --git a/include/deal.II/fe/fe_interface_values.h b/include/deal.II/fe/fe_interface_values.h
index 5231b202b0f90d9649803205345992358c30faeb..78727fb9c3687e8d4dbcfec4d85d35f4f3e85463 100644 (file)
--- a/include/deal.II/fe/fe_interface_values.h
+++ b/include/deal.II/fe/fe_interface_values.h
@@ -228,6 +228,10 @@ public:
 
   /**
    * Return the set of joint DoF indices. This includes indices from both cells.
+   * If reinit was called with an active cell iterator, the indices are based
+   * on the active indices (returned by `DoFCellAccessor::get_dof_indices()` ),
+   * in case of level cell (that is, if is_level_cell() return true )
+   * the mg dof indices are returned.
    *
    * @note This function is only available after a call to reinit() and can change
    * from one call to reinit() to the next.
@@ -542,10 +546,12 @@ FEInterfaceValues<dim, spacedim>::reinit(
     // Get dof indices first:
     std::vector<types::global_dof_index> v(
       fe_face_values->get_fe().n_dofs_per_cell());
-    cell->get_dof_indices(v);
+    cell->get_active_or_mg_dof_indices(v);
     std::vector<types::global_dof_index> v2(
       fe_face_values_neighbor->get_fe().n_dofs_per_cell());
-    cell_neighbor->get_dof_indices(v2);
+    cell_neighbor->get_active_or_mg_dof_indices(v2);
+
+
 
     // Fill a map from the global dof index to the left and right
     // local index.
@@ -594,7 +600,7 @@ FEInterfaceValues<dim, spacedim>::reinit(const CellIteratorType &cell,
   fe_face_values_neighbor = nullptr;
 
   interface_dof_indices.resize(fe_face_values->get_fe().n_dofs_per_cell());
-  cell->get_dof_indices(interface_dof_indices);
+  cell->get_active_or_mg_dof_indices(interface_dof_indices);
 
 
   dofmap.resize(interface_dof_indices.size());