data = np.full(pvec.shape, np.nan)