O(N) Iterative and O(NlogN) Fast Direct Volume Integral Equation Solvers with a Minimal-Rank H2-Representation for Large-Scale 3-D Electrodynamic Analysis

Abstract

Linear complexity iterative and log-linear complexity direct solvers are developed for the volume integral equation (VIE) based general large-scale electrodynamic analysis. The dense VIE system matrix is first represented by a new cluster-based multilevel low-rank representation. In this representation, all the admissible blocks associated with a single cluster are grouped together and represented by a single low-rank block, whose rank is minimized based on prescribed accuracy. From such an initial representation, an efficient algorithm is developed to generate a minimal-rank H2-matrix representation. This representation facilitates faster computation, and ensures the same minimal rank's growth rate with electrical size as evaluated from singular value decomposition. Taking into account the rank's growth with electrical size, we develop linear-complexity H2-matrix-based storage and matrix-vector multiplication, and thereby an O(N) iterative VIE solver regardless of electrical size. Moreover, we develop an O(NlogN) matrix inversion, and hence a fast O(NlogN) direct VIE solver for large-scale electrodynamic analysis. Both theoretical analysis and numerical simulations of large-scale 1-, 2- and 3-D structures on a single-core CPU, resulting in millions of unknowns, have demonstrated the low complexity and superior performance of the proposed VIE electrodynamic solvers. %The algorithms developed in this work are kernel-independent, and hence applicable to other IE operators as well.

0

Turn this paper into a lesson

ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…