diff --git a/include/convex_bodies/ballintersectconvex.h b/include/convex_bodies/ballintersectconvex.h
index b61ae5324..992aae1a1 100644
--- a/include/convex_bodies/ballintersectconvex.h
+++ b/include/convex_bodies/ballintersectconvex.h
@@ -50,6 +50,12 @@ class BallIntersectPolytope {
         return P.get_vec();
     }
 
+    bool is_normalized () {
+        return true;
+    }
+
+    void normalize() {}
+
     int is_in(PointType const& p) const
     {
         if (B.is_in(p)==-1)
diff --git a/include/convex_bodies/correlation_matrices/correlation_spectrahedron.hpp b/include/convex_bodies/correlation_matrices/correlation_spectrahedron.hpp
index a615b33c0..92c38dac6 100755
--- a/include/convex_bodies/correlation_matrices/correlation_spectrahedron.hpp
+++ b/include/convex_bodies/correlation_matrices/correlation_spectrahedron.hpp
@@ -261,6 +261,12 @@ class CorrelationSpectrahedron : public Spectrahedron<Point>{
     MT get_mat() const {
         return MT::Identity(this->d, this->d);
     }
+
+    bool is_normalized () {
+        return true;
+    }
+
+    void normalize() {}
 };
 
 #endif //VOLESTI_CONVEX_BODIES_CORRELATION_MATRICES_VOLESTI_CORRELATION_SPECTRAHEDRON_HPP
\ No newline at end of file
diff --git a/include/convex_bodies/orderpolytope.h b/include/convex_bodies/orderpolytope.h
index 79f41e161..12489f247 100644
--- a/include/convex_bodies/orderpolytope.h
+++ b/include/convex_bodies/orderpolytope.h
@@ -100,6 +100,10 @@ class OrderPolytope {
         return b;
     }
 
+    bool is_normalized ()
+    {
+        return _normalized;
+    }
 
     // print polytope in Ax <= b format
     void print() const
diff --git a/include/convex_bodies/spectrahedra/spectrahedron.h b/include/convex_bodies/spectrahedra/spectrahedron.h
index e8f341382..37b997ea1 100644
--- a/include/convex_bodies/spectrahedra/spectrahedron.h
+++ b/include/convex_bodies/spectrahedra/spectrahedron.h
@@ -341,6 +341,13 @@ class Spectrahedron {
     {
         return MT::Identity(lmi.dimension(), lmi.dimension());
     }
+    
+    bool is_normalized ()
+    {
+        return true;
+    }
+
+    void normalize() {}
 
     void resetFlags()
     {
diff --git a/include/convex_bodies/vpolyintersectvpoly.h b/include/convex_bodies/vpolyintersectvpoly.h
index a6dd97caf..ec2a7c5dd 100644
--- a/include/convex_bodies/vpolyintersectvpoly.h
+++ b/include/convex_bodies/vpolyintersectvpoly.h
@@ -92,6 +92,10 @@ class IntersectionOfVpoly {
         return P1.get_vec();
     }
 
+    bool is_normalized () {
+        return true;
+    }
+
     MT get_T() const {
         return P1.get_mat();
     }
diff --git a/include/convex_bodies/vpolytope.h b/include/convex_bodies/vpolytope.h
index 82ca66891..9100756ef 100644
--- a/include/convex_bodies/vpolytope.h
+++ b/include/convex_bodies/vpolytope.h
@@ -213,6 +213,10 @@ class VPolytope {
         return b;
     }
 
+    bool is_normalized () {
+        return true;
+    }
+
     // change the matrix V
     void set_mat(const MT &V2) {
         V = V2;
diff --git a/include/convex_bodies/zonoIntersecthpoly.h b/include/convex_bodies/zonoIntersecthpoly.h
index bbad43a68..e356db8db 100644
--- a/include/convex_bodies/zonoIntersecthpoly.h
+++ b/include/convex_bodies/zonoIntersecthpoly.h
@@ -65,6 +65,12 @@ class ZonoIntersectHPoly {
         return HP.get_vec();
     }
 
+    bool is_normalized () {
+        return true;
+    }
+
+    void normalize() {}
+
     std::pair<PointType,NT> InnerBall() const
     {
         return Z.InnerBall();
diff --git a/include/convex_bodies/zpolytope.h b/include/convex_bodies/zpolytope.h
index 1ad8868ea..68ce0cf0e 100644
--- a/include/convex_bodies/zpolytope.h
+++ b/include/convex_bodies/zpolytope.h
@@ -299,6 +299,11 @@ class Zonotope {
         return b;
     }
 
+    bool is_normalized ()
+    {
+        return true;
+    }
+
     // change the matrix V
     void set_mat(MT const& V2)
     {