Related papers: An FFT-accelerated direct solver for electromagnet…
Poisson's equation is the canonical elliptic partial differential equation. While there exist fast Poisson solvers for finite difference and finite element methods, fast Poisson solvers for spectral methods have remained elusive. Here, we…
An efficient direct solver for solving the Lippmann-Schwinger integral equation modeling acoustic scattering in the plane is presented. For a problem with $N$ degrees of freedom, the solver constructs an approximate inverse in…
Miniaturized optical resonators with spatial dimensions of the order of the wavelength of the trapped light offer prospects for a variety of new applications like quantum processing or construction of meta-materials. Light propagation in…
This paper introduces a novel matrix-free approach for full waveform inversion in anisotropic elastic media, incorporating density variation through the utilization of the distorted Born iterative method. This study aims to overcome the…
The convolution potential arises in a wide variety of application areas, and its efficient and accurate evaluation encounters three challenges: singularity, nonlocality and anisotropy. We introduce a fast algorithm based on a far-field…
An iterative algorithm is adopted to construct approximate representations of matrices describing the scattering properties of arbitrary objects. The method is based on the implicit evaluation of scattering responses from iteratively…
We apply boundary integral equations for the first time to the two-dimensional scattering of time-harmonic waves from a smooth obstacle embedded in a continuously-graded unbounded medium. In the case we solve the square of the wavenumber…
Methods for solving Maxwell's equations are integral part of optical metrology and computational lithography setups. Applications require accurate geometrical resolution, high numerical accuracy and/or low computation times. We present a…
The paper describes a numerical method for solving acoustic multibody scattering problems in two and three dimensions. The idea is to compute a highly accurate approximation to the scattering operator for each body through a local…
This paper introduces a new method for discretizing and solving integral equation formulations of Maxwell's equations which achieves spectral accuracy for smooth surfaces. The approach is based on a hybrid Nystr\"om-collocation method using…
We study an inverse scattering problem for Maxwell's equations in terminating waveguides, where localized reflectors are to be imaged using a remote array of sensors. The array probes the waveguide with waves and measures the scattered…
Maxwell equations are solved in a layer comprising a finite number of homogeneous isotropic dielectric regions ended by anisotropic perfectly matched layers (PMLs). The boundary-value problem is solved and the dispersion relation inside the…
We analyze electromagnetic scattering of TM polarized waves from a diffraction grating consisting of a periodic, anisotropic, and possibly negative-index dielectric material. Such scattering problems are important for the modelization of,…
The interface problem describing the scattering of time-harmonic electromagnetic waves by a dielectric body is often formulated as a pair of coupled boundary integral equations for the electric and magnetic current densities on the…
This paper is concerned with a numerical solution to the scattering of a time-harmonic electromagnetic wave by a bounded and impenetrable obstacle in three dimensions. The electromagnetic wave propagation is modeled by a boundary value…
An implicit causal space-time Galerkin scheme applied to the contrast current density volume integral equation gives rise to a marching-on-in-time scheme known as the MOT-JVIE, which is accelerated and stabilized via a fully embedded FIR…
A qubit lattice algorithm (QLA) is developed for Maxwell equations in a two-dimensional Cartesian geometry. In particular, the initial value problem of electromagnetic pulse scattering off a localized 2D dielectric object is considered. A…
In this paper, we present a numerical algorithm for the accurate and efficient computation of the convolution of the frequency domain layered media Green's function with a given density function. Instead of compressing the convolution…
A fast method for the computation of layer potentials that arise in acoustic scattering is introduced. The principal idea is to split the singular kernel into a smooth and a local part. The potential due to the smooth part is computed…
The advent of massively parallel supercomputers, with their distributed-memory technology using many processing units, has favored the development of highly-scalable local low-order solvers at the expense of harder-to-scale global very…