Xu X M, Liu Q H
Department of Electrical and Computer Engineering, Duke University, Durham, NC 27708-0291, USA.
IEEE Trans Ultrason Ferroelectr Freq Control. 2001 Mar;48(2):522-9. doi: 10.1109/58.911735.
This paper presents the application of the conjugate-gradient (CG) fast Fourier transform (FFT) (CG-FFT) method and the CG nonuniform FFT (CG-NUFFT) method for the integral equation arising from acoustic scattering problems. In the conventional method of moments (MoM) for integral equations, the CPU and memory requirements are O(N3) and O(N2), respectively, where N is the number of unknowns in the problem. The CG-FFT method, which combines the iterative conjugate-gradient method with FFT, reduces these requirements to O(KN log2 N) and O(N), respectively, where K is the number of CG iterations. The CG-NUFFT method differs from the CG-FFT method in that it makes use of nonuniform FFT algorithms instead of FFT to allow a nonuniform discretization. Therefore, the CG-NUFFT method can solve the integral equation with both uniform and nonuniform grid while retaining the efficiency of the CG-FFT method. These two methods are applied to solve for two-dimensional constant density acoustic scattering problems. Numerical results demonstrate that they can solve much larger problems than the MoM.
