Andreas Rosén

Jakob Jonsson, Andreas Rosén, Emil Timlin

We establish the theoretical foundation for a variant of the method of fundamental solutions (MFS), where the source points $\{q_j\}_{j=1}^\infty$ accumulate towards the domain in a Whitney fashion, meaning that their separation is proportional to the distance to the domain. We prove that the normalized Lusin wavelets $ψ_j(w) = b_j(w-q_j)^{-2}$ constitute a generalized basis, known as a frame, for the Hardy subspace of $L_2$-traces of holomorphic functions on the domain. Consequently, our method, where $ψ_j$ are used as basis functions in the MFS, enables a numerically stable approximation of solutions to Laplace boundary value problems, even when the solutions lack analytic continuation across the boundary. Despite the source points accumulating towards the domain, our computations achieve at least 12 digits of accuracy uniformly up to the boundary, including cases when the solution lacks analytic continuation or when the boundary has corners.

Andreas Rosén

We present a new integral equation for solving the Maxwell scattering problem against a perfect conductor. The very same algorithm also applies to sound-soft as well as sound-hard Helmholtz scattering, and in fact the latter two can be solved in parallel in three dimensions. Our integral equation does not break down at interior spurious resonances, and uses spaces of functions without any algebraic or differential constraints. The operator to invert at the boundary involves a singular integral operator closely related to the three dimensional Cauchy singular integral, and is bounded on natural function spaces and depend analytically on the wave number. Our operators act on functions with pairs of complex two by two matrices as values, using a spin representation of the fields.

Andreas Rosén

We construct new spin singular integral equations for solving scattering problems for Maxwell's equations, both against perfect conductors and in media with piecewise constant permittivity, permeability and conductivity, improving and extending earlier formulations by the author. These differ in a fundamental way from classical integral equations, which use double layer potential operators, and have the advantage of having a better condition number, in particular in Fredholm sense and on Lipschitz regular interfaces, and do not suffer from spurious resonances. The construction of the integral equations builds on the observation that the double layer potential factorises into a boundary value problem and an ansatz. We modify the ansatz, inspired by a non-selfadjoint local elliptic boundary condition for Dirac equations.