[x, y] = np.meshgrid(xp, xp)