A Butterfly-Based Direct Integral-Equation Solver Using Hierarchical LU Factorization for Analyzing Scattering From Electrically Large Conducting Objects

A butterfly-based direct combined-field integral-equation (CFIE) solver for analyzing scattering from electrically large, perfect electrically conducting objects is presented. The proposed solver leverages the butterfly scheme to compress blocks of the hierarchical LU-factorized discretized CFIE operator and uses randomized butterfly reconstruction schemes to expedite the factorization. The memory requirement and computational cost of the direct butterfly-CFIE solver scale as <inline-formula> <tex-math notation="LaTeX">$O(N\log ^{2}N)$ </tex-math></inline-formula> and <inline-formula> <tex-math notation="LaTeX">$O(N^{1.5}\log N)$ </tex-math></inline-formula>, respectively. These scaling estimates permit significant memory and CPU savings when compared to those realized by low-rank decomposition-based solvers. The efficacy and accuracy of the proposed solver are demonstrated through its application to the analysis of scattering from canonical and realistic objects involving up to 14 million unknowns.

[1]  D. Wilton,et al.  Electromagnetic scattering by surfaces of arbitrary shape , 1980 .

[2]  Mark Tygert,et al.  Fast algorithms for spherical harmonic expansions, III , 2009, J. Comput. Phys..

[3]  G. Franceschetti,et al.  On the degrees of freedom of scattered fields , 1989 .

[4]  C. Chui,et al.  Article in Press Applied and Computational Harmonic Analysis a Randomized Algorithm for the Decomposition of Matrices , 2022 .

[5]  Per-Gunnar Martinsson,et al.  A fast direct solver for scattering problems involving elongated structures , 2007, J. Comput. Phys..

[6]  Dan Jiao,et al.  An �-Matrix-Based Integral-Equation Solver of Reduced Complexity and Controlled Accuracy for Solving Electrodynamic Problems , 2009 .

[7]  Zaiping Nie,et al.  An MPI-OpenMP Hybrid Parallel -LU Direct Solver for Electromagnetic Integral Equations , 2015 .

[8]  Eric Michielssen,et al.  Investigation on the butterfly reconstruction methods for MLMDA-based direct integral equation solver , 2014, 2014 USNC-URSI Radio Science Meeting (Joint with AP-S Symposium).

[9]  Eric Michielssen,et al.  An HSS Matrix-Inspired Butterfly-Based Direct Solver for Analyzing Scattering From Two-Dimensional Objects , 2016, IEEE Antennas and Wireless Propagation Letters.

[10]  Jin-Fa Lee,et al.  A fast direct matrix solver for surface integral equation methods for electromagnetic wave scattering from non-penetrable targets , 2012 .

[11]  Per-Gunnar Martinsson,et al.  Randomized algorithms for the low-rank approximation of matrices , 2007, Proceedings of the National Academy of Sciences.

[12]  J. Parrón,et al.  Multiscale Compressed Block Decomposition for Fast Direct Solution of Method of Moments Linear System , 2011, IEEE Transactions on Antennas and Propagation.

[13]  Mark Tygert,et al.  Fast algorithms for spherical harmonic expansions, II , 2008, J. Comput. Phys..

[14]  W. Hackbusch,et al.  Introduction to Hierarchical Matrices with Applications , 2003 .

[15]  Jun Hu,et al.  Hierarchical Matrices Method and Its Application in Electromagnetic Integral Equations , 2012 .

[16]  Per-Gunnar Martinsson,et al.  An O(N) Direct Solver for Integral Equations on the Plane , 2013, 1303.5466.

[17]  Per-Gunnar Martinsson,et al.  Fast direct solvers for integral equations in complex three-dimensional domains , 2009, Acta Numerica.

[18]  Jun Hu,et al.  On MLMDA/Butterfly Compressibility of Inverse Integral Operators , 2013, IEEE Antennas and Wireless Propagation Letters.

[19]  E. Michielssen,et al.  A multilevel matrix decomposition algorithm for analyzing scattering from large structures , 1996 .

[20]  Weng Cho Chew,et al.  Scattering from elongated objects: Direct solution in O(N log2 N) operations , 1996 .

[21]  Yingzhou Li,et al.  Interpolative Butterfly Factorization , 2016, SIAM J. Sci. Comput..

[22]  Yaniv Brick,et al.  Fast Direct Solver for Essentially Convex Scatterers Using Multilevel Non-Uniform Grids , 2014, IEEE Transactions on Antennas and Propagation.

[23]  Mario Bebendorf,et al.  F ¨ Ur Mathematik in Den Naturwissenschaften Leipzig Hierarchical Lu Decomposition Based Preconditioners for Bem Hierarchical Lu Decomposition Based Preconditioners for Bem , 2022 .

[24]  W. Chai,et al.  An ${\cal H}^{2}$-Matrix-Based Integral-Equation Solver of Reduced Complexity and Controlled Accuracy for Solving Electrodynamic Problems , 2009, IEEE Transactions on Antennas and Propagation.

[25]  J. Shaeffer,et al.  Direct Solve of Electrically Large Integral Equations for Problem Sizes to 1 M Unknowns , 2008, IEEE Transactions on Antennas and Propagation.

[26]  A. Boag,et al.  A multilevel fast direct solver for EM scattering from quasi-planar objects , 2009, 2009 International Conference on Electromagnetics in Advanced Applications.

[27]  Laurent Demanet,et al.  A Fast Butterfly Algorithm for the Computation of Fourier Integral Operators , 2008, Multiscale Model. Simul..

[28]  Jiming Song,et al.  Multilevel fast multipole algorithm for electromagnetic scattering by large complex objects , 1997 .

[29]  Amir Boag,et al.  Multilevel evaluation of electromagnetic fields for the rapid solution of scattering problems , 1994 .

[30]  Lexing Ying,et al.  Butterfly Factorization , 2015, Multiscale Model. Simul..