diff --git a/Solver/STensor43.h b/Solver/STensor43.h
index 4f75313505c04b159979f05615c3c865c841a9eb..2ccec282b67cc3a8aab6e199ef73e51c980ec336 100644
--- a/Solver/STensor43.h
+++ b/Solver/STensor43.h
@@ -78,6 +78,11 @@ class STensor43 {
     for (int i = 0; i < 81; i++) _val[i] += other._val[i];
     return *this;
   }
+  STensor43& operator -= (const STensor43 &other)
+  {
+    for (int i = 0; i < 81; i++) _val[i] -= other._val[i];
+    return *this;
+  }
   STensor43& operator *= (const double &other)
   {
     for (int i = 0; i < 81; i++) _val[i] *= other;