SUPERCOMPUTER AWARE APPROACH FOR THE SOLUTION OF CHALLENGING ELECTROMAGNETIC PROBLEMS

Abstract|It is a proven fact that The Fast Fourier Transform(FFT) extension of the conventional Fast Multipole Method (FMM)reduces the matrix vector product (MVP) complexity and preservesthe propensity for parallel scaling of the single level FMM. In thispaper, an e–cient parallel strategy of a nested variation of the FMM-FFT algorithm that reduces the memory requirements is presented.The solution provided by this parallel implementation for a challengingproblem with more than 0.5 billion unknowns has constituted the worldrecord in computational electromagnetics (CEM) at the beginning of2009.1. INTRODUCTIONRecent years have seen an increasing efiort in the development of fastand e–cient electromagnetic solutions with a reduced computationalcost regarding the conventional Method of Moments. Among others,the Fast Multipole Method (FMM) [1] and its multilevel version, theMLFMA [2,3] have constituted one of the most important advances inthat context.This development of fast electromagnetic solvers has gone handin hand with the constant advances in computer technology. Dueto this simultaneous growth, overcoming the limits in the scalabilityof the available codes became a priority in order to take advantageof the large amount of computational resources and capabilities thatare available in modern High Performance Computer (HPC) systems.For this reason, works focused on the parallelization improvement ofthe Multilevel Fast Multipole Algorithm (MLFMA) [4{13] have gainedinterest in last years.Besides, the FMM-Fast Fourier Transform (FMM-FFT) deservesbe taken into account as an alternative to beneflt from massivelyparallel distributed computers. This variation of the single-level FMMwas flrst proposed in [14] as an acceleration technique applied to almostplanar surfaces. Later on, a parallelized implementation was applied togeneral three-dimensional geometries [15]. The method uses the FFTto speedup the translation stage resulting in a dramatic reduction ofthe matrix-vector product (MVP) time requirement with respect tothe FMM. Although in general the FMM-FFT is not algorithmically ase–cient as the MLFMA, it has the advantage of preserving the naturalparallel scaling propensity of the single-level FMM in the spectral (