Globally Adaptive Control Variate for Robust Numerical Integration

Many methods in computer graphics require the integration of functions on low-to-middle--dimensional spaces. However, no available method can handle all the possible integrands accurately and rapidly. This paper presents a robust numerical integration method, able to handle arbitrary nonsingular sca...

Full description

Saved in:
Bibliographic Details
Published inSIAM journal on scientific computing Vol. 36; no. 4; pp. A1708 - A1730
Main Authors Pajot, Anthony, Barthe, Loïc, Paulin, Mathias
Format Journal Article
LanguageEnglish
Published Society for Industrial and Applied Mathematics 01.01.2014
Subjects
Online AccessGet full text

Cover

Loading…
More Information
Summary:Many methods in computer graphics require the integration of functions on low-to-middle--dimensional spaces. However, no available method can handle all the possible integrands accurately and rapidly. This paper presents a robust numerical integration method, able to handle arbitrary nonsingular scalar or vector-valued functions defined on low-to-middle--dimensional spaces. Our method combines control variate, globally adaptive subdivision and Monte-Carlo estimation to achieve fast and accurate computations of any nonsingular integral. The runtime is linear with respect to standard deviation while standard Monte-Carlo methods are quadratic. We additionally show through numerical tests that our method is extremely stable from a computation time and memory footprint point of view, assessing its robustness. We demonstrate our method on a participating media voxelization application, which requires the computation of several millions integrals for complex media.
Bibliography:ObjectType-Article-2
SourceType-Scholarly Journals-1
ObjectType-Feature-1
content type line 23
ISSN:1064-8275
1095-7197
DOI:10.1137/130937846