Thank you for all the help so far, now I'm trying to check for convergence of my results.
First, is it enough to just check how the truncation error changes with sweeps (at a certain MAXSTATE value) or I need to both check the iteration history of the ground state energy and the truncation error?
Second, checking the supplementary materials provided (Computer Physics Communications 185 (2014) 3430–3440) for the Hubbard ladder, why are there 381 values of energy and truncation weights for every iteration? Why is the total truncation error for that iteration is the sum of all the 381 truncation weight values (based on iteration.py)? How is truncation error defined (e.g. is it 1 - sum of all reduced density eigenvalues corresponding to all the picked states?) ?
As far as I understand the DMRG algorithm every iteration/sweep has a corresponding bond dimension in which as you progress through iterations (until you reach the "SWEEPS"th iteration), bond dimension increases until you reach "MAXSTATES" at the "SWEEPS"th iteration.
Regards, Robertson Esperanza