Gfhathat = np.real(fft(Gfhat[:n + 1, :], 2 * (n + 1), axis=1))