This repository accompanies the paper "Singularity swap quadrature for nearly singular line integrals on closed curves in two dimensions" by L. af Klinteberg, arXiv:2304.11865 [math.NA].

The code is a Julia implementation of the quadrature method introduced in the paper. All plots reported in the paper can be generated using the scripts in scripts/.

Start Julia from the repo root with julia --project=.

Install all dependencies using

using Pkg; Pkg.instantiate() # First time setup
Pkg.test("TrapzSSQ")         # Test that it works

Note: The repository contains the Manifest.toml for maximum reproducibility. This was generated using Julia version 1.8.5. Instantiation may fail with other Julia versions, which can usually be resolved by running Pkg.resolve().

Run the demo scripts:

include("scripts/demo_laplace.jl")
include("scripts/show_decay.jl")
include("scripts/show_convergence.jl")

This will output the plot files in the current directory as png and pdf files.