Truncated multigrid versus pre-corrected FFT/AIM for bioelectromagnetics: When is O(N) better than O(NlogN)?

The effectiveness of multigrid and fast Fourier transform (FFT) based methods are investigated for accelerating the solution of volume integral equations encountered in bioelectromagnetics (BIOEM) analysis. The typical BIOEM simulation is in the mixed-frequency regime of analysis because the field variations in the simulation domain are dictated by a combination of the free space wavelength, geometrical features, and the wavelengths/skin depths in tissues. In this case, multigrid-based methods (when appropriately truncated at high-frequency levels) can achieve O(N) complexity that is asymptotically superior to the O(NlogN) complexity of FFT-based ones. Nevertheless, the constant in front of their asymptotic complexity estimate is larger and their accuracy-efficiency tradeoffs are different. Numerical experiments are performed to compare these methods and the results show that multigrid-based methods begin to outperform FFT-based ones for N∼103.

[1]  C Gabriel,et al.  The dielectric properties of biological tissues: I. Literature survey. , 1996, Physics in medicine and biology.

[2]  Niels Kuster,et al.  The Virtual Family—development of surface-based anatomical models of two adults and two children for dosimetric simulations , 2010, Physics in medicine and biology.

[3]  Kai Yang,et al.  Erratum: Comparison of precorrected FFT/adaptive integral method matching schemes (Microwave and Optical Technology Letters (2011) 53 (1368-1372) DOI: 10.1002/mop.26230) , 2011 .

[4]  A. Boag,et al.  Multilevel nonuniform grid algorithm for acceleration of integral equation-based solvers for acoustic scattering , 2010, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[5]  Ali E. Yilmaz,et al.  A 2-D decomposition based parallelization of AIM for 3-D BIOEM problems , 2011, 2011 IEEE International Symposium on Antennas and Propagation (APSURSI).

[6]  T. Eibert,et al.  A multilevel interpolating fast integral solver with fast fourier transform acceleration , 2010, 2010 URSI International Symposium on Electromagnetic Theory.

[7]  P. M. Berg,et al.  The three dimensional weak form of the conjugate gradient FFT method for solving scattering problems , 1992 .

[8]  D. Wilton,et al.  A tetrahedral modeling method for electromagnetic scattering by arbitrarily shaped inhomogeneous dielectric bodies , 1984 .

[9]  Guneet Kaur,et al.  A Multiple-Grid Adaptive Integral Method for Multi-Region Problems , 2010, IEEE Transactions on Antennas and Propagation.

[10]  Jacob K. White,et al.  A precorrected-FFT method for electrostatic analysis of complicated 3-D structures , 1997, IEEE Trans. Comput. Aided Des. Integr. Circuits Syst..

[11]  Chen Zhang,et al.  Chinese adult anatomical models and the application in evaluation of RF exposures , 2011, Physics in medicine and biology.

[12]  Leung Tsang,et al.  A new multilevel Green's function interpolation method for large-scale low-frequency EM simulations , 2005, IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems.

[13]  A. Brandt.  Multilevel computations of integral transforms and particle interactions with oscillatory kernels , 1991 .

[14]  M. Bleszynski,et al.  AIM: Adaptive integral method for solving large‐scale electromagnetic scattering and radiation problems , 1996 .