M = 100
N = 101
x = (np.arange(1, M*N + 1).reshape(M, N, order='F')).astype('float64')
h = (np.arange(1, M*N + 1).reshape(M, N, order='F')).astype('float64') + 1j

h = ifftshift(np.conj(h), axes=1)
out = np.zeros(N, dtype='complex128')
for i in range(N):
    out[i] = np.sum(x * np.roll(h, i, axis=1))