diff --git a/cpp/GaussianFactor.h b/cpp/GaussianFactor.h
index a530c3f1a..05aa803e3 100644
--- a/cpp/GaussianFactor.h
+++ b/cpp/GaussianFactor.h
@@ -12,6 +12,7 @@
 #include <boost/shared_ptr.hpp>
 #include <boost/tuple/tuple.hpp>
 #include <boost/serialization/map.hpp>
+#include <boost/foreach.hpp>
 #include <list>
 #include <set>
 
@@ -90,6 +91,14 @@ public:
 	    As_.insert(terms[i]);
 	}
 
+	GaussianFactor(const std::list<std::pair<Symbol, Matrix> > &terms,
+	    const Vector &b, const SharedDiagonal& model) :
+		model_(model), b_(b)  {
+		std::pair<Symbol, Matrix> pair;
+		BOOST_FOREACH(pair, terms)
+			As_.insert(pair);
+	}
+
 	/** Construct from Conditional Gaussian */
 	GaussianFactor(const boost::shared_ptr<GaussianConditional>& cg);