leftvecs = [np.array(1, ndmin=3)] + [None] * (nr_sites - var_sites)