A Preconditioned COCG Method for Solving Complex Symmetric Linear Systems Arising from Scattering Problems

In this paper, we study a class of modified incomplete LDLT preconditioners followed by the conjugate orthogonal conjugate gradient (COCG) method for solving the large complex symmetric system with linear equations Ax = b, arising from the scattering problems. The preconditioner is derived from the incomplete Cholesky factorization with diagonal compensation. With an ordering scheme, our factorized preconditioners can reduce the memory requirement as well as the COCG iteration numbers. Numerical tests on harmonic analysis for plane wave scattering from a metallic plate and a metallic sphere coated by a lossy dielectric layer show the efficiency of our method.

[1]  T. Sarkar,et al.  ANALYSIS OF TRANSIENT ELECTROMAGNETIC SCATTERING WITH PLANE WAVE INCIDENCE USING MOD-FDM , 2007 .

[2]  J. Meijerink,et al.  An iterative solution method for linear systems of which the coefficient matrix is a symmetric -matrix , 1977 .

[3]  Yousef Saad,et al.  ILUT: A dual threshold incomplete LU factorization , 1994, Numer. Linear Algebra Appl..

[4]  H. Oraizi,et al.  COMBINATION OF MLS, GA & CG FOR THE REDUCTION OF RCS OF MULTILAYERED CYLINDRICAL STRUCTURES COMPOSED OF DISPERSIVE METAMATERIALS , 2008 .

[5]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[6]  Z. Nie,et al.  Application of diagonally perturbed incomplete factorization preconditioned conjugate gradient algorithms for edge finite-element analysis of Helmholtz equations , 2006, IEEE Transactions on Antennas and Propagation.

[7]  RCS computation of airplane using parabolic equation , 2005, 2005 Asia-Pacific Microwave Conference Proceedings.

[8]  T. Manteuffel.  An incomplete factorization technique for positive definite linear systems , 1980 .

[9]  M. Hestenes,et al.  Methods of conjugate gradients for solving linear systems , 1952 .

[10]  J. Nédélec.  Mixed finite elements in ℝ3 , 1980 .

[11]  Roland W. Freund,et al.  Conjugate Gradient-Type Methods for Linear Systems with Complex Symmetric Coefficient Matrices , 1992, SIAM J. Sci. Comput..

[12]  Patrick R. Amestoy,et al.  An Approximate Minimum Degree Ordering Algorithm , 1996, SIAM J. Matrix Anal. Appl..

[13]  Electromagnetic Scattering by a Set of Objects: an Integral Method Based on Scattering Operator , 2006 .

[14]  B. Carpentieri.  Fast iterative solution methods in electromagnetic scattering , 2008 .

[15]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[16]  H. V. D. Vorst,et al.  A Petrov-Galerkin type method for solving Axk=b, where A is symmetric complex , 1990 .

[17]  Iain S. Duff,et al.  Users' guide for the Harwell-Boeing sparse matrix collection (Release 1) , 1992 .

[18]  Mohammed F. Hadi,et al.  Integral PML Absorbing Boundary Conditions for the High-Order M24 FDTD Algorithm , 2007 .

[19]  Carretera de Valencia,et al.  The finite element method in electromagnetics , 2000 .

[20]  Henk A. van der Vorst,et al.  Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..

[21]  J. Pasciak,et al.  Computer solution of large sparse positive definite systems , 1982 .

[22]  A. Bunse-Gerstner,et al.  On a conjugate gradient-type method for solving complex symmetric linear systems , 1999 .

[23]  X. Hou,et al.  On the Condition Number of Matrices from Various Hybrid Vector FEM-BEM Formulations for 3-D Scattering , 2006 .

[24]  R. Freund,et al.  QMR: a quasi-minimal residual method for non-Hermitian linear systems , 1991 .

[25]  Y.-L. Pan.  FINITE ELEMENT METHOD SIMULATION OF PHOTOINDUCTIVE IMAGING FOR CRACKS , 2008 .

[26]  Chih-Wei Huang,et al.  Radar target recognition by projected features of frequency-diversity RCS , 2008 .

[27]  Andreas C. Cangellaris,et al.  GT-PML: generalized theory of perfectly matched layers and its application to the reflectionless truncation of finite-difference time-domain grids , 1996, IMS 1996.

[28]  Andreas C. Cangellaris,et al.  GT-PML: generalized theory of perfectly matched layers and its application to the reflectionless truncation of finite-difference time-domain grids , 1996, 1996 IEEE MTT-S International Microwave Symposium Digest.

[29]  Xing Wang,et al.  SCATTERING ANALYSIS OF A PRINTED DIPOLE ANTENNA USING PBG STRUCTURES , 2008 .

[30]  Homayoon Oraizi,et al.  ULTRA WIDE BAND RCS OPTIMIZATION OF MULTILAYERD CYLINDRICAL STRUCTURES FOR ARBITRARILY POLARIZED INCIDENT PLANE WAVES , 2008 .

[31]  M. Benzi.  Preconditioning techniques for large linear systems: a survey , 2002 .

[32]  Weng Cho Chew,et al.  Applying divergence-free condition in solving the volume integral equation , 2006 .

[33]  Khalid Fawzi Ahmed Hussein.  Fast Computational Algoritm for EFIE Applied to Arbitrarily-Shaped Conducting Surfaces , 2007 .

[34]  J. Volakis,et al.  Finite element method for electromagnetics : antennas, microwave circuits, and scattering applications , 1998 .

[35]  Owe Axelsson,et al.  Diagonally compensated reduction and related preconditioning methods , 1994, Numer. Linear Algebra Appl..

[36]  P. L. Rui,et al.  Implicitly Restarted Gmres Fast Fourier Transform Method for Electromagnetic Scattering , 2007 .