DocumentCode :
2959034
Title :
High Performance Non-uniform FFT on Modern X86-based Multi-core Systems
Author :
Kalamkar, Dhiraj D. ; Trzasko, Joshua D. ; Sridharan, Srinivas ; Smelyanskiy, Mikhail ; Kim, Daehyun ; Manduca, Armando ; Shu, Yunhong ; Bernstein, Matt A. ; Kaul, Bharat ; Dubey, Pradeep
Author_Institution :
Parallel Comput. Lab., Intel Labs., Bangalore, India
fYear :
2012
fDate :
21-25 May 2012
Firstpage :
449
Lastpage :
460
Abstract :
The Non-Uniform Fast Fourier Transform (NUFFT) is a generalization of FFT to non-equidistant samples. It has many applications which vary from medical imaging to radio astronomy to the numerical solution of partial differential equations. Despite recent advances in speeding up NUFFT on various platforms, its practical applications are still limited, due to its high computational cost, which is significantly dominated by the convolution of a signal between a non-uniform and uniform grids. The computational cost of the NUFFT is particularly detrimental in cases which require fast reconstruction times, such as iterative 3D non-Cartesian MRI reconstruction. We propose novel and highly scalable parallel algorithm for performing NUFFT on x86-based multi-core CPUs. The high performance of our algorithm relies on good SIMD utilization and high parallel efficiency. On convolution, we demonstrate on average 90% SIMD efficiency using SSE, as well up to linear scalability using a quad-socket 40-core Intel(R) Xeon(R) E7-4870 Processors based system. As a result, on dual socket Intel(R) Xeon(R) X5670 based server, our NUFFT implementation is more than 4x faster compared to the best available NUFFT3D implementation, when run on the same hardware. On Intel(R) Xeon(R) E5-2670 processor based server, our NUFFT implementation is 1.5X faster than any published NUFFT implementation today. Such speed improvement opens new usages for NUFFT. For example, iterative multi channel reconstruction of a 240×240×240 image could execute in just over 3 minutes, which is on the same order as contemporary non-iterative (and thus less-accurate) 3D NUFFT-based MRI reconstructions.
Keywords :
fast Fourier transforms; iterative methods; microprocessor chips; multiprocessing systems; parallel algorithms; performance evaluation; Intel(R) Xeon(R) E5-2670 processor based server; NUFFT; SIMD utilization; X86-based multicore systems; dual socket Intel(R) Xeon(R) X5670 based server; fast Fourier transform; high performance nonuniform FFT; iterative 3D nonCartesian MRI reconstruction; iterative multichannel reconstruction; medical imaging; nonequidistant samples; parallel algorithm; partial differential equations; quad-socket 40-core Intel(R) Xeon(R) E7-4870 processors based system; radio astronomy; Convolution; Hardware; Image reconstruction; Interpolation; Kernel; Magnetic resonance imaging; Transforms; Non-uniform FFT; Parallelization; Scalability; Vectorization;
fLanguage :
English
Publisher :
ieee
Conference_Titel :
Parallel & Distributed Processing Symposium (IPDPS), 2012 IEEE 26th International
Conference_Location :
Shanghai
ISSN :
1530-2075
Print_ISBN :
978-1-4673-0975-2
Type :
conf
DOI :
10.1109/IPDPS.2012.49
Filename :
6267881
Link To Document :
بازگشت