diff --git a/Geo/STensor3.h b/Geo/STensor3.h
index 06b0ea5890e09d71210bc3f43aae7995c1c661f2..f310a23ed552bcdc5d7ae6226d60457b233b73ea 100644
--- a/Geo/STensor3.h
+++ b/Geo/STensor3.h
@@ -364,6 +364,13 @@ class STensor3 {
     de(2,2) -= p;
     return de;
   }
+  void eig(fullMatrix<double> &V, fullVector<double> &S, bool s=false) const
+  {
+    fullMatrix<double> me(3, 3), left(3, 3);
+    fullVector<double> im(3);
+    this->getMat(me);
+    me.eig(S, im, left, V, s);
+  }
 };
 
 // tensor product