Numerical integration over implicitly defined domains with topological guarantee

Numerical integration over the implicitly defined domains is challenging due to topological variances of implicit functions. In this paper, we use interval arithmetic to identify the boundary of the integration domain exactly, thus getting the correct topology of the domain. Furthermore, a geometry-based local error estimate is explored to guide the hierarchical subdivision and save the computation cost. Numerical experiments are presented to demonstrate the accuracy and the potential of the proposed method.