coord/ops: add a modulos function

coord/ops.hpp

@@ -883,6 +883,28 @@ namespace cruft {
         return k;
     }
 
+
+    //-------------------------------------------------------------------------
+    template <
+        typename A,
+        typename B,
+        typename = std::enable_if_t<
+            is_coord_v<std::decay_t<A>> && is_coord_v<std::decay_t<B>>
+        >
+    >
+    constexpr auto
+    mod (A &&a, B &&b)
+    {
+        return arithmetic<
+            std::decay_t<A>,
+            std::decay_t<B>
+        >::template eval (
+            std::modulus{},
+            std::forward<A> (a),
+            std::forward<B> (b)
+        );
+    }
+
     ///////////////////////////////////////////////////////////////////////////
     // trigonometric functions
     template <