Computing weakly singular and near-singular integrals over curved boundary elements

We present algorithms for computing weakly singular and near-singular integrals arising when solving the 3D Helmholtz equation with curved boundary elements. These are based on the computation of the preimage of the singularity in the reference element’s space using Newton’s method, singularity subtraction, the continuation approach, and transplanted Gauss quadrature. We demonstrate the accuracy of our method for quadratic basis functions and quadratic triangles with several numerical experiments, including the scattering by two half-spheres.