27 #ifndef MAQUIS_DMRG_IDENTITY_MPS_H
28 #define MAQUIS_DMRG_IDENTITY_MPS_H
32 template <
class Matrix,
class InSymm>
37 t.
set(0, 0, identity_matrix<Matrix>(phys_psi));
void set(index_type li, index_type ri, op_t const &op, value_type scale_=1.0)
MPS< Matrix, typename grouped_symmetry< InSymm >::type > mpo_to_smps_group(MPO< Matrix, InSymm > const &mpo, Index< InSymm > const &phys_i, std::vector< Index< typename grouped_symmetry< InSymm >::type > > const &allowed)
MPS< Matrix, typename grouped_symmetry< InSymm >::type > identity_dm_mps(std::size_t L, Index< InSymm > const &phys_psi, std::vector< Index< typename grouped_symmetry< InSymm >::type > > const &allowed)