  A high-order accurate, explicit kernel-split, panel-based, Fourier-Nystr\"om
discretization scheme is developed for integral equations associated with the
Helmholtz equation in axially symmetric domains. Extensive incorporation of
analytic information about singular integral kernels and on-the-fly computation
of nearly singular quadrature rules allow for very high achievable accuracy,
also in the evaluation of fields close to the boundary of the computational
domain.
