对一个巨大的多维数组进行n维FFT的资源

Question

我有一个50维数组，它的尺寸是255 x 255 x 255 x…（50倍）…x255。所以总共有50^255个浮点数。甚至连装在公羊身上的想法都超出了范围。此外，我需要对这个数组进行50维快速傅立叶变换（DFT）。我不能用python在普通PC上实现，我甚至无法想象在GPU上做。所以我想我得借助硬盘存储器，但那也太大了。我不需要这个实时的，我甚至可以花上几天的时间来运行它。我不知道我需要什么样的机器，或者说它有可能吗？感谢你的建议。超级电脑，电网，或者其他东西，即使太贵了，我也不担心投资。在

Answer 1

如果您found enough universes to save your data in，以下是您可以做的：

傅里叶变换是可分离的，这意味着一个接一个地计算每个轴的DFT将得到与计算n维DFT相同的结果：

for i in range(C.ndim):
    C[...] = numpy.fft.fft(C, axis=i)

使用2D张量再次检查值是否正确（因为我们有2D FFT numpy.fft.fft2进行比较）：

^{pr2}$