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.

Doubling Algorithm for the Discretized Bethe-Salpeter Eigenvalue Problem

Structure Preserving Parallel Algorithms for Solving the Bethe-Salpeter Eigenvalue Problem

Convergence Analysis of Structure-Preserving Doubling Algorithms for Riccati-Type Matrix Equations.

An Ulm-like Cayley Transform Method for Inverse Eigenvalue Problems with Multiple Eigenvalues

The Structure-preserving Doubling Algorithm for the Discrete Coupled Riccati Matrix Equations

A novel splitting strategy to accelerate solving generalized eigenvalue problem from Kohn--Sham density functional theory

Performance Enhancement of Doubling Algorithms for a Class of Complex Nonsymmetric Algebraic Riccati Equations

Convergence Analysis of the Doubling Algorithm for Several Nonlinear Matrix Equations in the Critical Case.

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

Highly accurate decoupled doubling algorithm for large-scale M-matrix algebraic Riccati equations

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

Properties of Definite Bethe-Salpeter Eigenvalue Problems

An Hp-Mixed Discontinuous Galerkin Method for the Biharmonic Eigenvalue Problem

A second-order structure-preserving discretization for the Cahn-Hilliard/Allen-Cahn system with cross-kinetic coupling

A Two-Scale Discretization Scheme for Mixed Variational Formulation of Eigenvalue Problems

Diagonal and normal with Toeplitz-block splitting iteration method for space fractional coupled nonlinear Schrödinger equations with repulsive nonlinearities

An iterative method for Helmholtz boundary value problems arising in wave propagation

Multigrid method for symmetric Toeplitz tridiagonal matrix: Convergence analysis and application to fractional Feynman-Kac equation

Generalization Error Estimates of Machine Learning Methods for Solving High Dimensional Schrödinger Eigenvalue Problems

Avoiding discretization issues for nonlinear eigenvalue problems

Local and Parallel Finite Element Discretizations for Eigenvalue Problems