      }
      
  private:
-@@ -546,7 +546,7 @@
+@@ -510,7 +511,7 @@
+ class basis_mps_init_generic : public mps_initializer<Matrix, SymmGroup>
+ {
+ public:
+-    typedef std::vector<boost::tuple<typename SymmGroup::charge, size_t> > state_type;
++    typedef std::vector<boost::tuple<typename SymmGroup::charge, std::size_t> > state_type;
+     basis_mps_init_generic(BaseParameters & params,
+                            std::vector<Index<SymmGroup> > const& phys_dims_,
+                            typename SymmGroup::charge right_end_,
+@@ -546,7 +547,7 @@
          if (mps.length() != state.size())
              throw std::runtime_error("Size of the initial state does not match the MPS size.");
          

Index: applications/dmrg/mps/framework/dmrg/mp_tensors/mps_initializers.h
===================================================================
--- applications/dmrg/mps/framework/dmrg/mp_tensors/mps_initializers.h	(revision 7814)
+++ applications/dmrg/mps/framework/dmrg/mp_tensors/mps_initializers.h	(revision 7819)
@@ -462,7 +462,7 @@
         std::vector<boost::tuple<charge, size_t> > state(mps.length());
         for (int i=0; i<mps.length(); ++i)
             state[i] = boost::make_tuple(C, occupation[i]);
-        mps = state_mps<Matrix>(state, phys_dims, site_type);

         if (mps[mps.length()-1].col_dim()[0].first != right_end)
             throw std::runtime_error("Initial state does not satisfy total quantum numbers.");
     }
Index: applications/dmrg/mps/framework/dmrg/mp_tensors/state_mps.h
===================================================================
--- applications/dmrg/mps/framework/dmrg/mp_tensors/state_mps.h	(revision 7814)
+++ applications/dmrg/mps/framework/dmrg/mp_tensors/state_mps.h	(working copy)
@@ -32,7 +32,7 @@
 #include <boost/tuple/tuple.hpp>
 
 template <class Matrix, class SymmGroup>
-MPS<Matrix, SymmGroup> state_mps(std::vector<boost::tuple<typename SymmGroup::charge, size_t> > const & state,
+MPS<Matrix, SymmGroup> state_mps(std::vector<boost::tuple<typename SymmGroup::charge, std::size_t> > const & state,
                                  std::vector<Index<SymmGroup> > const& phys_dims, std::vector<int> const& site_type)
 {
     typedef typename SymmGroup::charge charge;
