Numerical evaluation of surface integrals in three dimensions

In this paper, we consider the evaluation of surface integrals over piecewise smooth surfaces in three dimensions. The method consists in first replacing a parametrization for the surface and the integrand function by piecewise polynomial interpolants of them, and second, using a numerical integration scheme for the resulting integral. The order of convergence is higher than would be expected based on the underlying interpolation theory