diff --git a/gtsam/linear/HessianFactor.cpp b/gtsam/linear/HessianFactor.cpp
index 459c12f76..2614aea3f 100644
--- a/gtsam/linear/HessianFactor.cpp
+++ b/gtsam/linear/HessianFactor.cpp
@@ -503,8 +503,6 @@ GaussianFactor::shared_ptr HessianFactor::negate() const
 void HessianFactor::multiplyHessianAdd(double alpha, const VectorValues& x,
     VectorValues& y) {
 
-  std::cout << "size() " << size() << std::endl;
-
   for (size_t i = 0; i < size(); ++i) {
     pair<VectorValues::iterator, bool> it = y.tryInsert(keys_[i], Vector());
     Vector& yi = it.first->second;