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