diff --git a/gtsam/base/numericalDerivative.h b/gtsam/base/numericalDerivative.h
index 053a89da8..036f23f68 100644
--- a/gtsam/base/numericalDerivative.h
+++ b/gtsam/base/numericalDerivative.h
@@ -59,7 +59,7 @@ namespace gtsam {
 
   /** global functions for converting to a LieVector for use with numericalDerivative */
   inline LieVector makeLieVector(const Vector& v) { return LieVector(v); }
-  inline LieVector makeLieVectorD(double d) { return LieVector(Vector_(1, d)); }
+  inline LieVector makeLieVectorD(double d) { return LieVector((Vector) (Vector(1) << d)); }
 
   /**
    * Numerically compute gradient of scalar function
diff --git a/gtsam/base/tests/testLieMatrix.cpp b/gtsam/base/tests/testLieMatrix.cpp
index b21756dc4..88bae697a 100644
--- a/gtsam/base/tests/testLieMatrix.cpp
+++ b/gtsam/base/tests/testLieMatrix.cpp
@@ -50,7 +50,7 @@ TEST( LieMatrix, other_constructors ) {
 /* ************************************************************************* */
 TEST(LieMatrix, retract) {
   LieMatrix init(2,2, 1.0,2.0,3.0,4.0);
-  Vector update = Vector_(4, 3.0, 4.0, 6.0, 7.0);
+  Vector update = (Vector(4) << 3.0, 4.0, 6.0, 7.0);
 
   LieMatrix expected(2,2, 4.0, 6.0, 9.0, 11.0);
   LieMatrix actual = init.retract(update);
diff --git a/gtsam/base/tests/testLieScalar.cpp b/gtsam/base/tests/testLieScalar.cpp
index 11157a701..68655cc71 100644
--- a/gtsam/base/tests/testLieScalar.cpp
+++ b/gtsam/base/tests/testLieScalar.cpp
@@ -42,7 +42,7 @@ TEST( testLieScalar, construction ) {
 TEST( testLieScalar, localCoordinates ) {
   LieScalar lie1(1.), lie2(3.);
 
-  EXPECT(assert_equal(Vector_(1, 2.), lie1.localCoordinates(lie2)));
+  EXPECT(assert_equal((Vector)(Vector(1) << 2), lie1.localCoordinates(lie2)));
 }
 
 /* ************************************************************************* */
diff --git a/gtsam/geometry/Point2.h b/gtsam/geometry/Point2.h
index 0838650ea..dc9a1dac8 100644
--- a/gtsam/geometry/Point2.h
+++ b/gtsam/geometry/Point2.h
@@ -172,7 +172,7 @@ public:
   static inline Point2 Expmap(const Vector& v) { return Point2(v); }
 
   /// Log map around identity - just return the Point2 as a vector
-  static inline Vector Logmap(const Point2& dp) { return Vector_(2, dp.x(), dp.y()); }
+  static inline Vector Logmap(const Point2& dp) { return (Vector(2) << dp.x(), dp.y()); }
 
   /// @}
   /// @name Vector Space
diff --git a/gtsam/geometry/Pose3.cpp b/gtsam/geometry/Pose3.cpp
index adbe763b3..bfd2fcb9a 100644
--- a/gtsam/geometry/Pose3.cpp
+++ b/gtsam/geometry/Pose3.cpp
@@ -98,7 +98,7 @@ Vector Pose3::adjointTranspose(const Vector& xi, const Vector& y,
 /* ************************************************************************* */
 Matrix6 Pose3::dExpInv_exp(const Vector& xi) {
   // Bernoulli numbers, from Wikipedia
-  static const Vector B = Vector_(9, 1.0, -1.0 / 2.0, 1. / 6., 0.0, -1.0 / 30.0,
+  static const Vector B = (Vector(9) << 1.0, -1.0 / 2.0, 1. / 6., 0.0, -1.0 / 30.0,
       0.0, 1.0 / 42.0, 0.0, -1.0 / 30);
   static const int N = 5; // order of approximation
   Matrix res = I6;
diff --git a/gtsam/geometry/Rot2.h b/gtsam/geometry/Rot2.h
index 45012770f..d121beb12 100644
--- a/gtsam/geometry/Rot2.h
+++ b/gtsam/geometry/Rot2.h
@@ -170,7 +170,7 @@ namespace gtsam {
 
     ///Log map at identity - return the canonical coordinates of this rotation
     static inline Vector Logmap(const Rot2& r) {
-      return Vector_(1, r.theta());
+      return (Vector(1) << r.theta());
     }
 
     /// @}
diff --git a/gtsam/geometry/Rot3.h b/gtsam/geometry/Rot3.h
index d302a3502..847d31d65 100644
--- a/gtsam/geometry/Rot3.h
+++ b/gtsam/geometry/Rot3.h
@@ -194,7 +194,7 @@ namespace gtsam {
      * @return incremental rotation matrix
      */
     static Rot3 rodriguez(double wx, double wy, double wz)
-      { return rodriguez(Vector_(3,wx,wy,wz));}
+      { return rodriguez((Vector(3) << wx, wy, wz));}
 
     /// @}
     /// @name Testable
diff --git a/gtsam/geometry/Sphere2.cpp b/gtsam/geometry/Sphere2.cpp
index ce386ac17..3c7765bd2 100644
--- a/gtsam/geometry/Sphere2.cpp
+++ b/gtsam/geometry/Sphere2.cpp
@@ -170,7 +170,7 @@ Sphere2 Sphere2::retract(const Vector& v, Sphere2::CoordinatesMode mode) const {
     double alpha = p.transpose() * q;
     assert(alpha != 0.0);
     Matrix coeffs = (B.transpose() * q) / alpha;
-    Vector result = Vector_(2, coeffs(0, 0), coeffs(1, 0));
+    Vector result = (Vector(2) << coeffs(0, 0), coeffs(1, 0));
     return result;
   } else {
     assert (false);
diff --git a/gtsam/geometry/tests/testCal3DS2.cpp b/gtsam/geometry/tests/testCal3DS2.cpp
index c73ae1182..53a21a9fb 100644
--- a/gtsam/geometry/tests/testCal3DS2.cpp
+++ b/gtsam/geometry/tests/testCal3DS2.cpp
@@ -36,7 +36,7 @@ TEST( Cal3DS2, uncalibrate)
   double g = 1+k[0]*r+k[1]*r*r ;
   double tx = 2*k[2]*p.x()*p.y()       +   k[3]*(r+2*p.x()*p.x()) ;
   double ty =   k[2]*(r+2*p.y()*p.y()) + 2*k[3]*p.x()*p.y() ;
-  Vector v_hat = Vector_(3, g*p.x() + tx, g*p.y() + ty, 1.0) ;
+  Vector v_hat = (Vector(3) << g*p.x() + tx, g*p.y() + ty, 1.0) ;
   Vector v_i = K.K() * v_hat ;
   Point2 p_i(v_i(0)/v_i(2), v_i(1)/v_i(2)) ;
   Point2 q = K.uncalibrate(p);
diff --git a/gtsam/geometry/tests/testPoint3.cpp b/gtsam/geometry/tests/testPoint3.cpp
index ef4a3894c..5b49c92b2 100644
--- a/gtsam/geometry/tests/testPoint3.cpp
+++ b/gtsam/geometry/tests/testPoint3.cpp
@@ -40,7 +40,7 @@ TEST(Point3, Lie) {
   EXPECT(assert_equal(eye(3), H2));
 
   EXPECT(assert_equal(Point3(5,7,9), p1.retract((Vector(3) << 4., 5., 6.))));
-  EXPECT(assert_equal(Vector_(3, 3.,3.,3.), p1.localCoordinates(p2)));
+  EXPECT(assert_equal((Vector)(Vector(3) << 3.,3.,3.), p1.localCoordinates(p2)));
 }
 
 /* ************************************************************************* */
diff --git a/gtsam/geometry/tests/testRot2.cpp b/gtsam/geometry/tests/testRot2.cpp
index 3de661224..c67031a13 100644
--- a/gtsam/geometry/tests/testRot2.cpp
+++ b/gtsam/geometry/tests/testRot2.cpp
@@ -96,7 +96,7 @@ TEST(Rot2, logmap)
 {
   Rot2 rot0(Rot2::fromAngle(M_PI/2.0));
   Rot2 rot(Rot2::fromAngle(M_PI));
-  Vector expected = Vector_(1, M_PI/2.0);
+  Vector expected = (Vector(1) << M_PI/2.0);
   Vector actual = rot0.localCoordinates(rot);
   CHECK(assert_equal(expected, actual));
 }
diff --git a/gtsam/geometry/tests/testSphere2.cpp b/gtsam/geometry/tests/testSphere2.cpp
index 61dadb8eb..965f5389d 100644
--- a/gtsam/geometry/tests/testSphere2.cpp
+++ b/gtsam/geometry/tests/testSphere2.cpp
@@ -212,9 +212,9 @@ inline static Vector randomVector(const Vector& minLimits,
 TEST(Sphere2, localCoordinates_retract) {
 
   size_t numIterations = 10000;
-  Vector minSphereLimit = Vector_(3, -1.0, -1.0, -1.0), maxSphereLimit =
-      Vector_(3, 1.0, 1.0, 1.0);
-  Vector minXiLimit = Vector_(2, -1.0, -1.0), maxXiLimit = Vector_(2, 1.0, 1.0);
+  Vector minSphereLimit = (Vector(3) << -1.0, -1.0, -1.0), maxSphereLimit =
+      (Vector(3) << 1.0, 1.0, 1.0);
+  Vector minXiLimit = (Vector(2) << -1.0, -1.0), maxXiLimit = (Vector(2) << 1.0, 1.0);
   for (size_t i = 0; i < numIterations; i++) {
 
     // Sleep for the random number generator (TODO?: Better create all of them first).
@@ -242,9 +242,9 @@ TEST(Sphere2, localCoordinates_retract) {
 TEST(Sphere2, localCoordinates_retract_expmap) {
   
   size_t numIterations = 10000;
-  Vector minSphereLimit = Vector_(3, -1.0, -1.0, -1.0), maxSphereLimit =
-      Vector_(3, 1.0, 1.0, 1.0);
-  Vector minXiLimit = Vector_(2, -M_PI, -M_PI), maxXiLimit = Vector_(2, M_PI, M_PI);
+  Vector minSphereLimit = (Vector(3) << -1.0, -1.0, -1.0), maxSphereLimit =
+      (Vector(3) << 1.0, 1.0, 1.0);
+  Vector minXiLimit = (Vector(2) << -M_PI, -M_PI), maxXiLimit = (Vector(2) << M_PI, M_PI);
   for (size_t i = 0; i < numIterations; i++) {
 
     // Sleep for the random number generator (TODO?: Better create all of them first).
diff --git a/gtsam/geometry/tests/timePose2.cpp b/gtsam/geometry/tests/timePose2.cpp
index a8f2f7137..28044068c 100644
--- a/gtsam/geometry/tests/timePose2.cpp
+++ b/gtsam/geometry/tests/timePose2.cpp
@@ -138,7 +138,7 @@ int main()
 
   // create a random pose:
   double x = 4.0, y = 2.0, r = 0.3;
-  Vector v = Vector_(3,x,y,r);
+  Vector v = (Vector(3) << x, y, r);
   Pose2 X = Pose2(3,2,0.4), X2 = X.retract(v), X3(5,6,0.3);
 
   TEST(Expmap, Pose2::Expmap(v));
diff --git a/gtsam/geometry/tests/timeRot2.cpp b/gtsam/geometry/tests/timeRot2.cpp
index 6d828ef9c..86dda2b8c 100644
--- a/gtsam/geometry/tests/timeRot2.cpp
+++ b/gtsam/geometry/tests/timeRot2.cpp
@@ -95,7 +95,7 @@ int main()
   // create a random direction:
   double norm=sqrt(16.0+4.0);
   double x=4.0/norm, y=2.0/norm;
-  Vector v = Vector_(2,x,y);
+  Vector v = (Vector(2) << x, y);
   Rot2 R = Rot2(0.4), R2 = R.retract(v), R3(0.6);
 
   TEST(Rot2_Expmap, Rot2::Expmap(v));
diff --git a/gtsam/geometry/tests/timeRot3.cpp b/gtsam/geometry/tests/timeRot3.cpp
index 64feb1909..dd98465ed 100644
--- a/gtsam/geometry/tests/timeRot3.cpp
+++ b/gtsam/geometry/tests/timeRot3.cpp
@@ -39,7 +39,7 @@ int main()
   // create a random direction:
   double norm=sqrt(1.0+16.0+4.0);
   double x=1.0/norm, y=4.0/norm, z=2.0/norm;
-  Vector v = Vector_(3,x,y,z);
+  Vector v = (Vector(3) << x, y, z);
   Rot3 R = Rot3::rodriguez(0.1, 0.4, 0.2), R2 = R.retract(v);
 
   TEST("Rodriguez formula given axis angle", Rot3::rodriguez(v,0.001))