zout[:size, size:-size] = band - np.abs(z[size:0:-1, :] - band)