Quasi-Monte Carlo for Bayesian design of experiment problems governed by parametric PDEs

Fuente: arXiv
Guardado en:
Detalles Bibliográficos
Autores principales: Kaarnioja, Vesa, Schillings, Claudia
Formato: Preprint
Publicado: 2024
Materias:
Acceso en línea:
Etiquetas: Agregar Etiqueta
Sin Etiquetas, Sea el primero en etiquetar este registro!
_version_ 1866910082552299520
author Kaarnioja, Vesa
Schillings, Claudia
author_facet Kaarnioja, Vesa
Schillings, Claudia
contents This paper contributes to the study of optimal experimental design for Bayesian inverse problems governed by partial differential equations (PDEs). We derive estimates for the parametric regularity of multivariate double integration problems over high-dimensional parameter and data domains arising in Bayesian optimal design problems. We provide a detailed analysis for these double integration problems using two approaches: a full tensor product and a sparse tensor product combination of quasi-Monte Carlo (QMC) cubature rules over the parameter and data domains. Specifically, we show that the latter approach significantly improves the convergence rate, exhibiting performance comparable to that of QMC integration of a single high-dimensional integral. Furthermore, we numerically verify the predicted convergence rates for an elliptic PDE problem with an unknown diffusion coefficient in two spatial dimensions, offering empirical evidence supporting the theoretical results and highlighting practical applicability.
format Preprint
id arxiv_https___arxiv_org_abs_2405_03529
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle Quasi-Monte Carlo for Bayesian design of experiment problems governed by parametric PDEs
Kaarnioja, Vesa
Schillings, Claudia
Numerical Analysis
65D30, 65D32, 65D40, 62K05, 62F15, 65N21
This paper contributes to the study of optimal experimental design for Bayesian inverse problems governed by partial differential equations (PDEs). We derive estimates for the parametric regularity of multivariate double integration problems over high-dimensional parameter and data domains arising in Bayesian optimal design problems. We provide a detailed analysis for these double integration problems using two approaches: a full tensor product and a sparse tensor product combination of quasi-Monte Carlo (QMC) cubature rules over the parameter and data domains. Specifically, we show that the latter approach significantly improves the convergence rate, exhibiting performance comparable to that of QMC integration of a single high-dimensional integral. Furthermore, we numerically verify the predicted convergence rates for an elliptic PDE problem with an unknown diffusion coefficient in two spatial dimensions, offering empirical evidence supporting the theoretical results and highlighting practical applicability.
title Quasi-Monte Carlo for Bayesian design of experiment problems governed by parametric PDEs
topic Numerical Analysis
65D30, 65D32, 65D40, 62K05, 62F15, 65N21
url https://arxiv.org/abs/2405.03529