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.
* 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.
*/
+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
// ---------------------------------------------------------------------
+#include <deal.II/boost_adaptors/bounding_box.h>
+
#include <deal.II/fe/mapping.h>
#include <deal.II/grid/tria.h>