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