A multiresolution approach to regularization of singular operators and fast summation

Singular and hypersingular operators are ubiquitous in problems of physics, and their use requires a careful numerical interpretation. Although analytical methods for their regularization have long been known, the classical approach does not provide numerical procedures for constructing or applying the regularized operator. We present a multiresolution definition of regularization for integral operators with convolutional kernels which are homogeneous or associated homogeneous functions. We show that our procedure yields the same operator as the classical definition. Moreover, due to the constructive nature of our definition, we provide concise numerical procedures for the construction and application of singular and hypersingular operators. As an application, we present an algorithm for fast computation of discrete sums and briefly discuss several other examples.