We present a parallel FFT algorithm for SIMD systems following the `Transpose Algorithm' approach. The method is based on the assignment of the data field onto a 1-dimensional ring of systolic cells. The systolic array can be universally mapped onto any parallel system. In particular for systems with next-neighbour connectivity our method has the potential to improve the efficiency of matrix transposition by use of hyper-systolic communication. We have realized a scalable parallel FFT on the APE100/Quadrics massively parallel computer, where our implementation is part of a 2-dimensional hydrodynamics code for turbulence studies. A possible generalization to 4-dimensional FFT is presented, having in mind QCD applications.
@article{arxiv.hep-lat/9710060,
title = {FFT for the APE Parallel Computer},
author = {Thomas Lippert and Klaus Schilling and Federico Toschi and Sven Trentmann and Raffaele Tripiccione},
journal= {arXiv preprint arXiv:hep-lat/9710060},
year = {2015}
}