From 8fd56db8e5a73e1981c04f573d3dd0a9d0510fc3 Mon Sep 17 00:00:00 2001
From: Luca Heltai <luca.heltai@sissa.it>
Date: Fri, 10 May 2019 18:06:11 +0200
Subject: [PATCH] Add constant get_boundary_points().

---
 include/deal.II/base/bounding_box.h | 6 ++++++
 include/deal.II/fe/mapping.h        | 2 +-
 source/base/bounding_box.cc         | 9 +++++++++
 source/fe/mapping.cc                | 2 ++
 4 files changed, 18 insertions(+), 1 deletion(-)

diff --git a/include/deal.II/base/bounding_box.h b/include/deal.II/base/bounding_box.h
index b6962409b1..d0ba528574 100644
--- a/include/deal.II/base/bounding_box.h
+++ b/include/deal.II/base/bounding_box.h
@@ -127,6 +127,12 @@ public:
   std::pair<Point<spacedim, Number>, Point<spacedim, Number>> &
   get_boundary_points();
 
+  /**
+   * Return a const reference to the boundary_points
+   */
+  const std::pair<Point<spacedim, Number>, Point<spacedim, Number>> &
+  get_boundary_points() const;
+
   /**
    * Check if the current object and @p other_bbox are neighbors, i.e. if the boxes
    * have dimension spacedim, check if their intersection is non empty.
diff --git a/include/deal.II/fe/mapping.h b/include/deal.II/fe/mapping.h
index acf799a63f..f28fb3237b 100644
--- a/include/deal.II/fe/mapping.h
+++ b/include/deal.II/fe/mapping.h
@@ -380,7 +380,7 @@ public:
    * true bounding box, since it does not take into account curved faces, and it
    * may be smaller than the true bounding box.
    *
-   * @param[in] cell The cell for which you want to compute the center
+   * @param[in] cell The cell for which you want to compute the bounding box
    *
    * @author Luca Heltai, 2019.
    */
diff --git a/source/base/bounding_box.cc b/source/base/bounding_box.cc
index 006721f78c..12123de5d4 100644
--- a/source/base/bounding_box.cc
+++ b/source/base/bounding_box.cc
@@ -153,6 +153,15 @@ BoundingBox<spacedim, Number>::get_boundary_points()
 
 
 
+template <int spacedim, typename Number>
+const std::pair<Point<spacedim, Number>, Point<spacedim, Number>> &
+BoundingBox<spacedim, Number>::get_boundary_points() const
+{
+  return this->boundary_points;
+}
+
+
+
 template <int spacedim, typename Number>
 double
 BoundingBox<spacedim, Number>::volume() const
diff --git a/source/fe/mapping.cc b/source/fe/mapping.cc
index 3716dac5d1..612d71c9fb 100644
--- a/source/fe/mapping.cc
+++ b/source/fe/mapping.cc
@@ -14,6 +14,8 @@
 // ---------------------------------------------------------------------
 
 
+#include <deal.II/boost_adaptors/bounding_box.h>
+
 #include <deal.II/fe/mapping.h>
 
 #include <deal.II/grid/tria.h>
-- 
2.39.5