1998 IEEE. Personal use of this material is permitted. However, permission to reprint/republish this material for advertising or promotional purposes or for creating new collective works for resale or redistribution to servers or lists, or to reuse any copyrighted component of this work in other works must be obtained from the IEEE.

IEEE Transactions on Antennas and Propagation
Volume 46 Number 5, May 1998

Wavelet Packets for Fast Solution of Electromagnetic Integral Equations

Page 618.

Abstract:

This paper considers the problem of wavelet sparsification of matrices arising in numerical solution of electromagnetic integral equations by the method of moments. Scattering of plane waves from two-dimensional (2-D) cylinders is computed numerically using a constant number of test functions per wavelength. Discrete wavelet packet (DWP) similarity transformations and thresholding are applied to system matrices to obtain sparsity. If thresholds are selected to keep relative residual error constant the matrix sparsity is of order O(N^p) with p < 2. This stands in contrast with O(N^2) sparsities obtained with standard wavelet transformations [1]. Numerical tests also show that the DWP method yields faster matrix-vector multiplication than some fast multipole algorithms.

References

1. R. L. Wagner and W. C. Chew, "A study of wavelets for the solution of electromagnetic integral equations," IEEE Trans. Antennas Propagat., vol. 43, pp. 802-810, Aug. 1995.
2. V. Rokhlin, "Rapid solution of integral equations of scattering theory in two dimensions," J. Comp. Phys., vol. 86, no. 2, pp. 414-439, 1990.
3. R. L. Wagner and W. C. Chew, "A ray-propagation fast multipole algorithm," Microwave Opt. Technol. Lett., vol. 7, no. 10, pp. 435-438, 1994.
4. E. Bleszynski, M. Bleszynski, and T. Jaroszewicz, "Fast integral-equation solver for electromagnetic scattering problems," in IEEE AP-S Int. Symp. Dig., Seattle, WA, June 1994, pp. 416-419.
5. F. X. Canning, "Sparse approximations for solving integral equations with oscillatory kernels," SIAM J. Sci. Statist. Comput., vol. 13, no. 1, pp. 71-87, 1992.
6. D. S. Wang and G. Welland, "Modeling of electromagnetic scattering using wavelet techniques for geometry modeling and expansion functions," in URSI Radio Sci. Meet. Dig., Detroit, MI, June 1992, p. 277.
7. H. Kim and H. Ling, "On the application of fast wavelet transform to integral equation solution of electromagnetic scattering problems," Microwave Opt. Technol. Lett., vol. 6, no. 3, pp. 168-173, 1993.
8. G. Wang, "A hybrid wavelet expansion and boundary element analysis of electromagnetic scattering from conducting objects," IEEE Trans. Antennas Propagat., vol. 43, pp. 170-178, Feb. 1995.
9. Z. Baharav and Y. Leviatan, "Impedance matrix compression using adaptively constructed basis functions," IEEE Trans. Antennas Propagat., vol. 44, pp. 1231-1238, Sept. 1996.
10. G. Beylkin, R. Coifman, and V. Rokhlin, "Fast wavelet transforms and numerical algorithms i," Communicat. Pure Appl. Math., vol. XLIV, pp. 141-183, 1991.
11. M. V. Wickerhauser, Adapted Wavelet Analysis from Theory to Software.Boston, MA: A. K. Peters, 1994.
12. I. Daubechies, "Ten lectures on wavelets," CBMS Lecture Notes, SIAM, Philadelphiak PA, vol. 61, 1992.
13. A. Cohen, I. Daubechies, and P. Vial, "Wavelets on the interval and fast wavelet transforms," Appl. Computat. Harmonic Anal., vol. 1, no. 1, pp. 54-81, 1993.