@@ -27,6 +27,9 @@ public:

     template <class Parser>

     RowMatrix(Parser &p, unsigned nrow, unsigned ncol);

+    template <class Parser>

+    RowMatrix(Parser &p, unsigned nrow, std::vector<unsigned> whichCols);

+

     unsigned nRow() const {return mNumRows;}

     unsigned nCol() const {return mNumCols;}

@@ -86,7 +89,6 @@ public:

     friend Archive& operator>>(Archive &ar, ColMatrix &mat);

 };

-

 // construct RowMatrix from file

 template <class Parser>

 RowMatrix::RowMatrix(Parser &p, unsigned nrow, unsigned ncol) : mNumRows(nrow),

@@ -106,6 +108,16 @@ mNumCols(ncol)

     }

 }

+// This should construct a matrix, only reading the columns in "whichCols"

+// from the file. The matrix should have "nrow" rows and "whichCols.size()"

+// columns.

+template <class Parser>

+RowMatrix(Parser &p, unsigned nrow, std::vector<unsigned> whichCols)

+{

+

+

+}

+

 // construct ColMatrix from file

 template <class Parser>

 ColMatrix::ColMatrix(Parser &p, unsigned nrow, unsigned ncol) : mNumRows(nrow),