diff --git a/lib/simd/Grid_vector_unops.h b/lib/simd/Grid_vector_unops.h
index e4942937..540a54f5 100644
--- a/lib/simd/Grid_vector_unops.h
+++ b/lib/simd/Grid_vector_unops.h
@@ -67,6 +67,14 @@ namespace Grid {
     }
   };
 
+  template<class scalar> struct DivIntFunctor {
+    Integer y;
+  DivIntFunctor(Integer _y) : y(_y) {};
+    scalar operator()(const scalar &a)  const {
+      return Integer(a)/y;
+    }
+  };
+
   template<class scalar> struct RealFunctor {
     scalar operator()(const scalar &a)  const {
       return real(a);
@@ -131,6 +139,10 @@ namespace Grid {
   inline Grid_simd<S,V> mod(const Grid_simd<S,V> &r,Integer y) {
     return SimdApply(ModIntFunctor<S>(y),r);
   }
+  template < class S, class V > 
+  inline Grid_simd<S,V> div(const Grid_simd<S,V> &r,Integer y) {
+    return SimdApply(DivIntFunctor<S>(y),r);
+  }
   ////////////////////////////////////////////////////////////////////////////
   // Allows us to assign into **conformable** real vectors from complex
   ////////////////////////////////////////////////////////////////////////////