diff --git a/src/include/operation/op_types.cc.h b/src/include/operation/op_types.cc.h
index 4e246cd..a3f91e4 100644
--- a/src/include/operation/op_types.cc.h
+++ b/src/include/operation/op_types.cc.h
@@ -237,6 +237,14 @@ namespace CNORXZ
         return *this;
     }
 
+    template <typename T, class IndexT>
+    template <class Op>
+    constexpr OpCont<T,IndexT>& OpCont<T,IndexT>::operator-=(const Op& o)
+    {
+	OI::a(mIndex, [](auto& a1, const auto& a2) { a1 -= a2; }, o);
+        return *this;
+    }
+
     template <typename T, class IndexT>
     constexpr OpCont<T,IndexT>& OpCont<T,IndexT>::operator=(const OpCont<T,IndexT>& o)
     {
@@ -328,6 +336,14 @@ namespace CNORXZ
         return *this;
     }
 
+    template <typename T, class IndexT>
+    template <class Op>
+    constexpr OpRoot<T,IndexT>& OpRoot<T,IndexT>::operator-=(const Op& o)
+    {
+	OI::a(mIndex, [](auto& a, const auto& b) { a -= b; }, o);
+        return *this;
+    }
+
     template <typename T, class IndexT>
     constexpr OpRoot<T,IndexT>& OpRoot<T,IndexT>::operator=(const OpRoot<T,IndexT>& o)
     {
diff --git a/src/include/operation/op_types.h b/src/include/operation/op_types.h
index d60ce71..bb516bb 100644
--- a/src/include/operation/op_types.h
+++ b/src/include/operation/op_types.h
@@ -159,6 +159,9 @@ namespace CNORXZ
         template <class Op>
         constexpr OpCont& operator+=(const Op& in);
 
+        template <class Op>
+        constexpr OpCont& operator-=(const Op& in);
+
         constexpr OpCont& operator=(const OpCont& in);
 
 	template <class PosT>
@@ -196,6 +199,9 @@ namespace CNORXZ
         template <class Op>
         constexpr OpRoot& operator+=(const Op& in);
 
+        template <class Op>
+        constexpr OpRoot& operator-=(const Op& in);
+
         constexpr OpRoot& operator=(const OpRoot& in);
         
 	template <class PosT>