diff --git a/gtsam/linear/GaussianFactorGraph.cpp b/gtsam/linear/GaussianFactorGraph.cpp
index 54e721cd7..ad43b675b 100644
--- a/gtsam/linear/GaussianFactorGraph.cpp
+++ b/gtsam/linear/GaussianFactorGraph.cpp
@@ -355,14 +355,13 @@ namespace gtsam {
      f->multiplyHessianAdd(alpha, x, y);
   }
 
-  /* ************************************************************************* */
-  void GaussianFactorGraph::multiplyHessianAdd(double alpha,
-      const double* x, double* y) const {
-  vector<size_t> FactorKeys = getkeydim();
-  BOOST_FOREACH(const GaussianFactor::shared_ptr& f, *this)
-      f->multiplyHessianAdd(alpha, x, y, FactorKeys);
-
-  }
+  ///* ************************************************************************* */
+  //void GaussianFactorGraph::multiplyHessianAdd(double alpha,
+  //    const double* x, double* y) const {
+  //vector<size_t> FactorKeys = getkeydim();
+  //BOOST_FOREACH(const GaussianFactor::shared_ptr& f, *this)
+  //    f->multiplyHessianAdd(alpha, x, y, FactorKeys);
+  //}
 
   /* ************************************************************************* */
   void GaussianFactorGraph::multiplyInPlace(const VectorValues& x, Errors& e) const {
diff --git a/gtsam/linear/GaussianFactorGraph.h b/gtsam/linear/GaussianFactorGraph.h
index 910b25d1e..1861e2607 100644
--- a/gtsam/linear/GaussianFactorGraph.h
+++ b/gtsam/linear/GaussianFactorGraph.h
@@ -311,8 +311,8 @@ namespace gtsam {
         VectorValues& y) const;
 
     ///** y += alpha*A'A*x */
-    void multiplyHessianAdd(double alpha, const double* x,
-        double* y) const;
+    //void multiplyHessianAdd(double alpha, const double* x,
+    //    double* y) const;
 
     ///** In-place version e <- A*x that overwrites e. */
     void multiplyInPlace(const VectorValues& x, Errors& e) const;