Multilevel quasi-Monte Carlo for optimization under uncertainty

From MaRDI portal
Publication:6135914

DOI10.1007/S00211-023-01364-WzbMATH Open1521.65007arXiv2109.14367OpenAlexW3203728798MaRDI QIDQ6135914FDOQ6135914


Authors: Philipp A. Guth, Andreas van Barel Edit this on Wikidata


Publication date: 28 August 2023

Published in: Numerische Mathematik (Search for Journal in Brave)

Abstract: This paper considers the problem of optimizing the average tracking error for an elliptic partial differential equation with an uncertain lognormal diffusion coefficient. In particular, the application of the multilevel quasi-Monte Carlo (MLQMC) method to the estimation of the gradient is investigated, with a circulant embedding method used to sample the stochastic field. A novel regularity analysis of the adjoint variable is essential for the MLQMC estimation of the gradient in combination with the samples generated using the CE method. A rigorous cost and error analysis shows that a randomly shifted quasi-Monte Carlo method leads to a faster rate of decay in the root mean square error of the gradient than the ordinary Monte Carlo method, while considering multiple levels substantially reduces the computational effort. Numerical experiments confirm the improved rate of convergence and show that the MLQMC method outperforms the multilevel Monte Carlo method and the single level quasi-Monte Carlo method.


Full work available at URL: https://arxiv.org/abs/2109.14367




Recommendations




Cites Work


Cited In (8)





This page was built for publication: Multilevel quasi-Monte Carlo for optimization under uncertainty

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6135914)