diff --git a/Geo/SVector3.h b/Geo/SVector3.h
index f58a069cf0c42467d2789d12631fd5714deddc72..8e36d65777ecf603ff6cbd433d4b21aafa0c6e3c 100644
--- a/Geo/SVector3.h
+++ b/Geo/SVector3.h
@@ -35,6 +35,7 @@ class SVector3 {
     double n = norm(); if(n){ P[0] /= n; P[1] /= n; P[2] /= n; }
     return n;
   }
+  SVector3 unit() const{ SVector3 y(*this); y.normalize(); return y; }
   void negate() { P[0] = -P[0]; P[1] = -P[1]; P[2] = -P[2]; }
   // why both [] and (), why not
   double &operator[](int i){ return P[i]; }