|
地球物理学进展 2010
A fast and accurate algorithm for integral equations in a 2D homogeneous medium
|
Abstract:
The integral equations in a 2D homogeneous medium are accurately calculated by the stabilized biconjugate-gradient fast Fourier transform (BCGS-FFT) algorithm. The new interpolating function is chosen as basis and testing functions to get the weak-form discretization of the integral equations. The discrete form of the integral equations is solved via the stabilized biconjugate-gradient iteration method, and the distribution of the electric field within the abnormal objects can be obtained. The product between the Green′s function and the electric field within the integral equations can be expressed in the form of convolution, which can be accelerated by fast Fourier transform. Numerical examples have shown the accuracy and efficiency of the algorithm.