Composite quadrature formulae for the approximation of wavelet coefficients of piecewise smooth and singular functions

We develop quadrature rules for those integrals that converge fast for piecewise smooth and singular functions. They do not require the evaluation of the scaling function, and the convergence does not depend on the smoothness of that function. The analysis and computation is based completely on the filter coefficients that define the scaling function. An application is presented from the field of electromagnetics, involving the inner product of a singular function with two-dimensional tensor-product wavelets.