Related papers: H-Matrix Accelerated Direct Matrix Solver for Maxw…
Several researchers have developed a rich toolbox of matrix compression techniques that exploit structure and redundancy in large matrices. Classical methods such as the block low-rank format and the Fast Multipole Method make it possible…
We study the elastic time-harmonic wave scattering problems on unbounded domains with boundaries composed of finite collections of disjoints finite open arcs (or cracks) in two dimensions. Specifically, we present a fast spectral Galerkin…
In this paper, we develop an accurate and efficient Nystr\"{o}m volume integral equation (VIE) method for the Maxwell equations for large number of 3-D scatterers. The Cauchy Principal Values that arise from the VIE are computed accurately…
In the present paper, we propose Krylov-based methods for solving large-scale differential Sylvester matrix equations having a low rank constant term. We present two new approaches for solving such differential matrix equations. The first…
For scattering problems of time-harmonic waves, the boundary integral equation (BIE) methods are highly competitive, since they are formulated on lower-dimension boundaries or interfaces, and can automatically satisfy outgoing radiation…
We develop a spectral method to solve the heat equation in a closed cylinder, achieving a quasi-optimal $\mathcal{O}(N\log N)$ complexity and high-order, spectral accuracy. The algorithm relies on a Chebyshev--Chebyshev--Fourier (CCF)…
The objective of this study is to present a novel, efficient, and fast direct method for solving linear systems of equations whose coefficient matrix is a tridiagonal Quasi-Toeplitz matrix. Such matrices are frequently encountered in the…
In this paper, we develop highly accurate Nystr\"{o}m methods for the volume integral equation (VIE) of the Maxwell equation for 3-D scatters. The method is based on a formulation of the VIE equation where the Cauchy principal value of the…
This survey describes a class of methods known as "fast direct solvers". These algorithms address the problem of solving a system of linear equations $\boldsymbol{Ax}=\boldsymbol{b}$ arising from the discretization of either an elliptic PDE…
In this paper we explore the possibility for solving the 3D Maxwell's equations in the presence of nonlinear and/or inhomogeneous material response. We propose using a hybrid approach which combines a bound- ary integral representation with…
A butterfly-based fast direct integral equation solver for analyzing high-frequency scattering from two-dimensional objects is presented. The solver leverages a randomized butterfly scheme to compress blocks corresponding to near- and…
We study convex-concave saddle point problems with bilinear coupling, covering linearly constrained convex optimization and more general nonsmooth or constrained models via a proximable term in the dual objective. In linearly convergent…
A variety of problems in device and materials design require the rapid forward modeling of Maxwell's equations in complex micro-structured materials. By combining high-order accurate integral equation methods with classical multiple…
Modeling via fractional partial differential equations or a L\'evy process has been an active area of research and has many applications. However, the lack of efficient numerical computation methods for general nonlocal operators impedes…
Recent developments in high peak-power table-top laser systems reaching highly relativistic light intensities have led to significant advances in laser-driven particle acceleration schemes (mainly the laser wakefield acceleration, LWFA)…
We propose an optimization algorithm to compute the optimal sensor locations in experimental design in the formulation of Bayesian inverse problems, where the parameter-to-observable mapping is described through an integral equation and its…
In this paper, we propose fast solvers for Maxwell's equations in rectangular domains. We first discretize the simplified Maxwell's eigenvalue problems by employing the lowest-order rectangular N\'ed\'elec elements and derive the discrete…
Hierarchical matrices provide a highly memory-efficient way of storing dense linear operators arising, for example, from boundary element methods, particularly when stored in the H^2 format. In such data-sparse representations, iterative…
We describe a fourth-order accurate finite-difference time-domain scheme for solving dispersive Maxwell's equations with nonlinear multi-level carrier kinetics models. The scheme is based on an efficient single-step three time-level…
The matrix Numerov method provides an efficient framework for solving the time-independent Schr\"odinger equation as a matrix eigenvalue problem. However, for singular potentials such as the Coulomb interaction, the expected fourth-order…