Multilevel fast multipole algorithm enhanced by GPU parallel technique for electromagnetic scattering problems

Along with the development of graphics processing Units (GPUS) in floating point operations and programmability, GPU has increasingly become an attractive alternative to the central processing unit (CPU) for some of compute-intensive and parallel tasks.In this article, the multilevel fast multipole algorithm (MLFMA) combined with graphics hardware acceleration technique is applied to analyze electromagnetic scattering from complex target. Although it is possible to perform scattering simulation of electrically large targets on a personal computer (PC) through the MLFMA, a large CPU time is required for the execution of aggregation, translation, and deaggregation operations. Thus GPU computing technique is used for the parallel processing of MLFMA and a significant speedup of matrix vector product (MVP) can be observed. Following the programming model of compute unified device architecture (CUDA), several kernel functions characterized by the single instruction multiple data (SIMD) mode are abstracted from components of the MLFMA and executed by multiple processors of the GPU. Numerical results demonstrate the efficiency of GPU accelerating technique for the MLFMA. © 2010 Wiley Periodicals, Inc. Microwave Opt Technol Lett 52: 502–507, 2010; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/mop.24963

[1]  M.J. Inman,et al.  Programming video cards for computational electromagnetics applications , 2005, IEEE Antennas and Propagation Magazine.

[2]  Ramani Duraiswami,et al.  Fast multipole methods on graphics processors , 2008, J. Comput. Phys..

[3]  M.M. Okoniewski,et al.  Acceleration of finite-difference time-domain (FDTD) using graphics processor units (GPU) , 2004, 2004 IEEE MTT-S International Microwave Symposium Digest (IEEE Cat. No.04CH37535).

[4]  S. Velamparambil,et al.  Analysis and performance of a distributed memory multilevel fast multipole algorithm , 2005, IEEE Transactions on Antennas and Propagation.

[5]  V. Rokhlin.  Rapid Solution of Integral Equations of Scattering Theory , 1990 .

[6]  Martin Nilsson,et al.  Iterative solution of Maxwell's equations in frequency domain , 2002 .

[7]  Z. Nie,et al.  Acceleration of the Method of Moments Calculations by Using Graphics Processing Units , 2008, IEEE Transactions on Antennas and Propagation.

[8]  Jiming Song,et al.  Multilevel fast‐multipole algorithm for solving combined field integral equations of electromagnetic scattering , 1995 .

[9]  John D. Owens,et al.  GPU Computing , 2008, Proceedings of the IEEE.

[10]  Laurence E. Turner,et al.  Graphics processor unit (GPU) acceleration of finite-difference time-domain (FDTD) algorithm , 2004, 2004 IEEE International Symposium on Circuits and Systems (IEEE Cat. No.04CH37512).

[11]  Weng Cho Chew,et al.  Fast algorithm for solving hybrid integral equations , 1993 .