diff --git a/gtsam/linear/Sampler.cpp b/gtsam/linear/Sampler.cpp
index 48972a953..5d58d165a 100644
--- a/gtsam/linear/Sampler.cpp
+++ b/gtsam/linear/Sampler.cpp
@@ -14,8 +14,7 @@
  * @author Alex Cunningham
  */
 
-#include <boost/random/normal_distribution.hpp>
-#include <boost/random/variate_generator.hpp>
+#include <boost/random.hpp>
 
 #include <gtsam/linear/Sampler.h>
 namespace gtsam {
@@ -51,7 +50,7 @@ Vector Sampler::sampleDiagonal(const Vector& sigmas) {
     } else {
       typedef boost::normal_distribution<double> Normal;
       Normal dist(0.0, sigma);
-      boost::variate_generator<boost::minstd_rand&, Normal> norm(generator_, dist);
+      boost::variate_generator<boost::mt19937_64&, Normal> norm(generator_, dist);
       result(i) = norm();
     }
   }
diff --git a/gtsam/linear/Sampler.h b/gtsam/linear/Sampler.h
index f2632308b..453658147 100644
--- a/gtsam/linear/Sampler.h
+++ b/gtsam/linear/Sampler.h
@@ -37,7 +37,7 @@ protected:
   noiseModel::Diagonal::shared_ptr model_;
 
   /** generator */
-  boost::minstd_rand generator_;
+  boost::mt19937_64 generator_;
 
 public:
   typedef boost::shared_ptr<Sampler> shared_ptr;