On two ways to use determinantal point processes for Monte Carlo integration
Organizations: †Univ. Lille, CNRS, Centrale Lille, UMR 9189 – CRIStAL, 59651 Villeneuve d’Ascq, France · *Inria Lille-Nord Europe, 40 avenue Halley 59650 Villeneuve d’Ascq, France · ‡DeepMind Paris, 14 Rue de Londres, 75009 Paris, France · ∗Inria Lille-Nord Europe, 40 avenue Halley 59650 Villeneuve d’Ascq, France
Abstract
The standard Monte Carlo estimator of relies on independent samples from and has variance of order . Replacing the samples with a determinantal point process (DPP), a repulsive distribution, makes the estimator consistent, with variance rates that depend on how the DPP is adapted to and . We examine two existing DPP-based estimators: one by Bardenet & Hardy (2020) with a rate of for smooth , but relying on a fixed DPP. The other, by Ermakov & Zolotukhin (1960), is unbiased with rate of order , like Monte Carlo, but its DPP is tailored to . We revisit these estimators, generalize them to continuous settings, and provide sampling algorithms.