diff --git a/Solver/STensor43.h b/Solver/STensor43.h
index 5e3ecfe7f0a0bfa4bea188856231750c9f3a6679..c764fde277f0aeba8b4fb61c0c364ec7e42a9fd6 100644
--- a/Solver/STensor43.h
+++ b/Solver/STensor43.h
@@ -37,7 +37,8 @@ class STensor43 {
             else
               _val[getIndex(i, j, k, l)]=0.0;
   }
-  inline void symmetricIdentity(const double v = 1.0)
+  // Symmetric identity tensor
+  STensor43(const double vik, const double vil)
   {
     for (int i = 0; i < 3; i++)
       for (int j = 0; j < 3; j++)
@@ -46,9 +47,9 @@ class STensor43 {
           {
             _val[getIndex(i, j, k, l)]= 0.;
             if ((i==k)&&(j==l))
-              _val[getIndex(i, j, k, l)]+=0.5*v;
+              _val[getIndex(i, j, k, l)]+=0.5*vik;
             if ((i==l)&&(j==k))
-              _val[getIndex(i, j, k, l)]+=0.5*v;
+              _val[getIndex(i, j, k, l)]+=0.5*vil;
           }
   }
   inline double &operator()(int i, int j,int k, int l)