How to calculate a two-dimensional FFT of a large matrix with CUDA?

Question

Now I have a large 16K*16K matrix,and the global memory is not enough.How to calculate the two-dimensional FFT of the matrix?