Abstract:In this paper, we present Gauss's law-preserving spectral methods and their efficient solution algorithms for curl-curl source and eigenvalue problems in two and three dimensions arising from Maxwell's equations. Arbitrary order $H(curl)$-conforming spectral basis functions in two and three dimensions are firstly proposed using compact combination of Legendre polynomials. A mixed formulation involving a Lagrange multiplier is then adopted to preserve the Gauss's law in the weak sense. To overcome the bottleneck of computational efficiency caused by the saddle-point nature of the mixed scheme, we present highly efficient solution algorithms based on reordering and decoupling of the resultant linear algebraic system and numerical eigen-decomposition of one dimensional mass matrix. The proposed solution algorithms are direct methods requiring only several matrix-matrix or matrix-tensor products of $N$-by-$N$ matrices, where $N$ is the highest polynomial order in each direction. Compared with other direct methods, the computational complexities are reduced from $O(N^6)$ and $O(N^9)$ to $O(N^3)$ and $O(N^4)$ with small and constant pre-factors for 2D and 3D cases, respectively, and can further be accelerated to $O(N^{2.807})$ and $O(N^{3.807})$, when boosted with the Strassen's matrix multiplication algorithm. Moreover, these algorithms strictly obey the Helmholtz-Hodge decomposition, thus totally eliminate the spurious eigen-modes of non-physical zero eigenvalues. Extensions of the proposed methods and algorithms to problems in complex geometries with variable coefficients and inhomogeneous boundary conditions are discussed to deal with more general situations. Ample numerical examples for solving Maxwell's source and eigenvalue problems are presented to demonstrate the accuracy and efficiency of the proposed methods.

Localized spectrum slicing

Localized subspace iteration methods for elliptic multiscale problems

Enhancing Scalability of a Matrix-Free Eigensolver for Studying Many-Body Localization

Enhancing Scalability of a Matrix-Free Eigensolver for Studying Many-Body Localization

Parallel Eigenvalue Calculation Based On Multiple Shift-Invert Lanczos And Contour Integral Based Spectral Projection Method

Spectral ACMS: A robust localized Approximated Component Mode Synthesis Method

Estimation of spectral gaps for sparse symmetric matrices

Efficient Variational Quantum Linear Solver for Structured Sparse Matrices

Spectral Decomposition of Real Symmetric Quadratic Λ-Matrices and Its Applications.

Fast computation of spectral densities for generalized eigenvalue problems

Spectral Theory of Sparse Non-Hermitian Random Matrices

A class of sparse Johnson--Lindenstrauss transforms and analysis of their extreme singular values

Spectral Transformation for the Dense Symmetric Semidefinite Generalized Eigenvalue Problem

Randomized estimation of spectral densities of large matrices made accurate

Hybrid Localized Spectral Decomposition for multiscale problems

Large-Scale Minimization of the Pseudospectral Abscissa

Highly efficient Gauss's law-preserving spectral algorithms for Maxwell's double-curl source and eigenvalue problems based on eigen-decomposition

A sparse decomposition of low rank symmetric positive semi-definite matrices

Inverse eigenvalue problem for Laplacian matrices of a graph

A Novel Symmetric Skew-Hamiltonian Isotropic Lanczos Algorithm for Spectral Conformal Parameterizations

The Hierarchical Subspace Iteration Method for Laplace--Beltrami Eigenproblems