zout[-size:, :size] = band - np.abs(zout[-size:, 2 * size:size:-1] - band)