From c24756d1c4343fef0a2dab59863c1c73b6a2fbb6 Mon Sep 17 00:00:00 2001
From: Wolfgang Bangerth <bangerth@colostate.edu>
Date: Tue, 27 Jun 2017 08:23:05 -0600
Subject: [PATCH] Provide the hp::DoFHandler class with a policy object.

---
 include/deal.II/dofs/dof_handler.h |  2 +-
 include/deal.II/hp/dof_handler.h   | 62 +++++++++++++++++++++++++-----
 source/hp/dof_handler.cc           |  1 +
 3 files changed, 55 insertions(+), 10 deletions(-)

diff --git a/include/deal.II/dofs/dof_handler.h b/include/deal.II/dofs/dof_handler.h
index b3862ebe3d..ae3fe3df73 100644
--- a/include/deal.II/dofs/dof_handler.h
+++ b/include/deal.II/dofs/dof_handler.h
@@ -1272,7 +1272,7 @@ DoFHandler<dim,spacedim>::n_boundary_dofs (const std::map<types::boundary_id, co
 namespace internal
 {
   /**
-   * returns a string representing the dynamic type of the given argument.
+   * Return a string representing the dynamic type of the given argument.
    * This is basically the same what typeid(...).name() does, but it turns out
    * this is broken on Intel 13+.
    *
diff --git a/include/deal.II/hp/dof_handler.h b/include/deal.II/hp/dof_handler.h
index 0a5b3bb5d0..679b2ee7f8 100644
--- a/include/deal.II/hp/dof_handler.h
+++ b/include/deal.II/hp/dof_handler.h
@@ -41,6 +41,17 @@ template <int dim, int spacedim> class Triangulation;
 
 namespace internal
 {
+  namespace DoFHandler
+  {
+    struct Implementation;
+
+    namespace Policy
+    {
+      template <int dim, int spacedim> class PolicyBase;
+      struct Implementation;
+    }
+  }
+
   namespace hp
   {
     class DoFLevel;
@@ -734,7 +745,12 @@ namespace hp
      */
     SmartPointer<const hp::FECollection<dim,spacedim>,hp::DoFHandler<dim,spacedim> > finite_elements;
 
-  private:
+    /**
+     * An object that describes how degrees of freedom should be distributed and
+     * renumbered.
+     */
+    std::unique_ptr<dealii::internal::DoFHandler::Policy::PolicyBase<dim,spacedim> > policy;
+
 
     /**
      * Free all used memory.
@@ -919,9 +935,25 @@ namespace hp
     // then just hand everything over to the other function that does the work
     return n_boundary_dofs(boundary_ids_only);
   }
+}
 
 
+namespace internal
+{
+  /**
+   * Return a string representing the dynamic type of the given argument.
+   * This is basically the same what typeid(...).name() does, but it turns out
+   * this is broken on Intel 13+.
+   *
+   * Defined in source/dofs/dof_handler.cc.
+   */
+  template <int dim, int spacedim>
+  std::string policy_to_string(const dealii::internal::DoFHandler::Policy::PolicyBase<dim,spacedim> &policy);
+}
+
 
+namespace hp
+{
   template <int dim, int spacedim>
   template <class Archive>
   void DoFHandler<dim, spacedim>::save(Archive &ar, unsigned int) const
@@ -934,12 +966,16 @@ namespace hp
     ar &has_children;
 
     // write out the number of triangulation cells and later check during
-    // loading that this number is indeed correct;
+    // loading that this number is indeed correct; same with something that
+    // identifies the policy
     unsigned int n_cells = tria->n_cells();
+    std::string  policy_name = dealii::internal::policy_to_string(*policy);
 
-    ar &n_cells;
+    ar &n_cells &policy_name;
   }
 
+
+
   template <int dim, int spacedim>
   template <class Archive>
   void DoFHandler<dim, spacedim>::load(Archive &ar, unsigned int)
@@ -968,12 +1004,20 @@ namespace hp
     // these are the checks that correspond to the last block in the save()
     // function
     unsigned int n_cells;
-
-    ar &n_cells;
-
-    AssertThrow(n_cells == tria->n_cells(),
-                ExcMessage("The object being loaded into does not match the triangulation "
-                           "that has been stored previously."));
+    std::string  policy_name;
+
+    ar &n_cells &policy_name;
+
+    AssertThrow (n_cells == tria->n_cells(),
+                 ExcMessage ("The object being loaded into does not match the triangulation "
+                             "that has been stored previously."));
+    AssertThrow (policy_name == dealii::internal::policy_to_string(*policy),
+                 ExcMessage (std::string ("The policy currently associated with this DoFHandler (")
+                             + dealii::internal::policy_to_string(*policy)
+                             +std::string(") does not match the one that was associated with the "
+                                          "DoFHandler previously stored (")
+                             + policy_name
+                             + ")."));
   }
 
   template <int dim, int spacedim>
diff --git a/source/hp/dof_handler.cc b/source/hp/dof_handler.cc
index 34a11d2df9..55b874de1e 100644
--- a/source/hp/dof_handler.cc
+++ b/source/hp/dof_handler.cc
@@ -20,6 +20,7 @@
 #include <deal.II/hp/dof_level.h>
 #include <deal.II/hp/dof_faces.h>
 #include <deal.II/dofs/dof_accessor.h>
+#include <deal.II/dofs/dof_handler_policy.h>
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/grid/tria_levels.h>
-- 
2.39.5