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.

Fast computation of spectral densities for generalized eigenvalue problems

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

Sublinear Time Spectral Density Estimation

Fast parallelizable methods for the Hermitian eigenvalue problem

Fast Calculation of Laplace-Beltrami Eigenproblems Via Subdivision Linear Subspace

Efficient Block Preconditioned Eigensolvers for Linear Response Time-Dependent Density Functional Theory.

The Convergence of Generalized Lanczos Methods for Large Unsymmetric Eigenproblems

Fast Non-Hermitian Toeplitz Eigenvalue Computations, Joining Matrixless Algorithms and FDE Approximation Matrices

A spectrum adaptive kernel polynomial method

Fast Spectrum Estimation of Some Kernel Matrices

Fast and backward stable computation of the eigenvalues and eigenvectors of matrix polynomials

Efficient Spectral and Spectral Element Methods for Eigenvalue Problems of Schrödinger Equations with an Inverse Square Potential

Solving Dense Generalized Eigenproblems on Multi-threaded Architectures

A projected preconditioned conjugate gradient algorithm for computing many extreme eigenpairs of a Hermitian matrix

Approximating Spectral Densities of Large Matrices

A Fast Algorithm for Computing the Spectral Radius and the Correspondent Eigenvector of the Pascal Matrices

Computing Generalized Eigenfunctions in Rigged Hilbert Spaces

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

Randomized estimation of spectral densities of large matrices made accurate

A Novel Lanczos-Type Procedure for Computing Eigenelements of Maxwell and Helmholtz Problems

Accelerated Algorithms for Eigen-Value Decomposition with Application to Spectral Clustering