Keep calculateBounds in PointView but redirect to EigenUtils for now

pdal/PointView.cpp

@@ -33,16 +33,13 @@
 ****************************************************************************/
 
 #include <iomanip>
-#include <numeric>
 
 #include <pdal/EigenUtils.hpp>
 #include <pdal/KDIndex.hpp>
 #include <pdal/PointView.hpp>
 #include <pdal/PointViewIter.hpp>
 #include <pdal/util/Algorithm.hpp>
 
-#include <Eigen/Dense>
-
 namespace pdal
 {
 
@@ -102,6 +99,18 @@ void PointView::setFieldInternal(Dimension::Id dim, PointId idx,
 }
 
 
+void PointView::calculateBounds(BOX2D& output) const
+{
+    pdal::calculateBounds(const_cast<PointView&>(*this), output);
+}
+
+
+void PointView::calculateBounds(BOX3D& output) const
+{
+    pdal::calculateBounds(const_cast<PointView&>(*this), output);
+}
+
+
 MetadataNode PointView::toMetadata() const
 {
     MetadataNode node;

pdal/PointView.hpp

@@ -182,6 +182,19 @@ class PDAL_DLL PointView : public PointContainer
         getFieldInternal(dim, idx, buf);
     }
 
+    /*! @return a cumulated bounds of all points in the PointView.
+        \verbatim embed:rst
+        .. note::
+
+            This method requires that an `X`, `Y`, and `Z` dimension be
+            available, and that it can be casted into a *double* data
+            type using the :cpp:func:`pdal::Dimension::applyScaling`
+            method. Otherwise, an exception will be thrown.
+        \endverbatim
+    */
+    void calculateBounds(BOX2D& box) const;
+    void calculateBounds(BOX3D& box) const;
+
     void dump(std::ostream& ostr) const;
     bool hasDim(Dimension::Id id) const
         { return layout()->hasDim(id); }