diff --git a/cpp/FactorGraph-inl.h b/cpp/FactorGraph-inl.h
index ae1af2d85..3cf6e7196 100644
--- a/cpp/FactorGraph-inl.h
+++ b/cpp/FactorGraph-inl.h
@@ -186,7 +186,8 @@ boost::shared_ptr<Ordering> colamd(int n_col, int n_row, int nrNonZeros, const m
 	// Convert elimination ordering in p to an ordering
 	boost::shared_ptr<Ordering> result(new Ordering);
 	for(int j = 0; j < n_col; j++)
-		result->push_back(initialOrder[j]);
+	  result->push_back(initialOrder[p[j]]);
+		//result->push_back(initialOrder[j]);
 	delete [] p; // delete colamd result vector
 
 	return result;