![Sample our Physical Sciences journals, sign in here to start your FREE access for 14 days]({"type":"image" , "src" : "/sda/110819/TOC-Subject-Sample-2021-Physical-Sciences.jpg"})
Abstract
The method of moments for the simplified scattering and radiation on thin wire models can lead to the large-scale linear systems with non-Hermitian Toeplitz-like coefficient matrices. Since such reduced linear systems are typically dense and ill conditioned, iterative solvers are rarely considered to deal with the discretized linear systems in these electromagnetic model simulations. However, utilizing the Toeplitz-like structure of coefficient matrices, we can use the fast Toeplitz matrix-vector product (MVP) to refresh iterative solvers for the discretized linear systems. Compared with the traditional solvers, the fast numerical method can reduce the computational cost from to
and the storage from
to
without using any lossy compression. Due to the Toeplitz-like structure of coefficient matrices, circulant preconditioners are considered to further accelerate the convergence of iterative solvers. Since the circulant matrices can be diagonalized via fast Fourier transforms (FFTs), the fast Toeplitz MVP is utilized to retain the total computational complexity of the proposed methods at
, where
is the number of spatial grid nodes. Theoretical analyses of the spectra of the preconditioned linear systems are also preliminarily investigated. Finally, extensive numerical experiments involving both scattering and radiation issues are employed to demonstrate the effectiveness of the proposed methods in aspects of the number of iterations and CPU time elapsed.
Acknowledgements
The authors thank Dr Bruno Carpentieri and Dr Ran Zhao for their valuable comments and suggestions, leading to a clearer exhibition of the present numerical experiments.
Notes
No potential conflict of interest was reported by the authors.
1 Recently, it has been found that the CORS method is more comparative than other popular iterative solvers for dealing with the (dense) non-Hermitian linear systems arising from the computational electromagnetic problems (e.g. refer to [Citation43–Citation45]), so here we consider to exploit the CORS method for solving the resulted systems of linear equations in this paper.