Laplacian Filters for Integral Equations: Further Developments and Fast Algorithms

This paper extends the concept of Laplacian filtered quasi-Helmholtz decompositions we have recently introduced, to the basis-free projector-based setting. This extension allows the discrete analyses of electromagnetic integral operators spectra without passing via an explicit Loop-Star decomposition as previously done. We also present a fast scheme for the evaluation of the filters in quasi linear complexity in the total number of unknowns. Together with the fact that only a logarithmic number of these filters are required for solving the h-refinement breakdown of electric field integral equation, this results in an effective preconditioner that rivals Calder\'on strategies in performance without relying on barycentric refinements. Numerical results confirm the theoretically predicted behavior and the effectiveness of the approach.

On a Low-Frequency and Contrast Stabilized Full-Wave Volume Integral Equation Solver for Lossy Media

In this paper we present a new regularized electric flux volume integral equation (D-VIE) for modeling high-contrast conductive dielectric objects in a broad frequency range. This new formulation is particularly suitable for modeling biological tissues at low frequencies, as it is required by brain epileptogenic area imaging, but also at higher ones, as it is required by several applications including, but not limited to, transcranial magnetic and deep brain stimulation (TMS and DBS, respectively). When modeling inhomogeneous objects with high complex permittivities at low frequencies, the traditional D-VIE is ill-conditioned and suffers from numerical instabilities that result in slower convergence and in less accurate solutions. In this work we address these shortcomings by leveraging a new set of volume quasi-Helmholtz projectors. Their scaling by the material permittivity matrix allows for the re-balancing of the equation when applied to inhomogeneous scatterers and thereby makes the proposed method accurate and stable even for high complex permittivity objects until arbitrarily low frequencies. Numerical results, canonical and realistic, corroborate the theory and confirm the stability and the accuracy of this new method both in the quasi-static regime and at higher frequencies.