Add the possibility to build bounding boxes with containers of points.

diff --git a/include/deal.II/base/bounding_box.h b/include/deal.II/base/bounding_box.h
index 8922fa6b367626a305b38af4c0e9caeb800611a6..b6962409b166c9ead86416e59290aa0fb665de43 100644 (file)
--- a/include/deal.II/base/bounding_box.h
+++ b/include/deal.II/base/bounding_box.h
@@ -67,7 +67,7 @@ enum class NeighborType
  *
  * Objects of this class are used to represent bounding boxes. They are,
  * among other uses, useful in parallel distributed meshes to give a general
- * description the owners of each portion of the mesh.
+ * description of the owners of each portion of the mesh.
  *
  * Bounding boxes are represented by two vertices (bottom left and top right).
  * Geometrically, a bounding box is:
@@ -89,7 +89,7 @@ enum class NeighborType
  * |      |/
  * V------.
  * @endcode
- * Notice the sides are always parallel to the respective axis.
+ * Notice that the sides are always parallel to the respective axis.
  *
  * @author Giovanni Alzetta, 2017.
  */
@@ -112,10 +112,20 @@ public:
                 &boundary_points);
 
   /**
-   * Return the boundary_points
+   * Construct the bounding box that encloses all the points in the given
+   * container.
+   *
+   * The constructor supports any Container that provides begin() and end()
+   * iterators to Point<spacedim, Number> elements.
+   */
+  template <class Container>
+  BoundingBox(const Container &points);
+
+  /**
+   * Return a reference to the boundary_points
    */
-  const std::pair<Point<spacedim, Number>, Point<spacedim, Number>> &
-  get_boundary_points() const;
+  std::pair<Point<spacedim, Number>, Point<spacedim, Number>> &
+  get_boundary_points();
 
   /**
    * Check if the current object and @p other_bbox are neighbors, i.e. if the boxes
@@ -189,6 +199,20 @@ inline BoundingBox<spacedim, Number>::BoundingBox(
   this->boundary_points = boundary_points;
 }
 
+
+
+template <int spacedim, typename Number>
+template <class Container>
+inline BoundingBox<spacedim, Number>::BoundingBox(const Container &points)
+{
+  boost::geometry::envelope(
+    boost::geometry::model::multi_point<Point<spacedim, Number>>(points.begin(),
+                                                                 points.end()),
+    *this);
+}
+
+
+
 template <int spacedim, typename Number>
 inline void
 BoundingBox<spacedim, Number>::extend(const Number &amount)

diff --git a/include/deal.II/boost_adaptors/bounding_box.h b/include/deal.II/boost_adaptors/bounding_box.h
index a85f573aa9435fc51b7c0c1481d5860ef7c10ef7..63f4be68658126434015ba1d82e969e810471a6d 100644 (file)
--- a/include/deal.II/boost_adaptors/bounding_box.h
+++ b/include/deal.II/boost_adaptors/bounding_box.h
@@ -21,6 +21,7 @@
 #include <deal.II/boost_adaptors/point.h>
 
 #include <boost/geometry.hpp>
+#include <boost/geometry/geometries/geometries.hpp>
 
 namespace boost
 {

diff --git a/source/base/bounding_box.cc b/source/base/bounding_box.cc
index 410d610d3590d2bcbb8343c9f3440deb7d911997..006721f78caea451e360bf20abb09e4c464b3d59 100644 (file)
--- a/source/base/bounding_box.cc
+++ b/source/base/bounding_box.cc
@@ -142,13 +142,17 @@ BoundingBox<spacedim, Number>::get_neighbor_type(
     }
 }
 
+
+
 template <int spacedim, typename Number>
-const std::pair<Point<spacedim, Number>, Point<spacedim, Number>> &
-BoundingBox<spacedim, Number>::get_boundary_points() const
+std::pair<Point<spacedim, Number>, Point<spacedim, Number>> &
+BoundingBox<spacedim, Number>::get_boundary_points()
 {
   return this->boundary_points;
 }
 
+
+
 template <int spacedim, typename Number>
 double
 BoundingBox<spacedim, Number>::volume() const