yp = np.zeros(y.shape, np.float_)