Related papers: An OSRC Preconditioner for the EFIE
Convolution-type integral equations arise from various fields, \textit{e.g.}, finite impulse response filters in signal processing and deblurring problems in image processing. When solving these equations, conventional numerical methods,…
By modeling a linear, anisotropic and inhomogeneous magnetodielectric medium with two independent set of harmonic oscillators, electromagnetic field is quantized in such a medium. The electric and magnetic polarizations of the medium are…
Electrical Impedance Tomography (EIT) is a non-invasive medical imaging method that reconstructs electrical conductivity mediums from boundary voltage-current measurements, but its severe ill-posedness renders direct operator learning with…
Unless special conditions apply, the attempt to solve ill-conditioned systems of linear equations with standard numerical methods leads to uncontrollably high numerical error. Often, such systems arise from the discretization of operator…
We propose a matrix-free finite element (FE) homogenization scheme that is considerably more efficient than generic FE implementations. The efficiency of our scheme follows from a preconditioned well-scaled reformulation allowing for the…
Integral-equation-based fast direct solvers for electromagnetic scattering can substantially reduce computational costs, especially in the presence of multiple excitations. We recently proposed a new high-frequency fast direct solver…
We present a message passing approach to Expected Free Energy (EFE) minimization on factor graphs, based on the theory introduced in arXiv:2504.14898. By reformulating EFE minimization as Variational Free Energy minimization with epistemic…
Based on the perfectly matched layer (PML) technique, this paper develops a high-accuracy boundary integral equation (BIE) solver for acoustic scattering problems in locally defected layered media in both two and three dimensions. The…
We develop the shape derivative analysis of solutions to the problem of scattering of time-harmonic electromagnetic waves by a bounded penetrable obstacle. Since boundary integral equations are a classical tool to solve electromagnetic…
It is of great interest to solve the inverse problem of stationary radiative transport equation (RTE) in optical tomography. The standard way is to formulate the inverse problem into an optimization problem, but the bottleneck is that one…
Moving horizon estimation (MHE) is a widely studied state estimation approach in several practical applications. In the MHE problem, the state estimates are obtained via the solution of an approximated nonlinear optimization problem.…
The approximation of $e^{tA}B$ where $A$ is a large sparse matrix and $B$ a rectangular matrix is the key ingredient in many scientific and engineering computations. A powerful tool to manage the matrix exponential function is to resort to…
This work aims to accelerate the convergence of proximal gradient methods used to solve regularized linear inverse problems. This is achieved by designing a polynomial-based preconditioner that targets the eigenvalue spectrum of the normal…
A path integral Monte Carlo method (PIMC) based on Feynman-Kac formula for mixed boundary conditions of elliptic equations is proposed to solve the forward problem of electrical impedance tomography (EIT) on the boundary to obtain…
In this article, we present a new preconditioner, MatExPre, for the high-frequency Helmholtz equation by leveraging the properties of matrix exponentials. Our approach begins by reformulating the Helmholtz equation into a…
The goal of this work is to study the electromagnetic scattering problem of time-domain Maxwell's equations in an unbounded structure. An exact transparent boundary condition is developed to reformulate the scattering problem into an…
ITER coil tolerances are re-evaluated using the modern understanding of coupling to least-stable plasma modes and an updated center-line-traced model of ITER's coil windings. This reassessment finds the tolerances to be conservative through…
This article deals with boundary integral equation preconditioning for the multiple scattering problem. The focus is put on the single scattering preconditioner, corresponding to the diagonal part of the integral operator, for which two…
This note describes an extended exercise on the finite-element (FE) simulation of an accelerator magnet. The students construct and simulate a magnet model using the FEMM freeware. They get the opportunity to exercise on the theory of FEs,…
We present a theoretical formulation of competition among electromagnetically induced transparency (EIT) and Raman processes. The latter become important when the medium can no longer be considered to be dilute. Unlike the standard…