diff --git a/point.hpp b/point.hpp
index 01e2be12..7929e7a7 100644
--- a/point.hpp
+++ b/point.hpp
@@ -203,6 +203,26 @@ namespace util {
         util::point<S,T>
     >
     furthest (util::view<const util::point<S,T>*>);
+
+
+    ///////////////////////////////////////////////////////////////////////////
+    /// computes the mean point across a view of points
+    template <typename InputT>
+    auto
+    center (util::view<InputT> points)
+    {
+        CHECK_NEZ (points.size ());
+
+        using point_type = typename std::iterator_traits<InputT>::value_type;
+        using value_type = typename point_type::value_type;
+
+        util::vector<point_type::elements,value_type> accum = 0;
+
+        for (auto const &i: points)
+            accum += i.template as<util::vector> ();
+
+        return (accum / points.size ()).template as<util::point> ();
+    }
 }
 
 #endif // __UTIL_POINT_HPP