` psi = PurificationMPS.from_infiniteT(M.lat.mps_sites(), bc=bc)`

, but it gives me a grand canonical ensemble result. I want to conserve some charges, for example, I want to conserve the total number of boson. But I dont know How to do it. I use Code: Select all

`psi.gauge_total_charge([10])`

Code: Select all

`psi.get_total_charge(only_physical_legs = True)`

Code: Select all

`psi.get_total_charge()`

Code: Select all

`for i in range(L): print(psi._B[i].qtotal)`