if pos < nr_sites - var_sites:
                # We always do this, because we don't do the last site again.
                eigvec.canonicalize(right=pos_end)
                leftvecs[pos + 1] = None
                rightvecs[pos] = _eig_rightvec_add(