  Many applications in magnetic confinement fusion require the efficient
calculation of surface integrals with singular integrands. The singularity
subtraction approaches typically used to handle such singularities are
complicated to implement and low order accurate. In contrast, we demonstrate
that the Kapur-Rokhlin quadrature scheme is well-suited for the logarithmically
singular integrals encountered for a toroidally axisymmetric confinement
system, is easy to implement, and is high order accurate. As an illustration,
we show how to apply this quadrature scheme for the efficient and accurate
calculation of the normal component of the magnetic field due to the plasma
current on the plasma boundary, via the virtual casing principle.
