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