coord: add product function

coord/ops.hpp

@@ -1098,6 +1098,21 @@ namespace util {
     }
 
 
+    //-------------------------------------------------------------------------
+    template <
+        typename K,
+        typename = std::enable_if_t<is_coord_v<K>>
+    >
+    auto
+    product (K const& k)
+    {
+        typename K::value_type accum = 1;
+        for (auto i: k)
+            accum *= i;
+        return accum;
+    }
+
+
     ///////////////////////////////////////////////////////////////////////////
 #define VECTOR_OP(OP)                                   \
     template <                                          \