From f75a11db8c520efc834cde768343939dfd70d72e Mon Sep 17 00:00:00 2001
From: Wolfgang Bangerth <bangerth@math.tamu.edu>
Date: Sun, 25 Jan 2015 12:22:42 -0600
Subject: [PATCH] Remove deprecated function
 MGTools::extract_inner_interface_dofs.

---
 doc/news/changes.h                            |  1 +
 include/deal.II/multigrid/mg_tools.h          |  9 --------
 .../deal.II/multigrid/mg_transfer_component.h |  8 +++----
 source/multigrid/mg_tools.cc                  | 23 -------------------
 source/multigrid/mg_tools.inst.in             |  4 ----
 source/multigrid/mg_transfer_component.cc     |  7 ++++--
 6 files changed, 10 insertions(+), 42 deletions(-)

diff --git a/doc/news/changes.h b/doc/news/changes.h
index 977ecd63df..92e96e6ff3 100644
--- a/doc/news/changes.h
+++ b/doc/news/changes.h
@@ -205,6 +205,7 @@ inconvenience this causes.
   - PreconditionBlock::size.
   - MGTools::count_dofs_per_component.
   - MGTools::apply_boundary_values.
+  - MGTools::extract_inner_interface_dofs.
   <br>
   This release also removes the deprecated class MGDoFHandler. The
   functionality of this class had previously been incorporated into
diff --git a/include/deal.II/multigrid/mg_tools.h b/include/deal.II/multigrid/mg_tools.h
index 85277d5bf2..a6a69dfe6f 100644
--- a/include/deal.II/multigrid/mg_tools.h
+++ b/include/deal.II/multigrid/mg_tools.h
@@ -213,15 +213,6 @@ namespace MGTools
   extract_inner_interface_dofs (const DoFHandler<dim,spacedim> &mg_dof_handler,
                                 std::vector<IndexSet>  &interface_dofs);
 
-  /**
-   * As above but with a deprecated data structure. This makes one additional
-   * copy.
-   */
-  template <int dim, int spacedim>
-  void
-  extract_inner_interface_dofs (const DoFHandler<dim,spacedim> &mg_dof_handler,
-                                std::vector<std::vector<bool> >  &interface_dofs) DEAL_II_DEPRECATED;
-
 
   template <int dim, int spacedim>
   void
diff --git a/include/deal.II/multigrid/mg_transfer_component.h b/include/deal.II/multigrid/mg_transfer_component.h
index b9d9d1e05b..9dfc70f5c0 100644
--- a/include/deal.II/multigrid/mg_transfer_component.h
+++ b/include/deal.II/multigrid/mg_transfer_component.h
@@ -346,12 +346,12 @@ private:
   unsigned int mg_selected_component;
 
   /**
-   * The degrees of freedom on the the refinement edges. For each level (outer
-   * vector) and each dof index (inner vector), this bool is true if the level
-   * degree of freedom is on the refinement edge towards the lower level
+   * The degrees of freedom on the refinement edges. For each level the index
+   * set denotes which level
+   * degrees of freedom are on the refinement edge towards the lower level,
    * excluding boundary dofs.
    */
-  std::vector<std::vector<bool> > interface_dofs;
+  std::vector<IndexSet> interface_dofs;
 
   /**
    * The constraints of the global system.
diff --git a/source/multigrid/mg_tools.cc b/source/multigrid/mg_tools.cc
index 68b1a2963c..70025fbb61 100644
--- a/source/multigrid/mg_tools.cc
+++ b/source/multigrid/mg_tools.cc
@@ -1405,29 +1405,6 @@ namespace MGTools
   }
 
 
-  template <int dim, int spacedim>
-  void
-  extract_inner_interface_dofs (const DoFHandler<dim,spacedim> &mg_dof_handler,
-                                std::vector<std::vector<bool> >  &interface_dofs)
-  {
-    std::vector<IndexSet> temp;
-    temp.resize(interface_dofs.size());
-    for (unsigned int l=0; l<interface_dofs.size(); ++l)
-      temp[l] = IndexSet(interface_dofs[l].size());
-
-    extract_inner_interface_dofs(mg_dof_handler, temp);
-
-    for (unsigned int l=0; l<interface_dofs.size(); ++l)
-      {
-        Assert (interface_dofs[l].size() == mg_dof_handler.n_dofs(l),
-                ExcDimensionMismatch (interface_dofs[l].size(),
-                                      mg_dof_handler.n_dofs(l)));
-
-        temp[l].fill_binary_vector(interface_dofs[l]);
-      }
-  }
-
-
   template <int dim, int spacedim>
   void
   extract_non_interface_dofs (const DoFHandler<dim,spacedim> &mg_dof_handler,
diff --git a/source/multigrid/mg_tools.inst.in b/source/multigrid/mg_tools.inst.in
index 8a5981891f..329b7ce229 100644
--- a/source/multigrid/mg_tools.inst.in
+++ b/source/multigrid/mg_tools.inst.in
@@ -106,10 +106,6 @@ for (deal_II_dimension : DIMENSIONS)
 	void
 	extract_inner_interface_dofs (const DoFHandler<deal_II_dimension> &mg_dof_handler,
 				      std::vector<IndexSet>  &interface_dofs);
-      template
-	void
-	extract_inner_interface_dofs (const DoFHandler<deal_II_dimension> &mg_dof_handler,
-				      std::vector<std::vector<bool> >  &interface_dofs);
 
       template
         void
diff --git a/source/multigrid/mg_transfer_component.cc b/source/multigrid/mg_transfer_component.cc
index af335e2247..185706e5d3 100644
--- a/source/multigrid/mg_transfer_component.cc
+++ b/source/multigrid/mg_transfer_component.cc
@@ -614,7 +614,10 @@ void MGTransferSelect<number>::build_matrices (
 
   interface_dofs.resize(mg_dof.get_tria().n_levels());
   for (unsigned int l=0; l<mg_dof.get_tria().n_levels(); ++l)
-    interface_dofs[l].resize(mg_dof.n_dofs(l));
+    {
+      interface_dofs[l].clear();
+      interface_dofs[l].set_size(mg_dof.n_dofs(l));
+    }
   MGTools::extract_inner_interface_dofs(mg_dof, interface_dofs);
 
   // use a temporary vector to create the
@@ -649,7 +652,7 @@ void MGTransferSelect<number>::build_matrices (
               const unsigned int component
                 = fe.system_to_component_index(i).first;
               if (component_mask[component] &&
-                  !interface_dofs[level][level_dof_indices[i]])
+                  !interface_dofs[level].is_element(level_dof_indices[i]))
                 {
                   const types::global_dof_index level_start
                     = mg_component_start[level][mg_target_component[component]];
-- 
2.39.5