Fast Update of Conditional Simulation Ensembles

Abstract : Gaussian random fields (GRF) conditional simulation is a key ingredient in many spatial statistics problems for computing Monte-Carlo estimators and quantifying uncertainties on non-linear functionals of GRFs conditional on data. Conditional simulations are known to often be computer intensive, especially when appealing to matrix decomposition approaches with a large number of simulation points. Here we study the settings where conditioning observations are assimilated batch-sequentially, i.e. one point or batch of points at each stage. Assuming that conditional simulations have been performed at a previous stage, we aim at taking advantage of already available sample paths and by-products in order to produce updated conditional simulations at minimal cost. We provide explicit formulas allowing to update an ensemble of sample paths conditioned on $n\geq 0$ observations to an ensemble conditioned on $n+q$ observations, for arbitrary $q\geq 1$. Compared to direct approaches, the proposed formulas prove to substantially reduce computational complexity. Moreover, these formulas enable explicitly exhibiting how the $q$ ''new'' observations are updating the ''old'' sample paths. Detailed complexity calculations highlighting the benefits of our approach with respect to state-of-the-art algorithms are provided and are complemented by numerical experiments.
Liste complète des métadonnées
Contributeur : David Ginsbourger <>
Soumis le : lundi 28 avril 2014 - 14:58:39
Dernière modification le : lundi 21 mars 2016 - 17:38:20
Document(s) archivé(s) le : lundi 28 juillet 2014 - 11:45:41


Fichiers produits par l'(les) auteur(s)


  • HAL Id : hal-00984515, version 1



Clément Chevalier, Xavier Emery, David Ginsbourger. Fast Update of Conditional Simulation Ensembles. 2014. <hal-00984515>



Consultations de
la notice


Téléchargements du document