diff --git a/Geo/STensor3.h b/Geo/STensor3.h
index cb427ab81125a16f1b15f309d78e0ef9b46c9372..228648c76862f1c254b6c3a32b00190687899ae0 100644
--- a/Geo/STensor3.h
+++ b/Geo/STensor3.h
@@ -201,22 +201,22 @@ class STensor3 {
     return _index[i][j];
   }
   inline void set_m11(double x){  _val[0] = x; }
-  inline void set_m21(double x){  _val[1] = x; }
-  inline void set_m31(double x){  _val[2] = x; }
-  inline void set_m12(double x){  _val[3] = x; }
+  inline void set_m21(double x){  _val[3] = x; }
+  inline void set_m31(double x){  _val[6] = x; }
+  inline void set_m12(double x){  _val[1] = x; }
   inline void set_m22(double x){  _val[4] = x; }
-  inline void set_m32(double x){  _val[5] = x; }
-  inline void set_m13(double x){  _val[6] = x; }
-  inline void set_m23(double x){  _val[7] = x; }
+  inline void set_m32(double x){  _val[7] = x; }
+  inline void set_m13(double x){  _val[2] = x; }
+  inline void set_m23(double x){  _val[5] = x; }
   inline void set_m33(double x){  _val[8] = x; }
   inline double get_m11(){ return _val[0]; }
-  inline double get_m21(){ return _val[1]; }
-  inline double get_m31(){ return _val[2]; }
-  inline double get_m12(){ return _val[3]; }
+  inline double get_m21(){ return _val[3]; }
+  inline double get_m31(){ return _val[6]; }
+  inline double get_m12(){ return _val[1]; }
   inline double get_m22(){ return _val[4]; }
-  inline double get_m32(){ return _val[5]; }
-  inline double get_m13(){ return _val[6]; }
-  inline double get_m23(){ return _val[7]; }
+  inline double get_m32(){ return _val[7]; }
+  inline double get_m13(){ return _val[2]; }
+  inline double get_m23(){ return _val[5]; }
   inline double get_m33(){ return _val[8]; }
   void getMat(fullMatrix<double> &mat) const
   {